Effects of Affiliative Human-Animal Interaction on Dog Salivary and Plasma Oxytocin and Vasopressin.

Oxytocin (OT) and vasopressin (AVP) are neuropeptides with diverse effects on social behavior, cognition and stress responses. Recent studies suggest that OT facilitates and responds to affiliative forms of human-animal interaction (HAI). However, previous studies measuring OT and AVP in dogs have been limited to measures from blood or urine, which present concerns related to the invasiveness of sample collection, the potential for matrix interference in immunoassays, and whether samples can be collected at precise time points to assess event-linked endocrine responses. Previous studies from our laboratory validated salivary measures of OT and AVP in dogs, however, it is currently unknown whether these measures respond dynamically to aspects of HAI. Here, we investigated the effects of affiliative forms of HAI on both plasma and salivary OT and AVP in dogs. We employed a within- and between-subjects design with a group of Labrador retrievers and Labrador retriever × golden retriever crosses (23 females, 15 males). Half of the dogs engaged in 10 min of free-form friendly interaction with a human experimenter (HAI condition), and the other half rested quietly in the same environment, without human interaction (control condition). We collected blood and saliva samples before, and immediately following both experimental conditions, and all samples were analyzed using enzyme-linked immunosorbent assays (ELISAs) following previously validated protocols. Dogs participating in HAI exhibited a significant increase in both salivary OT (+39%) and plasma OT (+5.7%) whereas dogs in the control group did not. Salivary AVP showed no change in the HAI group but increased significantly (+33%) in the control group. Plasma AVP decreased significantly following HAI (-13%) but did not change across time in the control condition. Within the dogs exposed to HAI, increases in salivary OT, and decreases in plasma AVP, were predicted by the extent of affiliative behavior between the dog and human (indexed by scores from a principal components analysis of social behaviors between the dog and human). Collectively our results suggest that measures of salivary OT and AVP provide useful biomarkers in studies of HAI, and afford a flexible and non-invasive toolkit than can be employed in diverse research contexts.

Introduction

Studies throughout the last three decades have explored the psychological and physiological effects of human–animal interaction (HAI). Often the aim of such studies is to characterize the mechanisms through which non-human animals affect human health and wellbeing, and in turn, how interaction with humans affects animal participants. Recently such studies have begun to focus on the neuropeptide oxytocin (OT), which, together with (structurally related) arginine vasopressin, is well known for its roles in facilitating selective social bonds, and regulating various aspects of social behavior and cognition in mammals (Carter et al., 2008). For example, central OT administration can facilitate maternal behavior in sexually naïve rodents (Pedersen and Prange, 1979; Caldwell and Young, 2006), and both OT and AVP are critical to the formation of partner preference in pair-bonded species (Carter et al., 1992). With regard to the establishment of selective social attachment, the OT receptor (OXTR) is highly expressed in the nucleus accumbens of some pair-bonded species, suggesting that OT may help to encode social reward through modulation of the mesolimbic dopamine pathway (Insel and Shapiro, 1992; Lim and Young, 2006).

Research on the role of OT in HAI has been conducted by measuring OT release during interactions between humans and dogs, and evaluating the effects of exogenous OT administration on dog behavior in this context. Studies of endogenous OT have revealed increases in blood or urinary OT concentrations in both humans and dogs following affiliative interaction in dog-human dyads (Miller et al., 2009; Nagasawa et al., 2009; Handlin et al., 2011; Rehn et al., 2014; Nagasawa et al., 2015). Additionally, preliminary evidence suggests that dogs bred for friendly and non-aggressive temperaments are characterized by high levels of plasma OT, relative to pet dogs (MacLean et al., 2017b). Studies using intranasal OT administration in dogs suggest that exogenous OT can facilitate social play (Romero et al., 2015), and promote bonding with both conspecifics and humans (Romero et al., 2014). Additionally, recent studies have demonstrated that exogenous OT can improve dogs’ sensitivity to human communication (Oliva et al., 2015; Macchitella et al., 2016). Thus, a rapidly growing body of research suggests that OT pathways may be centrally involved in affiliative forms of HAI (Beetz et al., 2012; MacLean and Hare, 2015; Carter and Porges, 2016).

