Face Patch Resting State Networks Link Face Processing to Social Cognition

Faces transmit a wealth of social information. How this information is exchanged between face-processing centers and brain areas supporting social cognition remains largely unclear. Here we identify these routes using resting state functional magnetic resonance imaging in macaque monkeys. We find that face areas functionally connect to specific regions within frontal, temporal, and parietal cortices, as well as subcortical structures supporting emotive, mnemonic, and cognitive functions. This establishes the existence of an extended face-recognition system in the macaque. Furthermore, the face patch resting state networks and the default mode network in monkeys show a pattern of overlap akin to that between the social brain and the default mode network in humans: this overlap specifically includes the posterior superior temporal sulcus, medial parietal, and dorsomedial prefrontal cortex, areas supporting high-level social cognition in humans. Together, these results reveal the embedding of face areas into larger brain networks and suggest that the resting state networks of the face patch system offer a new, easily accessible venue into the functional organization of the social brain and into the evolution of possibly uniquely human social skills.


Introduction
Primates are highly social animals who cope with the challenges posed by life in complex social groups through sophisticated mechanisms for the recognition, evaluation, and generation of social signals.To understand the neural circuits mediating primate social behavior, here we take a novel, bottom-up approach utilizing a particularly well-defined sensory circuit as our starting point: the neural machinery that processes faces.Faces transmit rich information about socially relevant dimensions such as personal identity, emotional expressions, and gaze direction [1].To extract this multidimensional facial information, primates have evolved specialized brain areas [2], which are tightly and specifically interconnected [3], to form a faceprocessing system.Yet how this face-processing system is embedded, and thus how it makes face information available to other systems supporting social and cognitive functions, is largely unknown.
The connections between face processing and cognition are important to understand because behavioral and developmental studies show that faces occupy a special status among other, socially less relevant objects.Faces selectively draw spatial attention [4] and attract saccades much faster than other objects do [5], indicating privileged routing of facial information into attentional and eye movement control systems.Faces also drive specific mnemonic, emotional, and communicative responses [6], again suggesting specialized circuitry linking face areas to recipient, non-face-selective areas elsewhere in the brain.
To reveal with which parts of the brain the face-processing system can exchange information, we used resting state functional magnetic resonance imaging (rsfMRI) seeded in functionally defined face areas of the macaque monkey, the main animal model for face processing.rsfMRI noninvasively measures functional connectivity (FC) between brain areas based on spontaneous low frequency activity correlations, at high spatial resolution, and with full brain coverage [7].We focus on FC because although FC generally shows good agreement with anatomical connectivity [8], the set of potential functional connections between brain areas is far greater than that of direct structural links, as it is not constrained to monosynaptic connections but also includes dynamic, polysynaptic connectivity [9].A previous study of face patch connectivity using electrical microstimulation, a method that reveals primarily monosynaptic connectivity, found face patches to form a structurally closed system with few output connections [3], thus begging the question how the face patch system interacts with other systems.rsfMRI provides this type of complementary information and-due to its wide use in basic and clinical research in humans-confers the additional advantage of being readily comparable between species [10], thus enabling insight into the evolution of face-recognition systems [2,11,12].
Face area FC maps are not only essential for understanding the neural mechanisms of face recognition, they also provide a unique inroad into understanding high-level social cognition and its evolutionary heritage.This is because face areas are thought to constitute a major input into the so-called "social brain", a set of brain areas devoted to the processing of social interactions [13].In fact, almost 60% of the variance in our attitudes towards others can be explained by facial information alone [14].In humans, one of the social brain's core regions is the temporoparietal junction (TPJ).The TPJ is thought to be critical for high-level social cognition, in particular theory of mind (TOM) [15], the capacity to attribute mental states to ourselves and others.Apes and monkeys display basic forms of TOM such as understanding what others see or know [16,17].However, the very existence of a TPJ homolog in monkeys has been debated since the days of Brodmann, in part because the high-level functions that are supported by human TPJ, such as understanding others' false beliefs, may not be present in the macaque [15,18].This uncertainty is in large part due to the difficulties in studying macaque social cognition in a controlled experimental setup [19,20].Our approach sidesteps this issue and enables us to assess whether the kind of processing architecture enabling social cognition in humans already exists in the macaque: if this architecture was entirely absent, this would indicate that certain aspects of social cognition are indeed uniquely human.Conversely, if we could uncover similar brain networks in the macaque, this would suggest that at least a minimal scaffolding for high-level social cognition is already present in a primate whose evolutionary lineage split from ours some 25 million years ago [21].
To achieve this goal, we make use of the fact that the human social brain overlaps with another large scale network, the so-called default mode network (DMN) [22][23][24].The DMN, readily identifiable with rsfMRI in both humans and monkeys [25,26], comprises a set of interconnected areas more active during rest than task performance and is thus thought to generate the brain's default activity [27].Importantly, the overlap between the human DMN and the social brain includes high-level social cognition areas like TPJ, medial posterior parietal cortex (PPC), and dorsomedial prefrontal cortex (dmPFC) [23,28].We thus assessed, using the face patch resting state networks as a proxy for the social brain, whether and where a similar overlap exists in the macaque brain.This approach allows us to identify candidate homologs of human high-level social cognition brain areas.

