|Home | About | Journals | Submit | Contact Us | Français|
Changes in the default mode network (DMN) have been linked to multiple neurological disorders including schizophrenia. The anticorrelated relationship the DMN shares with task-related networks permits the quantification of this network both during task (task-induced deactivations: TID) and during periods of passive mental activity (extended rest). However, the effects of different methodologies (TID vs. extended rest) for quantifying the DMN in the same clinical population are currently not well understood. Moreover, several different analytic techniques, including independent component analyses (ICA) and seed-based correlation analyses exist for examining functional connectivity during extended resting states. The current study compared both methodologies and analytic techniques in a group of patients with schizophrenia (SP) and matched healthy controls (HC). Results indicated that TID analyses, ICA, and seed-based correlation all consistently identified the midline (anterior and posterior cingulate gyrus) and lateral parietal cortex as core regions of the DMN, as well as more variable involvement of temporal lobe structures. In addition, SP exhibited increased deactivation during task, as well as decreased functional connectivity with frontal regions and increased connectivity with posterior and subcortical areas during periods of extended rest. The increased posterior and reduced anterior connectivity may partially explain some of the cognitive dysfunction and clinical symptoms that are frequently associated with schizophrenia.
The ability to define the neuronal basis of passive mental activity is of particular interest in disorders such as schizophrenia, in which the inability to differentiate between normal internal mentations and external stimuli (e.g., perceiving internally generated thoughts as externally generated auditory input) is a hallmark of the disorder (Linden 2008). Although early neuroimaging research focused on task-related differences in functional activation (Ragland, et al. 2007), more recent data suggests that patients with schizophrenia (SP) may also be characterized by abnormalities in resting state networks (Calhoun, et al. 2008a; Calhoun, et al. 2008b; Garrity, et al. 2007; Harrison, et al. 2007; Zhou, et al. 2007a; Zhou, et al. 2007b). Of all putative resting state networks (Damoiseaux, et al. 2006; De Luca, et al. 2006), the default mode network (DMN) has received the most attention due to its anti-correlated relationship with experimentally induced cognitive and sensory-motor states (Fox, et al. 2005).
In early neuroimaging studies, general linear model analyses were predominantly utilized to determine task-related activity (Binder, et al. 1999; Shulman, et al. 1997). These analyses consistently revealed a network of task-induced deactivations (TID) during the transition between baseline and task states, which was later formalized as the DMN (Raichle, et al. 2001). However, the magnitude of TID has been shown to vary both as a function of task difficulty and of the person’s ability to disengage from a task (Esposito, et al. 2005; Greicius and Menon 2004; McKiernan, et al. 2003), suggesting that differences in TID may be either obscured or accentuated based on task performance in clinical populations. More recently, studies have begun to examine the DMN during extended periods of passive mental activity without the potential confounds of a task. Several analytic techniques, such as independent component analyses (ICA: (Beckmann, et al. 2005; Calhoun, et al. 2008a; Jafri, et al. 2008; McKeown, et al. 2003; van de Ven, et al. 2004) and seed-based correlation techniques (Biswal, et al. 1995; De Luca, et al. 2006), have facilitated the analysis of extended resting state data. However, the strength of this methodological approach (i.e. the use of rest data) has also been viewed as a weakness due to the unstructured nature of the task and the possibility that clinical samples might engage in different internal mentations (Morcom and Fletcher 2007).
A recent review article (Buckner, et al. 2008) qualitatively compared the DMN network in healthy controls (HC) during TID and with seed-based (i.e., hippocampal seed) correlations. Results from both methods indicated a primary DMN sub-network, consisting of the ventral medial prefrontal cortex (vmPFC), the anterior cingulate cortex (ACC), the posterior cingulate cortex (PCC), the bilateral inferior parietal lobes (IPL) and the lateral temporal cortex, as well as a secondary DMN sub-network, consisting of the hippocampal formation (HF) and parahippocampal cortices (PHC). Network analyses suggested that while the ACC and PCC were the central hubs of this network, each hub displayed slightly different connectivity properties with other nodes, leading to speculation of different roles for the ACC and PCC in passive mental activity (Broyd, et al. 2009). In addition, a qualitative analysis indicated potentially greater volumes of functionally connected voxels in the IPL, prefrontal cortex, and HF in the seed-based correlation methods compared to TID, although these results were potentially biased by the placement of the seed within the hippocampus. To date, we are unaware of a study comparing these two main methodological approaches (i.e. TID vs. extended rest) for DMN identification or the two main analytic techniques (ICA vs. seed-based correlation) used to quantify the DMN during periods of extended rest in SP.
Studies examining DMN activity in SP using the TID methodology have produced mixed results. For example, one study found no significant differences in deactivation between SP and controls at lower cognitive loads during a working memory task (Pomarol-Clotet, et al. 2008). However, at higher cognitive loads (2-back vs. baseline), SP showed performance-related decreases in activation in frontal task-positive regions with concurrent decreases in deactivation of the medial frontal, ACC, anterior temporal, and right hippocampal regions. Another recent study, also utilizing a WM paradigm, showed similar decreases in task-related deactivation in the mPFC in SP and their first-degree relatives (Whitfield-Gabrieli, et al. 2009). However, these decreases in deactivation were present for both the 0- and 2-back conditions and were accompanied by increases in task-related activation of the right dorsolateral prefrontal cortex for SP compared to HC. These findings highlight the complex interaction which may exist between diagnosis and cognitive load in TID research. Other results suggest that while SP and controls have similar spatial deactivation maps during the performance of an auditory oddball task, SP took two seconds longer to return to baseline for the thalamus, PCC and precuneus (Kim, et al. 2009). In contrast to the above results, (Harrison, et al. 2007) reported an overall increase in midline deactivation for SP during a multi-source interference task in regions including the mPFC, ACC, and PCC. In addition, the magnitude of deactivation in the rostral ACC and medial prefrontal cortex was related to task performance in SP (i.e. greater deactivation indicated improved performance and increased emotional awareness) whereas controls showed a task-related relationship with the PCC.
The results from studies examining connectivity during periods of extended rest in SP have also yielded variable results (Broyd, et al. 2009). Some findings provide evidence for global decreases in functional connectivity in SP compared to controls (Liang, et al. 2006), the extent of which (Liu, et al. 2008) found to be positively correlated with illness duration. Similarly, (Bluhm, et al. 2007) found that SP had lower positive correlations between the PCC seed region and the ACC, the middle temporal gyrus and the cerebellum compared to controls, though no differences in the negative DMN correlations were found. This same study also identified a significant relationship between positive symptoms and PCC connectivity with speech-related areas in the temporal gyrus, pointing toward the possibility that aberrant DMN connectivity with these temporal regions may play a role in positive symptom severity in SP. Finally, other research suggests higher and more variable frequency fluctuations within the DMN for SP, which in turn may explain the decreased functional connectivity reported for this group (Calhoun, et al. 2008a; Garrity, et al. 2007).
In contrast to the above results, other research has found evidence of increased connectivity within the DMN for SP, specifically between the PCC, medial prefrontal cortex, and inferior temporal gyrus (Zhou, et al. 2007b). In a study comparing the temporal relationships between several functional networks, (Jafri, et al. 2008)found evidence of increased connectivity between the DMN and other task-positive resting-state networks in SP, implying that SP have increased competition and anti-correlation (negative correlations) between these networks.
The primary aim of the current study was to contrast methods for identifying the DMN (TID versus extended rest) in a population of SP and matched controls. Our second aim was to compare DMN activation during extended resting state using the two analytic techniques (ICA and seed-based correlational techniques) most commonly used in the current literature. Although both ICA and seed-based correlation analyses identify the DMN through temporal correlations, the concepts behind each are fundamentally distinct. The most central difference is that seed-based techniques are sensitive only to temporal correlations with the chosen seed and are thus inherently biased to that particular region or voxel. In contrast, ICA is a data-driven approach that identifies and separates temporally coherent networks based on all brain voxels. Our three primary hypotheses were 1) SP would exhibit a reduction in the magnitude of DMN deactivation relative to controls, suggesting a failure to “turn off” the DMN during periods of cognitive demand; 2) SP would demonstrate decreased connectivity during periods of extended rest (e.g., decreased magnitude of temporal correlations) between the PCC and ACC; and 3) that while both ICA and seed-based correlation analyses would be able to reliably identify the DMN, ICA would be more sensitive to group differences in the network’s functional (i.e. temporal) relationships.
All SP subjects were diagnosed by an experienced clinician (J.C.) or team member with the Structured Clinical Interview for DSM-IV Axis-I Disorders, Clinician Version (SCID-CV). Sixteen SP (15 male, 1 female) and 16 HC (15 male, 1 female) participated in the current experiment. One HC was identified as an outlier (three standard deviations above the mean) on several behavioral measures and was therefore excluded from further analyses along with the matched SP. SP with a history of other neurological disease, history of psychiatric hospitalizations within the previous six months or history of substance abuse within the past year were excluded from the study. SP were also required to be stable on an atypical, anti-psychotic medication (Aripiprazole 4; Ziprasidone: 1; Risperidone: 5; Quetiapine Fumarate: 4; Olanzapine: 2) for at least three months to be included in the current study. All HC were screened and excluded from the study based on a history of major medical conditions, neurological disease, major psychiatric disturbance, substance abuse or psychoactive prescriptive medications usage. There were no significant differences (p > 0.10) between SP and the remaining HC for all major demographic categories including age (SP: 40.2 +/− 8.2, HC: 40.1 +/− 8.8), education (SP: 12.6 +/− 2.4, HC: 13.0 +/− 1.4), or handedness (SP: 77.7 +/− 56.1, HC: 67.3 +/− 68.9) as assessed by the Edinburgh Handedness Inventory (Oldfield 1971). Informed consent was obtained from subjects according to institutional guidelines at the University of New Mexico and the New Mexico Department of Veterans Affairs.
Participants were scanned both at rest and while performing a task. The first data set consisted of three periods of extended rest (3 minutes each) in which participants simply maintained visual fixation on a white cross centered on a black background (visual angle = 1.54°). The second data set was collected as subjects performed a complex selective attention task (see (Mayer, et al. 2009) for complete details). Briefly, participants were instructed to simultaneously tap fingers from both hands to auditory (1000 Hz tone) and/or visual (flashing checkerboard; visual angle = 19.42° × 14.88°) stimuli presented for a period of 8 seconds. Each task was preceded by a cue signifying whether participants should attend and tap their fingers to the auditory stimuli while ignoring visual input (attend-auditory), attend and tap to visual stimuli while ignoring auditory input (attend-visual) or both (multimodal). In all conditions, auditory and visual stimuli were presented at .5, 1 or 2 Hz and could occur either in synchrony (multimodal condition) or out of phase (attend-visual and attend-auditory conditions). Between each block of the complex attention task, participants were asked to maintain visual fixation on a white cross (identical to that presented during extended rest, visual angle = 1.54°). The length of these baseline periods was pseudorandomly alternated between 10 and 14 seconds. For all subjects, the extended rest task was presented first followed by 2 runs of the complex selective attention task. This sequence was repeated 3 times for a total of 3 extended rest periods and 6 TID runs.
To reduce head motion during the scan, participants were required to practice the selective attention task twice outside of the scanning session. During the first practice session, participants performed the task until they were able to achieve basic task competency (approximately 5 minutes of practice). Participants then repeated the practice session with instructions to maintain their head in a fixed position while responding to task stimuli with both hands. There were no significant differences between SP and HC on any of the estimated motion parameter variables.
At the beginning of the scanning session, high resolution MPRAGE T1 [TE (echo time) = 4.76 ms, TR (repetition time) = 12 ms, 20° flip angle, number of excitations (NEX) = 1, slice thickness = 1.5 mm, FOV (field of view) = 256 mm, matrix size = 256 × 256] anatomic images were collected on a 1.5 Tesla Siemens Sonata scanner. For both the extended resting state and complex attention task a series of echo-planar images were collected using a single-shot, gradient-echo echoplanar pulse sequence [TR = 2000 ms; TE = 36 ms; flip angle = 90°; FOV = 256 mm; matrix size = 64 × 64]. A total of 28 sagittal 5-mm thick slices were selected to provide whole-brain coverage (voxel size: 4 × 4 × 5 mm). Two hundred and three echoplanar images were collected during each complex attention run and 92 echoplanar images were collected during each extended resting state run. For both tasks, the first three images of each run were eliminated to account for T1 equilibrium effects, leaving a total of 1200 images for the complex attention task and 267 images for the extended resting state data.
Functional images were generated and processed using a mixture of freeware and commercial packages including the Analysis of Functional NeuroImages (AFNI), GIFT, MATLAB (Mathworks, Inc., Sherborn, MA, USA) and FSL. For both tasks, data was spatially registered in both 2- and 3-dimensional space to minimize effects of head motion and temporally interpolated to correct for slice-time acquisition differences (AFNI). Data was then de-spiked, spatially blurred using an 8 mm3 Gaussian full-width half-maximum filter kernel, and normalized to a standard stereotaxic coordinate space (performed in FSL) and resampled to 3 mm isotropic voxels (Talairach and Tournoux 1988).
A stereotaxically-derived mask of DMN, based on the most commonly reported regions in the literature (Broyd, et al. 2009; Buckner, et al. 2008), was then constructed using the Wake Forest University (WFU) Pick atlas. The specific regions included in the mask were the bilateral PCC (BAs 23/31), inferior and superior parietal lobes (BAs 7/39/40), superior frontal gyrus (BAs 8/9/10), ACC (BAs 11/32), and parahippocampal gyrus. The ACC and PCC regions were weighted with a value of 3 given the central role of these regions in the DMN (Raichle, et al. 2001) and to improve correspondence with functional activity (Franco, et al. 2009). The resultant mask was blurred with an 8 mm3 Gaussian full-width half-maximum filter to match the input FMRI data.
The general steps used across methodologies and analytic techniques were extremely similar following the individual pre-processing steps described below. Voxel-wise t-tests were first conducted on unthresholded data to examine potential differences in the DMN between HC and SP (HC vs. SP). For all group comparisons, the resulting statistical parametric maps were corrected at p < .05 by applying a parametric threshold of p < .005 in conjunction with a minimal cluster size of 0.480 milliliters (Forman, et al. 1995). These two thresholds were derived from 10,000 Monte Carlo simulations that demonstrated the probability of obtaining a significant activation cluster for an entire volume (Type I error) was below p < 0.05 for both the task and the extended resting state data. In addition, the unthresholded data were also collapsed across both groups to obtain an overall DMN mask (again corrected at p < .05), which was then used for group-wise ROI analyses for each methodology and analytical technique. Parameters of interest in the group ROI analyses included the total number of activated voxels, the spatial goodness of fit (e.g., between individual DMN and WFU mask), average magnitude of deactivation (TID analysis only) and the magnitude of temporal correlation (extended resting state data only). The overall masks were also used to identify common DMN regions across methodologies and analytical techniques.
A deconvolution analysis was used to generate one hemodynamic response function (HRF) for each selective attention condition on a voxel-wise basis (AFNI). Each HRF was derived from the first 12 images (22 seconds) following the onset of the cue. The images corresponding to the peak of the HRF (eight to fourteen seconds post-stimulus onset) were summed and divided by the baseline to obtain an estimate of percent signal change. Voxel-wise t-tests were then conducted to determine differences in the magnitude of deactivation (i.e., percent signal change) across the two groups. Similar to the group-wise comparisons, each subject’s data was parametrically thresholded at p < .005 with a minimal cluster size of 0.480 ml. Next, the unthresholded percent signal change estimates (all selective attention conditions) were pooled across all subjects (SP + HC) and tested against the null hypothesis to produce the overall DMN TID mask. This mask was then used in conjunction with each subject’s thresholded data to determine the number of active voxels and the average magnitude of activation/deactivation in activated voxels on a per-subject basis.
The GIFT software package was used to perform a group ICA for all subjects (SP + HC) as a first step for identifying the DMN during extended periods of rest (Calhoun, et al. 2001). A modified minimum description length (MDL) criterion was used to establish the ideal number of independent components for the group ICA, which was estimated to be sixteen (Li, et al. 2007). The resulting 16 components were correlated with the DMN mask derived from the WFU atlas to objectively determine which group component exhibited the highest spatial correlation. Back-reconstruction was used to select each individual subject’s component that corresponded with the group DMN component. Components were then scaled to an empirically derived z-score by dividing by the standard deviation of the original time sequence. In the ICA framework, the empirical z-scores approximate the temporal correlation (magnitude of the z-score indicates strength of the linear relationship) between each voxel and the associated component and can therefore be either positive or negative. Voxel-wise t-tests were then conducted to investigate potential differences in the magnitude of correlation across the two groups (corrected at p < .05). In addition, the spatial correlation between each subject’s DMN component and the WFU mask was also determined. The data from each individual subject component image were then thresholded so that only the top 5% of voxels remained.
An overall ICA DMN mask was created by comparing all subjects’ components (collapsed across group) against a null hypothesis of zero and selecting significant voxels at a corrected alpha of p < .05. Similar to the TID analysis, the overall group DMN mask was multiplied by the thresholded individual components to determine the total number of active voxels within the mask and the average z-score of these voxels.
For the seed-based analysis, a regression analysis was first conducted on individual subjects’ time-series data to remove potential sources of physiological and machine noise (Fox, et al. 2005). First, individual anatomical images (i.e., T1) were segmented into maps of white matter, gray matter and cerebral spinal fluid (CSF) using FSL’s FAST algorithm (Zhang, et al. 2001). Second, the resultant CSF and white matter masks were used to obtain an average time-series for these tissues on an individual basis. All six movement parameters, the ROI based time-series for CSF and white matter, a constant term, and a linear term were then entered into a linear regression against the extended resting state time-series (Fox, et al. 2005). Gray matter was not utilized as a global signal regressor, given existing evidence that this may introduce spurious anti-correlations into the data (Murphy, et al. 2009). The resulting residualized time-series were then concatenated and used in the connectivity analyses.
Spheres with a fourteen millimetre diameter were created in the ACC (center of sphere: x = 0, y = 37, z = 20; Talairach and Tournoux atlas coordinates (Talairach and Tournoux 1988) and PCC (x = 0, y = −49, z = 24) to form the ROIs necessary for a seed-based analysis. The centroids for these spheres were based on spatial locations of anterior and posterior cingulate activity in the TID analysis (setting X = 0). Next, the average time-courses within the ACC and PCC were computed by averaging the time-series data for all voxels in the ROI spheres. Subsequently, each of these two averaged time-courses was then correlated with all brain voxels. The resultant r values were then transformed into Fisher’s z values. In the seed-based correlation framework, the individual Fisher’s z value approximates the temporal correlation between each voxel and the associated seed and therefore can be either positive or negative. Voxel-wise t-tests were then conducted to determine group-wise differences in the magnitude of the Fisher’s z values. In addition, the spatial correspondence between each subject’s correlational map and the WFU mask was also determined. Individual subject correlation maps were thresholded using similar parametric and spatial thresholds as the group-wise comparisons.
An overall DMN mask for the seed-based analyses was created by contrasting all subjects’ maps against zero and correcting for false positives, which was then used to determine ROI parameters identical to those used in the ICA analyses.
Two 2 × 3 × 3 [Group (SP vs. HNV) × Attention Condition (multimodal, attend-auditory, attend-visual) × Frequency (.5, 1, 2Hz)] mixed-measures ANOVAs were conducted with the coefficient of variation (COV) and error data as the dependent variables for the task data. In the context of the current task, COV is a measure of variability in performance that is most likely attributable to increased distractibility (Mayer, et al. 2009). Only the results pertaining to group differences will be reported for the purpose of the current study. The COV analysis showed a significant main effect of group (F1,28=6.45, p <.05) with SP (.18 +/−.03) displaying more variable response times across all conditions compared to HC (.15 +/−.02). Neither the group by condition nor group by frequency interactions was significant. In the error analysis, the main effect of group and the interactions were not significant (p > .10).
To reduce redundancy in the presentation of results, we first briefly describe the regions (Figure 1; Table I) that were shown to have significant overlap across both groups during the TID and the two extended resting state data analyses (ICA and seed-based correlations). Commonly identified DMN regions included the medial and superior frontal gyrus (BAs 8/9/10/11/6), ACC (BAs 24/32), PCC (BAs 29/30/23) extending to the precuneus (BA 31), and bilateral middle and superior temporal gyrus (BAs 39/19) continuing into the angular gyrus. In addition, all three methods also identified areas of significant overlap within the left inferior frontal gyrus (BA 47), the left medial temporal lobe (BAs 34/35) and the right inferior and middle temporal gyrus (BAs 20/21).
Both the seed-based and ICA analyses identified significant DMN connectivity within right inferior frontal gyrus (BA 47), right medial temporal lobe (BAs 34/35), bilateral middle/superior temporal gyrus (BA 21/22), thalamus, striatum and cerebellum. The seed-based and TID analyses identified the left inferior/middle temporal gyri (BAs 21). Functionally connected regions specific only to the ICA analyses included the right pre- and postcentral gyri (BAs 1/2/3/4), right fusiform gyrus (BA 37), and extensive areas within bilateral occipital lobes (BA 18) and cerebellum. Finally, seed-based correlation analyses uniquely identified DMN connectivity in the bilateral cerebellar tonsil.
Results from the independent samples t-test indicated that there were no significant voxel-wise differences in DMN activity across the two groups. Two ANOVAs were then conducted to determine if either the overall magnitude of deactivation or the number of deactivated voxels within the TID network (identified with pooled-group mask) differed as a function of diagnosis or rate of stimulation. Results from the first ANOVA indicated that the main effect of group was significant (F 1, 27 = 4.3, p < .05), with the magnitude of deactivation being significantly greater in the DMN network for SP (mean = −0.329 +/− 0.087) compared to HC (mean = −0.270 +/− 0.064). In spite of the apparent qualitative differences in the greater volume (e.g., number of voxels) of deactivation for SP (see Figure 2), the effect of group was not statistically significant (p > .10) for the ANOVA. This was likely a result of the large degree of variability in this measure for both the SP (mean = 2179.9 +/− 1068.3) and HC (mean = 1612.4 +/− 1031.2), which resulted in a moderate effect size (Cohen’s d = .54) but no significant group differences.
For the extended resting state analyses, group ICA was first performed on the collapsed group data followed by a spatial correlation with the WFU mask to objectively determine the component that most closely corresponded to the DMN (Franco, et al. 2009; Greicius and Menon 2004). Results indicated two components that were highly correlated with the WFU mask, and primarily differed by the relative contribution of either anterior (r = .52) or posterior (r = .49) clusters (mean correlation for remaining components = .13 +/− .09). The anterior component was comprised of large clusters of bilateral connectivity within the ACC, medial and superior frontal cortex, inferior frontal gyrus, and bilateral PCC extending into the para-central lobule. Additional regions of connectivity included small clusters within bilateral inferior/superior occipital gyrus and several clusters in the cerebellum. In contrast, for the posterior component the volume of connectivity was relatively small within the midline and lateral frontal structures. Instead, large clusters of connectivity were present within the bilateral PCC, precuneus, medial temporal lobe and ventral visual stream. Common areas of connectivity across both components included the middle/superior temporal gyrus, striatum, thalamus, and cerebellum.
Following component identification, voxel-wise t-tests were then conducted on the individual subject component z-score to investigate potential group differences for both the anterior and posterior component (Table II; Figure 3). Results are classified according to whether the regions were positively or negatively related to the component (direction of the z-score) and whether the magnitude of the z-score was greater for the SP or HC. Finally, several regions exhibited a variable pattern in which the data from one group was positively correlated while the data from the other group was negative. For the anterior component, results indicated greater positive z-scores for the SP group in the left middle temporal gyrus, angular gyrus and bilateral anterior nucleus of the thalamus. A more negative z-score for the SP group was seen in the left insula (BA 13) extending to the superior temporal gyrus, left middle temporal gyrus (BA 21) and the left cerebellar tonsil and culmen. There were no areas in which HC exhibited a significantly stronger relationship with the anterior component.
The results of the t-tests on the more posterior ICA component indicated that SP had a greater positive z-score than HC in the right precuneus and angular gyrus (BAs 19/39). A greater negative z-score was observed for HC in the right inferior frontal gyrus (BA 46) compared to SP. In contrast, SP showed a greater negative z-score in the right (BAs 21/38) and left (BA 22) middle and superior temporal gyri. Though SP and HC showed a similar magnitude of relationship between the PCC component and the right middle occipital/inferior temporal gyri (BA 19), the relationship was positive for HC and negative for SP. Alternately, SP exhibited a negative, and HC a positive, relationship with the left middle frontal gyrus (BA 11). There were three areas in which HC exhibited strong relationships with the posterior component while SP showed minimal to no relationship. These areas included the left middle frontal gyrus (BA 10) and the left medial frontal and anterior cingulate gyri (BA 6) where HC had higher negative z-scores; and the left PCC/precuneus (BA 31) where HC had higher positive z-scores. Finally, though HC exhibited almost no relationship between the PCC and the left middle occipital and inferior temporal gyri, SP shared a moderate negative relationship with those areas.
ACC and PCC seed-based correlational analyses were also conducted on the extended resting state data. Similar to previous reports (Uddin, et al. 2008), the network of regions that exhibited a positive correlation with the ACC and PCC seeds were very similar and closely corresponded to the common DMN network reported above. Likewise, this technique also resulted in a large network of anti-correlated (e.g., negative correlations) regions in the frontal, temporal, parietal and occipital lobes for both the ACC and PCC clusters, which have been previously discussed in other publications (Fox, et al. 2005; Uddin, et al. 2008).
Two voxel-wise t-tests examined differences in connectivity between SP and HC using either the ACC or the PCC as the seed region (Table III; Figure 4). SP exhibited a greater correlation between the ACC and the left caudate and bilateral thalamus. In contrast, HC exhibited a greater correlation between the ACC and the left middle/superior temporal gyri (BA 39). Similar to the ICA analyses there were regions where the groups exhibited opposing directional relationships. Specifically, HC showed a positive correlation between the ACC seed and the left insula and claustrum (BA 13), while SP showed a negative correlation between these regions.
The PCC connectivity analysis indicated significantly greater negative correlations for SP compared to HC within the right precentral and medial frontal gyrus (BA 6), the left posterior parietal lobules (BA 40) extending to the postcentral gyrus (BA 2/5), and the bilateral auditory cortical areas (superior/transverse temporal gyri, BA 41) and insula (BA 13). There were also regions where the groups exhibited opposing directional relationships. SP showed positive correlations between the PCC seed and the midline cerebellar structures and also subcortical regions (bilateral caudate and thalamus) while HC showed negative correlations with these same regions.
Three MANOVAs were then conducted to examine if the number of active voxels, temporal correlation, and goodness of fit (the three dependent measures derived from both ICA and seed-based correlation analyses) varied as a function of group across the entire DMN network. The MANOVA examining the number of significantly activated voxels indicated a significant multivariate effect of group (F 4,25 = 4.0, p < .05). Univariate tests indicated that SP (mean = 1847.5 +/− 160.1) had significantly more active voxels than HC (mean = 1629.4 +/− 172.9) in the posterior (F 1,28 = 12.9, p < .005), but not the anterior DMN ICA component (p > .10). Univariate tests for either of the seed correlation analyses were not significant (p > .10), and the effect sizes were small.
A similar pattern of results was obtained for the MANOVA examining the spatial goodness of fit (e.g., correlation) between each of the components and the DMN template constructed with the WFU mask. Specifically, the multivariate effect for group (F 4,25 = 3.7, p < .05) and the univariate tests for the ACC ICA component (F 1,28 = 4.7, p < .05) and PCC ICA component (F 1,28 = 4.6, p < .05) were all significant. Once again, the direction of effects for these two structures were reversed, with HC (mean = .50 +/− .04) exhibiting a better overall fit with the DMN template for the anterior component compared to SP (SP mean = .47 +/−.05), and SP (mean = .44 +/−.05) exhibiting a better overall fit with the posterior component compared to HC (mean = .39 +/−.06). Univariate tests for either of the seed correlation analyses were not significant (p > .10).
Finally, the multivariate effect of group was not significant for the MANOVA examining the magnitude of temporal correlation between groups (p > .10). However, there was a non-significant trend for both the ACC ICA component (F 1,28 = 3.7, p = .07) and the seed-based ACC component (F 1,28 = 2.9, p = .10) for univariate effects. In both methodologies, the magnitude of the temporal correlation was higher for HC compared to SP in the anterior network.
There is a growing evidence indicating that the DMN may be typified by more than one functional network (Buckner, et al. 2008; Calhoun, et al. 2008a; Damoiseaux, et al. 2006), including a more anterior and posterior network (Broyd, et al. 2009; Buckner, et al. 2008). As other studies have intimated (Broyd, et al. 2009; Foucher, et al. 2005; Garrity, et al. 2007; Uddin, et al. 2008) deficiencies in the DMN and its modulation between rest and task states may be the result of increased temporal coherence between some nodes of the DMN (the PCC and temporal, occipital and thalamic regions) and equivalent or decreased connectivity between other nodes (ACC, PCC, and frontal regions). The primary objective of this paper was to compare common methods and analytic techniques used to identify the DMN during both task and extended resting states, and to examine how the results from these methodologies might differ in a clinical sample. In general, current results indicated increased functional connectivity within more posterior (i.e., non-frontal) cortical and sub-cortical regions for SP, whereas HC tended to exhibit either similar (no statistically significant differences) or greater connectivity within frontal cortical regions. In addition, TID was also greater for SP compared to HC. Finally, the results from our pooled analyses (HC+SP) also revealed several important findings, which are briefly discussed first.
The first finding from the pooled analyses pertained to the reliable deactivation of the ACC, PCC/precuneus, and bilateral superior temporal lobe/IPL in the general linear tests (TID), and the significant functional connectivity displayed between these same regions in both the seed-based correlation analyses (extended rest) and ICA (extended rest). Although the spatial extent (i.e., volume) of connectivity within these structures was greater for the ICA and seed-based analysis than for TID, the consistency of findings across various methodologies and analytic techniques confirms that these structures form the key nodes of the DMN (Buckner, et al. 2008). In contrast, the identification and volume of DMN involvement within the temporal lobe and subcortical structures was much more variable (see Figure 1). Specifically, while all methods and analytic techniques identified the left medial temporal lobe (HF and PHC), only the extended resting state techniques (seed-based analysis and ICA) identified significant correlations within the right medial temporal lobe and bilateral middle temporal gyri. A finding of more variable involvement of temporal lobe structures in the DMN is consistent with previous reports (Buckner, et al. 2008). This is notable given that previous studies (cited in (Buckner, et al. 2008)) utilized hippocampal seeds, which would be expected to correlate more with medial temporal and temporal regions than the midline (ACC and PCC) seeds utilized in the current study. Collectively, current and previous results suggest that the medial temporal structures may form a sub-network of the DMN (Buckner, et al. 2008) that is more readily observed using ICA and seed-based correlation techniques than through TID. Deactivation of temporal structures during typical functional imaging studies may be dependent both on the type of task and the contrasting baseline state (Stark and Squire 2001). A similar caveat may also apply to the striatum and thalamus, which both showed substantial activations in the seed-based analyses and ICA, but were noticeably absent in the TID analyses.
In addition, our pooled ICA analyses revealed an anterior and a posterior DMN component that were both highly correlated with the template mask. The finding of a split DMN component is not unique to the current study (Calhoun, et al. 2008a; Damoiseaux, et al. 2006; Sambataro, et al. 2008; van de Ven, et al. 2004), and is of interest given the different putative roles that the ACC and PCC may play during passive mental activity (Broyd, et al. 2009; Uddin, et al. 2008). Previous work in healthy controls suggests that the PCC (e.g., central hub of the posterior DMN), is actively involved in the simultaneous monitoring of endogenous and exogenous states including performance error (Eichele, et al. 2008) and attentional drifting (Weissman, et al. 2006). In contrast, the anterior hub of the DMN (ACC) has been linked with appraisal of our own and others’ psychological states (Broyd, et al. 2009) as well as levels of emotional awareness and task performance in SP (Harrison, et al. 2007). Finally, the PCC has been shown to have greater connectivity with medial temporal structures (HF and PHC) compared to the ACC (Buckner, et al. 2008).
Different patterns of anterior-posterior DMN connectivity were also observed between SP and HC in the extended resting state ROI analyses. Specifically, across both analytic techniques (seed-based and ICA analyses) and experimental variables, controls exhibited increased or equal DMN connectivity in frontal cortical regions (HC ≥ SP) compared to increased posterior cortical and subcortical connectivity for SP (SP ≥ HC). The direction of these findings was consistent both in terms of the extent of connectivity (number of voxels) and spatial correlation, and showed a similar directional trend in terms of temporal correlation. ROI analyses of the anterior ICA component indicated that HC exhibited a higher spatial correlation than SP. Likewise, HC exhibited a trend for greater temporal connectivity within the anterior DMN for both the ICA and seed-based techniques. In contrast, ROI analyses for the posterior ICA component indicated that SP exhibited a greater number of active voxels and a higher spatial correlation than HC.
Voxel-wise comparisons of the extended resting-state data revealed a similar anterior-posterior distinction as ROI analyses. Consistent with previous studies (Bluhm, et al. 2007; Broyd, et al. 2009; Liang, et al. 2006), HC exhibited increased connectivity (larger z-scores) compared to SP within several regions of the medial and lateral frontal cortex for the posterior ICA analyses. In contrast, increased connectivity for SP compared to HC was observed in several posterior (e.g., non-frontal) cortical and sub-cortical regions during both extended resting state analyses. ICA results indicated greater connectivity for SP within posterior regions such as the middle temporal gyrus, thalamus, superior temporal gyrus and cerebellum for the anterior component. Likewise, SP had higher positive z-scores than HC for the posterior component in the right angular gyrus and precuneus and, interestingly, showed a much stronger negative relationship with bilateral middle/superior temporal and left middle occipital gyri, regions where HC exhibited almost no relationship. The PCC seed-based analyses resulted in greater anti-correlations for SP within the left posterior parietal cortex and bilateral auditory cortex.
The combination of reduced frontal connectivity and heightened posterior connectivity may partially explain why SP are both more prone to cognitive deficits and more likely to experience unusual sensory events. As suggested by current and previous work (Bluhm, et al. 2007), increased connectivity (positive or negative) between the posterior node of the DMN (centered around the PCC) and regions typically associated with auditory input may play a role in SP difficulty distinguishing between normal internal mentations and external stimuli. In contrast, reduced connectivity between frontal regions may partially explain the deficits in attentional, working memory and executive functioning that have come to characterize SP and their first-degree relatives (Bowie and Harvey 2005; Filbey, et al. 2008; Luck and Gold 2008), and were present during the current study (i.e., increased variability in performance during selective attention task). In addition, reduced frontal connectivity may play a role in SP inability to modulate the DMN efficiently during task (Ragland, et al. 2007) as was suggested by both our TID analyses.
Contrary to our original hypothesis, the magnitude of TID (i.e. percent signal change) was significantly greater for SP, suggesting that patients disengaged the DMN more than HC. Although a quantitative examination of the number of deactivated DMN voxels during task did not reach conventional levels of statistical significance, the overall volume of deactivation was larger for SP compared to HC while task-related activation appeared to be more widespread for HC (see Figure 2). However, impairments in performance (i.e. increased variability of response time) were observed in SP, making it difficult to determine whether disparities in TID were attributable to performance or whether performance was adversely affected by abnormal TID. As noted in the introduction, previous studies have reported increased TID, decreased TID or a mixture of both in different neuronal structures for SP compared to HC (Garrity, et al. 2007; Harrison, et al. 2007; Pomarol-Clotet, et al. 2008; Whitfield-Gabrieli, et al. 2009). Although these mixed findings render a simplistic explanation difficult, current and previous results suggest that SP may exhibit a heightened DMN state which proactively interferes with task performance (Garrity, et al. 2007; Pomarol-Clotet, et al. 2008; Whitfield-Gabrieli, et al. 2009). However, this relationship may be dependent on the cognitive load (McKiernan, et al. 2003; Pomarol-Clotet, et al. 2008) and type of task (Gusnard, et al. 2001; Tomasi, et al. 2006).
The primary limitation of the current study was the sample size. Several of the findings in the ROI analyses displayed moderate effect sizes, suggesting that these results would have likely reached conventional levels of statistical significance with a larger sample size. Additionally, the direction of the effect sizes were consistent with the observed pattern of increased anterior connectivity for HC and reduced posterior connectivity for SP, suggesting that increasing the sample size would only increase the robustness of current findings. Other limitations include the possibility that the current results were influenced by systematic differences in morphological brain structure often observed in SP (Honea, et al. 2005), or that these results may be secondarily related to the effects of anti-psychotic medications.
Finally, although several consistencies were found across all three DMN measures (similar spatial maps, increased posterior and subcortical connectivity for SP, and increased frontal connectivity for HC) some differences in the results were also present. Specifically, the greatest discrepancy between the two extended rest analytic techniques was the more variable location of voxel-wise differences in functional connectivity between the two groups, which have several potential explanations. Seed-based correlation involves choosing a particular region of interest and then correlating the remaining brain voxels with the averaged timecourse of that particular region (or subsection of that region). Thus, this method is biased not only to the region chosen (e.g. ACC, PCC or hippocampus) but to the location of the seed within that region (e.g. rostral versus caudal ACC). Thus the placement of the seed may partially explain why some of the subtle group differences in functional connectivity did not converge across the two analytic techniques in identical spatial locations. In contrast, ICA provides a measure of the strength of temporal relationships between voxels that is based on commonalities in signal properties across the entire brain, which allows this technique to identify and separate multiple independent networks simultaneously. The criterion of independence inherent to ICA may be ideal for more global (e.g. multi-network) analyses, but may not provide any advantage for the examination of region-specific hypotheses that are not as concerned with the autonomy of a given network.
In terms of method, both extended resting state analyses identified important group differences on the voxel-wise level which were not evident in the TID analyses. However, it is likely that information derived from TID compared to extended resting states may be fundamentally different and potentially complementary. TID analyses may be better suited for measuring the magnitude of deactivation in functional networks whereas extended rest analyses may excel at identifying intrinsically correlated (positive and negative) regions given that one does not need to account for the presence of a task. Further usage of TID and extended resting states conditions in tandem will provide important information about how abnormalities in baseline networks, within a specific sample, translate into functional and behavioral impairments when external cognitive demands are made.
This research was supported by grants from The Mind Research Network DOE Grant No. DE-FG02-99ER62764, from the National Institutes of Health Grant 1 R01 EB006841 and from NCRR Center on Neural Mechanisms of Schizophrenia; COBRE Grant No. P20 RR021938. Special thanks to Diana South for assistance with data collection.