We are not aware of any studies investigating AVP in the context of HAI, however, many of the effects of OT are also critically dependent on AVP (Carter, 1998; Landgraf and Neumann, 2004). Although AVP plays important roles in selective sociality, particularly among males (Caldwell et al., 2008), its effects are often antagonistic to those of OT. For example, many of the social effects of OT during HAI may be facilitated by OT’s attenuation of sympathetic arousal (Kis et al., 2014; Buttner, 2016) through actions in the hypothalamus (Dabrowska et al., 2011), and on the vagus nerve (Porges, 2003, 2007, 2011). In contrast, AVP activates the hypothalamic-pituitary-adrenal (HPA) axis, and is more strongly linked to anxiety and aggression (Coccaro et al., 1998; Neumann and Landgraf, 2012). We are aware of two studies investigating relationships between AVP and behavior in dogs, both of which revealed positive associations with anxiety or aggression (Hydbring-Sandberg et al., 2004; MacLean et al., 2017b). Given its structural and functional relationships with OT, we expect that AVP may also regulate or respond to aspects of HAI, although these effects are likely to differ from those for OT.

To date, all studies of OT and AVP in dogs engaged in HAI have quantified peptide concentrations in blood or urine samples. Although urine sampling can be performed non-invasively, it yields poor temporal resolution and characterizes long periods of peptidergic activity. In contrast, blood sampling can capture acute changes in peptide release, but is an invasive procedure that may induce stress and acute pain, and is not well suited to many of the contexts in which HAI studies are conducted. Additionally, OT and AVP rapidly bind to other molecules in blood (and likely urine) which can lead to matrix interference in assays, or erroneously low estimates of peptide concentrations (Martin and Carter, 2013; Brandtzaeg et al., 2016). We recently validated methods for quantifying OT and AVP in dog saliva samples, which can be collected non-invasively, and at precise time points during HAI. Enzyme-linked immunosorbent assays (ELISAs) yielded good parallelism and accuracy, and did not require an extraction procedure. Lastly, we previously measured salivary OT in nursing dams and detected an acute rise in OT associated with milk letdown, providing an initial biological validation of this measure (MacLean et al., 2017a).

In the current studies, we evaluated both plasma and salivary OT and AVP concentrations in dogs, before and after affiliative interaction with a human, or in a control condition. Because the time course of salivary OT/AVP release during HAI was unknown, we first conducted a short pilot study to identify time points associated with changes in salivary OT/AVP concentrations during HAI. We then conducted an experiment with dogs assigned to an HAI or control condition, and assessed OT/AVP changes over time in both groups, and as a function of behavior during the test period. Based on the studies described above, we hypothesized that dogs in the HAI condition would exhibit increases in OT across the study period, and that any such changes would be larger in the HAI group than the control group. We also hypothesized that within the HAI group, changes in OT would be predicted by the extent of affiliative behavior between the human and the dog during the study. Because no studies have investigated effects of HAI on AVP, we had no specific hypotheses for the nature of this response. However, given that AVP also plays important roles in social emotions and behavior, we expected that dogs in the HAI group would exhibit changes in AVP that differed from those of dogs in the control group.

Pilot Study

Prior to the main HAI study, we conducted a short pilot study with 10 dogs to identify the optimal sampling periods for detecting changes in salivary OT and AVP. Based on our studies with nursing dams (MacLean et al., 2017a), and recent studies measuring salivary OT in humans (de Jong et al., 2015), we expected changes in saliva to be rapid. Therefore, we assessed both plasma and salivary OT and AVP at 5 and 10 min following the start of HAI. In both the pilot study and Experiment 1, we opted to collect both blood and saliva, despite the possibility that blood draws may impose a mild stressor with potential to influence hormonal changes detected saliva. We adopted this strategy because no studies have measured salivary OT/AVP in dogs during HAI, and consequently any effects (or lack thereof) in saliva would be challenging to interpret without knowledge of how blood concentrations had changed across the same period.

Method

Subjects

