Abstract Stress responses within dyads are modulated by interactions such as mutual emotional support and conflict. We investigated dyadic psychobiological factors influencing intra-individual cortisol variability in response to different challenging situations by testing 132 owners and their dogs in a laboratory setting. Salivary cortisol was measured and questionnaires were used to assess owner and dog personality as well as owners' social attitudes towards the dog and towards other humans. We calculated the individual coefficient of variance of cortisol (iCV = sd/mean*100) over the different test situations as a parameter representing individual variability of cortisol concentration. We hypothesized that high cortisol variability indicates efficient and adaptive coping and a balanced individual and dyadic social performance. Female owners of male dogs had lower iCV than all other owner gender-dog sex combinations (F = 14.194, p<0.001), whereas owner Agreeableness (NEO-FFI) scaled positively with owner iCV (F = 4.981, p = 0.028). Dogs of owners high in Neuroticism (NEO-FFI) and of owners who were insecure-ambivalently attached to their dogs (FERT), had low iCV (F = 4.290, p = 0.041 and F = 5.948, p = 0.016), as had dogs of owners with human-directed separation anxiety (RSQ) or dogs of owners with a strong desire of independence (RSQ) (F = 7.661, p = 0.007 and F = 9.192, p = 0.003). We suggest that both owner and dog social characteristics influence dyadic cortisol variability, with the human partner being more influential than the dog. Our results support systemic approaches (i.e. considering the social context) in science and in counselling.

Citation: Schöberl I, Wedl M, Beetz A, Kotrschal K (2017) Psychobiological Factors Affecting Cortisol Variability in Human-Dog Dyads. PLoS ONE 12(2): e0170707. https://doi.org/10.1371/journal.pone.0170707 Editor: Urs M. Nater, University of Marburg, GERMANY Received: August 6, 2016; Accepted: January 9, 2017; Published: February 8, 2017 Copyright: © 2017 Schöberl et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting Information file. Funding: The study was funded by the Austrian Science Fund (FWF; https://www.fwf.ac.at/en/): Project P23345 to KK and MW, and the PhD program Cognition and Communication from the Austrian Science Fund (FWF): Project W1234-G17 to KK. The funding agencies were not involved in the collection, analysis and interpretation of data. Competing interests: The project was funded by the Austrian Science Fund (FWF); Manuela Wedl and Iris Schöberl were employed via this project. The funding agency was not involved in the collection, analysis, interpretation of data and writing the paper. The authors confirm that the funding does not compete with the PLOS ONE policies on sharing data and materials. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Introduction One of the core issues in social relationships is stress coping [1]. For example, socio-positive interactions between dyadic partners may dampen physiological stress responses and promote relaxation [2,3] as a result of mutual emotional support [4–6]. Consequently, social animals show better recovery from aversive experiences when in the company of supportive partners than when alone [7]. On the other hand, social conflicts are among the most powerful stressors [8–10]. It has been shown that emotional social support also exists between humans and dogs [4,11–13]. This may be explained by shared and evolutionary conservative social mechanisms, including brain and stress physiology [14–18]. Dogs have socially adapted to live with humans for more than 30 000 years [19–24] and therefore, are closely adapted to human social performance. For example, dogs recognize and integrate emotional information from humans [25–29] and adjust their behavior to human communication signals [13] and human emotions [30,31]. There may even be emotional contagion between owners and dogs [32,33]. Stress coping, i.e. the behavioral, physiological and mental effort to master challenging situations [34] differs between individuals [35,36]. This may also be because mental representations influence the individual validation of potential stressors [18]. In fact, inter-individual variability in coping styles seems to be adaptive [37]. We suggest that parallel to heart rate variability [38], intra-individual cortisol variability reflects the adaptive flexibility of the hypothalamic-pituitary-adrenal (HPA) axis to respond to a variable environment [39]. This is based on findings that decreased cortisol variability is related to a down-regulated HPA axis by chronic stress [39], to bad health status, moderate well-being [40,41] and limited executive functions [42]. Hence, intra-individual cortisol variability over standardized situations may be an informative measure of stress coping, with a high amplitude between arousal peaks and relaxation lows reflecting healthy regulation. Individual differences in glucocorticoid variability are modulated genetically and epigenetically as well as by social experience [36,37,43]. Comparatively little is known about psychobiological factors influencing intra-individual variability in glucocorticoid response to different challenging social situations. Men emphasize different coping strategies than women, focusing on “fight and flight” instead of “tend and befriend”, which is the strategy preferred by females (Kudielka & Kirschbaum 2005). Men and women also differ in their interactions with, and attitudes towards, other animals [44]. This may also affect interaction styles with animals and the behavioral responses of these animals. For example, dogs seem generally more relaxed when petted by females than males [45]. Furthermore, male dogs of male owners are more sociable and active than male dogs of female owners [46] and male dogs of male owners had the lowest cortisol concentration compared to all other owner gender-dog sex combinations during an experimental assessment of their attachment [47]. This may be related to different behavior in humans depending on the gender of the interaction partner [48–51]. For example, same-gender pairs perform better during a computer task and show more verbal interactions than mixed-gender pairs [49]. Trust seems to be greater within same-gender pairs, than mixed-gender pairs [50]. However, other studies show that competition is greater and cooperation less pronounced within same-gender pairs than mixed-gender pairs [48,51]. With the shared social mechanisms of humans and dogs in mind, we therefore suggest that owner gender-dog sex combination may affect dyadic human-dog interactions and thereby affect stress coping. Gender is also linked to personality. For example, women score higher in Neuroticism (NEO-FFI dimension 1) [52–54] and Agreeableness (NEO-FFI dimension 4) [53], whereas men score higher in Extraversion (NEO-FFI dimension 2) [54]. Personality is connected to stress coping as well [34]. Humans high in Conscientiousness (NEO-FFI dimension 5), Agreeableness and Extraversion cope better with challenges in their lives than humans low on these dimensions [55]. Furthermore, Neuroticism is related to low expectations towards being supported [55] and to increased anxiety [55,56]. Interestingly, owner Neuroticism was found to be linked to low cortisol concentration in dogs [47,57]. A considerable proportion of individual HPA variability may be linked to relationships with partners, with distinct roots in early life history. Securely attached children [58,59] respond behaviorally and physiologically adequately to challenging situations and will quickly down-regulate behavioral and physiological responses after successful coping [60]. Along the lines of the parent-offspring model of human-dog relationships [57,61], the same principles also seem to work between dogs and their owners: We found that dogs classified as securely attached to their owner had a smaller cortisol increase during an attachment assessment and play situations with their owner than insecurely attached dogs [47]. Based on the same data set, we presently asked whether psychobiological factors influence intra-individual cortisol variability in owners as well as in their dogs. We hypothesized that high intra-individual cortisol variability is an indicator of efficient and adaptive coping, because a flexible HPA axis supports adequate responses to a variable environment, e.g. high cortisol variability is related to good health status and well-being [39–41]. Therefore, we predict that dyadic parameters which facilitate socio-positive interactions are related to high cortisol variability, e.g. owners with secure relationship patterns towards their dogs, but also towards other humans, may have dogs with high cortisol variability. Agreeable or extraverted owners are expected to show high cortisol variability, because of their competent stress coping. As male dogs of female owners are less sociable and relaxed than male dogs of male owners [46] we predict low cortisol variability over different challenging situations for male dogs and their female owners compared to all other owner gender-dog sex combinations.

