Oxytocin Manipulation Alters Neural Activity in Response to Social Stimuli in Eusocial Naked Mole-Rats

The social decision-making network (SDMN) is a conserved neural circuit that modulates a range of social behaviors via context-specific patterns of activation that may be controlled in part by oxytocinergic signaling. We have previously characterized oxytocin’s (OT) influence on prosociality in the naked mole-rat, a eusocial mammalian species, and its altered neural distribution between animals of differing social status. Here, we asked two questions: (1) do patterns of activation in the SDMN vary by social context and (2) is functional connectivity of the SDMN altered by OT manipulation? Adult subordinate naked mole-rats were exposed to one of three types of stimuli (three behavioral paradigms: familiar adult conspecific, unfamiliar adult conspecific, or familiar pups) while manipulating OT (three manipulations: saline, OT, or OT antagonist). Immediate early gene c-Fos activity was quantified using immunohistochemistry across SDMN regions. Network analyses indicated that the SDMN is conserved in naked mole-rats and functions in a context-dependent manner. Specific brain regions were recruited with each behavioral paradigm suggesting a role for the nucleus accumbens in social valence and sociosexual interaction, the prefrontal cortex in assessing/establishing social dominance, and the hippocampus in pup recognition. Furthermore, while OT manipulation was generally disruptive to coordinated neural activity, the specific effects were context-dependent supporting the hypothesis that oxytocinergic signaling promotes context appropriate social behaviors by modulating co-ordinated activity of the SDMN.