Results
To determine face patch resting state networks (FPRSNs), we first identified face patches in six awake macaque monkeys using standard face localizers.We presented pictures of faces, bodies, and other object categories and contrasted activation during face presentations with activations during the presentation of nonface stimuli to reveal face areas.We identified one orbitofrontal (prefrontal orbital, PO) [29] and five temporal (middle lateral, ML; middle fundus, MF; anterior lateral, AL; anterior fundus, AF; anterior medial, AM) [2] face patches in all six animals.For subsequent analyses, whenever possible, bilateral homolog pairs were joined into one region of interest (ROI).Subsequently, the same monkeys underwent scans for rsfMRI during light isoflurane anesthesia.Aligning face localizer with rsfMRI scans allowed us to extract the time courses of spontaneous activity from each of the face patches and from regions outside of the face patch system.
As a first step towards characterizing FPRSNs, we determined connections within the face patch system.Our goals were to (i) reveal the hitherto unknown FC between temporal and orbitofrontal face patches, (ii) determine whether the established organizational principles of the temporal lobe face patch system, i.e., hierarchical and parallel organization, can be recovered from rsfMRI, and (iii) validate the intra-face patch connectivity pattern our noninvasive methodology reveals with data from an invasive approach, i.e., electrical microstimulation.To this end, we performed ROI-to-ROI correlation analyses after regressing out motion, heartbeat, and breathing artifacts from the data.We found FC between all face patches, with average Pearson correlation coefficients ranging between 0.34 for ML-MF and 0.02 for ML-AM (Fig 1, Wilcoxon signed rank tests, one-sided, corrected for multiple comparisons using the False Discovery Rate (FDR) at q = 0.05).The frontal patch PO, whose connectivity pattern was previously unknown, showed significant FC with the temporal patches MF, ML, and AM.AM, the face patch residing at the top of the temporal face-processing hierarchy [30], showed significant FC with PO, a likely output structure for AM, and with AF and AL, two input areas to AM located at the preceding level of the processing hierarchy, but not ML and MF, which are one level further removed.Thus, rsfMRI FC patterns recover the first organizational principle of the face-processing system, i.e., its hierarchical organization along the posterior-anterior axis of the superior temporal sulcus (STS).The second main organizational principle of the face patch system is parallelism: two processing streams reside in different cytoarchitectonic subdivisions of the STS, one in the fundus and the other in the lower bank of the STS.Hence, we tested whether the strength of FC followed known anatomical patterns of connectivity [31], i.e., whether face patches residing either in the fundus or on the lip of the STS (MF-AF and ML-AL) are more strongly correlated with each other, or whether their connectivity was equally strong across the fundus and the lip of the STS (AF-ML and AL-MF).Indeed, we found stronger FC within cytoarchitectonic subdivisions than across (median HL difference CI 95 = [0.0040.146]; Wilcoxon signed rank test, p = 0.03, two-sided), but this effect was mainly driven by the differential connectivity of ML.Anatomy also predicts a small but systematic bias for stronger interhemispheric connectivity between homolog than nonhomolog brain areas [32], and interhemispheric FC was in fact slightly higher between homolog face patches than between nonhomolog face patches (median HL difference CI 95 = [0.00370.2326]; Wilcoxon signed rank test, p < 0.02, one-sided).Finally, we compared the pattern of FC between the face patches from rsfMRI to previous results obtained with electrical microstimulation [3], which had revealed a highly specific set of connections between the STS face patches.As these specialized anatomical links provide the major scaffold for the functionally defined face processing network, we expected that rsfMRI should reveal a similar (although not identical) pattern of FC.A Spearman rank correlation between connection strengths indeed showed a high level of agreement between the two methods (r = 0.6314, p = 0.01; S1 Fig) .The same results were obtained using the robust correlation method Shepherd's Pi (r = 0.6352, p = 0.02) [33].This substantiates the noninvasive rsfMRI approach with results from a method known for its specificity in identifying direct neuronal connections [34].
The second main goal of this study was to determine the embedding of the face patches into the rest of the brain.We first focused on connections with other cortical areas.To this end, we aligned and brought all functional data into a common surface space, preserving high specificity for cortical grey matter and anatomical landmarks despite slight smoothing (1.25 mm kernel).We then performed a fixed effects (FFX) General Linear Model (GLM) group analysis for each of the face patches, respectively.After correction for multiple comparisons, these analyses revealed a highly convergent set of connected areas for AL, AF, MF, ML, and PO (and the posterior lateral face patch, PL, which was identified only in a subset of animals, S2 Face-selective neurons have been found in several subcortical structures such as the pulvinar [40] and the amygdala [41], which are not accessible to surface-based analyses.To determine FC of the cortical face patch system with subcortical structures, we performed a volumebased whole-brain analysis.After smoothing (2 mm Gaussian kernel), a FFX GLM revealed subcortical FC of the face patch system with a subregion of the claustrum, the amygdala, and the pulvinar, as had previously been shown using microstimulation [3], and additionally with the geniculate nucleus of the thalamus, the caudate nucleus, and the hippocampal formation (Fig 3).
It has been suggested that visual categories represented in spatially disjunct parts of ventral visual cortex are associated with unique patterns of connectivity [42].We thus tested which of the functional connections we observed were specific to the face patches.To this end, we isolated a patch in the anterior lip of the STS that responded to manmade objects during the localizer scans and then contrasted its FC to that of anatomically neighboring face patch AL (see Materials and Methods).AL showed stronger connectivity with the upper bank of the STS (including area TPO), the insula, lateral and medial parietal cortex (including area 7 and 23, respectively), lateral and medial prefrontal cortex (including areas F4, F5 and 6M, 9, respectively), as well as orbitofrontal cortex (area 13) (Fig 4).In contrast, the nearby object patch showed stronger connectivity with the inferotemporal cortex (including area TE) as well as the occipital cortex (including area VOT).Hence, faces and objects indeed display distinct patterns of FC.
Overall, the FC pattern of the face-patch system we identified includes areas like the lateral prefrontal cortex and the amygdala that contain face representations themselves, and others, like the premotor cortex, that likely do not.How strongly then does functional specialization for faces shape face patch FC? Specifically, we tested whether the strength of FC on the Cortical FC of the face patches.The central panel shows the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (corrected for multiple comparisons using cluster size thresholding at p < 0.05) on an inflated and flattened right hemisphere in F99 space.Highlighted are three broad networks of FC: areas in the prefrontal cortex (light blue), a premotor-parietal network (blue), and an occipitotemporal network including ventral stream areas (areal boundaries from Lewis & van Essen [35]).The surrounding panels show connectivity maps of the individual face patches (corrected for multiple comparisons using cluster size thresholding at p < 0.05), along with the areal boundaries of the conjunction analysis (light blue, blue, green).Representative locations of the face patches are outlined in black.See Table 1 for  whole-brain level depended on the selectivity of the target voxels, i.e., whether face patches were more connected to other face-selective voxels than to object-selective voxels.First, as a measure of selectivity, we computed d' between faces and objects for each voxel from the localizer data.Because rsfMRI connectivity falls off with distance, we also calculated the Euclidean distance from the voxel with peak selectivity within each respective face patch to the remaining voxels within the same hemisphere.We then matched voxels outside the face patch under consideration for selectivity and distance, weighing both factors equally (see Materials and Methods).Finally, we compared FC with matched face and object-selective voxels across 12 hemispheres, and found that for each face patch (AF, AL, MF, ML, PO), connection strength was higher for face than object voxels (mean differences: AF 0.04, AL 0.03, MF 0.04, ML 0.02, PO 0.03, paired t tests, all p < 0.03, one-sided).This shows that whole-brain rsfMRI can recover functional specificity within connection patterns, similar to what has been shown for structural connectivity in humans [43].
Can the neural systems that support the most advanced human sociocognitive skills be traced back to the macaque, a species with more limited sociocognitive abilities?In humans, social information processing networks, broadly defined, and the DMN overlap in the three cortical areas supporting the most high-level social-cognitive functions [22][23][24].Should the two networks intersect in macaques as well, and should this intersection occur at anatomical locations corresponding to those in humans, this would support a scenario of deep evolutionary heritage of these sociocognitive abilities.We first determined the DMN according to its original definition [25], i.e., by seeding rsfMRI in a bilateral ROI placed in medial PPC (areas 31/PGm, see Materials and Methods).As in previous studies [25,26,44], we observed a network comprising PPC, medial prefrontal, and lateral temporoparietal cortex (Fig 5a).We then computed conjunction maps between the FPRSNs and the DMN.Fig 5c shows that there is significant overlap in area TPO in the posterior STS where the human TPJ resides.This result replicates for all face patches (including AM, S10 Fig) .Furthermore, there is also overlap in areas 9M/10 in the dmPFC and, although less consistent in its precise location for each of the FPRSNs, areas PGm/23 in the medial PPC (Fig 5d), two further areas involved in high-level social cognition for which overlap with the DMN has been observed in humans [23].Because the overlap critically depends on the statistical threshold at which it is evaluated, we calculated Jaccard indices, which can be interpreted as percent overlap between two networks, over a wide range of uncorrected thresholds and compared the empirically observed degree of overlap with distributions of Jaccard indices obtained from the overlap of the FPRSNs with randomly generated maps that had the same spatial and statistical properties as the DMN.After correction for multiple comparisons, we found significant overlap between the DMN and all FPRSNs, including that of AM, until thresholds were so conservative that the likelihood of overlap was minimal (S11 Fig) .To assess whether the overlap with the DMN was specific to the FPRSNs, we also calculated the overlap between the object patch resting state network and the DMN.There was significantly more overlap between the FPRSN of AL and the DMN than between the object patch resting state network and the DMN at all thresholds (FDR-corrected, q = 0.01), and in fact, for most thresholds tested, there was no overlap between the object patch resting state network and the DMN at all (S12 Fig) .Furthermore, voxels which showed overlap between the FPRSN of AL and the DMN in posterior STS, dmPFC, and medial PPC were more strongly connected with AL than with the nearby object patch, while the opposite was the case for a region around the occipitotemporal sulcus, an area in which FPRSN and DMN also overlapped but that is not considered part of the monkey DMN (Fig 5, grey inset).Together, these In humans, the intersection between the DMN and the social brain isolates areas involved in high-level social cognition, such as TPJ in the posterior STS, dmPFC, and medial PPC [22][23][24].We tested, using the FPRSNs as a proxy for the social brain, whether and where a similar overlap exists in the macaque.(c) Voxels in area TPO in the dorsal bank of the posterior STS, located dorsally of FST, laterally to MST, and anterior and dorsal of MT, that show significant connectivity both with the PPC and with the respective face patch (AF, AL, MF, ML, and PO) at p < 10 −10 , uncorrected.The corresponding Jaccard Indices are 0.1358 for AF, 0.1635 for AL, 0.1786 for MF, 0.1366 for ML, and 0.1150 for PO.Also evident is overlap around the occipitotemporal sulcus, which is not part of the classical DMN.The grey inset shows that the strength of connectivity to voxels in which AL showed overlap with the DMN in the posterior STS was significantly higher for AL than for the object patch.In contrast, connectivity to the occipitotemporal sulcus was significantly higher for the object patch than for AL.results show that the macaque FPRSNs and DMN prominently and specifically overlap in brain areas that support high-level social cognition in the human brain.

Discussion
We find that the macaque face patches form a network linking highly face-selective regions in the temporal and orbitofrontal cortices.This face patch network is functionally embedded into a larger-scale, anatomically specific network of cortical and subcortical structures.This extended network significantly overlaps with the DMN, in particular in posterior STS, medial PPC, and dmPFC, which are involved in high-level social cognition in humans.The overlap is specific to the face patches, considering both the amount of overlap as well as the strength of FC in comparison to object resting state networks.Our results show the utility of a combined fMRI-rsfMRI approach in determining the embedding of functionally specific brain areas into larger-scale brain networks, and indicate that the face patch system may provide a unique inroad into understanding the complex organization of the social brain and a window into the evolution of primate social cognition.

An Extended Face-Processing System in the Macaque Monkey
Human face recognition has been proposed to rely on a "core system" consisting of interconnected face-selective areas and an "extended system" that utilizes inputs from the core system for cognitive, emotive, and mnemonic functions [6].Our whole-brain FC maps of the macaque face-processing system (Figs 2 & 3) show how the face patches are embedded into a largerscale network that shares, as we will discuss below, many of the properties of the proposed human extended system.In addition, these maps show how the face patches are nested into the general flow of information along the visual ventral stream.The ventral stream is organized along a main posterior-anterior axis, and along a dorsal-ventral axis with extensive lateral connectivity [45].This connectivity pattern can account for the coupling of face patches to occipital and temporal areas we observed and likely reflects the input-output relationships directly relevant for visual face-processing.A subcortical face-processing system has been proposed to exist, consisting of the superior colliculus, the pulvinar, and the amygdala [46].FC of the faceprocessing system with pulvinar and amygdala, as we found here, is compatible with this proposal of a separate, nonclassical set of subcortical inputs into the cortical face-processing network.The face patch system strongly interconnects with lateral prefrontal cortex, one of the main recipients of ventral stream output [45].Since connectivity from ventral stream face patches appears to include, but does not appear to be confined to, face specializations within the lateral prefrontal cortex [29,36] and is only partially specific to the face patches, facial information is likely made available for both face-domain-specific processing like face-specific working memory [36] and for domain-general cognitive processes in prefrontal cortex-like categorization [47] or attentional control [48].Further structures guiding spatial attention and eye movements are the supplementary eye field (SEF) in area F7, PPC, the pulvinar, the amygdala, and the caudate nucleus [49], all of which we found to be connected with the core face-processing system.These connections may aid in relaying information about the direction of attention of others extracted from a visual analysis of eyes and faces into the attentional system [50], and stronger connectivity of face patches than of nearby object patches with areas such as the SEF may underlie the behavioral advantages in directing saccades [5] and drawing spatial attention towards faces [4].Thus, the core face-processing system interfaces with attention and executive control systems through multiple functional routes.
Extensive further connectivity of the face patch system to executive systems beyond those for oculomotor control was evidenced by face patch-specific rsfMRI connectivity to several parts of premotor cortex.This came as a surprise to us, since it was not predicted by classical anatomy.Current anatomical evidence for direct connections between the parts of the STS that contain the face patches is limited to area F7 [51,52].Polysynaptic projections, however, from the STS to area F2 have recently been identified [53].The latter may arise from relay through the ventrolateral prefrontal cortex [53] or through well-established connectivity between parietal areas 7b and S2, which provide input to areas F4 and F5 [37].These premotor areas all contain visually responsive neurons and have been shown to be involved in the visual guidance of movements of face, eyes, and upper limbs; in particular, area F5 contains mirror neurons for socially relevant movements of the mouth such as lip smacks [54].Thus, connectivity between face and premotor areas may be stronger than previously thought and may support social communicative functions.
One of the main sets of functions proposed for the extended face-processing system lie within the emotional domain.We found the face patches to be connected with the amygdala, orbitofrontal cortex, and the insula, three structures implicated in the processing of emotions, which, in humans, have been shown to be involved in evaluating faces on social dimensions such as trustworthiness [55].The amygdala in particular also processes facial expression and gaze direction [56], two of the most important facial cues for social interactions.The orbitofrontal cortex, tightly interconnected with the amygdala, is thought to support the assignment of valence [57].Thus, functional links between the face patches and core structures of the emotional brain exist, which may serve the utilization of facial information for the generation of emotional responses.
The fourth main set of connections we observed linked the face patch system, and in particular AM, the area at the top of the face-processing hierarchy [30], to structures supporting long-term memory, notably entorhinal cortex, the hippocampal formation, and the claustrum.The latter has been shown to be heavily connected to the temporal lobe and has been hypothesized to act as a relay for sensory inputs to mediotemporal memory areas [58].Mediotemporal areas, including the hippocampus, contain face-selective cells thought to encode episodic memories [59].Thus, FC of the most anterior-ventral face area to mediotemporal lobe structures exists, possibly supporting the encoding and retrieval of memories of familiar individuals.Taken together, we find evidence for the core face-processing system to be functionally connected to areas that are known to support cognitive, emotional, communicative, and mnemonic functions.Macaques thus appear to possess an extended face-processing system as proposed for humans [6].

Comparison to rsfMRI Connectivity of the Human Face-Processing System
Neuroimaging studies investigating FC of human face-processing areas have found connectivity patterns that are broadly consistent with those we obtained in the macaque.In both species, correlations between the more posterior face-processing regions (the occipital face area (OFA) and the fusiform face area (FFA) in humans, and MF, ML, AF, and AL in the monkey) are stronger than with more anterior face-processing regions [60].FC outside the core face-processing system also displays a similar pattern in monkeys and humans, with overlap in the occipital, temporal, and frontal lobes [61,62], as well as subcortically, including the hippocampus, amygdala, caudate nucleus, and thalamus [63].Thus, internal and external FC of core face-processing areas is similar across the two primate species, consistent with the hypothesis of a deep evolutionary heritage of human face recognition abilities.The core face-processing systems in humans and monkeys are composed of multiple face areas.Establishing their homologies based on criteria like relative location and functional specialization [12] has proven difficult.Connectivity is a third, strong, and independent criterion for homology.The differences in face patch connectivity we found, in particular with mediotemporal lobe areas and the insula (S9 Fig), provide clear predictions for human combined fMRI/rsfMRI studies that localize all major face areas and map their FC.Currently available data on human face area connectivity, although not entirely consistent, points to a differentiation between dorsal and ventral face areas.Specifically, it has been found that the posterior STS displays stronger FC to premotor cortex than ventral face-processing areas, in particular the FFA [61,64] (but see [62]), in line with the finding that the STS face areas, but not the FFA, show structural connectivity to the ventrolateral prefrontal cortex [65].If this differentiation between dorsal and ventral face areas is confirmed, it would indicate homology of the entire macaque face-processing network with the human dorsal face areas, a quite radical view that has previously been put forward on other grounds [11].More data, in particular from the more variable human brain, will be needed to fully exploit the potential that these FC patterns hold in establishing homologies.Differential FC of the macaque face patches as we found hints at functional differentiations within the system and specialized roles these areas might play in social behavior.Comparing human and monkey rsfMRI connectivity networks to establish homologies as illustrated for the case of face areas is an approach that can be taken even further to understand the networks supporting complex cognitive functions.

Substrates of High-Level Social Cognition with a Deep Evolutionary Heritage
It has been a long-standing question whether certain high-level sociocognitive skills, e.g., the ability to reason about the contents of other persons' mental states, are uniquely human [15,18].For example, monkeys do not infer the belief that someone has about his/her own state of mind, something humans do routinely [66].There is behavioral evidence that monkeys display basic forms of TOM such as understanding what others, including human agents, see or know [17].However, even this interpretation remains contested, since such behavior may also arise from reasoning about the observable behavior of others without explicitly representing the others' mental state, i.e., without a TOM [20].Even more so, there is uncertainty about the neural basis that supports high-level social cognition.A case in point is human area TPJ, for which an old-world monkey homolog has been outrightly rejected [67], or proposed to reside either in parietal area 7a [68] or the posterior STS [69].While a whole battery of tasks is available to characterize high-level social cognition in humans, it has remained difficult to study the sociocognitive abilities of nonhuman primates experimentally.Here, we bypassed this issue and used a novel mapping strategy that utilizes the overlap of social brain areas and the DMN to identify putative homologs of human high-level social cognition systems.In humans, the main overlap between these two networks localizes to areas TPJ, PPC, and dmPFC.These areas have been strongly implicated in high-level social cognition, e.g., TOM in area TPJ and PPC, as well as the understanding of triadic interactions in the dmPFC [15,70,71].We now find that a similar overlap between FPRSNs and DMN exists in the macaque.The overlap is specific to FPRSNs and prominently includes area TPO in the dorsal posterior STS where the human TPJ resides.Little is known about the role of TPO in social cognition, but it shares several functional characteristics with human TPJ: Like TPJ, TPO is a polysensory area [72] that responds to biological motion [73] and action observation [74] and is involved in attention [75,76].Together with the connectivity overlap we observed, this suggests that a TPJ precursor or homolog exists in the macaque dorsal posterior STS and where it exists.Interestingly, this area is distinct from a more anterior STS region that has been shown to correlate with social network size [77].
Furthermore, we also find face-patch-specific overlap in areas 9M/10 in the dmPFC and areas PGm/23 in the medial PPC, which are well-established components of the human social brain and anatomically homologous in both species.Electrophysiological recordings during action-and error-monitoring of others [78,79] implicate the dmPFC in at least some aspects of social cognition in the monkey, which may form a precursor for the high-level sociocognitive functions that are supported by the dmPFC in humans [70].Furthermore, the human dmPFC is also involved in more basic forms of social processing that explicitly rely on facial information, such as gaze following [80], possibly a consequence of functional integration with face processing.The medial PPC, which has been linked to understanding social interactions [81], inferring other people's thoughts [82], and attributing mental states to others [83] in humans, has been shown to be active during action observation in monkeys [84], a basic ingredient for understanding the intentions of others.Thus, in addition to location and connectivity overlap, the functional properties of these areas are suggestive of a role in social cognition in the macaque.
Fig 5c also shows overlap around the occipitotemporal sulcus, including parts of areas V4, TFO, and TEO.In contrast to the other three regions, overlap in this area was not face-specific and connectivity to an STS object patch was even stronger than to a nearby face patch.This region has been found to be part of the human DMN [85] and is known to be anatomically connected both to the STS [45] and medial PCC [86,87], where our seed regions were located.Thus, connections exist that link this location to the DMN and face processing, rendering overlap in this region plausible.However, they do not suggest a specialized role of posteroventral cortex in social cognition.
Taken together, we find a pattern of overlap between the FPRSN and the DMN that includes the very areas that are selectively active in humans conducting the complex mental operations of TOM in macaque monkeys at rest.Hence, tapping into the social brain via an easily accessible sensory route uncovers similar regions as going through explicitly social cognition tasks in humans.This offers the exciting possibility to uncover the neural computations underlying sociocognitive operations.Since this overlap is present even under anesthesia, it is unlikely to reflect "mentalizing" as a default mode of processing [23], but it suggests a connectivity basis upon which a social default mode could have arisen.Our results thus point to a deep evolutionary heritage of a brain network composed of at least three areas for high-level social cognition.

Conclusions
We used fMRI and rsfMRI to noninvasively assess FC within the face patch system and the embedding of the face patches into larger brain networks.Our results demonstrate that faceprocessing areas interconnect with each other and with a set of nonface areas involved in cognitive, emotional, and mnemonic functions, forming an extended face-processing network.Importantly, we can also show that this extended network seeded in the face patches exhibits a similar overlap with the DMN as areas involved in high-level social cognition in humans, which allows us to localize a putative TPJ homolog in the macaque STS.This suggests that the face patch system offers an easily accessible, sensory venue into studying the social brain in monkeys, and thus into the evolution of possibly uniquely human social skills.

Subjects
All animal procedures met the National Institutes of Health Guide for Care and Use of Laboratory Animals, and were approved by the local Institutional Animal Care and Use Committees of The Rockefeller University (protocol number 12585-H) and Weill-Cornell Medical College (protocol number 2010-0029), where MR scanning was performed.Data were acquired in six male, pair-housed macaque monkeys (5 Macaca mulatta, 1 M. fascicularis, 5.4-7.3kg, age 3-5 yr).

Stimuli and Task
To localized face-selective ROIs, we used a standard face localizer [3].In short, subjects fixated on a white dot at the center of the screen while we presented images of human and/or monkey faces, human and/or monkey body parts and/or headless bodies, manmade objects, and fruits, intermixed with baseline periods in which only the fixation dot was shown in a block design.Each block lasted 24-30 s.Fluid reward was delivered after variable periods of time (2-4 s), during which the subject maintained fixation within 2 degrees of the fixation dot.Only runs in which the subjects reached at least 90% fixation stability were used for analyses.Visual stimulation and reward were controlled using in house software (Visiko, M. Borisov).Stimuli were projected on a back-projection screen using a video projector (NEC NP3250, refresh rate 60 Hz, resolution 1024 × 768 pixel) with a custom lens.Eye position was measured at 120 Hz using a commercial eye monitoring system (ISCAN).

Magnetic Resonance Imaging
Data were acquired on a 3 T scanner (Siemens TIM Trio).Functional data were acquired with an AC88 gradient insert (Siemens) and a custom 8-channel phased-array receive surface coil with a horizontally oriented single loop transmit coil (L.Wald, MGH/HST Martinos Center for Biomedical Imaging) while the monkeys were in sphinx position.Before scanning, the contrast agent ferumoxytol (8-10 mg of Fe per kg body weight) was injected into the femoral vein to increase the signal-to-noise ratio (SNR).For the face localizer experiments, we acquired between 16  ) acquired in interleaved order with phase partial Fourier 7/8, and two times generalized autocalibrating partially parallel acquisitions (GRAPPA) acceleration, covering the whole brain.Additionally, we obtained field maps that allowed subsequent EPI undistortion [89].For the resting state scans, we acquired 12 runs of 300 volumes of EPI data per animal, using the same sequences as in the localizer experiments.After induction with ketamine and dexmedetomidine hydrochloride, monkeys were lightly anesthetized with isoflurane (0.5%-0.6%) and placed in an MR-compatible monkey chair.The use of anesthesia follows the original definition of the monkey DMN [25] and conferred several technical advantages, including the elimination of motion artifacts and the ability to record cardiac and respiratory signals.Although anesthesia can affect systemic physiology, neural activity, vasoactive signal transmission and/or vascular reactivity [90], it has been shown that anesthesia preserves the correlation structure that is observed when the subjects are awake [91,92], that significant changes in correlation patterns occur only under much deeper levels of isoflurane anesthesia (>1.5%) [93] than the one we used (0.5%-0.6%), and that anesthetized monkey resting state networks, including the DMN, are strikingly similar to the same networks observed in awake humans [25,94].Electrocardiogram (sampling rate 400 Hz) and breathing rate (sampling rate 50 Hz) were acquired together with the imaging data.Anatomical images were obtained in a separate session using a T 1 -weighted magnetization-prepared rapid gradient echo (MPRAGE) sequence (FOV 128 mm, voxel size 0.5 × 0.5 × 0.5 mm, TR = 2.53 s, TE = 3.07 ms, ESP = 7.3 ms, BW = 190 Hz/Px, FA = 7°, 240 slices) and a custom 1-channel receive coil (L.Wald, MGH/HST Martinos Center for Biomedical Imaging) while the monkeys were anesthetized (isoflurane 1.5%-2%) and positioned in an MR-compatible stereotactic frame (Kopf Instruments).