We tested dog subjects from the breeding colony at Canine Companions for Independence (CCI) in Santa Rosa, CA, United States. The pilot sample included 10 dogs (eight female, two male, four Labrador retrievers, six Labrador retriever × golden retriever crosses, mean age = 1.7 years (range = 1.6–1.8 years)). All dogs were pair-housed in indoor–outdoor kennels with ad libitum access to water and daily access to large outdoor play yards. Subjects were tested in a quiet room inside a familiar building on CCI’s campus. Client consent was obtained for participation of all dogs and all animal procedures were approved by the Duke University IACUC (protocol #A138-11-06).

Procedure

Prior to the test, each dog was allowed to rest quietly in a crate (outside the test room) for 30 min. The dog was then taken to a nearby room and we collected a baseline blood and saliva sample. Saliva samples were collected using the Salimetrics Children’s swab as described in MacLean et al. (2017a). Blood samples were collected from the cephalic vein into vacutainers (3 mL) containing ethylenediaminetetraacetic acid (EDTA). Blood and saliva samples were collected concurrently to minimize the time required for these procedures. Specifically, at the start of the collection period, one experimenter placed the swabs between the dog’s cheek and mandibular teeth and gently held the dog’s mouth closed while the second experimenter performed the blood collection. All samples were immediately placed on ice after collection.

Following this initial sample, dogs were allowed to rest in a crate inside the test room for 5 min prior to the start of the behavioral interaction. After 5 min, dogs were released from the crate and allowed to interact freely with the experimenter. The experimenter attempted to engage the dog in friendly interaction, including gently petting the dog, and speaking to the dog in a friendly tone while making eye contact. However, the experimenter allowed the dog to lead these interactions, and dogs were always free to disengage and move away from the experimenter at any point during the interaction. If subjects attempted to engage the experimenter in play (e.g., performing a play bow, chasing, or nuzzling the experimenter), the experimenter engaged with the dog in these more active forms of interaction. After 5 min of HAI we collected a second blood and saliva sample from the dog, and immediately resumed HAI for another 5 min.

The final blood and saliva samples were collected 10 min after the start of HAI and dogs received a food reward at the conclusion of this period (following the final saliva collection). The second and third blood samples were collected from (1) the cephalic vein on the opposite forelimb from the initial draw, and (2) the jugular vein, in order to avoid repeated needle punctures in the same location. We did not use a catheter for repeated collections from the same site because pilot tests revealed that catheters did not reliably maintain access to the vein when dogs were allowed to move freely between collection periods, and the procedures required to position, wrap, and access the catheter repeatedly were deemed more likely to cause discomfort for the dog than single collections from different sites. All biological samples were immediately frozen at -20°C.

Hormone Analysis

All samples were analyzed by enzyme-linked immunoassay (ELISA) following protocols previously validated in our laboratory (MacLean et al., 2017a,b). OT samples were measured using a commercially available kit from Cayman Chemical (Item #500440) and AVP samples were measured using a commercially available kit from Enzo Life Sciences (ADI-900-017A). Saliva samples were not extracted based on the results of validation studies (MacLean et al., 2017a) but plasma samples were processed using solid phase extraction (to isolate free peptide concentrations) with the protocols described in MacLean et al. (2017b).

Statistical Analysis

All statistical analyses were performed in the R programming environment for statistical computing (R Core Team, 2017). For each matrix (saliva, plasma) and hormone (OT, AVP), we used linear mixed models to predict the log transformed hormone concentration as a function of a fixed effect of time (baseline, +5 min, +10 min) and a random effect of subject ID. We used planned Dunnett contrasts to assess mean differences between hormone concentrations at baseline and +5 min, and between baseline and +10 min.

Results and Discussion

Table shows the results of linear mixed models predicting changes in hormone concentrations as a function of time. We inspected the data for time points associated with the largest deviations from baseline values.

Salivary OT increased slightly from baseline to +5 min, with a further increase at +10 min that was marginally different from baseline (p = 0.05; Table ). Changes in plasma OT and both salivary and plasma AVP were minimal, with no clear deviations from baseline. Therefore, we identified the +10-min measure as the most likely to show an HAI related effect in salivary OT, and collected samples at this time point in Experiment 1.

Experiment 1

We tested 38 dogs (11 Labrador retrievers, 27 Labrador retriever × golden retriever crosses, 23 female, 15 male, mean age = 1.8 years, range = 1.6 – 2.2 years). All subjects were assistance dogs in training at Canine Companions for Independence and were housed and cared for as described above. Half of the subjects were assigned to the HAI condition, and half were assigned to the control condition (see below). Groups were matched closely based on sex and breed (Table ).

The HAI procedure was identical to that in the pilot study with the exception that we collected blood and saliva samples at only one time point (+10 min) after the baseline measure, based on preliminary data suggesting the greatest changes in salivary OT at this time. Therefore, the experimenter and dog had 10 continuous minutes of interaction prior to the post-test sample collection. In the control condition, dogs were placed inside an exercise pen (∼3 m × 3 m) in the test room. The experimenter remained in the room with the dog but did not interact, speak to, or make eye contact with the dog during the 10-min test period. We opted to use an exercise pen because this structure provided dogs with opportunities to locomote similarly to dogs in the HAI condition, but prevented dogs from physically accessing the experimenter. Additionally, because all subjects were accustomed to resting quietly in novel environments during their training for assistance work, it was unlikely that containment inside the exercise pen would impose a significant stressor. As in the HAI condition, after 10 min we collected a post-test blood and saliva sample. After the study, dogs in both conditions received a food reward.

Behavioral Coding

To assess whether specific behaviors or forms of social interaction (SI) were related to individual differences in hormonal response, we coded the duration (s) of several behaviors from video. Prior to coding, we reviewed all videos to determine which behaviors could be meaningfully coded, considering both the prevalence and observability for a range of theoretically relevant measures. Locomotion and postural variables were coded for subjects in both the control and HAI conditions, and consisted of the following: (1) locomotion: time walking, running, or jumping, with the onset of the behavior marked by three consecutive steps (to disambiguate locomotion from minor bodily repositioning), and of the offset of the behavior marked by a lack of movement for ≥ 1 s; (2) upright: time standing, walking, running, or jumping; (3) lying (prone): time lying with abdomen against the floor; (4) lying (supine): time lying on back or side with abdomen exposed; (5) sitting: time with rump on ground and forelegs extended.

For subjects in the HAI condition we also coded the following behaviors relating to interaction with the experimenter: (6) physical contact: time during which any part of the experimenter’s body was in physical contact with any part of the dog’s body; (7) licking: time in which the dog’s tongue was in contact with the experimenter’s body; (8) play: time in which the dog engaged in play bows, chasing, or gentle mouthing with the experimenter, with the onset marked by the first occurrence of any of these behaviors, and the offset marked by a period of ≥2 s without the dog engaging in any of these behaviors. A second independent rater scored all behaviors for random sample of ∼50% of observations, and inter-rater reliability was excellent for all measures (Pearson correlation, mean: R = 0.997) (min = 0.993, max = 1.0).

Hormonal Analysis

Plasma and salivary OT and AVP were assayed using the same ELISA kits and protocols employed in the pilot study. Inter-assay coefficients of variation were 11.1% (salivary OT), 4.2% (salivary AVP), 11.7% (plasma OT) and 11.2% (plasma AVP).

All statistical analyses were performed in the R programming environment for statistical computing (R Core Team, 2017). To compare changes in OT/AVP across time between conditions, we used linear mixed models with the log transformed OT/AVP concentrations predicted as a function of fixed effects for sex (male, female), time (pre, post), condition (control, HAI), the time × condition interaction, and a random effect for subject ID. This model was fitted separately with each endocrine measure as the dependent variable. Because we predicted differential effects over time between groups, we conducted planned contrasts from these models assessing the effect of time within condition, and differences between conditions at each time point. Due to the exploratory nature of the study, we used an alpha value of p ≤ 0.05 for significance testing, without correction for multiple comparisons.

Due to the number of behavioral variables, and substantial correlations between them, we conducted principal components analyses (PCA) to derive a smaller set of behavioral measures for statistical modeling. Because the postural and locomotion variables were coded for subjects in both the control and HAI conditions, we performed a PCA with these measures including data from all subjects. For subjects in the HAI condition, we conducted a second PCA including measures related to SI. The postural variable “lying (supine)” was observed predominantly in conjunction with abdomen petting in the HAI condition and was rarely observed in the control group (median duration in control group = 0). Thus, we omitted this measure from the PCA with variables related to locomotion, and included it in the PCA with variables related to SI. To limit the impact of outliers and skew, all behavioral variables were Yeo-Johnson transformed (Yeo and Johnson, 2000), centered, and scaled prior to fitting the PCA. To determine the number of components to retain, we conducted parallel analysis (Horn, 1965) comparing eigenvalues from the actual data to randomly resampled (with replacement, and dimensions equal to those of the original data) and simulated data (random data from a normal distribution). These analyses suggested retention of two components for both the locomotion and SI PCAs. To assess relationships between these behavioral variables and changes in OT/AVP across the experiment, we fitted linear models predicting the log transformed percent change in OT/AVP concentration as a function of fixed effects for sex (male, female) and scores from the first two principal components. To accommodate negative values in the dependent measure, we first added a constant to all values (the absolute value of the smallest observation +1) prior to log transformation. Because the control group did not have component scores for the SI variables, models for this group included only sex and principal component scores for the PCA on locomotion/postural variables. Coefficients for fixed effects were tested with likelihood ratio tests comparing the full model, to nested models with the removal of individual terms.

Results

Between Group OT/AVP Effects

Figure shows the mean OT/AVP concentrations across time in the HAI and control groups. Results from the full mixed model for each matrix are shown in Table , and the results of planned contrasts comparing OT/AVP concentrations between groups, and assessing within group changes over time, are shown in Table .

Group means for oxytocin (OT) and arginine vasopressin (AVP) measured in plasma and saliva before and after Experiment 1. Using between-subjects design, dogs were assigned to conditions consisting of 10 min of affiliative interaction with a human (HAI, N = 19) or 10 min of rest without social contact with the experimenter (Control, N = 19). Dogs in the HAI condition exhibited significant increases in both salivary and plasma OT, and a significant decrease in plasma AVP. Dogs in the control group exhibited no significant changes in salivary or plasma OT, but showed increases in salivary AVP across time. Error bars represent the within-subjects standard error (Cousineau, 2005) and should be interpreted with regard to within-group, but not between group differences. The results of planned between-group, and within-group contrasts are reported in Table . Non-transformed concentration means: Salivary OT – HAI group: pre = 630 pg/mL, post = 874 pg/mL; Control group: pre = 563 pg/mL, post = 642 pg/mL; Plasma OT – HAI group: pre = 18.75 pg/mL, post = 19.83 pg/mL; Control group: pre = 20.11 pg/mL, post = 19.84 pg/mL; Salivary AVP – HAI group: pre = 605 pg/mL, post = 546 pg/mL; Control group: pre = 546 pg/mL, post = 724 pg/mL; Plasma AVP – HAI group: pre = 6.45 pg/mL, post = 5.64 pg/mL; Control group: pre = 6.12 pg/mL, post = 5.89 pg/mL.

For salivary OT, planned comparisons revealed that whereas the control and experimental groups had comparable mean salivary OT levels at baseline, dogs in the experimental group had significantly higher post-test salivary OT than controls. Similarly, salivary OT exhibited a significant increase from baseline to post-test in the HAI, but not the control group (Figure and Table ). For salivary AVP, there were no between group differences at either time point, however, the control group exhibited a significant increase in salivary AVP whereas the HAI group did not (Figure and Table ).

Plasma OT did not differ between groups at either time point, however, only the HAI group exhibited a modest, but significant, increase in plasma OT over the course of the study (Figure and Table ). Lastly, there were no group differences in plasma AVP at either time point, however, the HAI group exhibited a significant decrease in plasma AVP across time, whereas the control group did not (Figure and Table ).

Behavioral Predictors of Changes in OT and AVP

We retained two components from the PCA including variables related to locomotion and posture (LP), which collectively explained 80% of variance in these behaviors. Variable loadings from this model are shown in Table . The first component (LP-PC1, 50% variance explained) had strong positive loadings for locomotion and upright posture and moderate negative loadings for lying in in a prone position. The second component (LP-PC2, 30% variance explained) had a strong positive loading for sitting and moderate negative loadings for locomotion, upright posture, and lying (prone).

For dogs in the HAI condition, we retained two components from the PCA with variables relating to SI, which collectively explained 81% of variance in these behaviors. Variable loadings from this model are shown in Table . The first component (SI-PC1, 53% variance explained) was loaded positively by play and licking, and negatively by physical contact and lying (supine). The second component (SI-PC2, 28% variance explained) was loaded positively by all four SI variables (physical contact, play, licking, and lying supine).

Associations between changes in OT/AVP concentrations and behavior during the test are shown in Table . In the HAI condition, the percent increase in salivary OT was positively associated with SI-PC2 scores, which reflect longer durations of physical contact, play, licking, and lying (supine) with the abdomen exposed to the experimenter. Despite this positive association with salivary OT, there were no associations between changes in plasma OT concentrations and any of the behavioral variables (Table ). There were also no associations between changes in salivary AVP and any of the behavioral variables, however, SI-PC2 scores were negatively related to the percent change in plasma AVP. On average, subjects in the HAI group exhibited a 10% decrease in plasma AVP across the study. However, subjects with SI-PC2 scores in the upper 50th percentile (high levels of affiliative behavior with experimenter) exhibited a larger decrease in plasma AVP (mean change = -18%, SEM = 6.89%) than subjects with SI-PC2 scores in the lower 50th percentile (mean change = -4%, SEM = 4.87). Within the HAI group, LP-PC2 scores were also negatively associated with the change in plasma AVP. Thus, subjects who engaged in more sitting, and less standing, active locomotion and lying (prone) exhibited larger decreases in plasma AVP across time. Within the control group, there were no significant associations between any of the behavioral measures and changes in peptide concentrations, in saliva or plasma (Table ).

Discussion

We assessed changes in dog salivary and plasma OT and AVP in response to HAI, or a control condition. Dogs in the HAI group exhibited increases in both salivary and plasma OT, whereas neither of these effects were observed in the control group. Dogs in the HAI group exhibited a decrease in plasma (but not salivary) AVP across time, whereas dogs in the control group exhibited an increase in salivary AVP, with no significant changes in plasma AVP. Within the HAI group, the extent of the increase in salivary OT, and the decrease in plasma AVP were predicted by the degree of affiliative contact between the dog and experimenter. Therefore, both the between and within group differences are most likely attributable to the SIs between the experimenter and dog, and are unlikely to be accounted for by alternative explanations such as a stress response following the initial sample collection.

Although we observed a significant increase in both plasma and salivary OT, the effect was much more pronounced in saliva, echoing the findings of our pilot study. There are at least three reasonable explanations for this finding. First, plasma OT responses can occur extremely rapidly, and may be most evident within 90 s of a triggering event (Jonas et al., 2009). In contrast, time series analyses suggest a delay in the transfer of hormones from plasma to saliva, and hormonal peaks in saliva have been documented to occur ∼10 min after those in blood (Hernandez et al., 2014). Although the timing, and specific mechanisms through which OT and AVP reach saliva in dogs are unknown, it is probable that changes in salivary concentrations lag behind those in plasma, despite the fact that changes in salivary OT occur quickly relative to other salivary hormones (de Jong et al., 2015). Because this study was designed predominantly for salivary measures, we collected samples at the time point characterized by the highest salivary OT levels in our pilot study. Therefore, is it possible that the largest changes in plasma OT occurred quickly, as has been documented in previous studies (Handlin et al., 2011). Nonetheless, we detected a small but significant increase in plasma OT within the HAI group, suggesting similar effects in both blood and saliva.

Second, whereas OT can be measured without interference in non-extracted dog saliva, other components of dog plasma interfere with OT ELISAs, necessitating additional preparatory procedures such as solid phase extraction for the measurement of free OT (MacLean et al., 2017b). Although extraction procedures can eliminate interfering substances (e.g., albumin proteins), they are often characterized by poor recovery of the target analyte, and thus may eliminate ‘signal’ as well ‘noise.’ Unpublished data from our laboratory suggest that common extraction procedures fail to recover substantial amounts of free OT in plasma, despite performing well with kit standards. Therefore, in addition to benefits related to welfare, by virtue of not requiring extraction procedures, salivary measures of OT also have methodological advantages relative to blood sampling.

Third, it has recently been demonstrated that OT rapidly binds to other molecules in plasma, and bound OT is likely to evade detection (Brandtzaeg et al., 2016). Although the binding properties of OT in saliva are unknown, it is possible that OT remains free – and consequently detectable – in saliva, more so than in plasma. This possibility is supported by the high concentrations of OT that we observed in saliva relative to plasma, an effect that persists even with solid phase extraction of dog saliva samples (MacLean et al., 2017a).

In addition to these effects on OT, dogs in the HAI condition also exhibited a decrease in plasma AVP, whereas dogs in the control group exhibited an increase in salivary AVP. Because AVP can activate the HPA axis (Scott and Dinan, 1998; Aguilera and Rabadan-Diehl, 2000), one plausible explanation is that these differences reflect a short-term reduction in stress reactivity in the HAI group, and increased stress reactivity in the control group. This possibility is consistent with many other studies documenting that positive SIs can buffer stress responses (DeVries et al., 2003) – including interspecies interactions (Schöberl et al., 2016) – whereas social isolation can have the opposite effect. Although dogs in the control group remained in the same room with the experimenter, they were physically separated from him, and this separation may have imposed a mild psychological stressor. Despite this possibility, we observed minimal behavioral indications of stress or anxiety in this population of dogs, who are accustomed to resting in novel environments without physical contact from humans as a part of their training for assistance work. A second possibility is that the initial blood and saliva collection imposed an acute stressor in both groups, but that this event was socially buffered by OT in the HAI condition. The latter possibility is consistent with many studies documenting OT’s ability to supress HPA activity in response to a stressor (Carter and Altemus, 1997). Both possibilities invoke differential stress responses between the two groups, but differ regarding the hypothesized cause of this stress. In the former scenario, physical separation from the experimenter acts as a mild stressor for subjects in the control, but not the HAI condition. In the latter scenario, both groups are exposed to the same mild stressor (i.e., initial blood draw), but stress is buffered more effectively in the HAI group. Although we cannot distinguish between these, and other explanations at present, future studies will benefit by exploring the dynamics between OT, AVP and HPA activity in the context of HAI. Additionally, the differential changes in AVP between groups were observed in different matrices (increased salivary AVP in the control group, decreased plasma AVP in the HAI group), and the explanation for this phenomenon remains unknown. Given that AVP is released in response to both acute and chronic stress, it will be important to investigate the significance of both plasma and salivary vasopressin in future studies explicitly focusing on dogs’ biological responses to potentially stressful events.

Lastly, in addition to the between group differences that we observed, within the HAI group the extent of the increase in salivary OT, and decrease in plasma AVP, were predicted the nature of interactions between the dog and experimenter. Specifically, dogs who engaged in higher levels of physical contact, play, licking and lying in a supine position with the experimenter exhibited the largest increases in OT, and decreases in AVP. Thus, as in previous studies, changes in OT and AVP concentrations depended not only on human contact, but also on the nature of these interactions (Rehn et al., 2014; Nagasawa et al., 2015).

Collectively, our results corroborate previous findings suggesting that OT responds dynamically to HAI, and provide the first data on HAI’s effects on AVP in dogs. Notably, although HAI-related effects on OT have also been observed in humans, we are not aware of any studies examining AVP in humans in the context of interaction with animals. Given that we observed significant decreases in dogs’ plasma AVP concentrations following HAI, we expect that future studies will benefit by incorporating this measure in humans as well. Collectively, these studies suggest that salivary measures of OT and AVP provide non-invasive biomarkers which respond to aspects of affiliative social behavior, and provide researchers with a new set of tools for exploring the roles of OT and AVP in the biology of HAI.

Author Contributions

EM, NG, KL, and CSC designed the experiments. EM and KL collected the data. EM performed immunoassays and analyzed the data. EM, LG, NG, KL, WLM, and CSC wrote the paper.

Conflict of Interest Statement

The authors declare that this research was funded by the WALTHAM® Centre for Pet Nutrition (WCPN) which employed author NG, whose contribution is outlined in the author contribution statement. WLM was also employed by Martin-Protean LLC. All other authors have no conflicts of interest to declare.

Table 1

Pilot study results.

	T5–T0			T10–T0
	β	SE	p	β	SE	p
Salivary OT	0.14	0.17	0.40	0.33	0.17	0.05
Plasma OT	0.03	0.04	0.48	0.01	0.04	0.83
Salivary AVP	-0.19	0.15	0.22	-0.05	0.15	0.73
Plasma AVP	-0.05	0.12	0.67	-0.07	0.12	0.56

Table 2

Subject demographics by condition for Experiment 1.

Condition	Breed	Female	Male
HAI	LAB	5	1
	LGX	6	7
Control	LAB	4	1
	LGX	8	6

Table 3

Results of linear mixed models predicting hormone concentrations as a function of time, condition, sex, and the time by condition interaction.

	Oxytocin						Vasopressin
	Saliva			Plasma			Saliva			Plasma
	β	X2	p	β	X2	P	β	X2	p	β	X2	p
Time	0.13	1.95	0.16	-0.01	0.22	0.64	0.28	4.19	0.04	-0.01	0.22	0.64
Condition	0.11	0.73	0.39	-0.07	0.88	0.35	0.12	0.31	0.58	-0.07	0.88	0.35
Sex	0.01	0.00	0.95	0.06	0.70	0.40	0.19	0.95	0.33	0.06	0.70	0.40
Time × Condition	0.20	2.28	0.13	0.07	2.93	0.09	-0.39	4.01	0.05	0.07	2.93	0.09

Table 4

Results of planned contrasts between the HAI and control groups at each time point, and across time within groups.

	Group comparison (HAI vs. Control)						Within group change (Pre vs. Post)
	Pre-test			Post-test			HAI group			Control group
	β	SE	p	β	SE	p	β	SE	p	β	SE	p
Salivary OT	0.11	0.13	0.39	0.31	0.13	0.02	0.33	0.09	<0.01	0.13	0.09	0.16
Salivary AVP	0.12	0.21	0.58	-0.27	0.21	0.20	-0.10	0.13	0.44	0.28	0.13	0.04
Plasma OT	-0.07	0.07	0.35	0.00	0.07	0.97	0.06	0.03	0.05	-0.01	0.03	0.64
Plasma AVP	0.05	0.07	0.49	-0.05	0.07	0.47	-0.13	0.05	<0.01	-0.04	0.05	0.39

Table 5

Variable loadings from principal component analyses including variables related to locomotion and posture (all subjects) and social interaction (HAI condition only).

		Loadings
	Variable	PC1	PC2
Locomotion/Posture PCA (all subjects)	locomotion	0.62	-0.17
	upright	0.64	-0.16
	sitting	-0.02	0.83
	lying (prone)	-0.44	-0.50
Social interaction PCA (HAI subjects)	physical contact	-0.62	0.23
	play	0.55	0.32
	licking	0.32	0.71
	lying (supine)	-0.46	0.58

Table 6

Associations between changes in hormone concentrations in plasma and saliva and behavior during the test.

		Social interaction			Social interaction			Locomotion and Posture			Locomotion and Posture
		-PCI			-PC2			-PCI			-PC2
	Matrix / hormone	β	x2	p	β	x2	p	β	x2	p	β	x2	p
HAI group	Salivary OT	-0.60	1.55	0.21	0.61	4.37	0.04	0.72	2.70	0.10	0.22	0.60	0.44
	Salivary AVP	-0.02	0.00	0.97	0.16	0.22	0.64	0.10	0.04	0.84	0.18	0.26	0.61
	Plasma OT	-0.26	0.44	0.51	0.04	0.04	0.84	0.25	0.51	0.48	0.18	0.60	0.44
	Plasma AVP	-0.15	0.22	0.64	-0.49	6.15	0.01	0.21	0.52	0.47	-0.40	3.96	0.05
Control group	Salivary OT	-	-	-	-	-	-	-0.06	0.07	0.80	-0.05	0.03	0.85
	Salivary AVP	-	-	-	-	-	-	-0.36	1.80	0.18	0.20	0.35	0.56
	Plasma OT	-	-	-	-	-	-	0.12	0.89	0.35	0.13	0.58	0.45
	Plasma AVP	-	-	-	-	-	-	0.13	2.49	0.12	-0.04	0.15	0.70