Materials and Methods To investigate the influence of psychobiological factors on owner and dog intra-individual cortisol modulation, we tested human-dog dyads in a laboratory setting and measured their cortisol concentration before and after different challenging social situations. The current study follows a pilot study with only male dogs with their male or female owners [46,57,62] and formed part of an interdisciplinary project that investigated human-dog relationship and attachment. In previous papers [47,63] part of the methods we used were already described in detail and so we only report the overall design of the study and methods not sufficiently described elsewhere, and ask readers to refer to our previous papers for more details where applicable. Ethical approval Participation in our study was voluntary; dog owners were informed that they could stop the test situation at any time and also signed a written information and consent form. Data collection was conducted according to the standards of the Code of Ethics of the World Medical Association (Declaration of Helsinki), the EU Directive 2010/63/EU for animal experiments and Uniform Requirements for manuscripts submitted to Biomedical journals and the German Society of Psychology (Ethische Richtlinien der DGPs und des BDP). Ethical review for this study was conducted by the animal-welfare committee of the Faculty of Life Sciences, University of Vienna (approval number: 2014–015). Ethical approvement for the human participation in our study was granted by the Ethics Commission of the DGPs (Deutsche Gesellschaft für Psychologie, approval number AB 07_2011). Subjects and general procedure We collected data from 132 human-dog dyads. Intact pet dogs (mean age ± SD: 3.95 ± 1.83 years; mean weight ± SD: 29.53 ± 13.19 kg) and their primary attachment figure, the “owner” (mean age ± SD: 43.76 ± 10.71 years) were tested. All dogs had been adopted by their owners as puppies (mean age ± SD: 9.77 ± 3.50 weeks). Owner gender-dog sex combination was counterbalanced (35 female owners with female dogs, 35 male owners with male dogs, 31 female owners with male dogs and 31 male owners with female dogs). Different dog breeds and mixed breeds participated, but were limited to no more than three dogs of the same pure breed per owner gender-dog sex combination. Owner-dog dyads participated in two test sessions. A third session was scheduled with a subsample of 59 dyads [47], but is not considered in the present paper. The two sessions took about 90 minutes each and were scheduled in a test room at the University of Vienna. At the beginning of both sessions the general procedure was explained, which took about ten minutes. A camcorder (Canon Legria-HF-G10) with a wide-angle conversion lens (Canon WD-H58W) was fixed on the wall of the test room to videotape the two sessions. An experimenter observed the procedure from outside the room via a monitor, which was connected to the camcorder. During the first session, the focus was on the behavior of dog and owner in response to a novel room, separation from each other and reunion afterwards as well as to a challenging task. The aim of the second session was the measurement of heart rate and observation of behavior during a play situation and during two threatening situations. All tests are explained below. Furthermore, the owners completed questionnaires and took saliva samples from themselves and their dog during the sessions (Table 1). In addition, owners collected saliva samples during one non-test control day from themselves and their dog. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Table 1. Time schedule of the first and second session. https://doi.org/10.1371/journal.pone.0170707.t001 Questionnaires Owners completed questionnaires to assess human and dog personality, owner-dog relationship, as well as owner-dog interaction style and demographics. The questionnaires and the scales, which were revealed via principal component analysis (PCA), are the same as used in a previous study [47]. The NEO Five-Factor Inventory is a 60-item questionnaire designed to measure personality in adult humans in five domains: Neuroticism, Extraversion, Openness, Agreeableness, and Conscientiousness [52,64].

