Abstract
Despite the prominence of parietal activity in human neuroimaging investigations of sensorimotor and cognitive processes, there remains uncertainty about basic aspects of parietal cortical anatomical organization. Descriptions of human parietal cortex draw heavily on anatomical schemes developed in other primate species, but the validity of such comparisons has been questioned by claims that there are fundamental differences between the parietal cortex in humans and other primates. A scheme is presented for parcellation of human lateral parietal cortex into component regions on the basis of anatomical connectivity and the functional interactions of the resulting clusters with other brain regions. Anatomical connectivity was estimated using diffusion-weighted magnetic resonance image (MRI)-based tractography, and functional interactions were assessed by correlations in activity measured with functional MRI at rest. Resting-state functional connectivity was also assessed directly in the rhesus macaque lateral parietal cortex in an additional experiment, and the patterns found reflected known neuroanatomical connections. Cross-correlation in the tractography-based connectivity patterns of parietal voxels reliably parcellated human lateral parietal cortex into 10 component clusters. The resting-state functional connectivity of human superior parietal and intraparietal clusters with frontal and extrastriate cortex suggested correspondences with areas in macaque superior and intraparietal sulcus. Functional connectivity patterns with parahippocampal cortex and premotor cortex again suggested fundamental correspondences between inferior parietal cortex in humans and macaques. In contrast, the human parietal cortex differs in the strength of its interactions between the central inferior parietal lobule region and the anterior prefrontal cortex.
Introduction
Research with both human and nonhuman primates shows that the parietal cortex is involved in numerous cognitive functions, from perception to control of action and higher-order processes, such as numerosity and memory. It is well established that parietal cortex consists of a mosaic of subregions, which are distinct both functionally and anatomically (Grefkes and Fink, 2005; Caminiti et al., 2010; Nelson et al., 2010). Our knowledge of the different subregions and their functional properties is primarily based on the macaque model. However, the suitability of this model has been called into question (Passingham, 2009). The parietal cortex has expanded substantially in the human brain compared with other primates and even other hominids, which might have resulted in an increase or a rearrangement of its regions (Vanduffel et al., 2002b; Grefkes and Fink, 2005). Particular emphasis has been given to the possibility that human inferior parietal lobule is quite distinct to that of other species (Husain and Nachev, 2007).
The function of any brain area is constrained by its extrinsic connections. Indeed, patterns of anatomical connectivity can inform the segregation of functionally distinct areas (Passingham et al., 2002; Averbeck et al., 2009). Here, we aim to identify and characterize key component regions of human parietal cortex using a combination of non-invasive magnetic resonance imaging (MRI) techniques: probabilistic diffusion tractography and resting-state functional connectivity. The primary aim is to provide a parcellation of human parietal cortex that can inform non-invasive neuroimaging studies. An additional aim is to compare the parietal regions in the human brain with macaque parietal regions whose neurophysiology and connectional anatomy are well studied.
First, we determine the connectivity profiles of all human parietal voxels using diffusion-weighted magnetic resonance imaging (DW-MRI) and tractography and then parcellate the parietal cortex into distinct clusters based on the cross-correlations between these connectivity profiles (Johansen-Berg et al., 2004; Anwander et al., 2007; Tomassini et al., 2007; Beckmann et al., 2009; Schubotz et al., 2010). Second, we examine the correlations between spontaneous activity in voxels in the parietal cortex and in other brain areas, a technique known as “resting-state functional connectivity” (O'Reilly et al., 2010; Xiang et al., 2010). The concept behind this resting-state functional MRI (fMRI) is that when the brain is “idle,” i.e., not involved in any explicit task, correlations in slowly fluctuating spontaneous brain activity tend to reflect the intrinsic functional networks of the brain (Fox and Raichle, 2007). We focused on a limited set of eight areas known, in macaques, to project to circumscribed regions of parietal cortex. By comparing the correlations between the eight target areas with the tractography-based parietal clusters, we can characterize the distinct parietal subregions involved in various large-scale corticocortical interactions. Given the novelty of this approach, we first investigate resting-state interactions between parietal cortex and the same target areas in the macaque brain, to assess the degree to which the resting-state functional connectivity reflects known anatomical connections and thus to validate this method.
Materials and Methods
The research described in this paper consists of three parts. First, we parcellate the human parietal cortex on the basis of DW-MRI tractography. Second, we explore the potential of resting-state functional connectivity to describe interactions between the parietal cortex and other brain regions in both rhesus macaques and humans. Performing these analyses in macaques allow us to assess the degree to which functional connectivity reflects the prominent known anatomical connections in the same species. It is necessary to use macaques because the connections of macaque parietal cortex are well documented, but very little is known of the connections of human parietal cortex. Finally, having established the potential of resting-state functional interactions as indices of connectivity, we describe the patterns of functional interaction associated with each of the clusters found in the human data tractography-based parcellation of the human parietal cortex (Fig. 1A).
Diffusion-weighted data and tractography-based parcellation.
Diffusion-weighted images were acquired in eight healthy subjects (four female; age range, 20–36 years; mean ± SD age, 26.9 ± 6.0 years) on a 1.5 T Siemens Sonata MR scanner, with maximum gradient strength of 40 mT × m−1. This number of participants has been established previously as sufficient for obtaining reliable parcellation results (Klein et al., 2007). All participants gave informed written consent in accordance with ethical approval from the local ethics committee. Participants lay supine in the scanner, and cushions were used to reduce head motion. Diffusion-weighted data were acquired using echo planar imaging (72 × 2 mm thick axial slices; matrix size, 128 × 104; field of view, 256 × 208 mm2; giving a voxel size of 2 × 2 × 2 mm). Diffusion weighting was isotropically distributed along 60 directions using a b value of 1000 s × mm−2. For each set of diffusion-weighted data, five volumes with no diffusion weighting were acquired throughout the acquisition. Three sets of diffusion-weighted data were acquired for subsequent averaging to improve the signal-to-noise ratio. The total scan time for the diffusion-weighted imaging protocol was 45 min. A structural scan was acquired for each participant in the same session, using a T1-weighted three-dimensional fast, low-angle shot (3D FLASH) sequence [repetition time (TR), 12 ms; echo time (TE), 5.65 ms; flip angle, 19°; with elliptical sampling of k space, giving voxel size of 1 × 1 × 1 mm].
Analyses were performed using tools from FreeSurfer (Athinoula A. Martinos Center for Biomedical Imaging), FSL (for Functional MRI of the Brain Software Library) (Smith et al., 2004), Caret (Van Essen et al., 2001), and custom software written in Matlab (MathWorks). A parietal cortex region of interest (ROI) was drawn by hand in each participant's right hemisphere directly on the cortical surface model produced by FreeSurfer (Fig. 1B–D). In the absence of detailed cytoarchitectonic measurements, it is necessary to use macroscopic boundaries that can be reliably identified in all individuals as the boundaries of the ROI. The fundus of the inferior and superior postcentral sulci formed the rostral boundary of the ROI in its ventral and dorsal aspects, respectively. The ventral boundary was drawn between the ventral tip of the inferior postcentral sulcus and the conjunction of the horizontal and vertical branches of the Sylvian fissure, in the ventral supramarginal gyrus, and then to the conjunction of the ascending posterior segment of the superior temporal sulcus and horizontal posterior segment of the superior temporal sulcus. The boundary then followed the horizontal posterior segment of the superior temporal sulcus and continued to the nearest point on the intraparietal sulcus (IPS). The crown of the hemisphere constituted the medial boundary of the ROI. The point on the crown of the hemisphere nearest the parieto-occipital sulcus was the posterior boundary for the medial aspect of the mask, and it extended from there to the nearest point on the IPS. The lateral and medial banks of the full length of the entire IPS were included in the mask. The expanded fundus of the descending segment of the IPS, which is thought to correspond to extrastriate visual area V3A (Swisher et al., 2007), was excluded from the mask. Although there is individual variation in some aspects of the sulcal anatomy of the parietal cortex, these boundaries identified a similar location in the brains of all subjects after affine registration into Montreal Neurological Institute (MNI) space. It should be noted that, although the presence of both functional and anatomical differences between the hemispheres is well established in the human brain, especially in relation to language, the presence of hemispheric differences in the monkey brain is still a topic of debate (for review, see Passingham, 2008). Therefore, this study focuses on the right hemisphere, which is not dominant for language in humans.
DW-MRI data were preprocessed using tools from FDT (for Functional MRI of the Brain Diffusion Toolbox; part of FSL 4.1). Eddy-current distortions were corrected using affine registration of all volumes to a target volume with no diffusion weighting. Voxelwise estimates of the fiber orientation distribution were calculated using Bedpostx, limited to estimating two fiber orientations at each voxel, because of the b value and number of gradient orientations in the diffusion data (Behrens et al., 2007). For each participant, probabilistic tractography was run from vertices at the gray matter/white matter boundary surface within the parietal ROI to assess connectivity with every brain voxel (downsampled to 5 mm isotropic voxels), using a model accounting for multiple fiber orientations in each voxel (Behrens et al., 2007). Crucially, tractography was seeded from surface vertices (as opposed to voxels). From each seed vertex, the normal to the cortical surface was given by the FreeSurfer cortical model, and this information was used to track toward the brain instead of tracking in both directions as is usually done in voxel-based tractography. This approach dramatically reduces the risks of creating spurious sample tracts that cross gyral walls (for details, see www.fmrib.ox.ac.uk/fsl/fdt/fdt_surface.html). It also provides quantitative estimates of connectivity from every surface vertex of the parietal mask to every voxel in the rest of the brain.
A connectivity matrix between parietal vertices and each brain voxel was derived as described previously (Johansen-Berg et al., 2004) and used to generate a symmetric cross-correlation matrix of dimensions (number of seeds × number of seeds) in which the (i,j)th element value is the correlation between the connectivity profile of seed i and the connectivity profile of seed j. The rows or this cross-correlation matrix were then permuted using k-means segmentation for automated clustering to define different clusters. The goal of clustering the cross-correlation matrix is to group together regions that share the same connectivity with the rest of the brain. To increase the chances of obtaining continuous (non-scattered) clusters, we included a distance constraint (Tomassini et al., 2007). The parietal cortex contains a deep sulcus, the IPS, and the Euclidean distance between its banks is small although they are far apart in cortical distance. The flattening of the brain in FreeSurfer space, however, made it possible to implement the distance constraint along a geodesic that followed the cortical surface. The resulting clusters are then constrained to consist of vertices that are spatially contiguous, although the border between clusters is still guided by remote connectivity information.
The number of clusters in the k-means clustering must be set by the experimenter. To determine the optimal number of clusters resulting in consistency across participants, we used an iterative procedure (Beckmann et al., 2009). First, we looked for clusters to separate the parietal cortex mask into three separate partitions: a superior parietal lobule (SPL) mask, an inferior parietal lobule (IPL) mask, and a posterior IPS mask (Fig. 1D). Following Beckmann et al. (2009), each of these masks was then parcellated into the maximum number of areas that led to a consistent parcellation in all eight participants. Each of the participant's individual clusters was then transformed from FreeSurfer space to the MNI template brain. The center of gravity of each cluster in MNI space and its overlap with the FSL implementation of the Harvard–Oxford probabilistic cortical atlas, an atlas based on manual definition of sulcal/gyral landmarks, was then established.
Resting-state fMRI data acquisition and preprocessing.
Human resting-state fMRI data and T1-weighted images were collected for a different group of 12 healthy volunteers (nine females; age range, 31–61 years; mean ± SD age, 43.08 ± 9.17 years). This group of participants did not overlap with the group of participants in the diffusion-weighted imaging experiment. All participants gave informed written consent in accordance with ethical approval from the local ethics committee. Participants lay supine in a 1.5 T Siemens Sonata MR scanner. They were instructed to close their eyes and lie still. Cushions were used to reduce head motion. Whole-brain blood oxygen level-dependent (BOLD) fMRI data was collected for 11 min from each participant, using the following parameters: 45 axial slices; in-plane resolution, 3 × 3 mm; slice thickness, 3 mm; no slice gap; TR, 3400 ms; TE, 41 ms; 200 volumes. A structural scan was acquired for each participant in the same session, using a T1-weighted 3D FLASH sequence (voxel size of 1 × 1 × 1 mm).
Data were analyzed using tools from FSL (Smith et al., 2004) and SPM (for Statistical Parametric Mapping) (www.fil.ion.ucl.ac.uk/spm). The first six volumes of each functional dataset were discarded, after which the preprocessing was performed: motion correction, non-brain removal, spatial smoothing [using Gaussian 5 mm full-width at half-maximum (FWHM) kernel], grand-mean intensity normalization of the entire four-dimensional dataset by a single multiplicative factor, and high-pass temporal filtering (Gaussian-weighted least-squares straight line fitting, with σ = 50.0 s). Registration of functional images to the skull-stripped structural and MNI template was done using FLIRT (Jenkinson et al., 2002).
Macaque resting-state fMRI and anatomical scans were collected for 10 healthy macaques (Macaca mulatta; five females; age range, 3 years and 7 months to 4 years and 10 months; mean age, 4 years and 5 months; weight range, 4.34–8.53 kg). Protocols for animal care, magnetic resonance imaging, and anesthesia were performed under authority of personal and project licenses in accordance with the United Kingdom Animals (Scientific Procedures) Act (1986). Anesthesia was induced using intramuscular injection of ketamine (10 mg/kg), xylazine (0.125–0.25 mg/kg), and midazolam (0.1 mg/kg). Macaques also received injections of atropine (0.05 mg/kg, i.m.), meloxicam (0.2 mg/kg, i.v.), and ranitidine (0.05 mg/kg, i.v.). Local anesthetic (5% lidocaine/prilocaine cream and 2.5% bupivacaine injected subcutaneously around the ears to block peripheral nerve stimulation) was also used at least 15 min before placing the macaque in the stereotaxic frame. The anesthetized animals were placed in an MRI-compatible stereotactic frame (Crist Instruments) in a sphinx position and placed in a horizontal 3 T MRI scanner with a full-size bore. Scanning commenced ∼2 h after induction, when the ketamine was unlikely still to be present in the system. Anesthesia was maintained using the lowest possible concentration of isoflurane to ensure that macaques were lightly anesthetized. The depth of anesthesia was assessed using physiological parameters (heart rate and blood pressure, as well as clinical checks before the scan for muscle relaxation). During the acquisition of the functional data, the inspired isoflurane concentration was in the range 1.0–1.8% (mean, 1.45%), and the expired isoflurane concentration was in the range 0.9–1.7% (mean, 1.38%). Isoflurane was selected for the scans as resting-state networks have been demonstrated previously to be present using this agent (Vincent et al., 2007). Macaques were maintained with intermittent positive pressure ventilation to ensure a constant respiration rate during the functional scan, and respiration rate, inspired and expired CO2, and inspired and expired isoflurane concentration were monitored and recorded using VitalMonitor software (Vetronic Services Ltd.). In addition to these parameters, core temperature and SpO2 were monitored throughout the scan. A four-channel phased-array coil was used for data acquisition (H. Kolster, MRI Coil Laboratory, Laboratory for Neuro- and Psychophysiology, Katholieke Universiteit Leuven, Leuven, Belgium). Whole-brain BOLD fMRI data was collected for 53 min and 26 s from each animal, using the following parameters: 36 axial slices; in-plane resolution, 2 × 2 mm; slice thickness, 2 mm; no slice gap; TR, 2000 ms; TE, 19 ms; 1600 volumes. A structural scan (three averages) was acquired for each macaque in the same session, using a T1-weighted magnetization-prepared rapid-acquisition gradient echo sequence (either 0.5 × 0.5 × 0.5 or 0.5 × 0.5 × 1.0 mm voxel resolution).
The first six volumes of each functional dataset were discarded, and the following preprocessing was performed: non-brain removal, 0.1 Hz low-pass filtering to remove respiratory artifacts, motion correction, spatial smoothing (using Gaussian 3 mm FWHM kernel), grand-mean intensity normalization of the entire four-dimensional dataset by a single multiplicative factor, and high-pass temporal filtering (Gaussian-weighted least-squares straight line fitting, with σ = 50.0 s). Registration of functional images to the skull-stripped structural and a macaque template was done using FLIRT (Jenkinson et al., 2002).
Resting-state fMRI data analysis.
To establish the functional connectivity of the parietal subregions yielded by the tractography-based parcellation, we created a number of target masks in MNI space from the probable homologs of brain areas that are known to be connected to the parietal cortex in the macaque. Functional interactions were determined with target areas in the ventral premotor cortex (PMv) centered at [56, 6, 22] (Binkofski et al., 1999), dorsolateral prefrontal cortex (dlPFC) centered at [40, 40, 24] (Rajkowska and Goldman-Rakic, 1995), the parahippocampal gyrus (PH) centered at [28, −34, −16] (Suzuki et al., 2005), the anterior prefrontal cortex (aPFC) centered at [42, 50, −2] (Mars et al., 2008; Boorman et al., 2009), V5/middle temporal area (MT) centered at [44, −67, 0] (Dumoulin et al., 2000), inferior dorsal premotor cortex (PMd) centered at MNI coordinates [36, −2, 50] (Hinkley et al., 2009), PMd centered at MNI coordinates [30, −17, 64], (Amiez et al., 2006), and the frontal eye fields (FEF) centered at [24, −6, 44] (see Figs. 4, 9). All masks were of identical size (6 mm3), so as not to bias the analysis. Equivalent masks were created on the macaque template brain. The equivalent coordinates were determined by transferring the masks, using affine transformation, to the macaque template of McLaren et al. (2009), which uses the coordinates system of the atlas by Saleem and Logothetis (2006). The coordinates of the target areas were PMv [16, 24, 24], dlPFC [17, 37, 22], PH [11, 0, 7], aPFC [10, 44, 17], V5 [20, −1, 18], inferior PMd [16, 22, 29], PMd [6, 20, 36], and FEF [15, 27, 25] (see Figs. 3A, 8A).
We used a voxelwise approach to map resting-state functional connectivity between the parietal cortex and characteristic time series associated with specific target regions. This analysis was performed in both human and macaques to (1) establish the reliability of the resting-state functional connectivity method as an index of known structural connectivity in the macaque brain and (2) allow a preliminary comparison of the functional interactions between the macaque and human brain. First, we calculated the major Eigen time series representing activity in each of the target region. The major Eigen time series is the single time series that best reflects coherent activity across the mask in that it represents the largest amount of variance across the set of voxels within the region. Then, we calculated the correlation between each voxel in the parietal cortex and the first Eigen time series of specific target areas for each participant. This was done using the seed-based correlation analysis tool, which is part of FSL (fsl_sbca) and is described by O'Reilly et al. (2010). Correlations between the parietal voxels and each target area were calculated separately. The average time series of the whole brain and the six movement parameters were used as confound regressors. The resulting correlation masks were then entered into a second-level group general linear model analysis. The resulting images were thresholded at p < 0.05 for the parietal cortex region of interest (cf. O'Reilly et al., 2010).
Having established the reliability of the resting-state fMRI method and compared the patterns of correlation between the macaque and human brains, we then sought to formally compare the pattern of resting-state functional connectivity in parietal cortex associated with each of the target masks (PMv, dlPFC, PH, aPFC, V5, inferior PMd, PMd, and FEF) with the location of each of the parietal clusters established by tractography-based parcellation. Functional correlation maps of the target areas with the parietal cortex were overlaid with the thresholded clusters resulting from the tractography-based parcellation to reveal the resting-state connectivity fingerprint (Passingham et al., 2002), illustrating the interactions of each target area with the parietal clusters. For this analysis, the thresholded resting-state functional correlation map of each target area was binarized and multiplied by the thresholded (voxels belonging to a particular cluster in at least five participants), binarized tractography-based cluster masks. For each of the resulting overlay images resulting from this multiplication, the number of nonzero voxels indicates the strength of overlap between the area of functional interaction with a given target area and the tractography-based parietal cluster. The number of voxels in the overlay images were then first corrected for the size of each cluster and then normalized with respect to the total number of voxels correlating with the target mask across all clusters.
Results
Inferior parietal lobule
Tractography-based parcellation of the IPL subdivision (Fig. 1D, blue) identified five clusters, arranged consistently along the anteroposterior axis in all participants (Fig. 2A,B). The most anterior cluster (Fig. 2, red) with a center of gravity at [49, −25, 30] was most likely to be located in the parietal operculum according to Harvard–Oxford probabilistic atlas based and the Eickhoff et al. Anatomy Toolbox (Eickhoff et al., 2005). This region was similar in shape and location to the parietal opercular region (PFop) defined on cytoarchitectonic grounds by Caspers et al. (2006, 2008). This region, at least in the left hemisphere, is activated when people, but not macaques, observe tools being used, even if the macaques have experience of using the tools (Peeters et al., 2009). The tractography-based parcellation makes clear that this region is distinct from a more posterior region that is more commonly activated during grasping movements.
The two clusters (Fig. 2, blue and green) located directly posterior to PFop have centers of gravity at [53, −32, 44] and [50, −44, 43] and were most likely to be located in anterior and posterior parts of the supramarginal gyrus. Again the areas correspond with ones defined on the basis of cytoarchitecture by Caspers et al. (2008). The more anterior of the two clusters (Fig. 2, blue) resembles the conjunction of the PF and PFt regions defined by Caspers and colleagues. It is near the medial border of this region that activity is most often seen when human subjects grasp objects under visual guidance or manipulate them under visual guidance (Binkofski et al., 1999; Grefkes et al., 2002; Culham et al., 2003; Frey et al., 2005; Grol et al., 2007). As such, it might be a homolog of the monkey anterior intraparietal area (AIP). The more posterior of the two clusters (Fig. 2, green) resembles the PFm region defined by Caspers et al. (2008). Activity is recorded in this region under quite distinct circumstances. For example, when people make decisions, activity here increases in proportion with the evidence that the alternative choice would have been the better one to have taken and is especially prominent at the point that subjects switch between choices (Daw et al., 2006; Boorman et al., 2009).
The next most posterior cluster (Fig. 2, magenta), with a center of gravity of [46, −55, 45], corresponds most to the angular gyrus. Activity in this region is recorded when people redirect visuospatial attention from one location to another (Mort et al., 2003; Thiel et al., 2004). The most posterior cluster (Fig. 2, yellow/orange), with a center of gravity at [37, −67, 39], covers the most posterior and ventral parts of the angular gyrus and the most anterior part of the lateral occipital cortex. Activity in the most posterior IPL region has been noted when people successfully retrieve memories (Henson et al., 1999). Once again, the clusters correspond to cytoarchitectonic regions reported by Caspers et al. (2008). The anterior angular gyrus cluster resembles PGa and the posterior angular gyrus cluster resembles PGp.
We performed an explorative resting-state functional connectivity analysis looking at the correlation of IPL voxels with four target regions: PMv, PH, dlPFC, and aPFC. In the macaque, PMv is known to be strongly connected to the anterior parts of the IPL, including areas AIP, PF, and PFG (Matelli et al., 1986; Tanne-Gariepy et al., 2002). Resting-state functional connectivity of PMv with the IPL indeed was confined to the anterior IPL in macaques, suggesting that it reflects the known anatomical connections (Fig. 3). A similar pattern, with predominant interactions between anterior IPL and PMv, was seen in the human (Fig. 4A). In the macaque, the ventral part of dlPFC, where we placed our target, is also known to be interconnected with anterior IPL (Pandya and Yeterian, 1996), and indeed resting-state functional connectivity with dlPFC was most prominent in the same region of the IPL (Fig. 3). In contrast, PH is most strongly connected to posterior IPL regions in the macaque (Seltzer and Pandya, 1984; Suzuki and Amaral, 1994; Blatt et al., 2003). Once again, this pattern was also evident in the correlations between PH and IPL in both the macaque and the human (Fig. 4C). These results illustrate that resting-state functional connectivity reflects some of the most prominent anatomical connections found in macaque tracing studies and provides preliminary evidence of similarities between human and macaque IPL.
It has been suggested that lateral parts of the parietal cortex (Denys et al., 2004a) or, perhaps more specifically, the IPL may have particularly expanded in humans compared with other primates (Simon et al., 2004; Husain and Nachev, 2007). Because there are a number of similarities between posterior IPL and anterior IPL in both humans and macaques, the central part of the IPL is the part of the IPL that is particularly likely to have expanded in humans compared with macaques (Rushworth et al., 2009). This region is engaged in complex cognitive tasks, such as numerical processing, task exploration, and task change (Simon et al., 2002; Daw et al., 2006). We therefore investigated the functional connectivity between the IPL and an area that is particularly extended in the homonids compared with the macaque brain, namely the aPFC (Semendeferi et al., 2001). In macaque, parietal cortex and aPFC are not strongly connected; there are no connections between parietal cortex the frontal polar area 10 (Petrides and Pandya, 2007), although there are some between the PFG region in the central part of the IPL and anterior parts of lateral prefrontal cortex (Rozzi et al., 2006). In humans, aPFC showed strong interactions with the central parts of the IPL (Fig. 4D). In contrast, there was very little evidence for interactions between IPL and the aPFC regions used in the current investigation in the macaque even when the statistical threshold was reduced (Figs. 3, 4D).
These results indicate that IPL interactions with lateral PFC in the macaque are not as prominent as in humans or do not extend as far forward as in humans. Two additional analyses were performed to further investigate these results. First, to test whether the absence of interactions between aPFC and the central IPL in the macaque was simply a consequence of a difference in signal strength in aPFC in the two species, we looked for evidence of interactions between aPFC and the superior temporal cortex, a posterior brain area near to, but outside, the parietal cortex; connections between these two structures are known to exist in the macaque (Petrides and Pandya, 2007). Evidence for superior temporal/aPFC interactions was apparent, indicating that the absence of IPL/aPFC interactions in the macaque is not attributable to lack of signal in the aPFC. Second, we further investigated the pattern of IPL/lateral PFC interactions in the macaque by placing a seed in macaque area PFG, in the medial IPL (coordinates [22, 5, 29] in the atlas of Saleem and Logothetis, 2006) and establishing the location of functional interactions between this seed and the PFC. As shown in Figure 5A, lateral PFC interactions were evident near the location of our dlPFC seed but did not extend farther forward. In addition, prominent interactions were present between IPL and medial frontal cortex, particularly the cingulate motor areas (Fig. 5B).
Finally, the human resting-state functional connectivity was formally compared with the tractography-based parcellation described above. Resting-state functional connectivity was determined between the five clusters reported in the parcellation and the PMv, dlPFC, PH, and aPFC target masks in second-level random-effects analyses, which allow inferences at the population level. These four target masks correlated with different IPL clusters (Fig. 6). The most anterior IPL clusters (particularly red and blue clusters in Fig. 2A) showed strong interactions with PMv. This result is consistent with earlier diffusion-tractography work, showing that this region was most likely to connect to the ventral frontal cortex via the third branch of the superior longitudinal fasciculus (Croxson et al., 2005; Rushworth et al., 2006). It is also consistent with the fact that PMv, like the second most anterior IPL cluster, is active during reaching and grasping task (Tomassini et al., 2007). The second and third most anterior IPL regions (Fig. 2A, blue and green clusters) showed the strongest interactions with dlPFC. The most posterior IPL cluster (Fig. 2, yellow), covering the posterior part of the angular gyrus, showed interactions with PH, again mimicking results obtained from diffusion tractography (Rushworth et al., 2006). As such, the results from the most anterior and most posterior portions of the IPL resemble those found in macaques. Resting-state correlation with the aPFC, which was much stronger in the human compared with the macaque, was strongest in the second and third most posterior clusters (Fig. 2A, magenta and green). Both this IPL region and aPFC are active when participants are implementing a decision to change their behavior based on growing evidence in favor of an alternative course of action (Fig. 2C) (Daw et al., 2006; Boorman et al., 2009).
Superior parietal lobule and intraparietal sulcus
The SPL and adjacent posterior IPS were consistently parcellated into four and two subdivisions, respectively. One cluster within the posterior IPS regions was located around [32, −76, 24] and overlapped with visual area V7 (Swisher et al., 2007) outside the parietal cortex. The remaining five regions, however, have consistently been designated parietal cortex, and so our analyses focused on them (Fig. 7). One area (Fig. 7, red) was located in an anterior position on the medial bank of the IPS with a center of gravity at [30, −41, 53]. This area overlaps with the region identified as the ventral intraparietal area (VIP) on the basis of its activation by both visual and somatosensory stimulation adjacent to and on the face by Sereno and Huang (2006) and area 7PC on the basis of cytoarchitecture (Scheperjans et al., 2008b). This region is also close to an area identified by Bremmer et al. (2001) as involved in motion processing using visual, tactile, and auditory stimuli and subsequently suggested as the equivalent of macaque VIP. Adjacent to this area, but located more medially, was an area (Fig. 7, blue) with a center of gravity at [12, −50, 63] that overlaps with the cytoartchitectonically defined area 5L (Scheperjans et al., 2008b). Activity in this area is present during reaching movements and when reaching movements are adjusted when following a moving target (Diedrichsen et al., 2005).
Behind these two anterior SPL regions, we located a farther cluster in the more posterior part of the medial bank of the IPS (Fig. 7, green), with a center of gravity at [28, −55, 55] that overlapped with the hIP3 region identified on the basis of cytoarchitecture by Scheperjans et al. (2008b). Activity is also seen in this area during visually guided reaching (Hinkley et al., 2009) and pointing (Astafiev et al., 2003), other visually guided hand movements (Grefkes et al., 2004) and when visuomotor contingencies are changed and updated (Rushworth et al., 2001). A retinopically organized area, IPS4, responding to visual stimulation has also been reported at this location (Swisher et al., 2007). A very similar region has been referred to as DIPSA (Vanduffel et al., 2002b; Denys et al., 2004a; Durand et al., 2009; Georgieva et al., 2009), aIPS (Astafiev et al., 2003), and SPL ant (Koyama et al., 2004), and a region referred to as ant IPs (Corbetta et al., 2000) lies slightly more ventrally but close by. Because activity in this region is particularly strong when hand and arm movements are made or hand and arm intentions are updated (Rushworth et al., 2001; Grefkes et al., 2004; Hinkley et al., 2009), it is reminiscent of the medial intraparietal region (MIP) in the macaque (Snyder et al., 1997). An alternative account, based partly on the shape sensitivity of the area, has suggested that it may resemble a posterior subregion of the AIP area found in the macaque (Orban et al., 2006; Durand et al., 2007, 2009).
A fourth region was situated in a relatively medial and posterior position (Fig. 7, magenta) with a center of gravity at [19, −63, 53]. Its position therefore overlapped with the region that Scheperjans et al. (2008b) referred to as human 7A [and which bears no relation to the posterior IPL region in the macaque, which unfortunately has the same name (Cavada and Goldman-Rakic, 1989)]. The slightly more lateral position of the center of gravity of the area in the current study may be a consequence of not having included as much cortex on the medial wall in the present analysis. A retinotopic map, now often referred to as IPS3, has been reported in this area (Swisher et al., 2007) or close by (Sereno et al., 2001), and activity is more prominent in this region when the rules for guiding the allocation of visual attention are changed as opposed to the rules for guiding hand movements are changed (Rushworth et al., 2001). It has therefore been pointed out that this region bears some resemblance to the lateral intraparietal area (LIP) in the macaque. A very similar region has also been called DIPSM (Vanduffel et al., 2002b; Denys et al., 2004a; Durand et al., 2009; Georgieva et al., 2009), pIPS (Astafiev et al., 2003), and SPL post (Koyama et al., 2004).
Another cluster (Fig. 7, yellow) occupied a more ventral portion of the medial bank of the posterior IPS with a center of gravity at [21, −78, 43]. Its position corresponded closely to that of the region referred to as POIPS (Vanduffel et al., 2002b; Denys et al., 2004a; Durand et al., 2009; Georgieva et al., 2009). The two most posterior retinopically mapped divisions of parietal cortex, IPS1 and IPS2, fall within this region (Silver et al., 2005; Swisher et al., 2007). Activity in this region is modulated by the allocation of visual attention and during delays while subjects memorize targets for intended eye movements (Schluppeck et al., 2005, 2006; Szczepanski et al., 2010). Again, similarities between IPS1 and IPS2 and the macaque LIP region have been noted. Attentional modulation of visual responses in anterior and dorsal parts of this cluster, together with the previously discussed cluster (Fig. 7, magenta), may possibly have been referred to as pos IPS, whereas attentional modulation of visual responses in more ventral and posterior parts, together with those in V7, may have been referred to as vIPS (Corbetta et al., 2000; Kincade et al., 2005). Despite the resemblances between response patterns reported in the human in the yellow and magenta clusters (Fig. 7) and the LIP region in the macaque, it is important to note that there are also important differences. For example, on the basis of a comparison of fMRI data from both species, it has been argued that the representation in monkey LIP is more biased to the contralateral hemifield than is the case in the human (Patel et al., 2010).
We used a resting-state functional connectivity analysis to determine the brain areas with which the SPL and IPS regions were interacting. An initial explorative resting-state functional connectivity analysis examined the correlation of parietal cortex voxels with four target regions in both macaque and human brains: the extrastriate visual motion complex (V5/MT), PMd, an inferior PMd region, and FEF. The first region was at the approximate position of area V5/MT in the posterior superior temporal sulcus (Dubner and Zeki, 1971), which, within the macaque parietal cortex, is known to be interconnected with LIP and VIP (Blatt et al., 1990; Lewis and Van Essen, 2000). Within the macaque, parietal cortex functional connectivity with V5 is particularly strong in the posterior lateral bank of the IPS and fundus of the IPS at the approximate positions of LIP and VIP (Saleem and Logothetis, 2006) (Fig. 8). In the human subjects, however, prominent functional connectivity with V5 was seen on the medial bank of the IPS (Fig. 9A). Two regions of correlation of activity in IPS were prominent: a small region in an anterior position near the first cluster [identified with 7PC and VIP (Fig. 7, red)] and a large one in a posterior position near the fourth and fifth clusters [identified with 7A and IPS3 (Fig. 7, magenta) and with IPS1/2 (Fig. 7, yellow)].
In the macaque, the pattern of functional connectivity with the inferior PMd region was complementary to the regions identified on the basis of their functional connectivity with V5/MT. Rather than predominantly implicating the LIP region on the lateral bank of the IPS, functional connectivity with inferior PMd was present in the medial bank of the IPS in the MIP region and adjacent SPL (Saleem and Logothetis, 2006). Such distinct patterns of functional interaction are reminiscent of the differences in structural connections of these regions; whereas LIP and VIP but not MIP are interconnected with V5, inferior PMd is more strongly interconnected with MIP than it is with LIP (Caminiti et al., 1996; Johnson et al., 1996; Matelli et al., 1998; Lewis and Van Essen, 2000). Similarly, there were apparent regional differences in the functional connectivity of V5/MT and the inferior PMd region in the human brain (Fig. 9).
We also compared the functional connectivity patterns associated with a more dorsal part of PMd and the FEF. In the macaque, there was evidence of functional interactions between PMd and a wide swath of parietal cortex that was centered on area 5 and MIP on the medial bank of the IPS but which extended into LIP and VIP on the lateral bank and fundus of the IPS (Saleem and Logothetis, 2006). In contrast, the functional interactions of FEF were most prominent in LIP. These results are comparable with the structural connections of PMd. PMd is known to be connected with parts of area 5 and extending into the adjacent parietal reach region including MIP (Caminiti et al., 1996; Johnson et al., 1996; Matelli et al., 1998; Lewis and Van Essen, 2000), but there is also evidence for interconnections between LIP and PMd (Lewis and Van Essen, 2000). In contrast, FEF is interconnected with LIP and VIP (Lewis and Van Essen, 2000). Within human parietal cortex functional connectivity with PMd was prominent with the rostral SPL and rostral medial bank of the IPS (Fig. 9C). Functional connectivity with the human FEF was focused on a large region in the anterior medial bank of the IPS [the area identified with 7PC and VIP (Fig. 7, red)] and a smaller region in the posterior medial bank of the IPS [the area identified with 7A and IPS3 (Fig. 7, magenta)].
As with the IPL, these results show that resting-state functional connectivity in SPL/IPS is found in areas that are known to be structurally interconnected and illustrate that patterns of functional connectivity are broadly comparable in the human and macaque. In the next stage of the analysis, the human resting-state functional connectivity was formally compared with the tractography-based parcellation described above. Resting-state functional connectivity was determined between the parietal clusters and V5, inferior PMd, PMd, and FEF (Fig. 10). Evidence for functional connectivity with V5/MT was found in the most posterior SPL cluster on the posterior medial bank of the IPS [identified with IPS1/2 (Fig. 7, yellow)], the adjacent fourth IPS region [identified with 7A and IPS3 (Fig. 7, magenta)], and the first IPS region [identified with 7PC and VIP (Fig. 7, red)]. Evidence for functional connectivity with PMd was found predominantly in three SPL/IPS clusters [identified with 5L (Fig. 7, blue), hIP3/IPS4 (Fig. 7, green), and 7PC/VIP (Fig. 7, red)], with some limited connectivity also with the IPL (Wise et al., 1997). Evidence for functional connectivity within inferior PMd also included the same three SPL/IPS clusters, but interactions with the hIP3/IPS4 cluster (Fig. 7, green) were more prominent than had been the case with the more superior PMd region. The difference was particularly apparent at higher thresholds. This pattern resembles the pattern in the macaque in which PMd is interconnected with a number of SPL regions but a more inferior PMd region is particularly interconnected with MIP (Matelli et al., 1998).
Functional connectivity with the FEF (Fig. 10) was most prominent in the anterior lateral SPL/IPS region identified with 7PC and VIP (Fig. 7, red). Evidence for functional connectivity between the more posterior lateral SPL/IPS region identified with 7A and IPS3 (Fig. 7, magenta) and the SPL region identified with 5L (Fig. 7, blue) was also present. In addition, it should be noted that there is some uncertainty regarding the location of FEF in the human brain. In the macaque, FEF is located on the rostral bank of the arcuate sulcus. In the human brain, fMRI studies often activate a region more caudal than would be expected on the basis of macaque data, whereas stimulation studies position human FEF just anterior to the superior precentral sulcus (Amiez and Petrides, 2009). The target area used in the present study is located at the junction of the superior frontal sulcus and the precentral sulcus.
Discussion
The present study characterized the anatomical and functional organization of human lateral parietal cortex using non-invasive MRI techniques and compared it with the organization of the macaque parietal cortex. First, we used diffusion tractography to parcellate the human parietal cortex into 10 subregions. Then resting-state fMRI was used to examine the functional connectivity patterns of these subregions. Finally, comparing the resting-state fMRI data in humans and macaques revealed fundamental similarities in lateral parietal organization and demonstrated that the resting-state functional connectivity of lateral parietal cortex reflected major known anatomical connections. There are fundamental similarities in lateral parietal organization in humans and macaques, but some features distinguished human parietal cortex. An anterior IPL region active during tool observation was identified (Fig. 2, red cluster), and a central IPL region was characterized by strong interactions with aPFC.
Diffusion-tractography approaches have identified component areas within the pre-supplementary motor area (pre-SMA)/SMA region (Johansen-Berg et al., 2004), lateral premotor and ventrolateral prefrontal cortex (Anwander et al., 2007; Klein et al., 2007; Tomassini et al., 2007; Schubotz et al., 2010), and cingulate cortex (Beckmann et al., 2009), but this is the first time such a large expanse of cortex has been parcellated. Previous attempts to establish relationships between human and macaque parietal cortex have focused on retinotopic mapping, identification of similarities in functional activation profiles, and postmortem cytoarchitectonic analysis (Table 1). Because it has different strengths and weaknesses to these methods, the diffusion-tractography parcellation and functional connectivity approach provides independent and convergent evidence about the key component subdivisions of human parietal cortex.
The IPL parcellation yielded five clusters, and in most cases it was possible to identify their relationships with cytoartchitectonically defined areas PFop, PFm, PGa, and PGp (Caspers et al., 2008). One cluster (Fig. 2, blue), however, corresponded to several cytoarchitectonically defined areas, PF, PFt, and hIP2 (Caspers et al., 2006; Choi et al., 2006), which were presumably grouped together because of similarities in their connections. The differences in anatomical connections that gave rise to parcellation should mean that each IPL cluster differs in the access it has to the information represented in other cortical areas. Activity related to grasping, effecting strategy and response changes, attentional shifts, and memory retrieval, were preferentially associated with different clusters (Fig. 2C). Notably, the location of the most anterior IPL cluster was similar to that of activity recorded during tool use observation that has only been observed in human IPL (Peeters et al., 2009).
It has been claimed that human IPL may contain unique regions not found in other primates (Husain and Nachev, 2007; Peeters et al., 2009). Resting-state functional connectivity, however, in the posterior and anterior IPL, with PH, dlPFC, and PMv respectively, was similar in both humans and macaques (Fig. 4) and probably reflects similar anatomical connections in both species. In contrast, the human central IPL region was notable in that, unlike macaque central IPL, its activity was correlated with that of aPFC. It is in this region of the IPL that calculation- and language-related activity is found (Simon et al., 2002, 2004), and activity changes occur as people weigh up and change response strategy (Daw et al., 2006; Boorman et al., 2009). The aPFC region is close to the frontal pole in which activity has been related to the highest levels of cognitive control when “branching” between different subtasks is required (Koechlin et al., 1999). The two areas are often co-active (Daw et al., 2006; Boorman et al., 2009). There is evidence that the frontal pole is especially large in humans (Semendeferi et al., 2001) and that its increase in size is attributable to its lateral expansion, in hominoids, into the approximate region we found to be connected to mid-IPL. Consistent with our results, Petrides and Pandya (2007) have reported no connections between frontal polar area 10 and parietal cortex in the macaque, although Rozzi et al. (2006) note some between mid-IPL and anterior lateral prefrontal cortex. Brain differences emerge during speciation in a number of ways, including as a result of new connections invading an area or as specialized divisions of an area becoming spatially separate (Krubitzer, 1995, 2007), and such changes may underlie the present findings. The central IPL region is an important one to be considered in any attempt to account for the distinctive expansion of parietal cortex seen in humans even when comparison with other hominids (Bruner, 2010).
Parcellation of SPL and posterior IPS resulted in five clusters. Once again, it was possible to identify correspondences between the clusters and cytoarchitectonically defined regions (Fig. 7C), including 7PC, 5L, hIP3, and 7A (Scheperjans et al., 2008b). Retinotopic and other mapping procedures have been used to identify sensorimotor regions in IPS and SPL; the diffusion-tractography approach provides convergent evidence for the anatomical separation of many of these areas, including VIP (Bremmer et al., 2001; Sereno and Huang, 2006), also referred to as IPS5 (Konen and Kastner, 2008), IPS3, and IPS4 (Sereno et al., 2001; Swisher et al., 2007). One ventromedial posterior IPS cluster (Fig. 7, yellow), however, covered a region associated with two retinotopic maps, IPS1 and IPS2 (Silver et al., 2005), which have activity related to visuospatial attention and oculomotor intention (Schluppeck et al., 2005, 2006). IPS1, IPS2, and IPS3 have all been proposed as possible homologs of macaque LIP. The assignment of IPS1 and IPS2 to one cluster and the assignment of IPS3 to another cluster suggests that IPS1 and IPS2 share a number of connections and predicts that their functions may be especially similar. There is some evidence that both IPS1 and IPS2 are distinguished by the lateralization of their activity (Szczepanski et al., 2010). A modern, detailed study of the cytoartchitecture of this region has yet to be published, but again the diffusion-tractography results suggest that it may differ from the adjacent 7A region that encompasses much of IPS3.
The patterns of functional connectivity associated with the different SPL/IPS parietal clusters clarify the relationship between the parietal cortex in the human and other primates and are suggestive of function. Macaque VIP is known to be connected with V5, PMd, PMv, and FEF (Luppino et al., 1999; Lewis and Van Essen, 2000), and the anterior lateral SPL/IPS putative 7PC/VIP/IPS5 region (Fig. 7, red) was found to be functionally connected with the same regions in the present study (Fig. 10). The anterior medial, putative 5L, region (Fig. 7, blue) interacted strongly with PMd. Connections between area 5 and PMd have also been documented in the macaque (Caminiti et al., 1996; Johnson et al., 1996; Matelli et al., 1998). The mid-lateral SPL/IPS region near the IPS4 map (Fig. 7, green) was found to be functionally connected to PMd, particularly inferior PMd, in a manner reminiscent of the connections known to exist between MIP and inferior PMd in macaque (Matelli et al., 1998). As in macaque MIP, such a pattern would be expected if the region were to be involved in the directional control of movement and complements evidence that a similar region is active during visually guided hand movements even in the absence of eye movements (Rushworth et al., 2001; Grefkes et al., 2004; Hinkley et al., 2009). The functional connectivity between the most posterior IPS region (Fig. 7, yellow), near IPS1 and IPS2 maps, and V5 suggested that it resembles the ventral division of macaque LIP (Lewis and Van Essen, 2000). There was evidence of functional connectivity between the fourth, posterior lateral SPL/IPS cluster (magenta, IPS3 location) and FEF, suggesting that it also bore resemblances to LIP. Human SPL/IPS organization differed from macaque parietal organization in that the areas engaged in functional interactions with both FEF and V5, important interactions for any area concerned with visuospatial attention, were most prominent on the medial rather than the lateral bank of the IPS (Grefkes and Fink, 2005).
In summary, we show the feasibility of using non-invasive MRI methods to establish the anatomical and functional subregions of the human parietal cortex and to compare these with the macaque. Although only tracer injection studies in animal models provide definitive evidence for existence of synaptic connections between brain regions, comparative MRI-based studies are important because they facilitate translation to the human brain. Similarly, although studies such as the present one cannot replace detailed cytoarchitectonic descriptions of brain areas, they indicate ways of exploiting such descriptions in in vivo investigations of human brain anatomy and function in both health and illness. The possibility that functional connectivity is mediated by polysynaptic pathways cannot be ruled out. Nevertheless, similar, previous functional connectivity studies have recently identified correspondences between medial parietal and adjacent posterior cingulate areas in humans and macaques (Margulies et al., 2009). The qualitative similarity in human and macaque IPL is notable, suggesting that the neurophysiology and anatomy of macaque IPL is likely to be an important guide for understanding human IPL (Caminiti et al., 2010). Nevertheless, within this framework, human specializations, especially interactions between mid-IPL and aPFC, can be discerned.
Footnotes
This work was supported by a Marie Curie Intra-European Fellowship within the 6th European Community Framework Programme (R.B.M.), the Wellcome Trust (M.P.N., M.G.B., H.J.-B.), the Medical Research Council (United Kingdom) (R.B.M., S.J., J.X.O., A.S.M., and M.F.S.R.), the Multiple Sclerosis Society (United Kingdom) (V.T.), the Italian Foundation for Multiple Sclerosis (V.T.), the Royal Academy of Engineering (K.L.M.), and the Economic and Social Research Council (United Kingdom) (K.L.M.). We thank Urs Schüffelgen for help with Caret analyses.
- Correspondence should be addressed to Rogier B. Mars, Department of Experimental Psychology, University of Oxford, Tinbergen Building, 9 South Parks Road, Oxford OX1 3UD, UK. rogier.mars{at}psy.ox.ac.uk