INTRODUCTION
The social decision-making network (SDMN) is a highly conserved interconnected group of brain regions controlling behaviors related to sex, social dominance, parenting, and affiliation across vertebrates Hofmann, 2011, 2012). Originally, the network was described by linking the social behavior network (adjacent tegmentum, AH, BNST, MeA, LS, PG, PO, and VMH) and the mesolimbic reward circuit (BLA, BNST, caudate, hippocampus, LS, NAcc, VP, and VTA) (Newman, 1999;O'Connell and Hofmann, 2011). A case can also be made for affiliated nodes such as the olfactory regions (AON, MOB, and Tu) and the medial prefrontal cortex (ACC, CG, IL, and PrL) due to their role in social cognition (Brennan and Kendrick, 2006;Tobin et al., 2010;Lee and Harris, 2013;Nakajima et al., 2014;Bicks et al., 2015;Toor et al., 2015;Oettl et al., 2016;Williamson et al., 2018). Collectively, these regions interact to coordinate incoming social information with context appropriate social responses. Indeed, context dependent behavioral plasticity is likely attributed to changes in coordinated neural activity between nodes of the network, rather than to differential activity of individual brain regions per se (Goodson and Kabelik, 2009;Teles et al., 2015;Johnson et al., 2016).
Function of the SDMN is influenced by oxytocin (OT). OT is a neuropeptide implicated in both social (e.g., maternal care, affiliation, and stress) and sexual (e.g., arousal, ejaculation, and motivation) behaviors across vertebrate species (O'Connell and Hofmann, 2012;Anacker and Beery, 2013). For example, manipulation of OT signaling alters multiple social behaviors including prosociality/aggression, social recognition, short-term social memory, alloparenting, and pup care (Ferguson et al., 2000(Ferguson et al., , 2001Francis et al., 2000;Champagne et al., 2001;Consiglio et al., 2005;Olazábal and Young, 2006a,b;Choleris et al., 2007;Beery et al., 2008;Reddon et al., 2014;Chang et al., 2015). Furthermore, variability in distribution patterns of OT receptors in SDMN regions suggests that oxytocinergic signaling contributes to species-specific adaptations in social behavior (reviewed in Anacker and Beery, 2013). Finally, both central and sitespecific (NAcc) manipulation of OT receptor signaling disrupts coordinated activity among SDMN regions (Johnson et al., 2016(Johnson et al., , 2017. Thus, OT is a key mechanism for sculpting social behavior within and between species, prospectively through facilitating context-specific changes in coordinated activity between nodes of the SDMN. The naked mole-rat (Heterocephalus glaber) exhibits the most extreme form of sociality known in mammals: eusociality. These small, approximately mouse-sized, rodents live in large colonies of up to ∼300 animals with reproduction restricted to one breeding female, the queen, and one to three males (Brett, 1991b;Lacey and Sherman, 1991). The rest of the colony consists of non-breeding subordinates of varying age, which are highly social toward members of their own colony. In contrast, naked mole-rats can be very xenophobic and highly aggressive toward intruding members of another colony . Colony members engage in diverse behaviors such as foraging and food-sharing, cooperative care of pups, vocal communication, communal huddling, and colony maintenance and defense (Withers and Jarvis, 1980;Jarvis, 1981;Brett, 1991a;Pepper et al., 1991). There is individual variability in performance of these behaviors, resulting in stable yet plastic task specialization (Jarvis, 1981;Lacey and Sherman, 1991;Mooney et al., 2015b). Subordinates have been further split into two subcastes: workers and soldiers, responsible for colony maintenance and colony defense, respectively (Jarvis, 1981;Lacey and Sherman, 1991;Mooney et al., 2015b).
The oxytocinergic system contributes to the remarkable sociality found in naked mole-rats. Subordinates of both sexes have more OT neurons than breeders in the PVN (Mooney and Holmes, 2013) and peripheral administration of OT to subordinates increases prosocial behaviors in-colony (huddling) and during unfamiliar conspecific interaction tests (proximity, investigation) (Mooney et al., 2014). Among subordinates, workers have higher levels of c-Fos/OT immunoreactive neurons than soldiers in the PVN, accompanied with lower levels of aggression, further suggesting OT promotes prosocial behavior in workers (Hathaway et al., 2016). Naked mole-rats also express OT receptors in the SDMN. They have more OT receptors in the CeA, MeA, BNST, and NAcc in comparison to solitary cape mole-rats, suggesting that OT action in these regions is associated with colonial living (Kalamatianos et al., 2010). Furthermore, OT receptor density varies within naked mole-rats with breeding males showing higher binding than breeding females in the NAcc and males overall showing higher binding than females in the MeA (Kalamatianos et al., 2010;Mooney et al., 2015a).
We hypothesize that OT signaling mediates coordinated neural activity within nodes of the SDMN to promote contextspecific social behavior in naked mole-rats. To test this, we had two experimental goals: (1) to determine if coordinated neural activity in the SDMN varies according to social context and (2) to examine if activity of SDMN regions is altered by central OT manipulation. To achieve these goals, we exposed subordinate naked mole-rats of both sexes to three unique social stimuli: a familiar adult conspecific, an unfamiliar adult conspecific, or 1-week-old pups from their home colony while treating with either saline, OT, or an OT antagonist (OTA). Following social exposure, we used immunohistochemistry to stain for the immediate early gene c-Fos to assess activation of brain regions in the SDMN. We predicted that coordinated c-Fos expression within the SDMN regions would differ by social context, and that manipulation of central OT receptor signaling would disrupt these patterns of connectivity.

Animals and Housing
A total of 75 adult subordinate naked mole-rats were used (32 females and 43 males). Animals were considered adults if they were both over 20 g in weight and over 1-year old. Animals weighed between 22 and 60 g were housed on a 12:12 light/dark cycle at 28-30 • C and given ad libitum access to sweet potato and a wet 19% protein mash (Harlan Laboratories, Inc.). Animals lived in colonies comprised of large (45.75 cm L × 24 cm W × 15.25 cm H) and small (30 cm L × 18 cm W × 13 cm H) polycarbonate cages connected by plastic tubes (25 cm L × 5 cm D). Animals were collected from one of nine colonies ranging in size from 19 to 49 individuals. All procedures adhered to federal and institutional guidelines and were approved by the University Animal Care Committee. A summary of the experimental workflow is presented in Figure 1.

Intracerebroventricular Cannulation
Experimental animals were implanted with intracerebroventricular cannulae targeting the lateral ventricle as previously described (Mooney and Holmes, 2015). Briefly, mole-rats were deeply anesthetized using isoflurane (induction: 3%, delivered at a rate of 1 L/min inhalation; maintenance: 2%, delivered at a rate of 1 L/min) and the surgical site was cleaned and sterilized with 70% EtOH and then 10% iodine solution (Betadine; repeated twice). Animals were positioned in a stereotaxic instrument (Benchmark TM , MyNeurolab.com, St. Louis, MO, United States) and a 1.5 cm incision was made on the top of the head. The skin and muscle were shifted to reveal the skull, which was then cleaned with sterile saline and dried. A small hole was drilled in the skull 0.9 mm lateral and 1 mm anterior to bregma on the right side of the animal and a 22-gauge stainless-steel guide cannula with a 2 mm pedestal (Plastics One, Roanoke, VA, United States) was lowered to 3 mm below the top of the skull. Cyanoacrylate gel was applied to the base of the pedestal to secure the guide cannula to the top of the skull. A dummy cannula was inserted in the guide cannula to prevent exposure, infection, or occlusion. The muscle was then laid back on the skull around the cannula and the skin was sutured together over the pedestal of the guide cannula. A small dab of cyanoacrylate gel was placed on the outside of the dummy cannula at the juncture where the pedestal of the guide cannula and the cap of the dummy cannula meet. Ten minutes prior to the completion of surgery, animals were injected with ketoprofen (Anafen R , Merial; 5 mg/kg BW). This injection was also administered once a day for 3 days post-surgery. Animals were placed on top of a heating pad in a cage with clean bedding for 8 h for recovery before being returned to their home colonies.

Familiar Conspecific Interaction Tests
Twenty-one experimental animals (FAM; 9 females and 12 males) were removed from their colony and individually placed in a clear polycarbonate cage (L 43 × W 22 × H 21 cm) lined with corncob bedding. After a 20-min habituation period, an adult conspecific from the experimental animal's home colony (same-sex pairing, N = 12, opposite-sex pairing, N = 10; sex matching randomly assigned) was placed in the cage in the corner farthest away from each experimental animal's current position. Each animal's behavior was recorded for 20 min with a Sony Handycam R . All animals were then returned to their home colony. Twentyfour hours after the baseline test, animals were again removed from their home colony and placed in a clear polycarbonate cage for 10 min before receiving intracerebroventricular infusions of either OT, an OTA, or saline. Seven animals received 0.25 µg of OT (OT acetate salt hydrate or α-hypophamine, No. O6379; Sigma), seven animals received 1 ng of the specific OTA (d(CH2)51,Tyr(Me)2,Thr4,Orn8,des-Gly-NH29)-Vasotocin, No. H-2908 BACHEM), and seven animals received 1 µl of the sterile saline. These doses were chosen as they produce behavioral effects in similarly sized mammals (Ferguson et al., 2000;Shahrokh et al., 2010;Samuelsen and Meredith, 2011). For infusion, animals were lightly anesthetized with isoflurane (2% delivered at a rate of 1 L/min) and the dummy cannula was replaced with an internal infusion cannula that sat 0.1 mm below the base of the guide cannula. One microliter of the drug or saline was infused via 500-µl Bas gas-tight syringes (MD-0050; Bio Analytical Systems) connected to the internal cannula with PE50 tubing. Infusions were automated at a rate of 1 µl/min with a Harvard infusion pump (Harvard Apparatus Inc. 22, Natick, MA, United States). Animals were then removed from anesthesia and recovered in ∼2-3 min. After 10 additional minutes, a novel stimulus animal was introduced following the same procedure as for the baseline test. Target behaviors for the baseline test and the manipulation test (Table 1) were scored by an observer blind to the experimental condition using Observer XT software (Noldus). Aggression, anogenital investigation, and investigation (flank/face) were scored in animals interacting with a familiar or unfamiliar conspecific. Pup carrying and pup investigation were scored in animals interacting with pups.

Unfamiliar Conspecific Interaction Tests
Thirty experimental animals (UNFAM; 15 females and 15 males) were used in this paradigm. All procedures and manipulations were identical to the familiar conspecific paradigm described above with the exception of the stimulus animal. In this case, an adult conspecific from an unfamiliar colony was used as the stimulus (same-sex pairing, N = 12, opposite-sex pairing, N = 10; sex matching randomly assigned). Ten animals received OT, 10 animals received OTA, and 10 animals received sterile saline.

Pup Interaction Tests
Twenty-four experimental animals (PUP; 8 females and 16 males) were used in this paradigm. Again, procedures and manipulations were identical to the FAM paradigm with the exception of the social stimulus. In this case, three pups (∼1 week old) from the experimental animal's home colony were placed in the cage in the corner farthest away from each animal's current position. Eight animals received OT, nine animals received OTA, and seven animals received sterile saline.

Tissue Collection and c-Fos Immunohistochemistry
One hundred minutes after the start of each interaction test (FAM, UNFAM, PUP) in which OT activity was manipulated, animals were overdosed with avertin (400 mg/kg) and rapidly decapitated.

c-Fos Quantification
The number of c-Fos immunoreactive cells was counted using either OpenCFU (Beta version 3.9.0) or ImageJ (Rasband, 2012;Geissmann, 2013). For a given brain region, counts for all animals were performed using the same program. As tissue tearing and damage was more likely on the side of the brain ipsilateral to the cannula implantation, counts were done unilaterally on the side contralateral to the implantation. On OpenCFU, this was done automatically with a threshold setting of 5 and a radius setting of 5. Counts were verified by an observer blind to the experimental condition for all sections. Areas counted using OpenCFU were the AH, BLA, BNST, CeA, LS, MeA, MS, NAcc, PG, PO, PVN, SON, and VMH. Using ImageJ, images were converted to 8-bit, a threshold of 180-190 applied, and particles larger than 10 pixels 2 counted. Areas counted using ImageJ were the ACC, AON, CA1, CA2, CA3, caudate, CG, dDG, IL, MOB, PIC, PrL, Tu, vDG, VP, and VTA. For each region of interest, a photomicrograph was taken on three consecutive slices of tissue using a microscope mounted camera at 200× magnification. The CA1, CA2, CA3, dDG, and vDG were taken at 400× magnification to prevent counting other surrounding regions. Photomicrographs contained the same area for all animals, corresponding to the given region; all c-Fos labeled cells in the photomicrograph were counted. The exceptions to this were the BNST, which only had pictures taken from two consecutive slices because only two sections reliably contained the BNST in each series, and the PAG and VTA, which only had one slice counted because more caudal sections were not collected during slicing.

Placement Confirmation
Before euthanizing animals, one saline-treated animal from each behavioral paradigm was also infused with 1 µl of india ink (10% v/v) in order to confirm that fluid was diffusing throughout the brain. For all animals, tissue was examined to ensure that the cannula clearly penetrated the lateral ventricle. Because we could not definitively confirm penetration of the ventricles in two FAM animals (one receiving saline and one receiving OTA), two UNFAM animals (both receiving OT), and three PUP animals (one receiving saline, one receiving OT, and one receiving OTA), these animals were excluded from analyses. Final sample sizes are reported in the section "Social behavior paradigms."

Statistical Analyses
All statistical analyses were performed on raw counts of c-Fos immunoreactive cells. First, a linear mixed-effects model was run using the nlme and lmerTest packages in R (R Development Core Team, 2011; Kuznetsova et al., 2016;Pinheiro et al., 2016). c-Fos immunoreactivity was the response variable and behavioral paradigm, OT manipulation, brain region, and sex were predictor variables, using animal ID and immunohistochemistry batch as random effects. The linear mixed-effects model was used to determine if any predictor variables could be collapsed for the rest of the analyses; due to non-significant effects, sex was dropped as a variable for subsequent analyses. We then performed brain region-specific linear mixed-effects models with c-Fos as the response variable, behavioral paradigm, and OT manipulation as predictor variables, and batch as a random effect, again using the nlme and lmerTest packages in R. A Bonferroni correction was used to adjust for multiple testing (30 tests = p-value set at 0.0017). Main effects reaching this criterion (p < 0.0017) were followed with Tukey's HSD post hoc tests. We next explored coordinated activity between brain regions by examining correlations in c-Fos expression across all brain regions measured for each behavioral paradigm-by-OT manipulation group. The Hmisc package in R was used to perform pair-wise Pearson correlations for all brain regions, followed by visualization using the corrplot package (Wei and Simko, 2016;Harrell, 2017). No thresholding was applied, but significant correlations are marked by white asterisks on the plots. To further explore the effect of behavioral paradigm and OT manipulation on how c-Fos is expressed between regions, networks were produced for each behavioral paradigm-by-OT manipulation group. For the networks, correlation p-values were corrected using the Benjamini-Hochberg procedure at a 5% false-discovery rate. These correlations were then extracted and plotted using Cytoscape (Shannon et al., 2003). Cytoscape's builtin Network Analyzer tool was used to plot node size according to degree (how many edges correspond to a given node) and edges according to correlation strength. On the networks, nodes are clustered and color-coded based on anatomy and literaturebased functions into the following: olfactory (AON, MOB, PIC, Tu), mPFC (ACC, CG, IL, PrL), social behavior network (AH, BNST, LS, MS, PG, PO, VMH), amygdala (BLA, CeA, MeA), reward (caudate, NAcc, VP, VTA), OT production (AntPVN, PostPVN, SON), and hippocampal (CA1, CA2, CA3, dDG, vDG).
Finally, to examine if OT manipulation affected behavior, we tested two linear mixed-effects models as described above. First, to test whether injection manipulation itself affected behavior (because animals were anesthetized for ICV injections), we modeled difference in behavior duration between baseline and test day as the response variable, treatment and sex as independent variables, and animal ID as the random effect variable. Second, to test whether behavior on the test day was altered by OT manipulation, we modeled duration of behavior on test day as the response variable, treatment and sex as independent variables, and animal ID as the random effect variable. Both models were repeated for each behavior tested per behavioral paradigm and corrected for multiple testing using the Bonferroni method (FAM/UNFAM: three comparisons = p-value < 0.017; PUP: eight comparisons = p-value < 0.0062).
Then, to determine if c-Fos expression is related to behavior, we clustered brain regions using factor analysis to reduce analyses performed: all saline-treated animals across paradigms were included using the dimension reduction function in SPSS (IBM Corp, 2016). Principal axis factoring with a Promax rotation was used for the unsupervised clustering of brain regions; four clusters (as listed in the section "Results") were produced. Next, for each paradigm (collapsed across OT manipulation), a Pearson correlation between the summed c-Fos counts for a given cluster and given behavior on test day was performed. For FAM and UNFAM animals, duration of anogenital investigation, face/flank investigation, and aggression was scored. For PUP animals, duration of pup carrying and pup interaction was scored. P-values were corrected for multiple testing using the Bonferroni method (FAM/UNFAM: 12 comparisons = p-value < 0.0042; PUP: 8 comparisons = p-value < 0.0062).
For brain region-specific analyses, only results that were statistically significant after correcting for multiple testing (Bonferroni method) are reported here (p < 0.0017). All effects with p < 0.05 are shown in Table 2 while raw data are available in the Supplementary Table 1. A main effect of behavioral paradigm [F(2,61) = 7.61, p = 0.0011; Figure 2A] revealed altered c-Fos expression in the MOB, with PUP animals having higher expression of c-Fos relative to UNFAM animals (p = 0.002). This same pattern was seen in the Tu [F(2,61 = 10.9, p < 0.0001; Figure 2B], with a significant difference between PUP and UNFAM animals (p = 0.00012) and a trend between FAM and UNFAM animals (p = 0.059). In the IL, a main effect of behavioral paradigm [F(2,61) = 7.3, p = 0.0015; Figure 2C] revealed reduced c-Fos in UNFAM animals relative to both FAM and PUP animals whereas in the VP [F(2,63), p = 0.0002; Figure 2D], c-Fos was higher in FAM animals relative to both UNFAM and PUP groups. Main effects of OT manipulation revealed higher c-Fos expression in OTA-treated animals relative to both saline-and OT-treated animals in the PIC [F(2,63) = 6.9, p = 0.00167; Figure 3A], ACC [F(2,61) = 7.4, p = 0.0014; Figure 3B], and PrL [F(2,61) = 8.3, p = 0.000862; Figure 3C]. UNFAM, animals exposed to unfamiliar conspecific; FAM, animals exposed to familiar conspecific; PUP, animals exposed to pups; OT, oxytocin; OTA, oxytocin antagonist. Photomicrographs showing the reduced c-Fos immunoreactivity in the IL of an UNFAM animal (E) compared to the IL of a (F) FAM animal. IL, infralimbic cortex; Tu, olfactory tubercle; VP, ventral pallidum; UNFAM, animals exposed to an unfamiliar conspecific; FAM, animals exposed to a familiar conspecific.

Correlation Matrix Plots and Network Visualization Reveal the Disruption of SDMN Pathways by OT Manipulation
The correlation matrices were performed within each behavioral paradigm-by-OT manipulation group (Figure 4). Regions were clustered and color-coded based on literature and known functional connections in order to better visualize relationships between related regions. The plots demonstrate that manipulation of OT signaling, either by OT or OTA treatment, alters correlated activity among SDMN regions. In saline-treated FAM animals (Figure 4A), significant correlations are largely restricted to core social behavior network and reward-related regions. While OT (Figure 4B) had modest effects on c-Fos expression patterns, OTA ( Figure 4C) caused a striking increase in positive correlations between almost all brain regions. In saline-treated UNFAM animals (Figure 4D), the SDMN-related regions are active and positively correlated. Treatment with OT ( Figure 4E) or OTA ( Figure 4F) disrupts these positive correlations. The saline-treated PUP animals ( Figure 4G) are similar to UNFAM animals in that brain regions are significantly, positively correlated to one another and this pattern is disrupted by both OT (Figure 4H) and OTA ( Figure 4I). The networks, like the correlation matrices, are plotted as behavioral paradigm-by-OT manipulation groups; however, unlike the correlation plots only relationships that were significant after applying Benjamini-Hochberg thresholding were included in the networks (5% false discovery rate). The networks confirm that OTA treatment alters correlated activity across the SDMN for all paradigms: by increasing connections in the FAM paradigm and disrupting connections in the UNFAM and PUP paradigms. OT had similar effects as OTA in the UNFAM and PUP paradigms. In the network for FAM salinetreated animals (Figure 5A), the NAcc, PIC, amygdalar regions, and some parts of the social behavior network are key nodes with highly correlated co-expression. OTA ( Figure 5C) results in most all regions except the prefrontal cortex regions and SON highly co-expressing, while OT ( Figure 5B) results in a more random expression pattern. In the network for UNFAM saline-treated animals (Figure 5D), the SBN, prefrontal cortex, amygdala, NAcc, SON, MOB, and dDG are important nodes. Almost all edges are removed by OT ( Figure 5E) treatment while OTA (Figure 5F) activates hippocampal regions but results in few correlations between other regions. Finally, in the network for PUP saline-treated (Figure 5G), there is high co-expression in most of the SBN, amygdala, rewardrelated regions, PIC, AntPVN, and dDG. OT ( Figure 5H) and OTA (Figure 5I) treatment eliminated most of these edges.

Factor Analysis Reveals Clusters of Related Regions Reflecting Known Patterns
There were no differences in behavior between baseline and test day, indicating that ICV injection did not in and of itself alter behavior. OT manipulation did not significantly affect behavior on test day. A summary of these results can be found in Supplementary Tables 1, 2. Factor analysis revealed four clusters of brain regions (Table 3). Cluster 1 contains most classic SDMN regions (AH, BNST, LS, MeA, MS, PG, PIC, PO, VMH), cluster 2 is prefrontal cortex and OT production-related regions (ACC, AntPVN, CG, IL, PostPVN, PrL, SON), cluster 3 is hippocampal (CA1, CA2, CA3, dDG, MOB, vDG), and cluster 4 is olfactory (AON, Tu). Correlations between the summed counts of c-Fos within each cluster and behaviors revealed no significant relationships between any particular behavior and brain region cluster.
that the SDMN is conserved in the highly social naked molerat and that it functions in a context-dependent manner. Furthermore, we report that OT manipulation, both via OT and OTA, is generally disruptive to coordinated neural activity, though specific effects are context-dependent. Finally, these data suggest that while structures such as the medial prefrontal cortex, hippocampus, and olfactory regions are important for processing social stimuli, they are not direct members of the SDMN.
By examining neural activity in SDMN regions following different behavioral paradigms, we have demonstrated that, perhaps not surprisingly, this network is functionally conserved in naked mole-rats. This is consistent with previous neurochemical analyses in SDMN regions in this species (Rosen et al., 2007(Rosen et al., , 2008O'Connell and Hofmann, 2012;Holmes et al., 2013;Mooney et al., 2015a;Beery et al., 2016). Factor analysis shows the examined brain regions from 4 clusters ( Table 3), the first of which includes all the SDMN brain regions together except the hippocampus. The other three clusters comprise anatomically and functionally related regions (i.e., medial prefrontal cortex, hippocampus, and olfactory regions, respectively). This pattern of clustering corroborates the existence of a SDMN in the naked mole-rat, and further suggests the other three clusters are perhaps related, but not directly a part of the SDMN. Interestingly, the SDMN showed correlated c-Fos immunoreactivity in the FAM saline-treated group. Given the colony living of naked mole-rats, their baseline experience is consistent interaction with multiple familiar conspecifics. When Newman proposed the social behavior network, one limitation was that without social stimuli, there is no baseline activity for the associated regions (Newman, 1999;O'Connell and Hofmann, 2011). In the case of eusocial naked mole-rats, it is likely the SDMN is always active and integral for navigating colony life.
The SDMN is thought to modulate behaviors necessary for navigating varied social environments (O'Connell and Hofmann, 2011). The behavioral paradigm-specific patterns of neural activity that we report support the hypothesis that coordinated activation of brain regions within the SDMN is contextdependent. Examination of the correlation matrices (Figure 4) reveals that while saline-treated animals in all three behavioral paradigms show coordinated activation of SDMN regions, the network incorporates the medial prefrontal cortex regions in UNFAM animals and hippocampal regions in PUP animals. The medial prefrontal cortex has been implicated in social approach, social ascent, facilitating aggressive behavior and in establishing social rank in mice, hence activation could be associated with need for assessing status and establishing hierarchy between the two unfamiliar animals (Wang et al., 2011;Takahashi et al., 2014;Lee et al., 2016;Williamson et al., 2018). It is interesting to speculate that the coordinated activation of hippocampal regions in PUP animals reflects learning of new colony members and is consistent with the increase in cell proliferation observed following pup exposure in prairie voles and fatherhood in male California mice, both species that exhibit biparental care (Ruscio et al., 2008;Hyer et al., 2016). Alternatively, this could be related to spatial memory required for returning pups to their nest as seen in female rats (Kinsley et al., 1999). In addition to paradigm-specific patterns of correlated activity, we also report paradigm-specific effects on c-Fos expression in individual brain regions. For example, the VP had significantly higher c-Fos immunoreactivity in FAM 3 | Factor analysis results following application of principal axis factoring to saline-treated animals (collapsed across task), identifying four clusters of brain regions.

Cluster
Brain regions animals ( Figure 2C). This region is involved in formation of social attachment in socially monogamous prairie voles (Lim and Young, 2004;Barrett et al., 2013;Zheng et al., 2013).
In UNFAM animals, who are meeting a novel animal and presumably assessing relative status, there was significantly lower activation of the MOB, Tu, and IL (Figures 2A-C). The olfactory system is essential to pair bonding in prairie voles (Curtis et al., 2001), and, while naked mole-rat subordinates are reproductively inactive, they are highly affiliative with colony members. Naked mole-rats show preference for familiar olfactory cues over unfamiliar ones (Toor et al., 2015), likely contributing to in-colony social recognition and bonding. Interestingly, the IL (in addition to the PrL) is also activated in a shifting social environment in which social ascent is possible, as in the case of removing an alpha male in a mouse hierarchy (Williamson et al., 2018). Administration of a histone deacetylase inhibitor to the IL increases stress behaviors in a Syrian hamster model of social defeat (McCann et al., 2017), further implicating the IL in social interactions involving dominance.
The data reported here also support the hypothesis that OT mediates coordinated neural activity within the SDMN to promote context-specific social behavior. OTA treatment altered coordinated c-Fos expression in all tasks but while it increased connectivity in FAM animals, it reduced connectivity in UNFAM and PUP animals (Figures 4, 5). Brain region-specific analyses revealed that in all behavioral paradigms, OTA significantly increased c-Fos expression in the PIC, ACC, and PrL (Figure 3), suggesting these regions are particularly sensitive to perturbations in OT signaling. While examination of OT receptor distribution in naked mole-rats did not evaluate anterior regions like the ACC and PrL, subordinate NMRs have negligible OT receptor binding in the PIC (Kalamatianos et al., 2010). Thus, OTA could be acting elsewhere to indirectly influence activity in these regions. Similar to OTA, we also found evidence that OT treatment alters connectivity in the SDMN with reduced coordinated activity in animals from the UNFAM and PUP paradigms (Figure 5). This is not necessarily surprising as exogenous OT is unlikely to reflect endogenous levels and, furthermore, our central delivery ensures OT would act at multiple targets, interfering with "normal" OT signaling. Importantly, we cannot infer whether more or fewer connections are optimal as we have yet to directly relate coordinated activity to specific behaviors. Rather, we can conclude that central OT manipulation affects connectivity in the SDMN in various social contexts and future studies will identify how specific regions are involved.
Indeed, the NAcc is a likely hub modulating the role of OT between other SDMN regions. OT acts in the NAcc to modulate sociosexual interactions and mating in prairie voles (Johnson et al., 2017): administration of an OT receptor antagonist directly to the NAcc changes connectivity between NAcc and other SDMN regions. The NAcc is also involved in reproductive behavior and social dominance in naked molerats. For example, OT receptor binding is higher in breeding male NAcc compared to breeding females (Mooney et al., 2015a) and expression of genes involved in social suppression of reproduction (Kiss1, Npvf, Gpr147, Tac3r) varies by sex and status (Faykoo-Martinez et al., 2018). The current data demonstrate that the NAcc is a key node in all three behavioral paradigms (Figure 5). In the FAM network for saline-treated animals, NAcc coordinated activity is enhanced by OTA but knocked out by OT treatment. In the UNFAM network for saline-treated animals, the NAcc is the key active node out of the mesolimbic reward circuit and is knocked out by both OT and OTA. Finally, in the PUP network for salinetreated animals, all mesolimbic reward circuit regions are key nodes. This includes the NAcc, which is then knocked out to varying extents by OT and OTA treatment. Given its role in assigning valence to social stimuli, it will be important to evaluate OT effects directly in the NAcc to better understand how this region alters coordinated neural activity, similar to Johnson et al. (2017), across different social contexts.
We have taken a broad approach in our first attempt at teasing apart function in the naked mole-rat SDMN. For this reason, we did not break down groups by subcaste or focus on a single sex for stimulus animal. Subordinates are prepubertal, do not exhibit sex differences in behavior , and take a week to begin demonstrating behaviors characteristic of reproductive maturation (Margulis et al., 1995;Clarke and Faulkes, 1997;Mooney et al., 2015b;Swift-Gallant et al., 2015;Toor et al., 2015). Thus, it is likely that the sex of stimulus animals would cause a negligible effect given the very short period of removal from colony. It, of course, also warrants mention that our experimental design was likely insufficiently powered to robustly examine OT manipulation by behavioral paradigm interactions. Our sample size was limited by the challenges of rearing naked mole-rats in captivity (taking 1 year to reach adulthood coupled with unpredictable breeding). Low statistical power might also explain why we did not detect any significant effects of OT manipulation on behavior or significant relationships between c-Fos activity and behavior, though the lack of behavioral results following central administration of OT is not necessarily unexpected. While we have previously shown that peripheral OT administration changes in-colony behavior (e.g., time spent huddling) and proximity to novel conspecifics, we did not measure in-colony behavior in the current study and also used a different testing apparatus for outpairing (single chamber vs. double chamber) (Mooney et al., 2014). To address these limitations, we provide the raw data in Supplementary Information with the hope that future work by us and others will build on this sample.
Here we report the first formal investigation of activation of the SDMN in the highly social naked mole-rat. We have demonstrated that coordinated neural activity within this network and with related regions varies according to social context, and that this coordinated activation is altered by manipulation of the OT system. The pattern of connectivity associated with each behavioral paradigm suggests a role for the NAcc in social valence and sociosexual interaction, the mPFC in assessing/establishing social dominance, and the hippocampus in pup recognition.

AUTHOR CONTRIBUTIONS
MH and SM contributed conception and design of the study. SM and MF-M collected the data. MF-M performed statistical analyses and wrote the first draft of the manuscript. SM performed the experiments and wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

FUNDING
This work was funded by a NSERC CGS M and NSERC PGS D to MF-M, QEII-GSST OGS to SM, and NSERC grants (RGPIN 2011-402633, RGPIN 2018-04780, and RGPAS 2018 and an Ontario Early Researcher Award to MH.