The Monash Canine Personality Questionnaire (MCPQ-R) [65] was translated to German in cooperation with a bilingual expert. A principal component analysis (PCA) from the translated questionnaire revealed five reliable axes (reliability via Cronbach-Alpha—given in the same order as the axes: 0.809, 0.837, 0.794, 0.777, 0.712): “Active-excitable”, “Obedient-reliable”, “Insistent-goal directed”, “Nervous-anxious” and “Cool-friendly”.

The Relationship Scales Questionnaire (RSQ) [66,67] captures different aspects of attachment via the scales “Separation anxiety”, “Closeness anxiety”, “Lack of trust” and “Wish to be independent”, which are associated with adult interpersonal relationship patterns [68].

The FERT (Fragebogen zu Erfahrungen mit Tieren) is a questionnaire in German on the relationship with companion animals [69]. The questionnaire was adapted for dog owners. PCA revealed four reliable axes (reliability via Cronbach-Alpha—given in the same order as the axes: 0.817, 0.812, 0.745, 0.624): “Dog as social supporter” and “Communication”, “Secure-caregiving” and “Insecure-ambivalent”.

Our self-developed interaction questionnaire included questions on owner interaction style towards the dog during daily life and revealed two reliable axes via PCA (reliability via Cronbach-Alpha—given in the same order as the axes: 0.773, 0.681): “Affirmative interaction style” (e.g. positive reinforcement such as giving treats, play, and verbal praise) and “Aversive interaction style” (e.g. positive punishment such as leash-jerk, telling the dog off, shaking the dog’s neck). Saliva samples Saliva samples were taken from owners and dogs to measure cortisol concentration. It is known that blood cortisol concentrations peak approximately 20 minutes after a dog encounters a stressor [45] and 20–40 minutes after a human encounters a stressor, depending on the duration and stressfulness of the task [70]. Thus, saliva samples from both owners and their dogs were taken before the tests and 15 minutes after the end of each test. Both sessions were scheduled at similar times in the afternoon within each dyad, mainly to avoid potential morning peaks in cortisol concentration [71]. Additionally, owners were asked to take saliva samples in the morning (8 am), at noon (2pm) and in the evening (8pm) during one resting day without testing (the calmest day of the week), between the first and second session to obtain baseline values. The experimenter instructed the owner how to take saliva samples from themselves and their dog at the beginning of the first session; all samples were taken by the owner himself/herself. Salimetrics oral swabs were used to obtain saliva samples from the owners and Salimetrics children's swabs for the dogs`samples. First the owner put the oral swab under his/her tongue and then he/she took the sample from the dog at the same time; the children's swab was inserted in the dog’s cheek pouch for about 90 seconds. Before taking the saliva sample the dog’s salivation was stimulated by smelling cheese. Although it is known that cheese does not interfere with measuring cortisol [72] the dog received a piece of cheese only after the saliva sample was taken. Each sample was frozen at -20°C (Sarstedt tubes) until analysis. Cortisol concentration was analyzed by enzyme immunoassay (EIA) [73]. Tests Picture viewing test. This test aims to collect information about the dog-human relationship and orientation toward each other. To activate the attachment system, the owner was asked to leave the dog alone in the novel room for three minutes immediately after entering the room. During this time the owner and experimenter (IS) observed the dog via the monitor outside the room. Then the owner entered the room again and looked at 14 dog pictures placed on the walls and windows. The dog was allowed to move freely within the room; the owner was asked to behave towards the dog as he/she normally would. During this test the experimenter was outside the room observing owner and dog via the monitor. The experimenter came back into the room as soon as the owner finished the task, or after eight minutes. Challenge task. The Challenge task aims to collect information on how owner and dog cope with practical tasks. The owner was asked to lead the dog as efficiently and safely as possible over a wire mesh bridge (height 0.6 m, length 5 m). Then the owner was asked to lead the dog onto a wobbly podium where the dog should remain for at least five seconds. The dyads were given four minutes to complete each of these tasks. Finally, the owner was asked to give the dog the commands “sit” and “lie down” with a maximum of two repetitions for each command. The owners were asked to interact with their dog as they usually do during daily life; owners were free to give the commands verbally, via gesture or both and were also allowed to use food if this was also part of the daily routine. Adaptation to heart rate monitor. The owner was asked to attach heart rate monitors to him/herself and his/her dog. This took about five minutes. Then the owner was asked to play with the dog for five minutes to adapt to the heart rate monitor. Toys as well as treats could be used. Afterwards the dog was on leash for five minutes and was asked by the owner to sit or lie down, but not to move. Heart rate and heart rate variability were measured from owners and dogs, but will be published separately and the technical details are described in another paper [63]. Mild threat. Two counterbalanced randomly ordered threat situations, one test with and one without the owner present, were conducted. For safety reasons, the dog was tethered to the wall with a leash during both threats. For the “mild threat with the owner present” the owner was instructed to behave in the same way he/she would do in similar situations in daily life. While the “mild threat without the owner present” was conducted, the owner could observe the situation from outside via a monitor. During both tasks the experimenter (IS) observed from outside to interrupt if necessary. The threat was conducted by another female experimenter, the “threatening person”, who wore a long black coat with a hood and a ski mask; just the eyes of the threatening person were visible. After entering the room, the threatening person knocked at the door from inside to get the dog’s attention. Then the threatening person stared at the dog and moved towards it. After each step the threatening person stopped for three seconds and kept staring at the dog. After the third step the threatening person turned away and left the room at the end of the first threat. The same procedure was followed for the second threat, but this time the situation was resolved at the end: the threatening person turned away from the dog, moved to the corner opposite to the dog, removed the ski mask and coat and approached the dog in a friendly manner while talking and offering it pieces of cheese. Then the formerly threatening person left the room and the dog was let off the leash. This test was described in detail in a former paper [47]. Data analysis Cortisol data were checked for outliers across the entire sample size of 132 dyads. Differences were found in absolute cortisol concentration as well as in standard deviations of cortisol concentration in male owners taking medication (n = 8) and male owners without medication (n = 58). Therefore, we excluded cortisol concentration values of those eight owners under medication from further analysis. Neither hormonal contraception nor medication use by women influenced cortisol concentrations. Cortisol concentration values of one dog were excluded from further analysis because of pseudo-pregnancy. Afterwards the remaining cortisol data were corrected for outliers by taking the mean cortisol value of all dog samples and separately, of all human samples, except the morning samples during the resting day. Error correction was done separately for the morning samples during the resting day, because humans exhibit a circadian pattern of cortisol with a cortisol peak in the morning [74, 75] and also in dogs, circadian patterns are found with differing morning cortisol concentrations [76–79]. All values higher than three standard deviations of the mean were excluded from further analysis [80], which was the case in 1.6% of all cortisol concentration values for dogs and in 2.3% of all cortisol concentration values for owners. The mean salivary cortisol concentration value ± standard deviation for session one and two over all owners was 6.85 ± 5.01 ng/ml; over all dogs 2.33 ± 2.71 ng/ml. The mean cortisol maximum ± standard deviation for session one and two over all owners was 10.75 ± 6.18 ng/ml; over all dogs 4.66 ± 4.05 ng/ml. The coefficient of variation of cortisol for each individual (iCV = sd/mean*100) was taken as an indicator of intra-individual cortisol variability in dogs and owners. Therefore, for each owner and each dog the mean value and standard deviation was calculated over all cortisol samples from the first and second session as a base for iCV. Only individuals with at least four applicable cortisol concentration values after error correction were taken into account for the calculation of the iCV. This resulted in 116 dogs (87.9%) and 115 owners (87.1%) being included in the analysis of iCV. A maximum of seven values was possible for owners and for dogs to calculate the iCV (mean number of values for iCV ± SD: Owners = 6.86 ± 0.51, Dogs = 6.55 ± 0.77). Statistical analysis was performed with SPSS 21 software. The Shapiro-Wilk Test was used to test for normal distribution. Data were not normally distributed; thus, nonparametric tests (Friedman, Mann-Whitney U, Wilcoxon and Spearman’s rank correlation) were used for univariate statistics. General Linear Models (GLM) were calculated for owner iCV and dog iCV. As owner and dog iCV was not normally distributed a log transformation was conducted to gain normally distributed data. Owner personality (NEO-FFI), dog personality (MCPQ-R), owner to dog attachment (FERT), owner relationship to other humans (RSQ), owner interaction style towards the dog (interaction questionnaire), as well as owner and dog age, were tested as covariates. Owner gender-dog sex combination was included as a factor. We selected these explanatory variables as main effects and removed them in order of decreasing significance (if p>0.1). Terms with p<0.1 were retained in the final model; to ensure that excluded terms did not explain the variance they were re-entered one by one into the final model [81]. Only terms with p<0.05 were considered a significant influence on the dependent variable. All significances (p<0.05) are given two-tailed. As described in a previous paper [47], alpha correction for multiple comparisons was not considered here because this generally increases the risk of type-II error at a comparatively low potential of decreasing type-I error [82]. Instead, Cohen’s effect size was calculated for univariate statistics, with d = 0.2 indicating a small effect, d = 0.5 a medium effect and d = 0.8 a large effect [83]. Because non-significant results do not mean that there is no effect [84] we also include tendencies (p<0.1) within this paper.

