Abstract
Human neuroimaging has revealed a specific network of brain regions—the default-mode network (DMN)—that reduces its activity during goal-directed behavior. So far, evidence for a similar network in monkeys is mainly indirect, since, except for one positron emission tomography study, it is all based on functional connectivity analysis rather than activity increases during passive task states. Here, we tested whether a consistent DMN exists in monkeys using its defining property. We performed a meta-analysis of functional magnetic resonance imaging data collected in 10 awake monkeys to reveal areas in which activity consistently decreases when task demands shift from passive tasks to externally oriented processing. We observed task-related spatially specific deactivations across 15 experiments, implying in the monkey a functional equivalent of the human DMN. We revealed by resting-state connectivity that prefrontal and medial parietal regions, including areas 9/46d and 31, respectively, constitute the DMN core, being functionally connected to all other DMN areas. We also detected two distinct subsystems composed of DMN areas with stronger functional connections between each other. These clusters included areas 24/32, 8b, and TPOC and areas 23, v23, and PGm, respectively. Such a pattern of functional connectivity largely fits, but is not completely consistent with anatomical tract tracing data in monkeys. Also, analysis of afferent and efferent connections between DMN areas suggests a multisynaptic network structure. Like humans, monkeys increase activity during passive epochs in heteromodal and limbic association regions, suggesting that they also default to internal modes of processing when not actively interacting with the environment.
Introduction
Human neuroimaging studies revealed robust increases in cortical activity during passive tasks in contrast to epochs when individuals are engaged in externally oriented decisions (Shulman et al., 1997; Mazoyer et al., 2001; Laird et al., 2009). The ensemble of brain regions active during passive tasks, called the default-mode network (DMN), consistently includes the anterior and posterior cingulate cortex, medial and lateral parietal cortex, and medial prefrontal cortex (Raichle et al., 2001; Buckner et al., 2008). These regions form nexuses of distributed connectivity throughout heteromodal cortex (Greicius et al., 2003; Buckner et al., 2009) and possess metabolic properties that set them apart from other regions (Raichle et al., 2001; Vaishnavi et al., 2010). In recent years, the DMN has attracted considerable interest, as well as controversy, mainly in relation to its functions.
To test specific hypotheses about the functional role of the DMN, it is important to clarify whether default mode activity is present in other primate species. Also, the definition of the DMN in monkeys would permit the investigation of its single-cell properties and the use of invasive interventional protocols to probe its functional role—which is less feasible in humans (Tsao et al., 2006; Ekstrom et al., 2008; Moeller et al., 2008; Caggiano et al., 2009). Despite the absence of linguistic abilities and presumably certain forms of self-referential processing (Passingham, 2009), many protoforms of human skills can be found in monkeys, thus allowing for comparative monkey–human studies for a wide range of tasks (Nakahara et al., 2002; Vanduffel et al., 2002; Orban et al., 2004; Tsao et al., 2008). Cytoarchitectonic maps in the vicinity of the human DMN have revealed anatomical similarities between humans and monkeys (Parvizi et al., 2006; Vogt et al., 2006; Buckner et al., 2008; Binder et al., 2009), supporting the notion that elements of the DMN may be conserved across primate species, although analogies for regions including the inferior parietal lobule are less certain (Caminiti et al., 2010). Notably, activity reductions under specific tasks were recently documented in monkeys by single-cell recordings in the posterior cingulate (Hayden et al., 2009) and by positron emission tomography (PET) in medial prefrontal and parietal areas (Kojima et al., 2009). Furthermore, functional connectivity investigations in anesthetized monkeys provide evidence for a network of regions in the monkey that may be a functional equivalent of the human DMN (Vincent et al., 2007, 2010; Margulies et al., 2009).
However, the defining functional property of the DMN—that it increases activity in passive states relative to active periods of engagement with the environment—has not been explored in a manner parallel to the work in humans. Here we conducted an extensive meta-analysis of monkey fMRI data collected under many experimental conditions to investigate whether a consistent network of regions—a default mode network—is present in macaques. Following identification of a DMN in the monkey, we analyzed its network structure in a task-independent manner using resting-state connectivity. Furthermore, we compared the resting-state connectivity results with published anatomical connectivity data from tracing studies in the monkey.
Materials and Methods
Task experiments and data collection.
Animal care met the Belgian and European guidelines and was approved by the K.U. Leuven Medical School. Experimental procedures were approved by the local Ethics Committees of the K.U. Leuven. Task-related fMRI data were acquired between 2003 and 2007 during 15 experiments conducted in our laboratory (Nelissen et al., 2005, 2006, 2009; Sawamura et al., 2005; Durand et al., 2007; Peeters et al., 2009; Wardak et al., 2010; Joly et al., 2011) (Table 1). They were collected in one female and nine male monkeys with a 1.5 T MR scanner (Siemens Sonata). All experiments followed a block design with fixation as control condition. The active task varied across experiments and belonged to one of the following groups of experiments: visual, attention, action-observation, and auditory. The resting-state fMRI data were collected in four fixating monkeys (three male and one female) during the years 2008–2009 with a 3.0 T MR scanner (Siemens Trio). Before scanning, a contrast agent (MION or Sinerem) was injected into the animal's femoral or saphenous vein (Vanduffel et al., 2001). Eye position was monitored using a pupil–corneal reflection tracking system (RK-726PCI, Iscan) at 60 Hz. Monkeys were rewarded for fixation on a small red target (0.30 × 0.35 visual degrees) in the middle of the screen (Vanduffel et al., 2001, 2002). The dimension of the fixation window was 2 × 3 visual degrees. The functional volumes were gradient-echo-planar images (GE-EPI) covering the whole brain.
Data preprocessing.
fMRI time series were preprocessed and analyzed using the SPM5.0 toolbox (Wellcome Trust Centre for Neuroimaging). Spatial preprocessing consisted of realignment and rigid coregistration with a template anatomy [monkey M12 = MM1 from Ekstrom et al. (2008), 0.35 mm isotropic voxels]. To compensate for EPI distortions and interindividual anatomical differences, functional images were warped by a nonrigid matching technique (Nelissen et al., 2005, 2006). The monkey F99 image was used to define the target stereotaxic space. The images were resliced to 1 mm isotropic voxels and smoothed with an isotropic Gaussian kernel at 2 mm FWHM. The resting-state fMRI data were additionally preprocessed in preparation for functional connectivity analysis, according to the procedure described by Vincent et al. (2007). The preprocessing steps included the following: (1) bandpass filtering between 0.0025 and 0.05 Hz; (2) regression of white matter and ventricle signals, and their first derivatives; and (3) regression of three-dimensional motion parameters, and their first derivatives.
Meta-analysis of task fMRI data.
A mask image capturing the full extent of the template brain was used for General Linear Model (GLM) estimation. Head-realignment and eye-movement parameters were included in the model as covariates of no interest. Experiment-specific deactivation maps were delineated by contrasting the β-value map in the fixation condition with the β-value map in all the active-task conditions. Statistical t-score maps were obtained from the β-value maps by fixed-effect analysis. A consensus map for deactivations across experiments was created by binarizing the deactivation maps from single experiments (p < 0.01, uncorrected), summing these binarized masks, and thresholding by number of subjects (N ≥ 8).
Definition of DMN areas.
Since the monkey DMN was defined though a meta-analysis across multiple subjects and experiments, uncertainty in the exact definition of included areas is possible, due to the alignment of multiple datasets to a template space and because fMRI responses often include regions that span multiple areas. The monkey DMN regions were determined by selecting local maxima (identified by their stereotaxic coordinates) in the consensus map. Each peak was projected on the anatomical image of the template (monkey F99) using the MRIcron software, and classified by comparing its position with respect to the monkey atlas of Paxinos et al. (2008). We opted to use this atlas because it already has an extended group of users (Frey et al., 2011). Subsequently, the peak and the monkey DMN map were projected together on the anatomical template image to evaluate the contribution of regions adjacent to the DMN peak. In cases where the extended region of deactivation clearly invaded neighboring monkey areas in the atlas, the classification was adapted accordingly and uncertainty acknowledged.
Functional connectivity analysis of resting-state fMRI data.
Regions of interest (ROIs) in the monkey DMN were generated using the classified DMN peaks. They were defined as spheres of 2 mm radius, centered on the respective stereotaxic coordinates. For each ROI, a representative time course was obtained by averaging the signals across voxels. Whole-brain connectivity maps were created and converted to Z-scores by the Fisher's r-to-z transformation (Zar, 1998), correcting the degrees of freedom for the autocorrelation in the time-series (Shumway and Stoffer, 2006). A fixed-effect analysis was used to create group-level correlation maps (Genest, 1992). Finally, the group Z-score maps were back-converted to correlation maps by the Fisher's z-to-r transformation (Zar, 1998).
Following the same statistical approach, a correlation matrix between the ROI time courses was calculated. ROIs without any significant correlation to other ROIs (p < 0.001) were excluded from subsequent analyses. A residual (unthresholded) correlation matrix between correlated ROIs was used for a network analysis across DMN nodes (Andrews-Hanna et al., 2010). The goal of this analysis was to determine from the pairwise regional correlations whether the DMN areas clustered into coherent subsystems (Buckner et al., 2008). Graph analysis of the correlation matrix was implemented using the Kamada–Kawai algorithm, a spring-embedded algorithm that pulls connected nodes together and pushes disconnected nodes apart in a manner such that the total energy of the system is minimized (Kamada and Kawai, 1989). Network hubs, defined as nodes connected to most of the other nodes, were identified by means of the betweenness-centrality. This is a quantitative measure of how connected a particular node is to other nodes (Freeman, 1977). The nonhub nodes, connected by definition to only part of the network, were analyzed to detect the presence of network subsystems. To this end, the hub nodes were excluded from the correlation matrix, and the remaining correlations were processed by hierarchical cluster analysis using the average linkage algorithm (Everitt et al., 2001). In addition, the network structure emerging from resting-state connectivity data was qualitatively compared to anatomical connectivity data in monkeys (efferent/afferent connections), partially available through the CoCoMac database (Kötter, 2004).
Spatial comparison with the human DMN.
We performed a spatial comparison between the monkey and the human DMN, on the basis of the monkey DMN nodes detected by our meta-analysis and of previous results in humans. As no meta-analysis on rewarded fixation in humans is available for a direct comparison with monkey data, we referred for the human DMN to the meta-analysis conducted on fMRI datasets published by Laird et al. (2009). For the definition of the areas in the network, we applied the same procedure as for the monkey DMN. In this case, we used the Colin image registered to the MNI space (as available in MRIcron) for anatomical reference. The cytoarchitectonic classification of human areas was conducted using the SPM Anatomy Toolbox (Eickhoff et al., 2005, 2010; Toga et al., 2006); for the areas left unclassified, we referred to comparative cytoarchitectonic studies (Vogt et al., 1987, 2006; Petrides and Pandya, 1999; Ongür et al., 2003). These procedures have multiple sources of uncertainty, so the correspondences should be taken as best estimates.
Results
We delineated brain areas potentially belonging to the monkey DMN by performing a meta-analysis of monkey fMRI data (Table 1). Specifically, we calculated a consensus map to define areas with consistent task-induced deactivation across different experiments (single-experiment threshold at p < 0.01, N ≥ 8). This analysis revealed a widespread network including the anterior and the posterior cingulate, and the posterior parietal cortices, confirming the existence of a monkey DMN. By overlaying the peaks in the consensus map onto the monkey atlas of Paxinos et al. (2008), we tentatively defined areas 9/46d, 24/32, 8b, caudal temporal parietal occipital (TPOC), 23, 31, v23, PGm, S2i, and PEa to belong to the monkey DMN (Fig. 1, Table 2).
After the definition of the monkey DMN, we examined whether this can be fractionated in functional–anatomical subdivisions, as suggested by human DMN studies (Andrews-Hanna et al., 2010; Sestieri et al., 2011). Based on a methodology already applied on human data (Andrews-Hanna et al., 2010), we measured resting-state connectivity between the monkey DMN nodes to identify in an independent manner the presence of hubs and subsystems. To this end, we used fMRI data collected in four monkeys during fixation (see Table 1). Seed-based functional connectivity analysis showed areas 9/46d and 31 to be connected to the largest part of the remaining DMN areas (Fig. 2A,B). The connectivity measures from other DMN nodes, as for example from areas 8b (Fig. 3A) and v23 (Fig. 3B), showed only parts of the DMN, suggesting the presence of subsystems.
To test specifically for the existence of subsystems in the DMN, we conducted a connectivity-based network analysis, which allowed us to quantitatively delineate relations among regions. The betweenness-centrality measure (Freeman, 1977) confirmed regions centered at or near areas 9/46d and 31 to be network hubs and to constitute the DMN core (Fig. 4A). A subset of areas, such as area S2i and PEa, were not found to be significantly correlated with any of the other nodes (p < 0.001). By applying cluster analysis on the connected and nonhub nodes, we defined two subsystems. These included areas 24/32, 8b, and TPOC and areas 23, v23, and PGm, respectively (Fig. 4B). Task-related activity profiles revealed the functional selectivity of the network core and its subsystems (see the significant differences in percentage signal change across tasks, Fig. 4C).
Discussion
A previous PET study in monkeys showed regional blood flow increases in the medial prefrontal and posterior cingulate cortex during rest as compared to an active working memory task, providing support for the existence of the monkey DMN (Kojima et al., 2009). Moreover, task-induced suppression of neuronal activity has been shown in the monkey posterior cingulate cortex during attention and working memory tasks (Hayden et al., 2009). By means of our meta-analysis on monkey fMRI data, we defined a set of brain regions that are consistently deactivated across many experimental conditions (Shulman et al., 1997; Gusnard and Raichle, 2001; Mazoyer et al., 2001; Buckner et al., 2008; Laird et al., 2009). This set of areas putatively constitutes a monkey correspondence of the human DMN. Our results are valuable for several reasons: (1) they allow an accurate description of the DMN functional architecture, relying on the large number of anatomical and neurophysiologic studies conducted in monkeys; (2) they show that the DMN is relatively well preserved during ∼25 million years of separate evolution between monkeys and humans; (3) they suggest novel putative functional correspondences between monkey and human brain areas, based on the DMN foci in the two species; and (4) they show that many, though not all, subunits of the DMN are monosynaptically connected with each other.
Monkey DMN components
In this study, the functional organization of the brain areas involved in the DMN was characterized in the monkey using multiple tasks. It is worth mentioning that deactivations from a single experiment may not be sufficient to consistently delineate the DMN, for which instead a meta-analysis across multiple experiments is more reliable (Shulman et al., 1997; Mazoyer et al., 2001; Buckner et al., 2008; Laird et al., 2009). The anatomical definition of the areas in the monkey DMN, classified using a recently published parcellation of the monkey cortex (Paxinos et al., 2008), may permit comparison of our results with the large body of available anatomical, connectional, and electrophysiological data in the monkey. However, caution is needed with regard to the heterogeneity in areal definitions across studies and inherent uncertainties of the data.
Our results of resting-state connectivity from the monkey DMN areas largely reflect known anatomical connections in the monkey brain, although the match is not perfect (Table 3). The monkey DMN, like its human equivalent, consists of a set of interconnected subsystems that converge on hubs including the posterior cingulate cortex (PCC) and regions within prefrontal cortex (PFC) (Buckner et al., 2008; Andrews-Hanna et al., 2010; Sestieri et al., 2011). These hubs correspond with areas 9/46d and 31, which are estimated to be anatomically connected to many of the other monkey DMN areas (Barbas and Pandya, 1989; Kobayashi and Amaral, 2003; Morecraft et al., 2004; Parvizi et al., 2006; Petrides and Pandya, 2006, 2007; Price, 2007; Seltzer and Pandya, 2009). Anatomical connectivity data fit well with functional connectivity based on resting-state data of two functionally distinct DMN subsystems, including respectively areas 24/32, 8b, and TPOC (first subsystem) and areas 23, v23, and PGm (second subsystem). Furthermore, not only resting-state based functional connectivity, but also anatomical connectivity showed that areas PEa and S2i are only weakly connected with the other DMN nodes. In a limited number of cases, we found resting-state connectivity in the absence of documented monosynaptic anatomical connectivity. This may imply either that functional connectivity exists between regions not monosynaptically connected or that anatomical connections exist but have not been described yet (Bullmore and Sporns, 2009), or alternatively may speak to the imperfections in our estimates.
The multisynaptic structure of the monkey DMN and the functional selectivity of its subsystems during goal-directed behavior suggest the network draws on multiple, specialized processing systems. This may explain why there is only partial overlap with the functional connectivity analyses. Functional connectivity analyses are based on correlated fluctuations of fMRI signals (Fox and Raichle, 2007), and therefore they cannot determine the presence of task-related activations or deactivations. Despite the similarity in the medial anterior and posterior areas, our monkey DMN map overlaps only in part with the functional connectivity map originating in the posterior cingulate cortex, as reported by Vincent et al. (2007). It was shown in later analyses of the same data that the choice of different seeds within the posteromedial cortex results in the definition of distinct correlated networks (Margulies et al., 2009). In this regard, our monkey DMN map fits best with the limbic–heteromodal network obtained by Margulies et al. (2009) by seeding specifically the region at or near areas 23/31. It will be interesting to explore further the relationship between network analysis of the monkey systems and their increased activity during passive epochs.
Comparison with the human DMN
The definition of the human DMN comes from multiple meta-analyses demonstrating consistent task-related deactivations across a network including the PFC, the PCC, and the inferior parietal lobule (Shulman et al., 1997; Mazoyer et al., 2001; Raichle et al., 2001; Buckner et al., 2008; Laird et al., 2009). Both anatomical (Parvizi et al., 2006) and functional (Vincent et al., 2007; Margulies et al., 2009) connectivity investigations focusing on the posterior cingulate cortex suggested correspondent systems the monkey.
To evaluate the potential similarity of the monkey DMN with its human equivalent, we performed a comparative analysis of the deactivation foci between species. To this end, we referred for the human DMN to the results of a previous meta-analysis conducted by activation likelihood estimation on 1711 neuroimaging publications and 7920 experimental contrasts (Laird et al., 2009). We did not carry out quantitative interspecies comparisons, since the tasks performed by monkeys and humans in the two meta-analyses were not the same. The qualitative comparison between spatial maps showed similarities in the broad pattern of deactivated areas (Fig. 5). Putative monkey–human correspondences were primarily found in the vicinity of the medial and dorsolateral prefrontal cortex, the insula/operculum, the lateral parietal cortex, and the posterior cingulate cortex. In both species, the posterior cingulate included multiple foci. It is possible that not only anatomical correspondence, but also functional similarities exist between areas 23/31 of monkeys and humans (Vogt et al., 1987, 2006); however, additional evidence is needed to substantiate such a claim.
We also detected possible interspecies functional similarity in prefrontal cortex for area 8b, as already predicted by comparative cytoarchitectonic studies (Petrides and Pandya, 1999). It is less clear whether human area OP4, for which a homology with monkey area PV has been proposed (Eickhoff et al., 2010), may instead correspond to monkey area S2i. The monkey DMN areas in PFC correspond only in part to human DMN areas according to classical principles of homology and cortical proximity (Orban et al., 2004). For example, monkey area 32 in medial PFC may correspond to two parts in the human brain, including regions 9/10/32ac and 10/24/32 pl, respectively. This result suggests that the prefrontal DMN areas underwent a functional reorganization, in line with the disproportional anatomical expansion of these regions throughout primate evolution (Ongür et al., 2003; Buckner et al., 2008). Interestingly, we also observed deactivations within the posterior-medial intraparietal sulcus in monkeys, but not in humans. Future investigations will be necessary to clarify whether this difference results from the different tasks performed by monkeys and humans or is more directly related to evolutionary changes in the functional architecture of the primate brain.
Putative species-general DMN functions
Several theories have proposed that the DMN is involved in broad external attention or monitoring. An alternative hypothesis is that the DMN participates in internal modes of cognition evoked spontaneously during passive moments (see Buckner et al., 2008 for a recent review). It is currently unclear which cognitive abilities are present in the monkey and which instead would be exclusively human. The present results suggest that, like humans, monkeys default to consistent sets of processes when left in undirected states. It is tantalizing to infer that monkeys engage in forms of spontaneous cognition detached from the external environment during idle moments. The observation that the anatomy of the DMN in monkeys includes heteromodal association areas and not sensory regions supports this speculative possibility.
In conclusion, our study robustly confirms the existence of a default mode of brain function in monkeys. The topological similarities of the monkey DMN with the human DMN imply a role of this system that generalizes across species, possibly related to internally directed cognition. The DMN may have acquired expanded functions throughout evolution, along with the development of linguistic, self-referential, and social abilities that characterize and distinguish human from nonhuman primates (Schilbach et al., 2008; Passingham, 2009). Further investigations of the monkey DMN areas, for example with intracranial recordings or with fMRI combined with excitatory/inhibitory interventions (e.g., Ekstrom et al., 2008), are needed to evaluate the electrophysiological and behavioral correlates of the default-mode activity and to better understand the functional role of the DMN.
Footnotes
D.M. and K.N. are postdoctoral fellows of the Fonds Wetenschappelijk Onderzoek–Vlaanderen (F.W.O.). This work received support from European Union Grant Brainsync FP7 FWP-200728, Inter University Attraction Pole 6/29, Programme Financing PFV/10/008, Geconcerteerde Onderzoeks Actie 10/19, Hercules II funds, Fonds Wetenschappelijk Onderzoek–Vlaanderen Grants G062208N10 and G083111N10, and National Science Foundation Grant BCS-0745436. The Martinos Center is supported by National Center for Research Resources Grant P41RR14075. We thank A. Coeman, C. Fransen, M. Depaep, W. Depuydt, P. Kayenbergh, G. Meulemans, and S. Verstraeten for technical assistance, Prof. A. Laird for the human data, and Prof. R. Tootell for his comments on the manuscript.
- Correspondence should be addressed to Dr. Wim Vanduffel, Athinoula A. Martinos Center for Biomedical Imaging, 13th Street, Building 149, Charlestown, MA 02129. wim{at}nmr.mgh.harvard.edu