Association of Working Memory With Distributed Executive Control Networks in Schizophrenia
Abstract
Objective:
Working memory impairments represent a core cognitive deficit in schizophrenia, predictive of patients’ daily functioning, and one that is unaffected by current treatments. To address this, working memory is included in the MATRICS Consensus Cognitive Battery (MCCB), a standardized cognitive battery designed to facilitate drug development targeting cognitive symptoms. However, the neurobiology underlying these deficits in MCCB working memory is currently unknown, mirroring the poor understanding in general of working memory deficits in schizophrenia.
Methods:
Twenty-eight participants with schizophrenia were administered working memory tests from the MCCB and examined with resting-state functional MRI. Intrinsic connectivity networks were estimated with independent component analysis. Each voxel’s time series was correlated with each network time series, creating a feature vector for voxel-level connectivity analysis. This feature vector was associated with working memory by using the distance covariance statistic.
Results:
The neurobiology of MCCB working memory tests largely followed the multicomponent model of working memory but revealed unexpected differences. The dorsolateral prefrontal cortex was not associated with working memory. The central executive system was instead associated with delocalized right and left executive control networks. The phonologic loop within the multicomponent model, a subsystem involved in storing linguistic information, was associated with connectivity to the left temporoparietal junction and inferior frontal gyrus. However, connections to the language network did not predict working memory test performance.
Conclusions:
These results provide supporting evidence for the multicomponent model of working memory in terms of the biology underlying MCCB findings.
Individuals with schizophrenia experience symptoms that include intrusive thoughts and sensory stimuli, psychosocial deficits in emotional expression, and cognitive deficits, generally referred to as positive, negative, and cognitive symptoms, respectively. Although antipsychotic medications successfully control symptoms related to intrusive sensory stimuli such as hallucinations, they are not able to treat the emotional and cognitive symptoms of schizophrenia. The success of these medications in controlling the overt symptoms of psychosis, contrasted with their frequent failure to restore patients to their level of functioning before the onset of illness, highlights the distress caused by other core symptoms of schizophrenia. For example, from the patient’s perspective, medical treatments that only control hallucinations are of limited value if other untreated symptoms prevent the patient from socializing and pursuing life and career goals. These other core symptoms include untreated and currently untreatable cognitive deficits.
Impaired working memory is a core cognitive deficit in schizophrenia. Problems in working memory, defined as the temporary storage and manipulation of information, are common in patients with the disease (1, 2). They predate the onset of psychotic symptoms and may underlie deficits in other cognitive domains (3, 4). Most important, they are strongly predictive of impairments in social and vocational outcomes such as employment status, educational achievement, and independent living (5, 6). Specifically, in adults with schizophrenia, working memory test performance predicts work or education functioning (6). In adolescents with schizophrenia, baseline working memory performance is significantly related to social, communication, personal, and community living skills at a follow-up visit one year later (5). No current treatment is able to ameliorate the distress caused by this core cognitive deficit (7). Consequently, working memory is included as a domain in the MATRICS Consensus Cognitive Battery (MCCB). The MCCB is a comprehensive battery of cognitive tests designed for and extensively validated in subjects with schizophrenia in order to facilitate drug development targeting cognitive symptoms (8). Although the neurobiology of many cognitive domains in the MCCB is currently unknown, the functional neuroanatomy of working memory in other contexts has been studied in detail in the field of cognitive neuroscience (9). However, the neurobiology of the specific tests in the MCCB and how this neurobiology relates to other tests of working memory such as the n-back or item recognition tests is currently unknown.
The neurobiology of working memory is associated with prefrontal brain regions, including the dorsolateral prefrontal cortex (dlPFC), the ventrolateral prefrontal cortex (vlPFC) encompassing the inferior frontal gyrus (IFG), as well as posterior parietal regions, including the intraparietal sulcus (IPS) and temporoparietal junction (TPJ) (10, 11). Within the prefrontal cortices, the primary locus of activity differs by task. The dlPFC seems more involved during complex cognitive tasks, including the n-back, in which subjects monitor a series of stimuli while deciding if the currently displayed item matches the stimulus presented n stimuli previously (12). In contrast, the vlPFC is more involved during tasks involving simpler mental manipulations, such as the Sternberg Item Recognition Paradigm, in which subjects memorize a set of stimuli and, after a short delay, determine if a probe stimulus was included in this memorized set (10, 11). The n-back is cognitively demanding, requiring continual updating of contents stored in memory. In contrast, item recognition tasks typically involve mental manipulations of static items stored in memory, with lower processing demands. These different processing demands may account for the differences in the neurobiological responses to each task.
Similarly, the neurobiology of working memory deficits in schizophrenia also differs depending on the specific task, with both increased or decreased prefrontal response observed during different tasks. For example, compared with healthy controls, the n-back test elicits decreased dlPFC activity in patients with schizophrenia, whereas the Sternberg Item Recognition Paradigm elicits increased dlPFC activity (13, 14). Similarly, the n-back test is associated with decreased vlPFC response in patients, including in the IFG, whereas the Sternberg Item Recognition Paradigm elicits increased IFG response (13, 15).
In summary, the neurobiology of working memory deficits in schizophrenia depends on the nature of the specific test, with both increased and decreased prefrontal activity in response to different tasks. Given these seemingly contradictory results when comparing other working memory tests, the neurobiology that underlies the working memory deficit that the MCCB is designed to measure is currently unknown, despite a wealth of information on working memory in general and working memory in schizophrenia in particular. More sophisticated analytic methods are required to better understand the neurobiology of the MCCB.
The intrinsic network architecture within the brain may provide insight into the neurobiology of working memory as assessed by the MCCB. In the emerging consensus of whole-brain neural connectivity, the brain is divided into segregated functional processing systems, or intrinsic connectivity networks (ICNs). These ICNs include networks dedicated to processing unimodal sensory information, such as the auditory network, as well as networks dedicated to higher-level processing, including the left executive central network (LECN) and right executive central network (RECN), respectively (16). Interactions between ICNs occur at discrete anatomical locations within the brain, frequently referred to as processing hubs. This intrinsic functional architecture is stable across processing states or psychological tasks, forming a baseline organization of connectivity within the brain (17). This fundamental network architecture facilitates specific task performance due to the efficient arrangement of its processing pathways (18). Task-related activity induces minor rearrangements in the intrinsic network architecture on the basis of processing demands (19). It is important to note that because this architecture is relatively constant during rest or active processing, investigations using resting-state functional MRI (fMRI) may be able to provide insights into the networks and regions where intrinsic network architecture predicts task performance.
This study investigates the neurobiology of MCCB-measured working memory deficits in schizophrenia by examining their relationship to intrinsic network architecture with resting-state fMRI and a novel network analysis technique, distance covariance (dCov) (20). This technique investigates multivariate associations between cognitive measures, neuroanatomy, and the delocalized processing systems that contribute to working memory. Results are interpreted using the framework provided by Baddeley’s multicomponent model of working memory (21).
The multicomponent model consists of separate storage buffers for different information types—the visuospatial sketchpad, the episodic buffer, and the phonologic loop—in combination with a central executive control system that acts on the information contained in the storage buffers (Figure 1). Each system corresponds to known functional neuroanatomy (21, 22). For example, verbal information is stored in language areas of the phonologic loop, such as Wernicke’s and Broca’s areas, located in the left TPJ and IFG, respectively, in most subjects. Similarly, the central executive is commonly localized to the bilateral dlPFC.
Recent investigations call into question the previously exclusive and one-to-one relationship between neuroanatomy and function in the multicomponent model. These studies suggest that the central executive system, in addition to the dlPFC, also includes the vlPFC as well as superior parietal regions, such as the IPS (10, 23). This suggests that the central executive system, rather than being localized to a single region, may instead be a distributed and delocalized processing system that extends beyond neuroanatomical boundaries to encompass multiple regions. Other large-scale ICNs encompass entire neural processing systems and are active during cognitive states ranging from unfocused rest to most if not all deliberative cognitive tasks. In fact, two prominent ICNs, the RECN and LECN, encompass the prefrontal and parietal regions attributed to the delocalized central executive system. The study of these large-scale networks, along with their effect on cognitive processing, is a current topic of active investigation in resting-state fMRI.
The present study examines MCCB working memory using resting-state fMRI, focusing on delocalized networks and the connections included in the multicomponent model (Figure 2). We tested connections between all ICNs and all voxels in order to comprehensively examine how intrinsic network architecture contributes to working memory performance in schizophrenia.
We hypothesized that the multicomponent model (Figure 1) would be reflected in the combined interactions of delocalized networks with anatomical regions (Figure 2). Specifically, we hypothesized that MCCB working memory in schizophrenia would be associated with connectivity to the left IFG, vlPFC, TPJ, and right IPS, and that the central executive system would be best represented by the RECN and LECN. These hypotheses were tested using a novel combination of multivariate statistical techniques.
Methods
Participants
Thirty-three subjects with schizophrenia or schizoaffective disorder were recruited via advertisements in the community and through outpatient clinics (Table 1). All patients met DSM-IV criteria for schizophrenia or schizoaffective disorder, confirmed by the Structured Clinical Interview for DSM-IV. Seven subjects were treated with first-generation antipsychotic medications, 23 with second-generation antipsychotics, and two were nonmedicated. All subjects were clinically stable outpatients and were treated with the same antipsychotic regimen for at least three months. Five subjects were excluded due to excessive movement (>3 mm) during scanning. Exclusion criteria included a current diagnosis of substance abuse, neurological disorders, or head trauma, as well as MRI exclusion factors (claustrophobia, weight >300 lb, metal in the body). Only decision-capable individuals were eligible for study participation. All subjects provided written informed consent after receiving a complete description of the study, and they received compensation for their participation. The study was approved by the Colorado Multiple Institutional Review Board.
Characteristic | N or M | % or SD |
---|---|---|
Sex (female) | 7 | 25.0 |
Ethnicity | ||
African American | 5 | 17.9 |
Hispanic | 2 | 7.1 |
Age (years) (M±SD) | 42.4 | 12.6 |
Education (years) (M±SD) | 13.4 | 1.81 |
Smoker | 14 | 50.0 |
Diagnosis | ||
Schizophrenia | 20 | 71.4 |
Schizoaffective | 8 | 28.6 |
Medications | ||
Typical antipsychotics | 7 | 25.0 |
Atypical antipsychotics | 23 | 82.1 |
No medication | 2 | 7.1 |
Brief Psychiatric Rating Scale total (M±SD) | 34.8 | 8.1 |
Scale for the Assessment of Negative Symptoms total (M±SD) | 4.36 | 2.9 |
Demographic and clinical characteristics of the study subjects (N=28)
Cognitive Measures
The MCCB has been described in detail elsewhere (8). Briefly, the battery consists of individual tests covering the cognitive domains of working memory (Wechsler Memory Scale–III, spatial span and letter-number span [LNS] subtests), speed of processing (Trail-Making Test, Part A, Brief Assessment of Cognition in Schizophrenia, symbol coding subtest, and Category Fluency Test), attention and vigilance (Continuous Performance Test–Identical Pairs Version), verbal learning (Hopkins Verbal Learning Test–Revised), visual learning (Brief Visuospatial Memory Test–Revised), reasoning and problem solving (Neuropsychological Assessment Battery, mazes subtest), and social cognition (Mayer-Salovey-Caruso Emotional Intelligence Test, managing emotions subtest). Raw scores were converted to normalized T scores, corrected for age and gender, and centered and scaled relative to a healthy population of subjects using the MCCB scoring program (24). The MCCB and fMRI scans were administered on the same day. Only scores from the MCCB working memory domain and the two domain subtests were used for the present analysis.
For comparison with the multicomponent model (Figure 1), the verbal subtest of MCCB working memory (LNS) was considered to primarily test functioning of the phonologic loop, whereas the nonverbal subtest (Wechsler Memory Scale–III, spatial span) was considered to test the visuospatial sketchpad. Additionally, both working memory subtests were expected to involve the central executive system. However, consistent with the multicomponent model, no subtest uniquely tested this system. Instead, the central executive system was identified on the basis of its association with the dlPFC and executive control networks. No MCCB working memory subtest was considered to adequately test the functioning of the episodic buffer, although a disputable association with the hippocampus has been suggested (25).
fMRI Parameters
Images were acquired on a 3-T scanner (General Electric, Milwaukee) using a standard quadrature head coil. An inversion-recovery echo planar image (TI=505 ms) was collected to improve coregistration of functional images. Functional scans were acquired with the following parameters: TR=2000 ms, TE=30 ms, field of view=240 mm2, matrix=64×64, voxel size=3.75×3.75 mm2, slice thickness=2.6 mm, gap=1.4 mm, interleaved, flip angle=70°. Resting fMRI scan duration was 10 minutes. Participants were instructed to rest with their eyes open.
Data Analysis
Preprocessing.
MRI data were preprocessed using SPM8 (http://www.fil.ion.ucl.ac.uk/spm/software/spm8) for each subject individually. The first four images were excluded for saturation effects. Echo planar images from each subject were realigned to the first volume, resliced to a 3-mm3 voxel size, normalized to the Montreal Neurological Institute template by using unified segmentation (26), and smoothed with an 8-mm full width at half maximum Gaussian kernel (phase 1) (Table 2).
Item | Description |
---|---|
Preprocessingb |
|
Network identificationc |
|
Voxel-to-network functional connectivityd |
|
Connectivity associated with working memorye |
|
Identification of networks that contribute to working memorye |
|
Independent component analysis.
Large-scale networks (ICNs) were identified by using independent component analysis (ICA) (phase 2) (Table 2). Spatial ICA was carried out on group-level data with GIFT v1.3i (http://icatb.sourceforge.net) (27). Twenty-one components were estimated on the basis of minimum description length criteria and extracted by using the infomax algorithm (28, 29). Voxel time series were whitened, variance normalized, and temporally concatenated with two principle component analysis data reduction steps of 70 and 21 components. Spatial maps were reconstructed with GICA3 and scaled to z scores (30). All spatial maps and time courses were visually inspected to identify noise components. Seven components classified as artifacts on the basis of spatial distributions in CSF or white matter or resulting from high-frequency oscillations were excluded from further analysis. To identify common ICNs, group mean ICA spatial maps were correlated with published ICA templates (16). Templates matching multiple ICA components were identified as subnetworks, whereas ICA components without template matches were classified on the basis of anatomy. For example, two ICA components with minimal spatial overlap matched the template for the dorsal default mode network (DMN). These were labeled as the anterior DMN and posterior DMN on the basis of anatomical differences in their spatial maps (see supplementary materials). Another ICA component did not clearly match any template but was strongly localized to the bilateral hippocampus and amygdala and was labeled as the bilateral medial temporal lobe network on the basis of the corresponding anatomy. Following ICA, a back-reconstructed subject-specific time series for each ICN was correlated with voxels’ time series in a functional connectivity analysis.
Functional connectivity.
CSF, white matter signals, and six movement parameters were regressed out of all time series individually for each subject. Time points with excessive movement were censored (31). Following nuisance signal removal and movement control precautions, time series for every gray matter voxel were correlated with each ICA component time series for each subject individually (phase 3) (Table 2). The resulting vector of correlations for each voxel efficiently and comprehensively summarizes a voxel’s connectivity to all networks within the brain; it was used as a feature vector after removing the influences of age, smoking status, and education level using linear regression residuals (32).
Distance covariance statistic.
dCov is a recently developed multivariate technique that tests the statistical independence of two vectors with arbitrary dimensions (20). The dCov statistic is zero if and only if the random vectors are independently distributed and increasingly positive otherwise. A complete description of dCov is available in the supplementary materials, including its statistical properties, details of its calculation, a sample use, and application in the present analysis.
dCov analyses were conducted using two methods. First, dCov was calculated for each gray matter voxel using subjects’ connectivity vector and MCCB scores as input (phase 4) (Table 2). This resulted in a nonthresholded whole-brain statistical map, showing associations between voxel-level connectivity and MCCB scores in the sample. Results were then corrected for multiple comparisons using a modification of Nichol’s cluster-level nonparametric test (33). Five thousand volumes were generated by permuting subjects’ MCCB scores for each voxel. Empirical p values were calculated as the proportion of times the permutation dCov statistic exceeded the original dCov statistic at that voxel. Empirical p values for all permutation voxels were calculated in an identical manner, and volumes were thresholded using a cluster-defining threshold at a p value <0.01, and resulting cluster sizes were calculated. A cluster was considered significant after correcting for multiple comparisons if the empirical p value of obtaining a cluster this size or larger was significant at a p value <0.05. Significant clusters were displayed as cortical surface projections, created using Caret (34).
Distance covariance was then applied to test connectivity using the elements of the individual vector, using an individual scalar correlation coefficient and the scalar MCCB score (phase 5) (Table 2).This procedure tested the individual elements of the connectivity vector that contributed to that voxel’s significance in the previous whole-brain analysis. Resulting dCov statistics for connectivity were thresholded and corrected for multiple comparisons using a cluster-level permutation procedure developed for network analysis, the network-based statistic (35). Connections were thresholded at a p value <0.005. An individual connection was considered significant after correcting for multiple comparisons if it was included in a cluster size with statistical significance at a p value <0.05. Significant connections were displayed as two-dimensional sagittal, coronal, and axial projections in Montreal Neurological Institute space, the “glass brain” format.
Results
MCCB Working Memory
All MCCB scores were corrected for age and gender, as well as centered and scaled relative to MCCB scores from a healthy population of subjects (mean=50 [SD=10]) (24). Subjects with schizophrenia demonstrated deficits in MCCB working memory domain (mean=43.43 [SD=11.56], one-sample t test with H0: mu=50: t=−3.01, df=27, p<0.006). Performance on the LNS test of verbal working memory was diminished (mean=41.07 [SD=12.01], one-sample t test with H0: mu=50: t=−3.93, df=27, p<0.006), whereas deficits on the test of nonverbal working memory (WMS-III spatial span) were less pronounced (mean=48.21 [SD=10.08], one-sample t test with H0: mu=50: t=−0.94, df=27, p>0.10). Full results for all MCCB cognitive domains and tests are reported in the supplementary materials.
Intrinsic Connectivity Networks
Fourteen ICNs matched templates for known networks to a relatively high degree (mean r=0.38). The RECN and LECN were present, encompassing the dlPFC and vlPFC along with right or left lateralized IPS (Figure 3). The language network encompassed the bilateral vlPFC (including the left IFG), TPJ, superior temporal gyri, and primary auditory cortices. An ICN consisting of bilateral medial temporal regions included the hippocampus.
Connectivity and MCCB Working Memory
Associations between working memory and connectivity were examined in a two-step procedure. First, we focused on associations to neuroanatomy, investigating where in the brain connectivity influences task performance. Second, after these regions were identified, we investigated which connections to and from the regions most contributed to task performance. This analysis highlighted both how different subregions within a larger network made different contributions to working memory and how each region’s interactions with multiple networks contributed to task performance.
Working memory and neuroanatomy.
Voxel-level whole-brain analysis of connectivity identified associations between MCCB working memory and several regions throughout the brain (Table 3, Figure 4). These included the vlPFC bilaterally as well as the left TPJ. Smaller clusters were located in the left hippocampus, dorsal anterior cingulate cortex, middle temporal gyrus, and the left anterior cerebellum extending into the inferior temporal gyrus. Smaller clusters in the right hemisphere were located along the right IPS, superior temporal sulcus, middle cingulate gyrus, presupplementary motor area, posterior hippocampus, and parietal operculum. Interestingly, dlPFC connectivity was not associated with MCCB working memory.
Montreal Neurological Institute coordinate | |||||
---|---|---|---|---|---|
Region | Statistic (distance covariance) | x | y | z | Cluster size (voxel) |
Left supramarginal gyrus | 2.08 | –51 | –49 | 28 | 107 |
Left angular gyrus | 1.76 | –57 | –55 | 28 | |
Left supramarginal gyrus | 1.64 | –45 | –43 | 25 | |
Left angular gyrus | 1.59 | –48 | –58 | 28 | |
Left supramarginal gyrus | 1.52 | –57 | –49 | 22 | |
Right inferior parietal lobule | 1.91 | 51 | –46 | 55 | 76 |
Right inferior parietal lobule | 1.73 | 45 | –52 | 55 | |
Right inferior parietal lobule | 1.61 | 51 | –52 | 49 | |
Right inferior parietal lobule | 1.59 | 45 | –46 | 49 | |
Right angular gyrus | 1.47 | 45 | –58 | 49 | |
Right inferior parietal lobule | 1.43 | 42 | –52 | 43 | |
Right orbitofrontal cortex | 1.83 | 45 | 47 | –2 | 122 |
Right middle frontal gyrus | 1.67 | 39 | 50 | 4 | |
Right superior frontal gyrus | 1.65 | 30 | 50 | 13 | |
Right orbitofrontal cortex | 1.55 | 48 | 44 | –11 | |
Right middle frontal gyrus | 1.55 | 39 | 53 | 13 | |
Right superior frontal gyrus | 1.46 | 27 | 56 | 19 | |
Right superior frontal gyrus | 1.41 | 21 | 53 | 13 | |
Right superior temporal gyrus | 1.78 | 63 | –49 | 22 | 42 |
Right middle temporal gyrus | 1.55 | 51 | –49 | 22 | |
Right superior temporal gyrus | 1.46 | 57 | –43 | 19 | |
Left orbitofrontal cortex | 1.71 | –36 | 38 | –14 | 32 |
Left orbitofrontal cortex | 1.64 | –39 | 44 | –8 | |
Left orbitofrontal cortex | 1.51 | –45 | 41 | –14 | |
Left orbitofrontal cortex | 1.37 | –33 | 29 | –14 | |
Left middle temporal gyrus | 1.71 | –66 | –34 | –11 | 18 |
Left middle temporal gyrus | 1.62 | –63 | –25 | –14 | |
Right middle cingulate cortex | 1.67 | 6 | –19 | 40 | 13 |
Right superior frontal gyrus | 1.66 | 27 | 38 | 46 | 13 |
Left inferior parietal lobule | 1.64 | –36 | –43 | 37 | 15 |
Left precentral gyrus | 1.63 | –33 | 11 | 34 | 22 |
Left precentral gyrus | 1.46 | –42 | 14 | 34 | |
Left orbitofrontal cortex | 1.61 | –54 | 32 | –2 | 23 |
Left orbitofrontal cortex | 1.44 | –51 | 35 | –11 | |
Right middle temporal gyrus | 1.6 | 54 | –58 | 16 | 16 |
Right middle temporal gyrus | 1.59 | 66 | –40 | 4 | 10 |
Right inferior frontal gyrus | 1.57 | 54 | 26 | 13 | 58 |
Right inferior frontal gyrus | 1.55 | 48 | 32 | 19 | |
Right inferior frontal gyrus | 1.49 | 45 | 41 | 28 | |
Right middle temporal gyrus | 1.55 | 66 | –19 | –8 | 11 |
Left cuneus | 1.55 | 0 | –79 | 34 | 13 |
Right angular gyrus | 1.54 | 54 | –64 | 37 | 5 |
Right medial prefrontal cortex | 1.53 | 12 | 56 | 31 | 16 |
Left cerebellar lobule IV, V | 1.52 | –15 | –52 | –14 | 10 |
Left fusiform gyrus | 1.45 | –21 | –46 | –14 | |
Right middle frontal gyrus | 1.51 | 45 | 17 | 52 | 6 |
Left medial prefrontal cortex | 1.45 | –6 | 35 | 40 | 8 |
Left medial prefrontal cortex | 1.43 | 3 | 35 | 43 |
Local peak multivariate associations in MATRICS Consensus Cognitive Battery working memory domains
Working memory and network connectivity.
In order to better understand how connectivity to the regions in Figure 5 contributed to working memory, we then examined the individual connections of these voxels.
Many of the systems of the multicomponent working memory model were significantly associated with MCCB working memory performance (Figure 5). These included the RECN and LECN for the central executive system, the right vlPFC and IPS for the visuospatial sketchpad, and the left vlPFC and TPJ for the phonologic loop. The small cluster within the hippocampus was not significantly connected to either the LECN or RECN (p>0.05 corrected).
In addition to connectivity predicted by the multicomponent model, other networks and regions contributed to MCCB working memory. Notably, and beyond the scope of the multicomponent model and the networks detailed above, the DMN influenced working memory through the IFG bilaterally. In contrast, language and bilateral medial temporal lobe ICNs minimally influenced working memory. See supplementary materials for full results, showing all voxels and connections significant after correcting for multiple comparisons.
Associations With Verbal and Nonverbal Working Memory
Results for the individual tests composing the MCCB working memory domain were largely similar to above, with minor differences between verbal and nonverbal subtests (see supplementary materials). For the verbal LNS subtest, anatomical results were strongly lateralized to the left hemisphere and largely limited to regions within the phonologic loop. For the nonverbal WMS-III spatial span, anatomical results were bilateral and included more extensive clusters throughout the cortex. In both cases, network results mirrored the composite MCCB working memory domain.
Discussion
Overall, associations between working memory task performance and functional connectivity were in good agreement with the predictions of the hypothesis and the multicomponent model (Figures 2 and 5). A result that was not expected, however, was that the dlPFC was not observed to directly influence MCCB working memory (Figure 4). Instead, dlPFC contributions to working memory were indirect, as parts of the RECN and LECN (Figure 3). While perhaps counterintuitive, these results are consistent with the concept of a distributed and delocalized central executive system, associated with but not limited to prefrontal regions (23, 36). Similarly, the RECN influenced working memory through the left TPJ as well (Figure 5), even though this region was not included in this network. Alternatively, working memory tasks in the MCCB may not require processing within the dlPFC. Many previous investigations of working memory were carried out with the n-back, a relatively difficult task requiring continual updating of the information contained within working memory (12). Tasks with low cognitive load often preferentially activate the vlPFC instead, with dlPFC active during high processing demands (10, 37). In this view, the absence of dlPFC may reflect the low cognitive demands of MCCB tests. Lastly, these results do not contradict previous associations of the central executive system with the dlPFC. Instead, they suggest a delocalized central executive system that includes the dlPFC but extends beyond any single cortical region.
The phonologic loop was strongly left lateralized, agreeing with the multicomponent model (Figure 2). Consistent with the predictions of the model, this was especially true for the verbal LNS test with associated large clusters in the left IFG and TPJ (for further details, see the online supplement). Interestingly, the language network itself was minimally associated with working memory. On the basis of these results, the phonologic loop is likely separate from the delocalized language network.
Regions of the visuospatial sketchpad, such as the right vlPFC and IPS, were present in the results, but the interpretation was not as straightforward. Processing during the nonverbal WMS-III spatial span subtest would be expected to rely on the visuospatial sketchpad. Clusters in the right vlPFC and IPS were associated with performance on this subtest (see supplementary materials). However, these associations were bilateral, a result consistent with meta-analysis of activation during working memory fMRI tasks (10).
It is unclear if the episodic buffer contributed to MCCB working memory tests. Although this system has been associated with the hippocampus, the supporting evidence is less well established (25). In the present investigation, an association to working memory was found to a small cluster within the hippocampus (Figure 4). However, these voxels were not influenced by connections from the RECN or LECN (Figure 5). Since none of the MCCB tests involved the constructions of mental scenes from episodic long-term memory, the episodic buffer is unlikely to play a prominent role and a lack of any clear associations is unsurprising, if uninformative.
Previous neuroimaging investigations into working memory in schizophrenia have primarily focused on activity resulting from the n-back fMRI task. This task is notable for its robust but nonspecific activations, and its dependence on dlPFC functioning (10, 38). In contrast, the association with the vlPFC but not the dlPFC may suggest that these tests are more likely to involve transformation of stored information and inhibition of irrelevant stimuli (10) (Figure 4). Alternatively, an absence of associations between dlPFC connectivity and MCCB working memory may also indicate that dlPFC dysfunction during the n-back in schizophrenia may extend these working memory tests as well (13). Future investigations will seek to clarify this distinction, as well as to investigate the relationship between connectivity and the impact of working memory dysfunction in schizophrenia.
The present study includes several limitations. Only subjects with a diagnosis of schizophrenia or schizoaffective disorder were included. This study design was chosen because the primary aim was to investigate performance on a cognitive test designed for patients with schizophrenia and validated in this population (8). The lack of a healthy comparison group does, however, limit study interpretation. Any extrapolation to the neurobiology of MCCB working memory in healthy subjects is premature, as results could be due to disease-specific neurobiology. The present analysis was limited to the MCCB working memory domain. Analysis of other MCCB cognitive domains and generalization of these results will be the focus of future investigations. An additional limitation is the modest sample size. Although likely adequately powered for many types of fMRI studies (39–41), the sample size may be small for the present analysis. The risk of false positive results cannot be ruled out; however, false negatives are a greater concern in this analysis, given the sample size (40). Lastly, the validity of the central executive system itself as a concept in psychology is not without controversy (36, 42, 43).
Conclusions
We have demonstrated that performance on the MCCB in subjects with schizophrenia can be predicted from functional connectivity, using resting-state fMRI. MCCB working memory tests largely agreed with Baddeley’s multicomponent model. However, the dlPFC was not directly associated with MCCB working memory. Instead, the central executive system was best represented by delocalized RECN and LECN. These results illuminate the complex neurobiology underlying working memory in terms of intrinsic network architecture of the brain at rest, in an unprecedented level of detail, encompassing every voxel. Lastly, with the increasing use of MCCB in pharmaceutical development, these results can help to provide insight into the neuronal pathways and connections mediating the potential therapeutic effects of novel therapeutic strategies designed to ameliorate the cognitive symptoms of schizophrenia.
1 : Schizophrenics show spatial working memory deficits. Arch Gen Psychiatry 1992; 49:975–982Crossref, Medline, Google Scholar
2 : Working memory impairments in schizophrenia: a meta-analysis. J Abnorm Psychol 2005; 114:599–611Crossref, Medline, Google Scholar
3 : Working memory deficit as a core neuropsychological dysfunction in schizophrenia. Am J Psychiatry 2003; 160:1809–1816Crossref, Medline, Google Scholar
4 : Evolution of neuropsychological dysfunction during the course of schizophrenia and bipolar disorder. Psychol Med 2011; 41:225–241Crossref, Medline, Google Scholar
5 : Neurocognitive deficits in adolescents with schizophrenia: longitudinal stability and predictive utility for short-term functional outcome. J Am Acad Child Adolesc Psychiatry 2007; 46:867–878Crossref, Medline, Google Scholar
6 : Cognitive and symptomatic predictors of functional disability in schizophrenia. Schizophr Res 2011; 126:257–264Crossref, Medline, Google Scholar
7 : Neurocognitive deficits and functional outcome in schizophrenia: are we measuring the “right stuff”? Schizophr Bull 2000; 26:119–136Crossref, Medline, Google Scholar
8 : The MATRICS Consensus Cognitive Battery, part 1: test selection, reliability, and validity. Am J Psychiatry 2008; 165:203–213Crossref, Medline, Google Scholar
9 Andrade J (ed): Working Memory in Perspective. Hove, United Kingdom, Psychology Press, 2001Google Scholar
10 : Neuroimaging studies of working memory: a meta-analysis. Cogn Affect Behav Neurosci 2003; 3:255–274Crossref, Medline, Google Scholar
11 : A meta-analysis of executive components of working memory. Cereb Cortex 2013; 23:264–282Crossref, Medline, Google Scholar
12 : N-back working memory paradigm: a meta-analysis of normative functional neuroimaging studies. Hum Brain Mapp 2005; 25:46–59Crossref, Medline, Google Scholar
13 : Beyond hypofrontality: a quantitative meta-analysis of functional neuroimaging studies of working memory in schizophrenia. Hum Brain Mapp 2005; 25:60–69Crossref, Medline, Google Scholar
14 :
15 : A functional magnetic resonance imaging study of working memory abnormalities in schizophrenia. Biol Psychiatry 2006; 60:11–21Crossref, Medline, Google Scholar
16 : Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb Cortex 2012; 22:158–165Crossref, Medline, Google Scholar
17 : Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci USA 2009; 106:13040–13045Crossref, Medline, Google Scholar
18 : The economy of brain network organization. Nat Rev Neurosci 2012; 13:336–349Crossref, Medline, Google Scholar
19 : Intrinsic and task-evoked network architectures of the human brain. Neuron 2014; 83:238–251Crossref, Medline, Google Scholar
20 : Measuring and testing dependence by correlation of distances. Ann Stat 2007; 35:2769–2794Crossref, Google Scholar
21 : Working memory: looking back and looking forward. Nat Rev Neurosci 2003; 4:829–839Crossref, Medline, Google Scholar
22 : Neural working memory; in Working Memory in Perspective. Edited by Andrade J. Hove, England, Psychology Press, 2001Google Scholar
23 : Frontal cortex as the central executive of working memory: time to revise our view. Cortex 2003; 39:871–895Crossref, Medline, Google Scholar
24 : The MATRICS Consensus Cognitive Battery, part 2: co-norming and standardization. Am J Psychiatry 2008; 165:214–220Crossref, Medline, Google Scholar
25 : Binding in visual working memory: the role of the episodic buffer. Neuropsychologia 2011; 49:1393–1400Crossref, Medline, Google Scholar
26 : Unified segmentation. Neuroimage 2005; 26:839–851Crossref, Medline, Google Scholar
27 : A method for making group inferences from functional MRI data using independent component analysis. Hum Brain Mapp 2001; 14:140–151Crossref, Medline, Google Scholar
28 : An information-maximization approach to blind separation and blind deconvolution. Neural Comput 1995; 7:1129–1159Crossref, Medline, Google Scholar
29 : Estimating the number of independent components for functional magnetic resonance imaging data. Hum Brain Mapp 2007; 28:1251–1266Crossref, Medline, Google Scholar
30 : Comparison of multi-subject ICA methods for analysis of fMRI data. Hum Brain Mapp 2011; 32:2075–2095Crossref, Medline, Google Scholar
31 : Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 2012; 59:2142–2154Crossref, Medline, Google Scholar
32 : Multiple comparison procedures for neuroimaging genomewide association studies. Biostatistics 2015; 16:17–30Crossref, Medline, Google Scholar
33 : Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp 2002; 15:1–25Crossref, Medline, Google Scholar
34 : An integrated software suite for surface-based analyses of cerebral cortex. J Am Med Inform Assoc 2001; 8:443–459Crossref, Medline, Google Scholar
35 : Network-based statistic: identifying differences in brain networks. Neuroimage 2010; 53:1197–1207Crossref, Medline, Google Scholar
36 : The central executive: a concept and some misconceptions. J Int Neuropsychol Soc 1998; 4:523–526Crossref, Medline, Google Scholar
37 : Working memory: a view from neuroimaging. Cognit Psychol 1997; 33:5–42Crossref, Medline, Google Scholar
38 : CNTRICS imaging biomarkers selection: working memory. Schizophr Bull 2012; 38:43–52Crossref, Medline, Google Scholar
39 : Estimating sample size in functional MRI (fMRI) neuroimaging studies: statistical power analyses. J Neurosci Methods 2002; 118:115–128Crossref, Medline, Google Scholar
40 : An empirical investigation into the number of subjects required for an event-related fMRI study. Neuroimage 2004; 22:879–885Crossref, Medline, Google Scholar
41 : Power calculation for group fMRI studies accounting for arbitrary design and temporal autocorrelation. Neuroimage 2008; 39:261–268Crossref, Medline, Google Scholar
42 : Executive control function: a review of its promise and challenges for clinical research. A report from the Committee on Research of the American Neuropsychiatric Association. J Neuropsychiatry Clin Neurosci 2002; 14:377–405Link, Google Scholar
43 : The elusive nature of executive functions: a review of our current understanding. Neuropsychol Rev 2007; 17:213–233Crossref, Medline, Google Scholar