Discussion We investigated cortisol variability in human-dog dyads and hypothesized that high cortisol variability may be an indicator of efficient and adequate coping strategies. As expected, high cortisol variability was related to low absolute cortisol concentration values in owners and partly also in dogs. Furthermore, high cortisol variability was related to cortisol decrease in owners during the PVT and in dogs during the second threat, which supports our hypothesis. Cortisol variability (iCV) seems to follow the same principle as heart rate variability (HRV); heart rate (HR) also correlates negatively with heart rate variability. High HR indicates arousal, whereas high HRV indicates relaxation, good well-being and health status. HRV reflects the ability of the heart to adapt to changing circumstances [38]. Decreased cortisol variability is, similar to low HRV, related to bad health status and well-being [39–41]. This is supported by the psychobiological factors affecting dyadic HPA modulation, notably dyadic sex/gender combination, personality and attachment style. Owners scoring high in Neuroticism had dogs with low cortisol variability. Neuroticism is linked to low expectations that social support will be forthcoming when needed [55], to major depression [85–87] and to anxiety [55,56]. Dogs are sensitive to their owners´ emotional states [25–29] and emotional contagion between owners and dogs is possible [32,33]. Thus, dogs may mirror the anxiety and negative expectations of neuroticistic owners in their cortisol variability. Interestingly, cortisol increase after stressful stimuli and also morning cortisol concentration are low in dogs from owners scaling high in Neuroticism [47,57]; but the opposite effect was found for cortisol variability in dogs. Cortisol reactivity may be more influenced by immediate interactions, e.g. dogs of neuroticistic owners approach them more often [62], which may dampen the immediate stress response. In contrast iCV, reflecting the regulatory capacity of the HPA axis, may be more integrative over time and contexts. Highly agreeable people are known to cope better with the challenges of daily life [55] than less agreeable ones. Agreeable owners are probably good at coping with social situations and hence at controlling social stress. This personality dimension is described as being cooperative, sympathetic and considerate. Thus, agreeable owners had higher cortisol variability then less agreeable ones, indicating efficient stress coping. Dog personality, as rated by owners, did not significantly explain cortisol variability in dogs and owners. In fact, owner characteristics seem to be more relevant in this respect [47,57,88] than just the characteristics of the dog itself. Still, perceiving one´s own dog as cool and friendly may feed back to the psychophysiological variability in owners, as it is probably easier to handle daily life with a relaxed dog, which in turn depends mainly on the performance of the human partner. In general, relaxed owners are likely to have relaxed and friendly dogs. As predicted, owners’ insecure-ambivalent attachment towards their dog was related to low cortisol variability in dogs. Attachment classification is known to predict sensitivity and reliability as a caregiver [89,90]. Sensitive caregivers interact thoughtfully with the child and are aware of the child’s needs [89,90]. Thus, owners low in sensitivity may be more likely to have dogs with an insecure attachment style than owners with a sensitive caregiving style [91]. Dogs of the former are more stressed during an attachment assessment than securely attached dogs, and their owners are more likely to be insecurely attached towards them [47]. Our findings confirm previous results that the owner’s relationship towards other humans is reflected in the owner-dog relationship [92]. Attachment-related interaction styles, such as contact-seeking, are displayed in interactions with the dog in the same way as they are expected to be in interaction with humans [93]. Thus, the attachment representation with the primary human caregiver may also be transferred to the dog’s cortisol variability. While gender is important, particularly the dyadic gender combination appears to influence social interactions in humans [48–51]. Men tend to show comparatively higher cortisol increases when their social prestige is challenged (e.g. by performing in front of an audience), whereas experiencing social rejection seems to be a more powerful trigger for HPA responses in women [94]. Furthermore, women report more fear and less happiness than men after a social stress test [95]. Men and women differ not just in their social emphases and coping strategies, but also in their interactions with, and attitudes towards, animals [44]. Women show stronger emotional relationships with companion animals than men do [46,96] and girls seek more contact with animals than boys [97]. During domestication, dogs became highly adapted to living with humans; therefore minor variations in the owners´ interaction styles may have distinct effects on the dogs´ physiological and behavioral responses (45). Dogs can probably discriminate human gender and may adapt their behavior according to the owner gender. For example, male dogs of female owners are less sociable and relaxed than male dogs of male owners [46]. This seems to be reflected in our result where female owners with male dogs were shown to have the lowest cortisol variability compared to all other dyadic gender combinations. As expected, we found considerable individual cortisol variability for owners and dogs over the different experimental challenge situations. Interestingly, a decrease of cortisol over both sessions could be observed. This may reflect the recovery from travel and apprehension of a novel situation with subsequent relaxation after the initial phase; pretest baseline samplings are sensitive to the novelty of the setting and an adaptation to this setting may lead to lower samples afterwards [71]. Furthermore, we did not find a clear increase in cortisol in owners and in dogs during most of the tests, probably because they were predictable and controllable. The owners’ and dogs’ emotional support for each other could also be an explanation for the decrease of cortisol in some individuals during both sessions. For logistic reasons we were not able to control how dyads reached the lab and whether this may influence cortisol concentrations. But we figure that this was in alignment with the routines of each dyad and we do not consider the lack of this kind of control as crucial for interpreting the outcome of our study. In future studies a longer initial phase to adapt to the lab and a control group of dogs being regularly at the lab could prevent this limitation. We also acknowledge that a still higher sample size per group would increase statistical power and may lead to more distinct results on group level. We conclude that individual cortisol variability over different challenging situations may be indicative of individual adaptation in humans as well as dogs. High intra-individual cortisol variability could be an indicator of efficient and adaptive coping, because a flexible HPA axis supports adequate responses to a variable social environment. We suggest that both owner and dog social characteristics influence dyadic cortisol variability, with the human partners being more influential than the dog. These findings underscore the importance of considering the human-dog dyad on a systemic level, i.e. including the social context in experimental science as well as in dog training and dog behavior therapy.

Acknowledgments We would like to thank all owners and their dogs for their participation in this study. Furthermore, we would like to thank all students who were involved in our project: Evi Myska, Sarah Prettner, Philipp Stöger, Sigrid Amon and Verena Ziemen. We also thank Matthias Nemeth and Teresa Schmidjell for their statistical support and Gesche Westphal-Fitch for critical revision and proof-reading.

Author Contributions Conceptualization: IS MW AB KK. Data curation: IS. Formal analysis: IS. Funding acquisition: KK MW. Investigation: IS. Methodology: IS MW AB KK. Project administration: IS KK MW. Resources: KK. Supervision: KK AB. Validation: KK. Writing – original draft: IS. Writing – review & editing: IS MW AB KK.