Analyses
Data were analyzed in Freesurfer (v5.1, https://surfer.nmr.mgh.harvard.edu)and Matlab (R2011b, The Mathworks) using custom code.The first five volumes of each functional run were excluded to prevent T 1 saturation effects.Preprocessing included slice scan time correction, motion correction, and geometric distortion correction by means of a field map.Outliers in the time courses were detected semiautomatically based on a threshold of median absolute deviation = 3.5 [95] in the mean whole-brain time course and later excluded from analyses.To create inflated cortical surface reconstructions, the gray-white matter boundary in the skullstripped anatomical scans was segmented, reconstructed, smoothed, and inflated separately for each hemisphere [96,97].
Localization of face patches.In each animal and each hemisphere, we identified five temporal (ML, MF, AL, AF, AM) face patches and one frontal (PO) face patch based on the functional localizer, following established procedures [2].Data were slightly smoothed with a Gaussian kernel (2 mm full-width at half-maximum), whitened using a first-order autoregressive model, and detrended for first and second order polynomials.For each animal, we calculated a GLM with the stimulation conditions as predictors as well as six orthogonalized nuisance regressors accounting for motion artifacts.As in previous studies, face patches were identified based on anatomical location and relative position [2] in uncorrected significance maps (ranging between p < 0.05 and p < 10 −110 ) resulting from the contrasts (faces versus objects and bodies).Whenever possible, homolog face patches were then joined into bilateral ROIs for further processing.This way, we were able to identify six face patch ROIs in all six animals.
Analysis of resting state data.To analyze the FPRSNs, we aligned the functional data from the localizer scans with the resting state scans.We then extracted seed time courses in the resting state data from the face patch ROIs in volume space as the first Eigenvariate (using the mean time course instead gave almost identical results).To remove the effects of heartbeat and respiration [98], we created nuisance regressors from electrocardiographic (EKG) and breathing rate measurements acquired during the scans using the RETROICOR algorithm [99] as implemented in the PhLEM toolbox [100].Runs for which the EKG and/or respiration data were incomplete or contained many artifacts (e.g., because the recording devices detached during the measurement) were removed from the analyses (<10% of the data).Additionally, MRI data were highpass filtered at 0.0025 Hz and detrended for first and second order polynomials.Global mean regression was not performed.For ROI analyses, seed time courses were extracted from the unsmoothed data after regressing out the effects of motion, heartbeat, and respiration.Correlations between residuals were Fisher z-transformed per run, averaged, and back transformed to yield one r value per subject; this r value was then z-transformed and tested against 0 over subjects, using one-sided Wilcoxon signed rank tests (we used one-sided tests whenever we had a hypothesis about the directionality of the effects, and two-sided tests otherwise).Average r values over subjects were obtained by averaging and back transforming the per-subject z-values [101].
For median differences, we report 95% confidence for the Hodges-Lehman estimator, a rank-based, unbiased estimator of the median [102].Whole-brain analyses were performed in surface and in volume space.For surface-based analyses, functional data were aligned to an average anatomical surface template based on five subjects in the study (see below).Data was smoothed on the surface with a 1.25 mm kernel, thus preserving anatomical specificity for tissue type and anatomical location.To identify the cortical areas that show resting state connectivity with the face patches, we ran a FFX GLM for each hemisphere on the surface template with the respective face patch time course as a predictor, as well as nuisance regressors for motion (6), heartbeat (4), and respiration (10).To attenuate the effects of individual subjects on the group results, we used the same number of runs from each animal (8, i.e., 80 min of data).We then contrasted the respective face patch predictor against 0 to obtain p-value maps of connected voxels (two-tailed t test).Cluster size thresholding based on 10,000 permutations of voxel locations assuming a z-distribution was used to correct for multiple comparisons [103].The cluster-forming threshold was p < 0.00001 for AF, AL, MF, ML, and PO, and p < 0.05 for AM.AM was analyzed at a more lenient statistical threshold because it had a lower SNR but nevertheless showed a very similar pattern of results to the other face patches, including FC to lateral prefrontal, temporal, parietal, and occipital areas.Conjunction maps [104] were calculated from uncorrected maps of AF, AL, MF, ML, and PO and then corrected for multiple comparisons using cluster size thresholding (p < 0.05, cluster-forming threshold p < 0.001, 10,000 permutations).Data were deposited in the Dryad repository: http://dx.doi.org/10.5061/dryad.80476[39].For volume-based analyses, data were nonlinearly aligned to an average volume template (see below) using JIP (http://www.nmr.mgh.harvard.edu/~jbm/jip/)and smoothed using a 2 mm Gaussian kernel.The further analysis steps were the same as in surface space.
To assess the specificity of face patch FC, we first identified voxels activated by manmade objects in the localizer scan from the contrast (objects versus faces, bodies, and fruits), analogously to how we identified the face patches.We isolated an object-selective patch in the anterior lip of the STS in each animal and hemisphere, approximately 4 mm posterior to the face patch AL, as previously described [105].Before extracting seed time courses in the resting state data from the object patch ROIs as the first Eigenvariate, we matched the number of voxels in the object patch to the number of voxels in AL for each animal and hemisphere, respectively, and then joined homolog object ROIs into one bilateral ROI per animal.Next, to avoid issues of multi-colinearity, we separately orthogonalized the time course of the object patch ROI and the time course of AL against the global brain-wide mean, using the Gram-Schmidt algorithm [106].This preserves variance that is specific to the respective ROI and effectively reduced the correlation between face patch and object patch predictors.We then ran a FFX GLM for each hemisphere in surface space with the AL and the object patch time courses as predictors, as well as nuisance regressors for motion (6), heartbeat (4), and respiration (10).Finally, we contrasted the respective face patch predictor against the object patch predictor.The resulting maps of differential face patch versus object patch connectivity were corrected for multiple comparisons using cluster size thresholding (p < 0.05, cluster-forming threshold p < 10 −2 , 10,000 permutations).To further assess the relationship between face selectivity and connectivity, we also compared the connectivity of each face patch to face-and object-selective voxels in volume space, at matched levels of category selectivity.To this end, we first calculated d' between faces and objects as a measure of selectivity for each voxel in each hemisphere from the localizer data sets.d' was defined as q where M(f) and M(nf) are the mean to the face and nonface category with the highest mean response, respectively, and σ 2 (f) and σ 2 (nf) are their variances, taken from the first level GLM.Because rsfMRI connectivity is known to fall off with distance [107,108], we also calculated the Euclidean distance from the voxel with peak selectivity within each respective face patch to the remaining voxels within the same hemisphere.Euclidean distance is a good approximation to more refined distance measures, e.g., measures based on tractography [109].Subsequently, for each face-selective voxel (d' ! 2) outside the face patch under consideration (AF, AL, MF, ML, PO), we searched for the object-selective voxel that was the closest match both in its degree of selectivity as well as in distance from the face patch.To this end, we first considered the differences between each face-selective voxel and all object-selective voxels in selectivity and distance, respectively, as XY coordinates; we then found the coordinates that were closest to the origin (0,0-no difference); the object-selective voxel with the shortest Euclidean distance (<0.3) in this hypothetical XY space was chosen for the pair.This procedure assured that pairs of face-and object-selective voxels were matched as closely as possible, weighting selectivity and distance equally.We confirmed that across hemispheres and face patches, this led to no significant difference between face and object selectivity for 50 out of 54 patch × hemisphere combinations, and to no significant difference in distance for 51 out of 54 patch × hemisphere combinations, using two-sided t tests at p < 0.05.Finally, we compared the connectivity (z-transformed correlation coefficients, corrected for multiple comparisons using the FDR [110] at q = 0.01) for matched face and object-selective voxels across 12 hemispheres for AF, AL, MF, ML, and six hemispheres for PO using paired, one-sided t tests.Thus, only differences in connectivity strength for matched pairs entered the analysis.These analyses were done in native volume space on slightly smoothed (2 mm Gaussian kernel) data.To assess the overlap with the DMN, we placed an additional seed in PPC (anatomically defined area 31/ PGm in the left and right hemisphere, 46 voxels) of each monkey, following the original definition of the monkey DMN [25] and previous studies assessing the overlap between the DMN and the social brain in humans (e.g., [28]).DMN resting state connectivity maps were then calculated as described above.To further quantify the overlap between each FPRSN and the DMN, we calculated the Jaccard Index (# of voxels intersection / # of voxels union) [111] for each pair of networks over a range of statistical thresholds (p < 0.05, 10 −3:-1:-90 ).Statistical significance of the overlap was assessed by generating z-distributed random noise fields which had the same smoothness as the DMN inside a brain mask; at each threshold, we created 5,000 such maps per face patch, thresholded them to the same number of significant voxels as the original DMN map, and computed the overlap of each noise map with the respective face patch resting state map.We then compared the empirically observed Jaccard Index to the thus obtained distribution of Jaccard Indices arising by chance, which yielded a p-value.To correct for multiple comparisons (face patches × thresholds), we used a FDR at q = 0.01.Finally, conjunction maps were computed for illustration purposes.To assess the specificity of overlap between FPRSN and the DMN, we determined the overlap between the resting state network of the object patch we identified in the anterior lip of the STS and the DMN.To assess whether there was a difference in the amount of overlap with the DMN, we then compared the ratio of overlapping to nonoverlapping DMN voxels between the AL FPRSN and the object patch resting state network using a χ 2 -test for the same range of p-value thresholds as before (p < 0.05, 10 −3:-1:-90 ).Finally, we also compared the strength of connectivity to the object patch in voxels for which AL showed significant overlap with the DMN, contrasting AL with the object patch in a FFX GLM.The resulting maps were corrected for multiple comparisons (the number of voxels in which AL showed overlap with the DMN) using a FDR at q = 0.05.All overlap analyses were done in volume space after alignment to the average volume template.
Generation of anatomical templates for group analyses.order to conduct group statistics over the whole brain, functional data had to be brought into a common space.To this end, we created an anatomical surface template for alignment and visualization of the results.This was done by iteratively aligning and averaging the inflated surfaces of the five M. mulatta in the study using Freesurfer [112].This average surface was also brought into F99 space [113] using landmark surface-based registration [114] in Caret (v5.65; http://www.nitrc.org/projects/caret/) for use with surface-based atlases.Additionally, we created a volume template in AFNI (http://afni.nimh.nih.gov) for volume-based group analyses.This template was aligned to the MNI-Paxinos macaque template brain [38,115].

Supporting Information
S1 Fig. Correlation between rsfMRI and microstimulation connectivity.Spearman rank correlations revealed a strong correlation (r = 0.5256, p = 0.0365, red) between the connection strengths obtained through noninvasive rsfMRI and electrical microstimulation [3].After the exclusion of one extreme data point (AL-AF), the rank correlation rose to r = 0.6314, p = 0.0116 (blue).The same results were obtained using the robust correlation method Shepherd's Pi (r = 0.6352, p = 0.0219) [33].To bring rsfMRI and microstimulation connectivity onto the same scale, we first computed the median connectivity strength across animals per connection for the microstimulation data and then computed tied ranks within the columns of this and the resting state connectivity matrices.(PDF) S2 Fig. Cortical FC of the face patch PL.PL, which could be identified in four out of six animals in the study, showed a connectivity pattern (multiple comparisons corrected using cluster size thresholding at p < 0.05) that was similar to that of the other face patches, including connectivity to prefrontal, premotor, and occipitotemporal areas.Results are shown on inflated and flattened left and right hemispheres in F99 space.Representative locations of the face patches are outlined in black.For comparison, three broad networks of connectivity from the main conjunction analysis (including AF, AL, MF, ML, and PO) are highlighted: areas in prefrontal cortex (light blue), a premotor-parietal network (blue), and an occipitotemporal network including ventral stream areas.Areal boundaries are from Lewis & van Essen [35].Data shown here are publicly available at the Dryad Digital Repository [39].(PDF) S3 Fig. Cortical FC of the face patches.The central panel shows the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (multiple comparisons corrected using cluster size thresholding at p < 0.05) on an inflated and flattened left hemisphere in F99 space.Highlighted are three broad networks of connectivity: areas in prefrontal cortex (light blue), a premotor-parietal network (blue), and an occipitotemporal network including ventral stream areas (green).The surrounding panels show connectivity maps of the individual face patches (multiple comparisons corrected using cluster size thresholding at p < 0.05), along with the areal boundaries of the conjunction analysis (light blue, blue, green).Representative locations of the face patches are outlined in black.Areal boundaries are from Lewis & van Essen [35].See S1  [35].Shown are the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (multiple comparisons corrected using cluster size thresholding at p < 0.05) on an Jaccard Indices, over a wide range of statistical thresholds.The blue lines show the empirically observed overlap, while the green lines show the average degree of overlap between the respective face patch map and 5,000 randomly generated noise maps with spatial smoothness and number of significant voxels matched to those of the DMN map at each threshold.b) After correction for multiple comparisons (FDR, q = 0.01), there was significant overlap between the DMN and each of the FPRSNs until thresholds were so conservative that the likelihood of overlap was minimized.(PDF) S12 Fig. Specificity of overlap between the DMN and the FPRSNs.We compared the overlap between the resting state networks of the face patch AL and the DMN to the overlap between the resting state networks of a nearby object patch and the DMN.The ratio of Jaccard Indices for AL-DMN to object patch-DMN rapidly increases as statistical thresholds get more conservative.For most thresholds tested, there was no overlap between the object patch resting state network and the DMN (orange area).For the remaining thresholds, χ 2 -tests of the ratio of overlapping to nonoverlapping DMN voxels between AL and the object patch showed that face patch connectivity overlap always exceeded object patch connectivity overlap (corrected for multiple comparisons using a FDR at q = 0.01).(PDF) S1 Table .Areas left hemisphere.Area names are from Lewis & van Essen [35].(DOC)

Fig 1 .
Fig 1. FC within the face patch system.Most face patches were functionally connected to each other, with the exception of the frontal face patch PO, which was only connected to AM, MF, and ML and AM, which was only connected to AF, AL, and PO (corrected for multiple comparisons using FDR at q = 0.05).The inset shows the location of the face patches of an example subject on the average surface of the left hemispheres of the rhesus macaques used in this study.Nonsignificant correlations (ns) are set to 0. doi:10.1371/journal.pbio.1002245.g001 Fig), as well as connections that were unique to individual face patches (Fig 2, S3-S9 Figs).The FC pattern common to all face patches can be summarized by a conjunction analysis using the minimum statistic from each of the five maps (Fig 2, center).Both temporal and orbitofrontal face patches were connected to the (i) lateral prefrontal cortex (Fig 2, light blue), including areas 12 and 46, where face-selective neurons have been located [36]; (ii) regions of premotor cortex (Fig 2, blue), including areas F2, F4, F5, and F7, involved in the visual guidance of movements [37]; (iii) inferior parietal areas, including areas 7a and 7b (Fig 2, blue); (iv) areas of the temporal lobe, including the lower bank, fundus, and upper bank of the STS and parts of area TE (Fig 2, green); and (v) early visual cortex, especially areas V3 and V4 (Fig 2, green).This pattern of results was highly consistent across hemispheres (S3 Fig).Within this common pattern of connectivity shared across face patches, we also found significant variation across face patches (S9 Fig): (i) connectivity to the insula was more prominent for AL, ML, and PO than for MF and AF, (ii) AL and AF connectivity extended more posteriorly on the dorsolateral surface towards the central sulcus than any of the other face patches, and (iii) only AM showed connectivity to medial temporal lobe structures (entorhinal cortex, perirhinal cortical areas 35 and 36).

Fig 2 .
Fig 2.Cortical FC of the face patches.The central panel shows the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (corrected for multiple comparisons using cluster size thresholding at p < 0.05) on an inflated and flattened right hemisphere in F99 space.Highlighted are three broad networks of FC: areas in the prefrontal cortex (light blue), a premotor-parietal network (blue), and an occipitotemporal network including ventral stream areas (areal boundaries from Lewis & van Essen[35]).The surrounding panels show connectivity maps of the individual face patches (corrected for multiple comparisons using cluster size thresholding at p < 0.05), along with the areal boundaries of the conjunction analysis (light blue, blue, green).Representative locations of the face patches are outlined in black.See Table1 fora list of area names, S4 Fig for complete areal labeling of face patch connectivity on the conjunction map based on Lewis & van Essen [35], S5 Fig for areal labels of the premotor-parietal network based on Paxinos et al. [38], S6 Fig for a quantification of how many face patches connected with a given vertex, S7 Fig for a conjunction analysis restricted to the temporal face patches (AF, AL, MF, ML), and S8 Fig for results of a corresponding conjunction analysis in volume space.Data shown here are publicly available at the Dryad Digital Repository [39].
Fig 2.Cortical FC of the face patches.The central panel shows the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (corrected for multiple comparisons using cluster size thresholding at p < 0.05) on an inflated and flattened right hemisphere in F99 space.Highlighted are three broad networks of FC: areas in the prefrontal cortex (light blue), a premotor-parietal network (blue), and an occipitotemporal network including ventral stream areas (areal boundaries from Lewis & van Essen[35]).The surrounding panels show connectivity maps of the individual face patches (corrected for multiple comparisons using cluster size thresholding at p < 0.05), along with the areal boundaries of the conjunction analysis (light blue, blue, green).Representative locations of the face patches are outlined in black.See Table1 fora list of area names, S4 Fig for complete areal labeling of face patch connectivity on the conjunction map based on Lewis & van Essen [35], S5 Fig for areal labels of the premotor-parietal network based on Paxinos et al. [38], S6 Fig for a quantification of how many face patches connected with a given vertex, S7 Fig for a conjunction analysis restricted to the temporal face patches (AF, AL, MF, ML), and S8 Fig for results of a corresponding conjunction analysis in volume space.Data shown here are publicly available at the Dryad Digital Repository [39].doi:10.1371/journal.pbio.1002245.g002

Fig 3 .Fig 4 .
Fig 3. Subcortical FC of the face patches.A conjunction analysis of the rsfMRI maps of face patches AF, AL, MF, ML, and PO revealed FC with several subcortical areas.(a) Pulvinar and hippocampal formation (FDR corrected at q = 0.025); (b) amygdala (FDR corrected at q = 0.005); (c) claustrum (FDR corrected at q = 0.025); (d) caudate nucleus (FDR corrected at q = 0.025); (e) lateral geniculate nucleus (FDR corrected at q = 0.005).Results are overlaid on coronal (a-d) and sagittal (e) slices of the MNI-Paxinos template brain, in radiological convention (left is right).Coordinates are relative to the center of the anterior commissure.Cortical results are masked for display purposes.Note different scaling of correlation coefficients in (c) and (d).doi:10.1371/journal.pbio.1002245.g003

Fig 5 .
Fig 5. Overlap analysis of the DMN and FPRSNs.(a) The classical monkey DMN on coronal and sagittal slices.The arrow indicates the approximate location of the DMN seed.Note that the sagittal slices are the same as in c and d, respectively, at corresponding statistical thresholds.(b) Illustration of the overlap analysis.In humans, the intersection between the DMN and the social brain isolates areas involved in high-level social cognition, such as TPJ in the posterior STS, dmPFC, and medial PPC[22][23][24].We tested, using the FPRSNs as a proxy for the social brain, whether and where a similar overlap exists in the macaque.(c) Voxels in area TPO in the dorsal bank of the posterior STS, located dorsally of FST, laterally to MST, and anterior and dorsal of MT, that show significant connectivity both with the PPC and with the respective face patch (AF, AL, MF, ML, and PO) at p < 10 −10 , uncorrected.The corresponding Jaccard Indices are 0.1358 for AF, 0.1635 for AL, 0.1786 for MF, 0.1366 for ML, and 0.1150 for PO.Also evident is overlap around the occipitotemporal sulcus, which is not part of the classical DMN.The grey inset shows that the strength of connectivity to voxels in which AL showed overlap with the DMN in the posterior STS was significantly higher for AL than for the object patch.In contrast, connectivity to the occipitotemporal sulcus was significantly higher for the object patch than for AL.(d) Voxels in dmPFC (areas 9M/10) and medial PPC (areas PGm/23) that show significant connectivity both with the PPC and with the respective face patch (AF, AL, MF, ML, and PO) at p < 10 −7.5 , uncorrected.The corresponding Jaccard Indices are 0.2006 for AF, 0.2268 for AL, 0.2530 for MF, 0.2062 for ML, and 0.1821 for PO.The grey inset shows that the strength of connectivity to voxels in which AL showed overlap with the DMN in the dmPFC (area 9M) and medial PPC (area 23) was significantly higher for AL than for the object patch.Both inset maps are corrected for multiple comparisons with a FDR at q = 0.05, accounting for the number of voxels that show significant overlap between AL and the DMN at the same statistical thresholds as shown in panels c and d, respectively.All results are overlaid on the MNI-Paxinos template brain, in radiological convention (left is right).Coordinates are relative to the center of the anterior commissure.Area labels are based on Paxinos et al.[38].See S10 Fig for overlap between the DMN and AM connectivity, which replicates the main findings at a more lenient statistical threshold.
Fig 5. Overlap analysis of the DMN and FPRSNs.(a) The classical monkey DMN on coronal and sagittal slices.The arrow indicates the approximate location of the DMN seed.Note that the sagittal slices are the same as in c and d, respectively, at corresponding statistical thresholds.(b) Illustration of the overlap analysis.In humans, the intersection between the DMN and the social brain isolates areas involved in high-level social cognition, such as TPJ in the posterior STS, dmPFC, and medial PPC[22][23][24].We tested, using the FPRSNs as a proxy for the social brain, whether and where a similar overlap exists in the macaque.(c) Voxels in area TPO in the dorsal bank of the posterior STS, located dorsally of FST, laterally to MST, and anterior and dorsal of MT, that show significant connectivity both with the PPC and with the respective face patch (AF, AL, MF, ML, and PO) at p < 10 −10 , uncorrected.The corresponding Jaccard Indices are 0.1358 for AF, 0.1635 for AL, 0.1786 for MF, 0.1366 for ML, and 0.1150 for PO.Also evident is overlap around the occipitotemporal sulcus, which is not part of the classical DMN.The grey inset shows that the strength of connectivity to voxels in which AL showed overlap with the DMN in the posterior STS was significantly higher for AL than for the object patch.In contrast, connectivity to the occipitotemporal sulcus was significantly higher for the object patch than for AL.(d) Voxels in dmPFC (areas 9M/10) and medial PPC (areas PGm/23) that show significant connectivity both with the PPC and with the respective face patch (AF, AL, MF, ML, and PO) at p < 10 −7.5 , uncorrected.The corresponding Jaccard Indices are 0.2006 for AF, 0.2268 for AL, 0.2530 for MF, 0.2062 for ML, and 0.1821 for PO.The grey inset shows that the strength of connectivity to voxels in which AL showed overlap with the DMN in the dmPFC (area 9M) and medial PPC (area 23) was significantly higher for AL than for the object patch.Both inset maps are corrected for multiple comparisons with a FDR at q = 0.05, accounting for the number of voxels that show significant overlap between AL and the DMN at the same statistical thresholds as shown in panels c and d, respectively.All results are overlaid on the MNI-Paxinos template brain, in radiological convention (left is right).Coordinates are relative to the center of the anterior commissure.Area labels are based on Paxinos et al. [38].See S10 Fig for overlap between the DMN and AM connectivity, which replicates the main findings at a more lenient statistical threshold.doi:10.1371/journal.pbio.1002245.g005 and 51 runs of functional (T 2 Ã -weighted) gradient-echo echoplanar imaging (EPI) data per animal.Each run consisted of 196 volumes of 54 horizontally oriented slices (field of view [FOV] 96 mm, voxel size 1 × 1 × 1 mm, repetition time [TR] = 2 s, echo time [TE] = 16 ms, echo spacing [ESP] = 0.63 ms, bandwidth [BW] = 1,860 Hz/Px, flip angle [FA] = 80°, no gap inflated and flattened right hemisphere in F99 space.Representative locations of the face patches outlined in black.Data shown here are publicly available at the Dryad Digital Repository[39].(PDF) S5 Fig.Alternativepartitioning scheme for the premotor-parietal FC based on Paxinos et al.[38].Shown are the results of a conjunction analysis of the maps from AF, AL, MF, ML, and PO (multiple comparisons corrected using cluster size thresholding at p < 0.05) on an inflated and flattened left hemisphere in F99 space.Representative locations of the face patches are outlined in black.Data shown here are publicly available at the Dryad Digital Repository[39].(PDF) S6 Fig. Number of face patches connected to a given vertex.For each vertex, we quantified how many face patches were functionally connected after correction for multiple comparisons using cluster size thresholding at p < 0.05 (cf.outer panels of Fig 2 & S3 Fig).Results are shown on inflated and flattened left and right hemispheres in F99 space.Representative locations of the face patches are outlined in black/white.Data shown here are publicly available at the Dryad Digital Repository [39].(PDF) S7 Fig. Cortical FC of the temporal face patches.Shown are the results of a conjunction analysis of the maps from AF, AL, MF, and ML (multiple comparisons corrected using cluster size thresholding at p < 0.05) on inflated and flattened left and right hemispheres in F99 space.Highlighted are three broad networks of connectivity from the main conjunction analysis including the orbitofrontal face patch PO: areas in prefrontal cortex (light blue), a premotorparietal network (blue), and an occipitotemporal network including ventral stream areas (green).Areal boundaries are from Lewis & van Essen [35].Representative locations of the face patches are outlined in black.Data shown here are publicly available at the Dryad Digital Repository [39].(PDF) S8 Fig. FC of the FPRSN in volume space.Shown are the results of a conjunction analysis (uncorrected) of the rsfMRI maps of bilateral face patches AF, AL, MF, ML, and PO, overlaid on coronal slices of the MNI-Paxinos template brain, in radiological convention (left is right).Coordinates are relative to the center of the anterior commissure.Area labels are based on Paxinos et al. [38].(PDF) S9 Fig. Differences in FC between the face patches.(a) connectivity to the insula was more prominent for AL, ML, and PO than for MF and AF; (b) AL and AF connectivity extended more posteriorly on the dorsolateral surface towards the central sulcus than any of the other face patches (the premotor-parietal network from the main conjunction analysis is shown for reference); (c) only AM showed connectivity to medial temporal lobe structures (entorhinal cortex, perirhinal cortical areas 35 and 36).All contrasts were corrected for multiple comparisons using a FDR at q = 0.05.Results are shown on an inflated and flattened right hemisphere in F99 space, negative differences are truncated for display purposes.Representative locations of the face patches are outlined in black.Areal boundaries are from Lewis & van Essen [35].Data shown here are publicly available at the Dryad Digital Repository [39].(PDF) S10 Fig. Overlap analysis of DMN and the FPRSN of AM.(a) Voxels in area TPO in the dorsal bank of the posterior STS that show significant connectivity both with the PPC and AM at p < 10 −3.5 , uncorrected.(c) Voxels in dmPFC and medial PPC that show significant connectivity both with the PPC AM at p < 10 −1.6 , uncorrected.The overlap between the resting state network of AM and the DMN was highly consistent with that of the other face patches, localizing to three areas known to support high-level social cognition in humans.Results are overlaid on the MNI-Paxinos template brain.Coordinates are relative to the center of the anterior commissure.(PDF) S11 Fig. Overlap between the DMN and the FPRSNs, Jaccard Indices.a) Results of the permutation tests for significant overlap between the DMN and individual FPRSNs, quantified as

Table 1 .
[35]s right hemisphere.Area names are from Lewis & van Essen[35]. doi:10.1371/journal.pbio.1002245.t001 Table for a list of area names.Data shown here are publicly available at the Dryad Digital S4 Fig. Complete areal labeling of face patch FC in the right hemisphere based on Lewis & van Essen