Neural network of social interaction observation in marmosets
Abstract
A crucial component of social cognition is to observe and understand the social interactions of other individuals. A promising nonhuman primate model for investigating the neural basis of social interaction observation is the common marmoset (Callithrix jacchus), a small New World primate that shares a rich social repertoire with humans. Here, we used functional magnetic resonance imaging acquired at 9.4 T to map the brain areas activated by social interaction observation in awake marmosets. We discovered a network of subcortical and cortical areas, predominately in the anterior lateral frontal and medial frontal cortex, that was specifically activated by social interaction observation. This network resembled that recently identified in Old World macaque monkeys. Our findings suggest that this network is largely conserved between New and Old World primates and support the use of marmosets for studying the neural basis of social cognition.
Introduction
Observing and interpreting social interactions of other individuals is crucial for our everyday life. We share this remarkable ability with other primates that can also readily recognize social interactions such as play, fight, and grooming of other group members (Ambady et al., 2000; Ebenau et al., 2019; Spelke and Kinzler, 2007). Monkeys use this information to infer social hierarchies, adjust their behavior, and find allies (Bergman et al., 2003; Cheney et al., 1986). In Old World macaque monkeys, Sliwa and Freiwald, 2017 recently identified areas in parietal and frontal cortex that were exclusively activated by the observation of social interactions of other macaques using functional magnetic resonance imaging (fMRI). Interestingly, this network in macaques resembles the human default mode network (Mars et al., 2012) and the theory of mind network (Gallagher and Frith, 2003). As these networks share anatomical and functional characteristics, it suggests that the social interaction network may play a key role in these two human networks and be a conserved feature across primates (Sliwa and Freiwald, 2017).
Although macaque monkeys are widely used to study the neural basis of higher cognitive function, they may not be an ideal nonhuman primate model for studying the neural basis of social cognition as laboratory studies of social behavior are challenging due to their size and aggression toward conspecifics. Different species of macaques display varying type of prosociality but no cooperative breeding (Balasubramaniam et al., 2018; Ballesta and Duhamel, 2015; Beisner et al., 2016; de Waal, 1986; Joly et al., 2017; McCowan et al., 2011; Thierry et al., 2000) and are thus still an important model for studying some features and behaviors shared with humans (evolutionarily closer to humans than marmosets, larger size of societies, hierarchical societies, social behaviors including alliances, retaliation, redirected aggression). However, pairing and matching more than two macaque monkeys is often difficult under laboratory conditions. For these reasons, the common marmoset (Callithrix jacchus), a small New World nonhuman primate species that lives in family groups and shares certain similarities with humans such as prosocial behavior, social cognition (Burkart and Finkenwirth, 2015; Miller et al., 2016), and cooperative breeding (Burkart et al., 2009; Caselli et al., 2018), has garnered interest as a powerful primate model for human social behavior (Huang et al., 2020; Miller et al., 2016; Saito, 2015; Yokoyama et al., 2013; Yokoyama and Onoe, 2015).
Here, we aimed to identify the neural circuits involved in social interaction observation in marmosets using ultra-high-field (9.4 T) fMRI in fully awake animals. Similar to Sliwa and Freiwald’s study (2017), we acquired fMRI data from marmosets while they viewed videos of two marmosets interacting with each other (social condition) or of two marmosets engaged in independent goal-directed behavior (nonsocial condition).
Results
In this study, we employed a block-design task that consisted of videos showing social interactions of two marmosets (social condition), two side-by-side videos showing a marmoset in each video interacting with its own environment (nonsocial condition), and phase-scrambled versions of these videos. In the social conditions, the actions performed in the videos included eating, playing, fighting, and grooming between two marmosets. In the nonsocial conditions, the actions performed in the videos included eating, exploring objects, looking around, and scratching. A more detailed description of these actions is provided in Tables 1 and 2. Three common marmosets were scanned, and eight runs per marmosets were used for group analysis (24 runs in total).
Figure 1 shows group maps for four conditions (social, social scrambled, nonsocial, and nonsocial scrambled conditions) on fiducial brain surface (Figure 1, right hemisphere; see Figure 1—figure supplement 1 for the left hemisphere and Figure 1—figure supplement 2 for individual analyses, cluster-defining threshold. p<0.05, α = 0.05 for all contrasts). Both social (Figure 1a) and nonsocial (Figure 1c) videos activated a large network including visual (V1, V2, V3, V4, V4T, V6, V6A, 19M, MT), temporal (FST, TE1, TE2, TE3, TEO, IPa), parietal (MIP, VIP, LIP, AIP, PGM), and frontal areas (45, 47L, 46v, 10, 8av, 8ad, 8C, 6DR). The scrambled versions (Figure 1b, d) showed strong activation in visual (V1, V2, V3, V4, V4T, MT) and temporal areas (FST, TEO, IPa). Activations were also found in frontal cortex (47L, 45, 8av). The scrambled version of the nonsocial condition (Figure 1d) also showed activation in S2E, 1/2, 3a, 3b. Overall, the social and nonsocial conditions elicited a similar network, substantially stronger than those elicited by the scrambled versions.
To identify areas that were more active for the observation of social interactions, we compared the social with the nonsocial condition (Figure 2; see Figure 2—figure supplement 1 for individual analyses, cluster-defining threshold p<0.05, α = 0.05 for all contrasts). This analysis revealed stronger activations in frontal, somatosensory, and cingulate cortex for the social condition. Significant differences were present in frontal areas 10, 9, 14c, 14r, 45, 46d, 47L, 8ad, 8b, 8av, 8C, 6av, 6bv, and 6DC. More posteriorly, somatosensory areas S2I, S2PV, S2PR, 1/2, 3a, and 3b, and motor areas 4ab and ProM were more active for the social than nonsocial condition. We also observed activations in parietal areas (MIP, VIP), rostral temporal areas (STR, TPO, AuCPB, AuRPB), and in cingulate areas (25, 23b, and 31). The largest differences were present in motor/premotor cortex at the border of areas 4ab and 6DC, at the border of somatosensory area 1/2 and premotor areas 6av and 6bc, and in area 25 and orbitofrontal area 14. Differences between social and single interaction conditions are also found at the subcortical level in the SC, Tha, Amy, and CeB. Stronger activations for the nonsocial condition were mainly located in visual (V1, V2, V3, V4, V6), cingulate area 24b, and temporal areas (TF, TE2, TEO).
Overall numerous activations looked very similar in both hemispheres but are less strong in the left hemisphere (Figure 2; see Figure 2—figure supplement 2 for temporal signal-to-noise ratio maps).
To better control for low-level visual features, we also reported the interaction effect [social vs. social scrambled] minus [nonsocial vs. nonsocial scrambled] (Figure 3; see Figure 3—figure supplement 1 for individual analyses, cluster-defining threshold p<0.05, α = 0.05 for all contrasts). Overall, this analysis revealed most of the activations found by the previous contrast (Figure 2) but with less spread activations (Figure 3). Significant interaction effect was present in frontal areas 10, 9, 14c, 45, 46d, 47L, 8ad, 8b, 8C, 6bv, and 6DC, and more posteriorly, in somatosensory areas S2I, S2PR, 1/2, 3a, and 3b, motor area 4ab, and rostral temporal areas (STR, TPO, AuCPB). We still observed activations in parietal areas (VIP), in cingulate area 25, and at the subcortical level in the SC, Tha, Amy, and CeB. We did not observe any activation in area 8aV, indicating that the differences are not primarily related to saccades.
By extracting the time series from 12 regions of interest defined by the Paxinos et al., 2011 atlas and based on the group activation maps, we assessed the differences of response magnitude between the four conditions (Figure 4). Some areas showed stronger activation for the social condition than for the nonsocial condition (paired t-tests, p<0.5; right area 25 p=0.0435; right area 10 p=0.0075; left area 10L p=0.0012; right area 14c p=0.0016; right area VIP p=0.0052; right area 6va p=0.027; right area 6vb p=0.021; right area 46d p=0.0186; left area 46d p=0.005; and right area 45 p=0.0291). These areas also showed stronger activations for the social condition than for its scrambled version (paired t-tests, p<0.05; right area 25 p=0.0355; right area 10 p=0.049; left area 10L p=0.011; right area 14c p=0.002; right area VIP p=0.0257; right area 6va p=0.038; right area 46d p=0.0074; and right area 45 p=2.8 10−6) except for right area 6vb and left 46d.
Other areas did not show significant difference between social and nonsocial conditions but stronger activations for the social conditions than for its scrambled version (paired t-tests, p<0.05; left area 25 p=0.027; left area 6va p=0.012; right area 8ad p=0.007; left area 8ad p=0.014; right area 8av p=6.9 10−6; left area 8av p=0.004; and left area 45 p=0.004).
Discussion
In this study, we aimed to identify the neural circuitry underlying the observation of social interactions in New World marmosets. To that end, we employed a block-design task that consisted of videos showing social interactions of two marmosets (social condition), two side-by-side videos showing a marmoset in each video interacting with its own environment (eating, foraging) (nonsocial condition), and phase-scrambled versions of these videos. Both the social and nonsocial conditions showed strong occipitotemporal activation when compared to their scrambled versions, consistent with the location of face and body patches. In addition, these comparisons revealed activations in premotor and prefrontal areas. Like in macaques (Sliwa and Freiwald, 2017), the comparison of the social with the nonsocial condition revealed a network of mainly frontal cortical areas that is specifically activated by the observation of social interactions.
The occipitotemporal activations in the intact vs. scrambled movies resemble face-selective patches previously found in marmoset fMRI studies using photos of marmoset faces and body parts (Hung et al., 2015) or videos of faces (Schaeffer et al., 2020). They include the occipital (V2/V3), anterior dorsal (rostral TE), middle dorsal (caudal TE), posterior dorsal (FST), and posterior ventral (V4/TEO) face patches. This was to be expected as marmoset faces were visible for most of the time during these conditions. We also confirmed a recent report of a frontal face patch located in the lateral frontal cortex (45/47L) in both the social and nonsocial condition (Schaeffer et al., 2020). These areas are thus recruited for face processing independent of social interactions.
Recently, Sliwa and Freiwald, 2017 identified an exclusively social interaction brain network in macaques using fMRI. This network included predominately frontal and cingulate areas (14, 24c, 25, 32, 9, 10, 8B, 8D, 46D, 46V, 47L, 47O, 44, 6VR [F5], 6DC [F2], ProM). These overlap well with the frontal areas that we found here for a similar comparison in marmosets. In addition, we observed strong right hemisphere activation at the border of somatosensory area 1/2 and premotor areas 6av and 6bc and in the primary motor cortex 4ab in the marmoset.
We also found activation of area PG (homologous with area 7a in macaque) in both social, social scrambled, and nonsocial condition but not during nonsocial scrambled version. In macaques (Sliwa and Freiwald, 2017), area PG/7a is strongly activated during social interaction observation but also showed, albeit weaker, activation during nonsocial interaction observation. In macaques, this area contains motor neurons coding goal-directed motor acts in macaques (Rozzi et al., 2008), which were present both in our social and nonsocial condition. The role of area PG is still poorly understood in marmosets, except that it seems to be a major functional hub of the marmoset brain (Ghahremani et al., 2017). Further studies need to be performed to fully unravel its function, in particular for social scrambled stimulation.
Most of the areas activated by the social vs. nonsocial comparison (and present in the interaction effect [social vs. social scrambled] minus [nonsocial vs. nonsocial scrambled]) are already known to form a strong structural network in marmosets and many have been implicated in various aspects of social cognition in humans and macaques. It should be noted that none of the oculomotor areas were activated by this contrast, precluding any confounds with eye movements. Eye movements, and particularly visual saccades, mostly activate the lateral intraparietal area (LIP) and the frontal eye field in primates. Our group has mapped the saccade circuity in marmosets with multimodal approaches (electrophysiology and fMRI) (Ghahremani et al., 2017; Johnston et al., 2019; Schaeffer et al., 2019b; Selvanayagam et al., 2019). Together, this work demonstrates that saccades can be evoked by electrical stimulation of areas 45, 8aV, 8C, and 6DR. These regions are also activated in task-based fMRI with a visuomotor task and are functionally connected to the superior colliculus and other saccade-related regions as shown with RS-fMRI. In our contrast between the social and nonsocial conditions, only a small portion of these areas was activated and only in the right hemisphere. In addition, in the interaction effect contrast, no activation was present in area 8aV and only weak activations were observed in the dorsal part of 8aD and 6DR, suggesting that the activations we observed in our contrast are not directly related to oculomotor areas and so preclude most of the confound with eye movements. Furthermore, the PSC analysis of areas 8aV and 8aD showed no significant difference between social and nonsocial conditions. However, stimulation of area 45 evokes small saccades (Selvanayagam et al., 2019), and we observed here activation of the right area 45 in both social and nonsocial conditions as well as more activation for the social than nonsocial condition. The social condition may have elicited more small saccades related to the exploration of the facial expression and the interaction itself compared to the nonsocial interaction, which may involve more large explorations. This point raises some important limitations of our current setup. The animals were not rewarded for fixation or staying engaged during the task, and the quality of the video eye tracking did not allow us to confidently evaluate the animal’s engagement/fixation across conditions and runs.
Area 10 is densely connected with areas 14 and 9, to the temporal pole (STR) and parabelt auditory cortex (Buckner and Margulies, 2019; Burman et al., 2011). Area 10 also receives projections from parahippocampal, ventral somatosensory, and premotor areas. It is a key node of a social cognition circuitry in humans (Simmons et al., 2010) and linked to areas integrating information about social attention, intentionality, and other’s body actions in space in macaques (Jellema et al., 2000; Nummenmaa and Calder, 2009; Oram and Perrett, 1996).
Area 10 is also connected with area 25 in marmosets (Burman et al., 2011) and macaques (Joyce and Barbas, 2018). Connections with the amygdala, ventral striatum, hypothalamus, and periaqueductal gray, as well as with the anterior insula in macaques suggest a role of area 25 in emotion, visceromotor functions, and memory (Alexander et al., 2020; Alexander et al., 2019; Joyce and Barbas, 2018; Öngür et al., 2003). Our videos for the social condition showed marmosets playing or grooming each other, which may have triggered emotions and activated memories in the animals.
Tracer studies of area 14 in marmosets (Burman and Rosa, 2009), and some tracing data of the orbitofrontal cortex, encompassing area 14, have shown connections with amygdala, thalamus, and parahippocampal areas (Roberts et al., 2007). Area 14 seems to play a crucial role for appetitive value, reward processing, and decision making in humans (Elliott et al., 2010; FitzGerald et al., 2009) and in marmosets (Roberts and Wallis, 2000; Roberts and Clarke, 2019; Stawicka et al., 2020). Lesions in this area have been associated with social and emotional disturbance in humans (Bechara et al., 1996). Rostral area 14 also provides a major input to area 9 and send outputs to area 10 (Burman et al., 2011; Burman et al., 2006).
Sliwa and Freiwald, 2017 also identified a parieto-premotor network involved in social interaction analysis. This network included the posterior part of area AIP and premotor area F5 that have been shown to contain mirror neurons in macaques (Cléry et al., 2015; Cléry and Hamed, 2018; Rizzolatti et al., 1996; Rizzolatti and Fogassi, 2014; Rizzolatti and Matelli, 2003; Rizzolatti and Rozzi, 2018). This literature, in addition to the study of Sliwa and Freiwald, 2017, suggests a role for this mirror system in recognizing and understanding other agent’s intention and action. However, in this macaque study (Sliwa and Freiwald, 2017), this AIP/F5 network also responded to physical interactions between objects (not tested here), suggesting a broader role for this network. Here, we found that these areas were activated by the observation of nonsocial (i.e., ‘agent-object interaction’) and social interactions (i.e., ‘agent-agent interaction’) in marmosets. The mirror neuron system is still poorly understood in marmosets; however, neurons with characteristic ‘mirror’ properties have been identified in the superior temporal sulcus and in area 6V in marmosets (Suzuki et al., 2015).
Interestingly, the contrast between the social and nonsocial condition also revealed strong activations of parietal areas MIP and VIP in addition to premotor area 6V (corresponding to areas F4 and F5 in macaques, see Matelli et al., 1985). A parieto-frontal network related to VIP and F4 has been suggested to contribute to the construction of both the representation of one own’s body and of the body of others in macaques (Cléry et al., 2015). Both areas VIP and F4 are strongly linked to peripersonal space representation, defense, avoidance behavior, and social cognition (Brozzoli et al., 2013; Cléry et al., 2018; Cléry et al., 2017; Ishida et al., 2010; for review, see Cléry et al., 2015; Cléry and Hamed, 2018). Our current study suggests that areas VIP and 6V also contribute to social interaction analysis in marmosets.
Even though the monkeys were simply watching other marmosets interacting with other marmosets or objects, we also found strong activations of somatosensory (1/2, 3a, 3b, SII) and motor areas (4ab, PGM) that were absent in macaque monkeys (Sliwa and Freiwald, 2017). Activations were localized in areas that are also activated by tactile stimulation of the face and arm (Cléry et al., 2020), body parts where tactile stimulation would often occur during social interactions (grooming, playing). Particularly SII (Cléry et al., 2020; Huffman and Krubitzer, 2001; Krubitzer and Kaas, 1990; H-x et al., 2002) was activated in both hemispheres. A role of SII in social interaction processing is consistent with its involvement in the observation of others’ actions, as demonstrated both by human (Keysers et al., 2004; Blakemore et al., 2005; Lee Masson et al., 2018; for reviews, see Keysers et al., 2010; Bretas et al., 2020) and macaque studies (Hihara et al., 2015; Raos et al., 2014).
Our data showed stronger activations in the right hemisphere compared to the left hemisphere, even if the overall pattern looked similar. The temporal signal-to-noise ratio (tSNR) map showed a slightly lower tSNR in the left hemisphere (Figure 2—figure supplement 2), which may have led to the difference of activation we observed. However, this difference may also be due to a real lateralization effect during social behaviors. Indeed, it has been shown that the right hemisphere has a relative advantage over the left hemisphere, particularly in identifying social stimuli, understanding the intentions of others, social relationships, and social interactions in humans (see review, Hecht, 2014). Consequently, our data may reflect a similar lateralization for social processing in marmosets.
An interesting finding with our current study is the deactivation observed in both the social and nonsocial condition that were not present during the scrambled versions. This deactivation may have been triggered by an expectation of hearing vocalizations that were not met, but the actions during nonsocial interaction were mainly associated with eating so there is not necessarily vocalization associated with this kind of action. However, the fact that we observed this deactivation in both interaction conditions and not in the scrambled conditions may be the result of the marmosets' interest and engagement during the videos. Looking at conspecifics interacting with each other or with their environment is more entertaining than blurred videos. Further studies need to be performed to test a more complex paradigm such as videos showing social interactions with or without vocalizations.
In summary, we have identified the neural circuitry for social interaction observation in marmosets that resembles the one identified in Old World macaque monkeys (Sliwa and Freiwald, 2017). As Old and New World primates separated ~40 million years ago, this suggests that neural circuits involved in social interaction observation are conserved between macaques and marmosets or may be the result of convergent evolution (they may have underwent similarities in selection pressure but also show differences in their social behavior; Preuss, 2019). Social interaction observation recruited multiple brain networks, at both the cortical and subcortical levels, likely involved in emotion and affective processing, memory, peripersonal space representation, and the mirror neuron system. Overall, the identification of a preserved cortical and subcortical circuitry involved in social interaction observation in marmosets provides the map for studying the neural processes underlying primate social cognition.
Materials and methods
All experimental methods described were performed in accordance with the guidelines of the Canadian Council on Animal Care policy on the care and use of experimental animals and an ethics protocol #2017-114 approved by the Animal Care Committee of the University of Western Ontario.
Subjects and experimental setup
Request a detailed protocolThree common marmosets (C. jacchus) completed this study (one female, two males). Their ages and weights at the time of the experiment were 30 months and 363 g (M1), 38 months and 410 g (M2), and 35 months and 310 g (M3), respectively. The animals were prepared for awake fMRI experiments by implanting an MRI-compatible head restraint/recording chamber (for details, see Johnston et al., 2018). Each animal was initially habituated to the MRI environment over the course of 3 weeks. This included the acclimatization to the animal holder (tube closed by neck and tail plates), the head-fixation system, and periodic sounds produced by the MRI sequences (for details, see Schaeffer et al., 2019a). The behavior of the animals was monitored throughout the training to assess the tolerance and well-being of the marmosets (Silva et al., 2011). Their progression was based on their performance on the behavioral scale defined by Silva et al., 2011.
For MRI acquisition, the animal was placed in the sphinx position in an animal holder consisting of an MRI-compatible restraint system (Schaeffer et al., 2019a). The animal was first restrained using neck and tail plates, then the head was restrained by the fixation of the head chamber to the five-channel receive coil. An MRI-compatible camera (model 12M-i, 60 Hz sampling rate, MRC Systems GmbH, Heidelberg, Germany) was positioned in front of the animal holder to allow monitoring of the animals during the acquisition sessions by a veterinary technician. However, this camera and the setup did not allow a proper recording of eye movements as the pupil signal was frequently lost. In the MR scanner, the animal faced a translucent plastic screen placed at the front of the bore (119 cm from the animal’s head) where visual stimuli were displayed via back-projection after being reflected off of a first-surface mirror. Images were projected with a SONY VPL-FE40 projector. The maximum visual angle from the center to the side of the screen was 6.5°. Stimuli were presented via Keynote (version 10.1, Apple Incorporated, CA) and launched in synchronization with MRI TTL pulse triggers via a custom-written program running on a Raspberry Pi (model 3 B+, Raspberry Pi Foundation, Cambridge, UK).
Task and stimuli
Request a detailed protocolFour task conditions were presented: social condition videos where two marmosets directly interacted with each other (playing, grooming …) and nonsocial videos where two marmosets interacted with their own environment (eating, foraging) in two separate videos but displayed at the same time (Figure 1a, c, right panel). We chose real-world videos to maximize cognitive engagement (Finn and Bandettini, 2020; Meer et al., 2020). Videos were edited into 12 s clips using custom video-editing software (iMovie, Apple Incorporated). Two video clips were displayed within a block (6 s each), and 20 different videos clips were used in total. The videos within the same run were different between the two repetitions of the block condition but played following the same order. Tables 1 and 2 describe which actions were performed for each social (Table 1) and nonsocial (Table 2) blocks. Scrambled versions of social interaction and nonsocial interaction videos (Figure 1b, d, right panel) were created by random rotation of the phase information using a custom program (MATLAB, The MathWorks, Matick, MA). To preserve motion and luminance components, the same random rotation matrix was used for each frame in the scrambled conditions. Videos were presented in the center of the screen (7° height × 13° width) during acquisition sequences.
These four pseudo-randomized conditions were each displayed twice (social condition, nonsocial condition, and scrambled versions of each) following a block design (Figure 5). These eight task condition blocks (12 s each) were separated by baseline blocks (18 s each) where a 0.36° circular black cue was displayed in the center of a gray background. Three different stimulus sets were used and counterbalanced within the same animal and between animals. A sequence lasted 253 s (4.3 min). As motion and luminance varied across videos, mainly due to where the action took place, we counterbalanced these differences across the selected videos, conditions, and order of presentation to avoid any link with the condition block itself. No reward has been provided to the monkeys during the scan session.
Scanning
Request a detailed protocolWe performed data acquisition using a 9.4 T, 31 cm horizontal-bore magnet (Varian/Agilent, Yarnton, UK) and Bruker BioSpec Avance III HD console with the software package Paravision-6 (Bruker BioSpin Corp, Billerica, MA) at the Centre for Functional and Metabolic Mapping at the University of Western Ontario. We used an in-house, custom-built integrated receive coil with five channels (Schaeffer et al., 2019a) paired with a custom-built high-performance 15-cm-diameter gradient coil with 400 mT/m maximum gradient strength (xMR, London, Canada; Peterson et al., 2018). We used an in-house quadrature birdcage coil (12 cm inner diameter) for the transmit coil.
For each animal, a T2-weighted structural image was acquired, to allow an anatomical registration, with the following parameters: repetition time (TR) = 5500 ms; echo time (TE) = 53 ms; field of view (FOV) = 51.2 × 51.2 mm; voxel size = 0.133 × 0.133 × 0.5 mm; number of slices = 42 (axial); bandwidth = 50 kHz; GRAPPA acceleration factor = 2. For functional imaging, gradient-echo-based, single-shot echo-planar images covering the whole brain were acquired over multiple sessions (TR = 1500 ms; TE = 15 ms; flip angle = 40°; FOV = 64 × 64 mm; matrix size = 128 × 128; voxel size = 0.5 mm isotropic; number of slices = 42 [axial]; bandwidth = 500 kHz; GRAPPA acceleration factor = 2 [anterior-posterior]).
Analysis
Monkey 1 performed 8 runs in four sessions, monkey 2 performed 11 runs in five sessions, and Monkey 3 performed 9 runs in four sessions. We selected eight runs per monkeys to have the same number of runs across subjects for the group analysis (the selection was based on the best engagement of the animal as assessed by visual monitoring during the sessions). Time series were preprocessed using AFNI (Cox, 1996), FSL (Smith et al., 2004), and ANTS software (Advanced Normalization Tools, Avants et al., 2011). Raw functional images were converted to NifTI format using dcm2niix (Li et al., 2016) and reoriented from the sphinx position using FSL. The images were then despiked (AFNI’s 3dDespike) and volume registered to the middle volume of each time series (AFNI’s 3dvolreg). The motion parameters from volume registration were stored for later use with nuisance regression. Images were smoothed by a 1.5 mm full-width at half-maximum Gaussian kernel to reduce noise (AFNI’s 3dmerge). An average functional image was then calculated for each session and linearly registered (FSL’s FLIRT) to each animal’s T2-weighted image – the 4D time-series data was carried over using this transformation matrix. Anatomical images were manually skull-stripped, and this mask was applied to the functional images in anatomical space. The T2-weighted images were then nonlinearly registered to the NIH marmoset brain atlas (Liu et al., 2018) using Advanced Normalization Tools and the resultant transformation matrices stored for later transformation (see below). The olfactory bulb was manually removed from the T2-weighted images of each animal prior to registration as it was not included in the template image.
The task timing was convolved to the hemodynamic response (using AFNI’s ‘BLOCK’ convolution), and a regressor was generated for each condition to be used in a regression analysis (AFNI’s 3dDeconvolve) for each run. All four conditions were entered into the same model, along with polynomial (N = 5) detrending regressors, bandpass regressors, and the motion parameters derived from the volume registration described above (these nuisance regressors correspond to the baseline variance for the statistical analysis model). The resultant regression coefficient maps were then registered to template space using the transformation matrices described above. The T-value maps for each monkey were then compared at the individual level via paired-sample t-tests and fixed effect analyses across the eight runs for each marmoset (AFNI’s 3dttest++) and at the group level via a three-way analysis of variance and mixed effect analysis (AFNI’s 3danova3) using eight runs per marmosets (for a total of 24 runs). To protect against false positives, a cluster-size correction derived from Monte Carlo simulations was applied to the resultant t-test maps for α = 0.05 (using AFNI’s 3dClustsim).
The contrast between the social and nonsocial condition as well as the interaction effect [social vs. social scrambled] minus [nonsocial vs. nonsocial scrambled] was compared to previous fMRI and electrophysiological studies performed in marmosets (Ghahremani et al., 2017; Johnston et al., 2019; Schaeffer et al., 2019b; Selvanayagam et al., 2019) to check if the observed activations may be related to eye movements.
The resultant T-value maps are displayed on coronal sections or fiducial maps obtained with Caret (Van Essen et al., 2001; http://www.nitrc.org/projects/caret/) using the NIH marmoset brain template (Liu et al., 2018). Coordinates are presented with respect to the anterior commissure. The labeling refers to the histology-based atlas of Paxinos et al., 2011 for the cortical regions and to the atlas of Liu et al., 2018 for the subcortical regions (Figure 1, Figure 1—figure supplement 1, Figure 2, Figure 3).
Twelve regions of interest, areas 25, 10, 14c, VIP, PG, AIP, 6va, 6vb, 46d, 45, 8av and 8ad, were defined based on the Paxinos et al., 2011 atlas to evaluate the difference in response magnitude between conditions. First, masks of these regions were generated using AFNI’s 3dcalc in the right or left hemispheres of the ‘NIH_Marmoset_atlas_Paxinos_cortex’. Then, time series were extracted for each condition from the resultant regression coefficient maps using AFNI’s 3dmaskave. The mean and standard errors were computed and paired t-tests were calculated using MATLAB (version 2018a) after the data were tested for normality (Kolmogorov–Smirnov test).
Data availability
The datasets generated during this study are available at https://github.com/JClery/Social_interaction_paper and a copy archived at https://archive.softwareheritage.org/swh:1:rev:c9381727d3448d97d4f4263c1bb34c68ca296968/.
References
-
Toward a histology of social behavior: judgmental accuracy from thin slices of the behavioral streamAdvances in Experimental Social Psychology 32:201–271.https://doi.org/10.1016/S0065-2601(00)80006-4
-
The influence of phylogeny, social style, and sociodemographic factors on macaque social network structureAmerican Journal of Primatology 80:e22727.https://doi.org/10.1002/ajp.22727
-
Social power, conflict policing, and the role of subordination signals in rhesus macaque societyAmerican Journal of Physical Anthropology 160:102–112.https://doi.org/10.1002/ajpa.22945
-
Secondary somatosensory cortex of primates: beyond body maps, toward conscious self-in-the-world mapsExperimental Brain Research 238:259–272.https://doi.org/10.1007/s00221-020-05727-9
-
Cooperative breeding and human cognitive evolutionEvolutionary Anthropology: Issues, News, and Reviews 18:175–186.https://doi.org/10.1002/evan.20222
-
Marmosets as model species in neuroscience and evolutionary anthropologyNeuroscience Research 93:8–19.https://doi.org/10.1016/j.neures.2014.09.003
-
Cytoarchitectonic subdivisions of the dorsolateral frontal cortex of the marmoset monkey (Callithrix jacchus), and their projections to dorsal visual AreasThe Journal of Comparative Neurology 495:149–172.https://doi.org/10.1002/cne.20837
-
Cortical input to the frontal pole of the marmoset monkeyCerebral Cortex 21:1712–1737.https://doi.org/10.1093/cercor/bhq239
-
Architectural subdivisions of medial and orbital frontal cortices in the marmoset monkey (Callithrix jacchus)The Journal of Comparative Neurology 514:11–29.https://doi.org/10.1002/cne.21976
-
The prediction of impact of a looming stimulus onto the body is subserved by multisensory integration mechanismsThe Journal of Neuroscience 37:10656–10670.https://doi.org/10.1523/JNEUROSCI.0610-17.2017
-
Whole brain mapping of somatosensory responses in awake marmosets investigated with ultra-high-field fMRIJournal of Neurophysiology 124:1900–1913.https://doi.org/10.1152/jn.00480.2020
-
Frontier of self and impact predictionFrontiers in Psychology 9:73.https://doi.org/10.3389/fpsyg.2018.01073
-
AFNI: software for analysis and visualization of functional magnetic resonance neuroimagesComputers and Biomedical Research 29:162–173.https://doi.org/10.1006/cbmr.1996.0014
-
The integration of dominance and social bonding in primatesThe Quarterly Review of Biology 61:459–479.https://doi.org/10.1086/415144
-
The role of human orbitofrontal cortex in value comparison for incommensurable objectsJournal of Neuroscience 29:8388–8395.https://doi.org/10.1523/JNEUROSCI.0717-09.2009
-
Functional imaging of 'theory of mind'Trends in Cognitive Sciences 7:77–83.https://doi.org/10.1016/S1364-6613(02)00025-6
-
Frontoparietal functional connectivity in the common marmosetCerebral Cortex 27:3890–3905.https://doi.org/10.1093/cercor/bhw198
-
Cortical and thalamic connections of the parietal ventral somatosensory area in marmoset monkeys (Callithrix jacchus)Journal of Comparative Neurology 443:168–182.https://doi.org/10.1002/cne.10113
-
Cerebral lateralization of pro- and anti-social tendenciesExperimental Neurobiology 23:1–27.https://doi.org/10.5607/en.2014.23.1.1
-
Functional mapping of face-selective regions in the extrastriate visual cortex of the marmosetJournal of Neuroscience 35:1160–1172.https://doi.org/10.1523/JNEUROSCI.2659-14.2015
-
Shared mapping of own and others' Bodies in visuotactile bimodal area of monkey parietal cortexJournal of Cognitive Neuroscience 22:83–96.https://doi.org/10.1162/jocn.2009.21185
-
Neural representation for the perception of the intentionality of actionsBrain and Cognition 44:280–302.https://doi.org/10.1006/brcg.2000.1231
-
Methods for chair restraint and training of the common marmoset on oculomotor tasksJournal of Neurophysiology 119:1636–1646.https://doi.org/10.1152/jn.00866.2017
-
Alpha oscillations modulate preparatory activity in marmoset area 8adThe Journal of Neuroscience 39:1855–1866.https://doi.org/10.1523/JNEUROSCI.2703-18.2019
-
Comparing physical and social cognitive skills in macaque species with different degrees of social toleranceProceedings of the Royal Society B: Biological Sciences 284:20162738.https://doi.org/10.1098/rspb.2016.2738
-
Cortical connections position primate area 25 as a keystone for interoception, emotion, and memoryThe Journal of Neuroscience 38:1677–1698.https://doi.org/10.1523/JNEUROSCI.2363-17.2017
-
Somatosensation in social perceptionNature Reviews Neuroscience 11:417–428.https://doi.org/10.1038/nrn2833
-
The organization and connections of somatosensory cortex in marmosetsThe Journal of Neuroscience 10:952–974.https://doi.org/10.1523/JNEUROSCI.10-03-00952.1990
-
The first step for neuroimaging data analysis: dicom to NIfTI conversionJournal of Neuroscience Methods 264:47–56.https://doi.org/10.1016/j.jneumeth.2016.03.001
-
On the relationship between the "default mode network" and the "social brain"Frontiers in Human Neuroscience 6:189.https://doi.org/10.3389/fnhum.2012.00189
-
Patterns of cytochrome oxidase activity in the frontal agranular cortex of the macaque monkeyBehavioural Brain Research 18:125–136.https://doi.org/10.1016/0166-4328(85)90068-3
-
Movie viewing elicits rich and reliable brain state dynamicsNature Communications 11:5004.https://doi.org/10.1038/s41467-020-18717-w
-
Neural mechanisms of social attentionTrends in Cognitive Sciences 13:135–143.https://doi.org/10.1016/j.tics.2008.12.006
-
Architectonic subdivision of the human orbital and medial prefrontal cortexThe Journal of Comparative Neurology 460:425–449.https://doi.org/10.1002/cne.10609
-
ConferenceDevelopment of a gradient and shim insert system for marmoset imaging at 9.4T26th Annual Meeting of International Society for Magnetic Resonance in Medicine.
-
Critique of pure marmosetBrain, Behavior and Evolution 93:92–107.https://doi.org/10.1159/000500500
-
Premotor cortex and the recognition of motor actionsCognitive Brain Research 3:131–141.https://doi.org/10.1016/0926-6410(95)00038-0
-
The mirror mechanism: recent findings and perspectivesPhilosophical Transactions of the Royal Society B: Biological Sciences 369:20130420.https://doi.org/10.1098/rstb.2013.0420
-
Two different streams form the dorsal visual system: anatomy and functionsExperimental Brain Research 153:146–157.https://doi.org/10.1007/s00221-003-1588-0
-
BookChapter 28 - The mirror mechanism in the parietal lobeIn: Vallar G, Coslett H. B, editors. Handbook of Clinical Neurology. Elsevier. pp. 555–573.https://doi.org/10.1016/B978-0-444-63622-5.00028-0
-
Forebrain connectivity of the prefrontal cortex in the marmoset monkey (Callithrix jacchus): An anterograde and retrograde tract-tracing studyThe Journal of Comparative Neurology 502:86–112.https://doi.org/10.1002/cne.21300
-
Functional organization of inferior parietal lobule convexity in the macaque monkey: electrophysiological characterization of motor, sensory and mirror responses and their correlation with cytoarchitectonic AreasEuropean Journal of Neuroscience 28:1569–1588.https://doi.org/10.1111/j.1460-9568.2008.06395.x
-
The marmoset as a model for the study of primate parental behaviorNeuroscience Research 93:99–109.https://doi.org/10.1016/j.neures.2014.12.011
-
Face selective patches in marmoset frontal cortexNature Communications 11:4856.https://doi.org/10.1038/s41467-020-18692-2
-
Functional localization of the frontal eye fields in the common marmoset using microstimulationThe Journal of Neuroscience 39:9197–9206.https://doi.org/10.1523/JNEUROSCI.1786-19.2019
-
BookLongitudinal Functional Magnetic Resonance Imaging in Animal ModelsIn: Modo M, Bulte J. W. M, editors. Magnetic Resonance Neuroimaging: Methods and Protocols, Methods in Molecular Biology. Humana Press. pp. 281–302.https://doi.org/10.1007/978-1-61737-992-5
-
Mirror neurons in a new world monkey, common marmosetFrontiers in Neuroscience 9:459.https://doi.org/10.3389/fnins.2015.00459
-
An integrated software suite for surface-based analyses of cerebral cortexJournal of the American Medical Informatics Association 8:443–459.https://doi.org/10.1136/jamia.2001.0080443
-
Positron emission tomography imaging of the social brain of common marmosetsNeuroscience Research 93:82–90.https://doi.org/10.1016/j.neures.2014.12.006
Article and author information
Author details
Funding
Canadian Institutes of Health Research (FRN 148365)
- Stefan Everling
Canada First Research Excellence Fund (BrainsCAN)
- Stefan Everling
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Support was provided by the Canadian Institutes of Health Research (FRN 148365) and the Canada First Research Excellence Fund to BrainsCAN. We thank Miranda Bellyou, Cheryl Vander Tuin, and Hannah Pettypiece for animal preparation and care and Dr. Alex Li for scanning assistance.
Ethics
Animal experimentation: All experimental methods described were performed in accordance with the guidelines of the Canadian Council on Animal Care policy on the care and use of experimental animals and an ethics protocol #2017-114 approved by the Animal Care Committee of the University of Western Ontario. Animals were monitoring during the acquisition sessions by a veterinary technician.
Copyright
© 2021, Cléry et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,081
- views
-
- 307
- downloads
-
- 28
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Motivation depends on dopamine, but might be modulated by acetylcholine which influences dopamine release in the striatum, and amplifies motivation in animal studies. A corresponding effect in humans would be important clinically, since anticholinergic drugs are frequently used in Parkinson’s disease, a condition that can also disrupt motivation. Reward and dopamine make us more ready to respond, as indexed by reaction times (RT), and move faster, sometimes termed vigour. These effects may be controlled by preparatory processes that can be tracked using electroencephalography (EEG). We measured vigour in a placebo-controlled, double-blinded study of trihexyphenidyl (THP), a muscarinic antagonist, with an incentivised eye movement task and EEG. Participants responded faster and with greater vigour when incentives were high, but THP blunted these motivational effects, suggesting that muscarinic receptors facilitate invigoration by reward. Preparatory EEG build-up (contingent negative variation [CNV]) was strengthened by high incentives and by muscarinic blockade, although THP reduced the incentive effect. The amplitude of preparatory activity predicted both vigour and RT, although over distinct scalp regions; frontal activity predicted vigour, whereas a larger, earlier, central component predicted RT. The incentivisation of RT was partly mediated by the CNV, though vigour was not. Moreover, the CNV mediated the drug’s effect on dampening incentives, suggesting that muscarinic receptors underlie the motivational influence on this preparatory activity. Taken together, these findings show that a muscarinic blocker impairs motivated action in healthy people, and that medial frontal preparatory neural activity mediates this for RT.
-
- Medicine
- Neuroscience
The advent of midazolam holds profound implications for modern clinical practice. The hypnotic and sedative effects of midazolam afford it broad clinical applicability. However, the specific mechanisms underlying the modulation of altered consciousness by midazolam remain elusive. Herein, using pharmacology, optogenetics, chemogenetics, fiber photometry, and gene knockdown, this in vivo research revealed the role of locus coeruleus (LC)-ventrolateral preoptic nucleus noradrenergic neural circuit in regulating midazolam-induced altered consciousness. This effect was mediated by α1 adrenergic receptors. Moreover, gamma-aminobutyric acid receptor type A (GABAA-R) represents a mechanistically crucial binding site in the LC for midazolam. These findings will provide novel insights into the neural circuit mechanisms underlying the recovery of consciousness after midazolam administration and will help guide the timing of clinical dosing and propose effective intervention targets for timely recovery from midazolam-induced loss of consciousness.