Physiological stress reactivity and recovery related to behavioral traits in dogs (Canis familiaris)
Rian C. M. M. Lensen aff001; Christel P. H. Moons aff002; Claire Diederich aff001
Authors place of work:
Department of Veterinary Medicine (IVRU), University of Namur, Namur, Belgium
aff001; Department of Nutrition, Genetics and Ethology, Ghent University, Merelbeke, Belgium
Published in the journal:
PLoS ONE 14(9)
This study investigated whether stress responsiveness (in one context) can be used to predict dog behavior in daily life. On two occasions (NT1 = 32 puppies; NT2 = 16 young adults), dogs’ physiological stress response after a behavioral test at home was measured in terms of reactivity (10 min post-test) and recovery (40 min post-test) for three salivary markers: cortisol, chromogranin A (CgA) and secretory immunoglobulin A (sIgA). For each marker, it was determined whether dogs with a strong physiological response displayed different behavior in daily life compared to dogs with a weaker physiological response. The results revealed three main findings: first, for CgA and cortisol, different patterns were identified according to sample time. High reactivity related to desirable traits, whereas slow recovery after the behavioral test related to undesirable traits. The findings suggest that increased levels of CgA and cortisol 10 minutes after the behavioral test reflected an adaptive stress response, whereas elevated levels 40 minutes after the test reflected unsuccessful coping. Second, patterns for sIgA differed from CgA and cortisol: significant associations were only found with behavioral traits at T2, mostly considered desirable and related to Trainability. Possibly, the delayed reaction pattern of sIgA caused this difference between markers, as sIgA reflects the (secondary) immune response to stress, due to immunosuppressive effects of cortisol. Third, predictive capacity of puppies’ physiological stress response (T1) was inconclusive, and contrary relations were found with behavioral traits at T2, suggesting that developmental factors play an important role. This study provides new insights about the relation between stress physiology and behavioral traits, and methodological advice is given to study these patterns further. In conclusion, physiological markers could provide additional insights in dogs’ tendencies to display certain behaviors, especially at the young adult stage. Further studies are needed to confirm these patterns.
Biology and life sciences – Organisms – Eukaryota – Animals – Vertebrates – Amniotes – Mammals – Dogs – Animal types – Pets and companion animals – Biochemistry – Hormones – Lipid hormones – Cortisol – Steroid hormones – Psychology – Behavior – Animal behavior – Zoology – Anatomy – Body fluids – Saliva – Physiology – Social sciences – Medicine and health sciences – Immune physiology – Immunology – Immune response – People and places – Population groupings – Age groups – Young adults
Consistent individual-specific patterns in behavior and stress physiology have been reported in various animal species, but the relation between both types of parameters remains inconclusive (reviewed in [1,2]). It has been suggested that reactivity (i.e., an animal’s responsiveness to environmental stimuli [3,4]) is a fundamental factor shaping personality [5,6] and plays a role in various behavior problems in dogs (Canis familiaris) . In theory, a reactive animal responds more quickly and intensely to environmental stimuli, which can be expressed by different behavioral manifestations (e.g., excitability, fear or aggression [2,7]). Accordingly, one might expect that assessment of a dog’s stress responsiveness could be used to predict its behavior in daily life. If this were true, an objective, quantifiable assessment of stress responsiveness could be a valuable indicator for the potential to develop behavior problems. An early detection of individuals “at risk” allows for timely measures to be taken to prevent behavioral problems.
To test whether stress responsiveness can be used to predict general behavior in daily life, the physiological stress response of dogs in one context should be compared to their behavioral responses in another unrelated context (e.g., ). However, studies on dog personality often focus on behavioral parameters (e.g., [9,10]), whereas physiological measures are more commonly used for assessing levels of stress or stress responsiveness in specific contexts (e.g., [11,12]). Furthermore, physiological parameters that are included in personality studies (i.e., assessing consistent patterns over contexts and time ) are generally taken on the same occasion as the behavioral parameters (e.g., [14,15]). In the current study, we assessed whether phyiological measures of stress reactivity in one context related to behavioral traits in daily life, as assessed by dog owners.
Saliva sampling is a relatively easy, non-invasive method that has often been used in dogs to monitor short-term physiological changes . Cortisol is a well-known stress marker in dogs and reflects the activity of the hypothalamic-pituitary-adrenal (HPA) axis . Saliva concentrations of this marker increase in response to sudden non-social stressors , fear-inducing events , and novel environments . However, physiological markers are regulated through complex pathways [20,21]; and their expression could be influenced by factors like circadian fluctuations , physical health , activity , as well as stimulus intensity  and experience . Hence, the measurement of multiple phsyiological parameters is recommended for a more accurate estimation of stress levels [17,26].
A more immediate response (“fight or flight”) to a perceived or actual threat is regulated through the sympatho-adrenal-medullary (SAM) axis through the release of catecholamines (adrenalin and noradrenalin). However, these markers are difficult to measure in saliva due to their low concentrations, rapid degradation, and instability in the sample  and are found to be poor indicators of acute sympathetic activity in humans . Instead, Chromogranin A (CgA) can be measured, as it is co-released with adrenalin from the adrenal medulla and more stable in the circulatory system than catecholamines [27,29]. In dogs, plasma levels of CgA increased after physiological stress (insulin-induced hypoglycemia ); and increases in salivary CgA were found after psychogenic stress in humans (arithmetic task ), pigs (immobilization with nose snare ), and cows (social isolation ). However, little is known about this analyte as potential stress marker in dogs , so it was included in our study to explore its value for future studies.
Components of the immune system can also provide insights into the physiological stress response, due to the immunosuppressive effects of cortisol [23,27]. As such, markers of immune system activity provide an indirect (secondary) indication of an animal’s stress response. Decreased levels of secretory immonoglobulin A (sIgA) were found in dog saliva after acute stress  and 10 days after introduction into a novel kennel . Therefore, this analyte was included to assess its added value to the abovementioned markers for assessing the physiological stress response.
The aim of this exploratory study was to examine if physiological indicators of stress responsiveness are related to behavioral traits in dogs as reported by the owners. For this purpose, three endocrine and immune biomarkers (cortisol, CgA, sIgA) were measured in saliva collected before and at two time points (10 and 40 min) after a behavioral test and compared to owner ratings of their dog’s behavior in daily life (by means of the validated Canine Behavioral Assessment & Research Questionnaire: C-BARQ ). This was done at two test stages, for the same dogs: puppy (T1) and young adult (T2). We hypothesized that a strong physiological response would be indicative of undesirable behavioral traits, based on findings that saliva levels of cortisol increased [18,19] and levels of sIgA decreased  after dogs were exposed to (potentially) aversive stimuli. Furthermore, we hypothesized that dogs with a longer recovery time (i.e., time before levels of physiological markers returned to baseline) would be less successful in coping with stressful situations and thus might be more prone to display undesirable (“maladaptive”) behaviors [7,36].
Material and methods
Owners were recruited by participating veterinarians (N = 20) during their puppies’ first vaccination visit at the veterinary practice over a 21-month period (Table 1). The owners received two months of free dog food from a well-known brand upon completion of their participation, and the veterinarians received a gift voucher for an online general goods store.
At test stage 1 (T1), saliva was collected from 32 puppies (mean age ± SD: 17.57 ± 2.94 weeks) subjected to a behavioral test at home, and owner ratings were available for 31 of these puppies (procedures described below). At test stage 2 (T2), 9.91 (± 0.61) months later, 16 dogs were tested a second time at home (age: 60.85 ± 2.64 weeks), and owner ratings were available for all 16 dogs.
At both stages (T1 and T2), owners were asked to rate the behavior of their dog in daily life by means of an online questionnaire (C-BARQ: updated version of PennBARQ ) in Dutch or French. The C-BARQ describes 14 different behavior categories, each represented by the mean of a subset of questions (presented in [37,38]): Stranger-directed aggression (ten questions), Owner-directed aggression (eight questions), Dog-directed aggression (four questions), Stranger-directed fear (four questions), Dog directed fear (four questions), Non-social fear (six questions), Separation-related behavior (eight questions), Attachment and attention-seeking (six questions), Trainability (eight questions), Chasing (four questions), Excitability (six questions), Touch sensitivity (four questions), Energy level (two questions) and Dog rivalry (four questions). The higher the score, the more present a particular trait is in a dog according to the owner.
Owners completed the questionnaire before the test to ensure their dog’s behavior during the test would not influence their ratings. They were instructed to answer all questions but could skip a question pertaining to a situation if their dog had not experienced it yet. For any individual dog, if more than 25% of the questions for a category were missing, it was considered a missing value.
Saliva collection & analysis
At both test stages (T1 and T2), saliva samples were collected before (pre-test), as well as 10 (post10) and 40 minutes (post40) after a behavioral test (described in S1 Appendix). In short, the first 16 minutes of this test consisted of different social situations (e.g., entrance of a stranger, being ignored by persons in the room, being left alone), and in the second part (27 min), the dogs were exposed in a random order to several stimuli with varying visual and acoustic properties (e.g., sound of a barking dog, squeaky toy, weasel ball, vacuum cleaner). Two startling stimuli were presented at the end (umbrella opened suddenly and shaking of a metal cylinder containing keys), so that the overall test remained as little stressful as possible to the dogs.
The test was conducted in a familiar room at home to minimize influences of unfamiliar settings on the stress response [19,39]. To minimize influences from circadian fluctuations, testing only took place at one of two predetermined time points (10:00 or 14:00), and, on both occasions (T1 & T2), the dogs were tested on the same day (i.e., with a similar routine).
Saliva samples were collected according to a standardized protocol that was refined over several pilot trials . Owners habituated their dogs to the sample procedure, starting four days before the test (S1 Appendix), and collected the saliva samples to prevent influences of unfamiliar persons [19,40]. On the test day, owners withheld food from the dogs at least one hour before sampling and did not let the animal drink or engage in strenuous activities at least half an hour beforehand to prevent contamination of saliva samples and influence of physical activity on the physiological measures [41,42]. At each sampling occasion (pre-test, post10, post40), the owners put on clean latex gloves and collected saliva samples with four swabs. They put one swab in each cheek pouch and gently held the dog’s muzzle for one minute, after which they placed the two swabs back in their centrifugation tubes. Immediately after, they took samples with the two remaining swabs. The owners were instructed not to pet the dog during sampling, as this could influence the physiological response . The dogs were not restrained, and sampling was stopped if an animal appeared overly stressed by the procedure (e.g., low posture, avoidance, resistance). Sampling occurred within four minutes to prevent handling from influencing marker concentrations . Owners noted sampling time, sampling duration, last eating/drinking time, activity before sampling, and any additional remarks. They then placed the four centrifugation tubes, containing one swab each, in the freezer until the end of the behavioral test.
After the test, the samples were transported on ice to the tester’s home and stored within two hours (-20°C). Within three months , samples were transported on ice to the University of Namur, where they were centrifuged at 3000 rpm (1851.41 x g for 15 min at 10°C) and samples containing blood or other visible contamination were discarded (N = 61 or 11.30%; from an additional 32 samples [5.92%] no saliva could be extracted). The extracted saliva volume per swab was noted, then saliva was pooled per sampling occasion (four swabs), resulting in three samples per test day: pre-test, post10, post40. The samples were assayed for the following three analytes, given sufficient saliva was available (S1 Table): CgA, cortisol, sIgA (respective kits: Human Chromogranin A EIA kit, Yanaihara Institute, Shizuoka, Japan; Salivary cortisol enzyme immunoassay kit, Salimetrics LLC, State College, Pennsylvania; Dog IgA ELISA Quantitation set, Bethyl Laboratories, Montgomery, Texas). Samples had to be diluted to assay CgA (1:8) and sIgA (1:3000), so the required sample volumes for duplicate measurements were 8 μl, 50 μl and 3 μl for CgA, cortisol and sIgA, respectively .
To minimize assay influences, all samples from a dog (collected on one test day) were assayed on the same ELISA plate, and at both test stages the same dogs were assayed together on a plate. For example, if at T1 the ELISA plate contained samples from dogs 1–4 then the ELISA plate at T2 also contained samples from dogs 1–4. Intra-assay coefficients of variation (CV) were calculated by dividing the standard deviation of duplicate measurements of a sample by the mean of these measurements. Samples with an intra-assay CV > 15% were excluded from further analysis (N = 5, 0 and 16 for CgA, cortisol, and sIgA, respectively), so that analytical imprecision was less than half of the within-subject variation . Mean intra-assay CVs for the final samples were 4.56% (CgA), 5.16% (cortisol) and 4.58% (sIgA). Inter-assay CVs, based on repeated measurements of two control samples (aliquoted before freezing) were relatively high: 67.90% (CgA), 11.60% (cortisol), and 41.06% (sIgA). However, over the whole assay period (almost two years), the control samples were frozen longer than the recommended 6 months , so instability of the analytes might have contributed to this high variance. Acceptable inter-assay CVs of 9.31% (CgA), 4.98% (cortisol), and 6.38% (sIgA) were found for the same assay kits during pilot trials . Each dog served as its own control, and all samples from that dog (one test day) were assayed on the same plate, so inter-assay variation is expected to minimally influence the measured values for each dog (expressed as changes in value over time, see “statistics”).
The described procedures (behavioral test, saliva sampling) were approved by the Ethical Commission of the University of Namur, which did not consider them animal experiments based on the non-invasive nature (as attested by seven international researchers, experienced with these types of measurement). The protocol was agreed upon by the Belgian State authority for Directive 2010/63/EU.
The owners were informed about the study aims and procedures and signed a letter of informed consent. Their participation was voluntary, and they were told that they could interrupt the behavioral test at any moment. Also, endpoints were formulated for the behavioral test (see S1 Appendix) to protect the welfare of the dogs.
As high inter-individual variation in marker concentrations may confound interpretation of reactivity and recovery parameters , we analyzed differences between time points, rather than absolute concentration values. Pre-test values were subtracted from post10 and post40 values  to create two variables: Δ10, Δ40.
Given the low sample size (see S1 Table), and in consultation with a statistician, non-parametric analyses were used to identify potential patterns between physiological and behavioral parameters. Per marker and for each variable, the dogs were median-split to distinguish those with the largest physiological response from those with the smallest response. Individuals were considered to have a “strong response” if their saliva concentrations increased above median for cortisol and CgA or decreased below median for sIgA . Scatter plots were used to verify if the median split the sample into two distinctive groups. For three variables, a dog was moved from the “weak response” to the “strong response” group (T1: CgA Δ40; sIgA Δ10) or conversely (T2: CgA Δ10) because their values resembled more closely to the newly assigned group. The first author then used a Mann-Whitney U Test (IBM® SPSS® Statistics, Version 20) to determine whether these groups differed significantly (α < 0.050) for the different C-BARQ scores at T1, T2 and when comparing physiological responses at T1 with behavioral scores at T2.
The following sections describe significant associations. To minimize reporting bias , all tested associations and ranges in marker concentrations are reported in a supplemental table (S2 Table).
Associations between physiological and behavioral variables at the same test stage
For CgA and cortisol, significant associations were found with behavioral traits at both test stages (Table 2). Puppies with a strong physiological response 10 minutes after the behavioral test (Δ10) received lower scores for Separation-related behavior (CgA), whereas those with a slow recovery (Δ40) received lower scores for Trainability (CgA) and higher for Dog Rivalry and Stranger-directed fear (cortisol). These results suggest that dogs with a strong physiological response shortly after the test were less likely to display undesirable behaviors upon separation. This finding is in contrast with our hypothesis, as we expected more undesirable behaviors for dogs with a strong physiological response (both for Δ10 and Δ40). Indeed, we found that a strong physiological response 40 minutes after the test was associated with more undesirable behaviors in daily life, as reported by the dog owners.
At T2, similar patterns emerged for CgA and cortisol. High reactivity (Δ10) was associated with desirable behaviors in daily life: low Stranger-directed fear (CgA) and high Trainability (cortisol). In contrast, a slow recovery (Δ40) was associated with low Excitability (CgA) and high Energy level (cortisol), which might (arguably) be considered undesirable by owners (see discussion).
In contrast to the primary stress markers cortisol and CgA, significant associations between sIgA levels and behavior were only found at the young adult stage (T2), and only for desirable traits. A strong physiological change in sIgA levels, both at 10 and 40 minutes after the test, was associated with the desirable trait Trainability. In addition, slow recovery after the test (Δ40) was associated with low levels of Non-social fear and Touch sensitivity.
Associations between physiological variables at T1 and behavior at T2
Ambiguous associations were found between puppies’ stress response, measured by levels of cortisol and CgA at T1, and their future behavioral traits (T2). Puppies with high physiological reactivity (Δ10) at T1 received higher scores at T2 for Stranger-directed fear, and Excitability (CgA) and lower scores for Attachment and attention-seeking (cortisol; Table 3). Those puppies with sustaining levels of CgA after the test (Δ40) received higher scores for Excitability and Energy level later in their life.
For sIgA, again a positive association was found between large physiological reactivity following the test (Δ10) and the desirable trait Trainability. In contrast, a slow recovery (Δ40) at the puppy stage (T1) was associated with more Dog-directed fear later in life (T2).
This study investigated whether stress responsiveness, as measured via physiological markers in one context were related to dogs’ general behavior in daily life (as rated by their owners). At two stages (T1: puppy, T2: young adult), the dogs’ physiological stress response after a behavioral test at home was measured in terms of reactivity (10 min post-test) and recovery (40 min post-test) for three salivary markers: cortisol, chromogranin A (CgA) and secretory immunoglobulin A (sIgA). Based on these values, it was determined whether dogs with a strong physiological response (for each marker separately) displayed different behavior in daily life, compared to dogs with a weaker physiological response. We hypothesized that a strong physiological response 10 minutes after the test would be associated with undesirable behavioral traits, and we expected the same for a longer recovery time (i.e., time before physiological levels return to baseline). The results revealed three main findings: 1) for CgA and cortisol, different patterns were identified according to sample time post-test, 2) patterns for sIgA differed from CgA and cortisol, and 3) predictive capacity of puppies’ physiological stress response (T1) was inconclusive, and ambiguous relations were found with behavioral traits at T2.
First, for CgA and cortisol, high reactivity at 10 minutes after the behavioral test related to desirable traits, whereas slow recovery 40 minutes after the test related to undesirable traits. The former result suggests that dogs with a relatively strong physiological response following the test respond to stressors in an adaptive manner [23,49], reflected in daily life by less fear and higher trainability. In contrast, the latter result might indicate unsuccessful coping (slow recovery), as reflected by the associations with undesirable traits. This could also explain the high Energy (cortisol) and low Excitability (CgA) scores at T2. Increased activity can be a mechanism to cope with prolonged stress [24,50], and, if the dogs were in a negative emotional state, they may have responded less excitedly to the situations described in the C-BARQ score Excitability (e.g., going for walks, arrival of visitors, owner coming home).
The finding that associations of physiological responses with behavioral traits vary according to the sample times indicates the importance of measuring changes in marker concentrations, at multiple time points. Collection of baseline samples before the onset of a specific event allows for a correction of inter-individual variation in marker concentrations  and thus provides a more accurate indication of each dog’s physiological response to that event. This is supported by findings from Wormald et al. , who only found a significant difference between dogs that either passed or failed an aggression test, when considering the change in salivary cortisol levels (pre-post venipuncture). Likewise, Sherman et al.  found only significant associations between “emotional reactivity” (measured in a behavioral test) and cortisol levels in saliva and plasma when considering the change compared to baseline but not when using absolute post-test values.
Furthermore, the collection of multiple samples is necessary to monitor changes over time, thereby distinguishing dogs with a high (initial) reactivity from dogs with a more blunted (prolonged) stress response. In fact, we hypothesized that a strong physiological response (in itself) would be indicative of undesirable behavioral traits [2,19]. Instead, our findings suggest that the recovery time tells us more about a dog’s overall stress responsiveness. Possibly, dogs that are more prone to stress in daily life have frequent elevations of stress markers, which might have led to an attenuated stress response due to negative feedback . In line with this, Beerda et al.  found that dogs (N = 8) subjected to 5 weeks of social and spatial restriction, showed an increase over time in daily cortisol levels but had an attenuated cortisol response when exposed to a noise stressor. Dreschel and Granger  found that dogs living with other dogs had a lower increase in cortisol levels 20 and 40 minutes after exposure to a fearful event (thunderstorm recording), and their baseline cortisol levels tended to be higher. They argued these findings might reflect stress in daily life from living with another dog. Likewise, Rosado et al.  found higher plasma levels of cortisol in dogs diagnosed as aggressive, which they attributed to stress most likely resulting from an inconsistent and unpredictable living environment. Translating these findings to our study, it is possible that relatively small changes in marker concentration shortly after the behavioral test reflected an attenuated stress response to that specific event.
Second, this study revealed differences among physiological markers. Contrary to findings for cortisol and CgA, changes in sIgA levels only related to behavioral traits at T2, most of which were considered desirable. In particular, Trainability was positively associated with reactivity and recovery at T2 and also predicted by high reactivity at T1. Likewise, Kikkawa et al.  found sIgA levels to differ significantly between dogs (11–14 months old) selected (N = 25) or rejected (N = 44) for guide dog training. During a 2-week assessment period, the dogs were kenneled and their sIgA levels were measured on days 1, 2, 3, 7 and 14. The sIgA levels of subsequently selected guide dogs gradually increased, whereas sIgA levels of rejected dogs remained relatively low, despite similar housing and caring conditions. On day 14 the selected dogs had significantly higher sIgA levels than rejected dogs, which the authors interpreted as a difference in “adaptive ability”. These findings seem to reveal the same pattern as described above for cortisol and CgA (Δ40): consistently low sIgA levels suggest that the rejected dogs had low adaptive capacity, with slow recovery during a period of prolonged stress (kenneling). In our study, however, a strong response (i.e., relatively low levels of sIgA) was associated with high Trainability scores, which appears to be more in line with the Δ10 results for CgA and cortisol. These contrary findings are likely a result of different protocols: in our study the samples were collected after startling stimuli (two at the end of the test), whereas Kikkawa et al.  took samples at predetermined time points, not associated with a particular stimulus. Consequently, sIgA measures in our study appear to reflect the more immediate response to a stressor, in contrast to the long-term measures of Kikkawa et al. Indeed, in an extensive review on stress and secretory immunity in humans, Bosch et al.  identified an overall distinction between acute and chronic stress: chronic stress generally induces a decrease in sIgA levels, acute stressors mostly induce an increase. One would then expect the dogs from our “strong response” group to be more prone to chronic stress, and thus have a lower trainability (as found by Kikkawa et al. ). On the other hand, if we take into account a possible attenuation in cortisol response for dogs under chronic stress (as described above), it is possible that this results in less immunosuppression, hence relatively higher sIgA levels, in these dogs after an acute stressor. If this were true, then this would explain why the “weak response” group had relatively low Trainability scores. However, many factors affect the findings (including type of stressor, timing of sampling , and the balance of sympathetic-parasympathetic control [27,54]. Consequently these complex interactions need to be further explored and elucidated. Nevertheless, it is noteworthy that we found three associations with the Trainability trait that were consistent over time, and these findings along with those from Kikkawa et al.  suggest that sIgA holds potential as a marker of (long-term) coping ability. Further investigation is needed to unravel individual differences in sIgA levels, both in response to short-term and long-term stressors.
The differences between CgA and cortisol on one hand and sIgA on the other hand are in line with findings from other studies, and distinguishes them as primary and secondary markers of stress. Similar reaction patterns were found for CgA and cortisol in dog plasma , pig saliva  and human saliva . Svobodová et al.  reported no significant association between salivary cortisol and sIgA values after a stressful condition, though Skandakumar et al.  reported a negative (semi-logarithmic) relation between cortisol and sIgA levels in repeated saliva samples from six dogs. Possibly these contradictory findings reflect the delayed reaction pattern for sIgA, in comparison to cortisol: in response to a stressor cortisol is first released, which subsequently inhibits components of the immune system , thereby decreasing sIgA levels. So, strictly speaking sIgA is a marker of immune response, which responds to stress, instead of a primary stress marker. The delayed reaction pattern for sIgA could explain the correlations with desirable behavioral traits and sIgA measured 40 minutes after the behavioral test (as opposed to cortisol and CgA). After a stressor sIgA levels return to baseline only after 30–60 min . Consequently, 40 minutes after the test might have been too soon to detect differences between dogs in coping capacity (recovery), as we found for CgA and cortisol. Unfortunately, combined analyses with all markers (e.g., pairwise comparisons; Principal Component Analysis) were not possible due to incomplete datasets for the individual markers (S1 Table). So, more research is needed to further elucidate the interaction between different stress markers.
Third, predictive capacity of stress physiology at the puppy stage (T1) remains inconclusive. Apart from Trainability, sIgA responsiveness in puppyhood only related to an undesirable behavioral trait at T2 (Dog-Directed fear), in contrast with the findings described above. Only one association was found for cortisol, and contradictory associations were found for CgA: Stranger-directed fear at T2 was associated with high reactivity (Δ10) at T1 but with low reactivity at T2. Likewise, Excitability at T2 was associated with slow recovery (Δ40) at T1 but with fast recovery at T2. These contradictions could partly be explained by the different datasets at both test stages, as the findings describe relative patterns for each test population . For instance, a dog’s stress response may have been relatively strong compared to the other dogs at T1, but the same response may have been relatively weak at T2. Nevertheless, it is interesting that associations with the same behavioral traits were found at both test stages, even though inverted. Probably the relation with behavioral traits is shaped by developmental factors [43,58], and further studies should elucidate the predictive capacity of the described salivary markers.
One important limitation of this study is the small sample size. Data were inevitably lost during laboratory analyses (insufficient or contaminated saliva samples) and by excluding samples with an intra-assay CV > 15%. Effect sizes were relatively high, but this is partially caused by the small sample size (i.e., the denominator of the probabilistic index). Furthermore, potential influences like breed, sex, or neuter status [40,59] could not be analyzed, so the results should be interpreted with caution. Though the patterns identified in this study should be confirmed, they highlight important considerations for future studies.
Another limitation of this study is the high dependence on dog owners. Participation rate was low for this longitudinal and multidisciplinary study, despite proactive recruitment by veterinarians and offering incentives. Furthermore, not all owners adhered equally well to the specific test and sampling procedures. This may have caused inter-individual variation and possibly also some intra-individual variation. By median-splitting the group, the influence of such variation was expected to be minor: individuals with a large stress response would fall into the “strong response” group regardless of the exact change in marker concentration and vice versa. However, the use of the median was an arbitrary choice, and the created groups might not reflect ‘true’ differences in stress response. Therefore, the data were visually explored to ensure they were split into two distinctive groups.
Furthermore, dogs that would normally cope by means of owner interaction, activity or play had limited options for recovery, as interaction and activity were limited for standardization. Nevertheless, the saliva samples still provided valuable information about their coping ability when owner support or outlet options are limited. For interpretation of the results, standardization of the test situation was thus preferred. For future studies, it might be useful to note signs of stress at the time of sampling, instead of merely activity, and to question owners how their dog usually copes in a stressful situation. This would allow retrospective inspection whether or not the dogs recovered, and if not, whether this might have been caused by the absence of relevant maintenance stimuli . Alternatively, the potential value of hair samples to determine dogs’ general stress responsiveness could be further explored. Although hairs contain a mixture of analytes accumulated over weeks, it might give an estimation of dogs’ stress reactivity, when taking into account activity levels and coat color [61,62].
Another potential influence concerns the baseline levels before the test. The owners met the tester outside 20–30 min before the test but the dogs’ reactions to their owners’ return could not be controlled. Excited dogs might have had higher baseline levels of the measured markers, which would hypothetically result in lower differences between post-test and pre-test marker concentrations. However, no such patterns were identified in the dataset, and the findings described above appear to reflect plausible differences in reactivity. This suggests that the influence of different baseline levels was minor.
The findings of this exploratory study suggest that individual canine differences in physiological stress response relate to behavioral traits, mainly describing fearfulness, reactivity or responsiveness to training. Different patterns were found for physiological stress reactivity and recovery, suggesting that both are important for determining dog stress response and coping ability. High reactivity in response to a stressor is ambiguous to interpret, as it could indicate an adaptive or a maladaptive response. Slow recovery, on the other hand, reflects a maladaptive response. Differences in markers were also found, confirming that a combined analysis of multiple endocrine and immune biomarkers would give a more accurate indication of stress levels [17,26]. Finally, predictive capacity of puppies’ physiological stress response was inconclusive, and more research is needed to assess the influence of developmental factors. Though these tenuous conclusions are based on a relatively small dataset, they shed new light on the interpretation of physiological markers and their application in future studies. Hopefully, our findings will inspire other researchers to measure both reactivity and recovery and to use multiple markers for a more accurate indication of the physiological stress response. Further research is needed to confirm these patterns in a larger population.
S2 Table [pdf]
Tested associations (Mann-Whitney U test) between dogs’ physiological responses and C-BARQ scores.
1. Carere C, Caramaschi D, Fawcett TW. Covariation between personalities and individual differences in coping with stress: Converging evidence and hypotheses. Curr Zool. 2010;56: 728–740.
2. Rayment DJ, De Groef B, Peters RA, Marston LC. Applied personality assessment in domestic dogs: Limitations and caveats. Appl Anim Behav Sci. 2015;163: 1–18. doi: 10.1016/j.applanim.2014.11.020
3. Sforzini E, Michelazzi M, Spada E, Ricci C, Carenzi C, Milani S, et al. Evaluation of young and adult dogs’ reactivity. J Vet Behav. 2009;4: 3–10. doi: 10.1016/j.jveb.2008.09.035
4. Jones AC, Gosling SD. Temperament and personality in dogs (Canis familiaris): A review and evaluation of past research. Appl Anim Behav Sci. 2005;95: 1–53.
5. Wolf M, van Doorn GS, Weissing FJ. Evolutionary emergence of responsive and unresponsive personalities. Proc Natl Acad Sci U S A. 2008;105: 15825–15830. doi: 10.1073/pnas.0805473105 18838685
6. Koolhaas JM, de Boer SF, Coppens CM, Buwalda B. Neuroendocrinology of coping styles: Towards understanding the biology of individual variation. Front Neuroendocrinol. 2010;31: 307–321. doi: 10.1016/j.yfrne.2010.04.001 20382177
7. Overall KL. Manual of Clinical Behavioral Medicine for Dogs and Cats. St. Louis, Missouri: Elsevier; 2013.
8. Wormald D, Lawrence AJ, Carter G, Fisher AD. Physiological stress coping and anxiety in greyhounds displaying inter-dog aggression. Appl Anim Behav Sci. 2016;180: 93–99. doi: 10.1016/j.applanim.2016.04.007
9. Barnard S, Marshall-Pescini S, Passalacqua C, Beghelli V, Capra A, Normando S, et al. Does Subjective Rating Reflect Behavioural Coding? Personality in 2 Month-Old Dog Puppies: An Open-Field Test and Adjective-Based Questionnaire. PLoS One. 2016;11: e0149831. doi: 10.1371/journal.pone.0149831 26977588
10. Riemer S, Müller C, Virányi Z, Huber L, Range F. The predictive value of early behavioural assessments in pet dogs—A longitudinal study from neonates to adults. PLoS One. 2014;9: e101237. doi: 10.1371/journal.pone.0101237 25003341
11. Ley J, Coleman GJ, Holmes R, Hemsworth PH. Assessing fear of novel and startling stimuli in domestic dogs. Appl Anim Behav Sci. 2007;104: 71–84. doi: 10.1016/j.applanim.2006.03.021
12. Ng ZY, Pierce BJ, Otto CM, Buechner-Maxwell VA, Siracusa C, Werre SR. The effect of dog–human interaction on cortisol and behavior in registered animal-assisted activity dogs. Appl Anim Behav Sci. 2014;159: 69–81. doi: 10.1016/j.applanim.2014.07.009
13. Stamps J, Groothuis TGG. The development of animal personality: Relevance, concepts and perspectives. Biol Rev. 2010;85: 301–325. doi: 10.1111/j.1469-185X.2009.00103.x 19961473
14. Horváth Z, Igyártó B-Z, Magyar A, Miklósi Á. Three different coping styles in police dogs exposed to a short-term challenge. Horm Behav. 2007;52: 621–630. doi: 10.1016/j.yhbeh.2007.08.001 17900575
15. Sherman BL, Gruen ME, Case BC, Foster ML, Fish RE, Lazarowski L, et al. A test for the evaluation of emotional reactivity in Labrador retrievers used for explosives detection. J Vet Behav. 2015;10: 94–102. doi: 10.1016/j.jveb.2014.12.007
16. Lensen CMM, Moons CPH, Diederich C. Saliva sampling in dogs: How to select the most appropriate procedure for your study. J Vet Behav. 2015;10: 504–512. doi: 10.1016/j.jveb.2015.08.006
17. Hekman JP, Karas AZ, Sharp CR. Psychogenic stress in hospitalized dogs: Cross species comparisons, implications for health care, and the challenges of evaluation. Animals. 2014;4: 331–347. doi: 10.3390/ani4020331 25126422
18. Beerda B, Schilder MBH, van Hooff JARAM, de Vries HW, Mol JA. Behavioural, saliva cortisol and heart rate responses to different types of stimuli in dogs. Appl Anim Behav Sci. 1998;58: 365–381. doi: 10.1016/S0168-1591(97)00145-7
19. Dreschel NA, Granger DA. Physiological and behavioral reactivity to stress in thunderstorm-phobic dogs and their caregivers. Appl Anim Behav Sci. 2005;95: 153–168. doi: 10.1016/j.applanim.2005.04.009
20. Mormède P, Andanson S, Aupérin B, Beerda B, Guémené D, Malmkvist J, et al. Exploration of the hypothalamic-pituitary-adrenal function as a tool to evaluate animal welfare. Physiol Behav. 2007;92: 317–339. doi: 10.1016/j.physbeh.2006.12.003 17234221
21. Protopopova A. Effects of sheltering on physiology, immune function, behavior, and the welfare of dogs. Physiol Behav. 2016;159: 95–103. doi: 10.1016/j.physbeh.2016.03.020 26996275
22. Kikkawa A, Uchida Y, Nakade T, Taguchi K. Salivary secretory IgA concentrations in beagle dogs. J Vet Med Sci. 2003;65: 689–693. doi: 10.1292/jvms.65.689 12867728
23. Lucassen PJ, Pruessner J, Sousa N, Almeida OFX, Van Dam AM, Rajkowska G, et al. Neuropathology of stress. Acta Neuropathol. 2014;127: 109–135. doi: 10.1007/s00401-013-1223-5 24318124
24. Hiby EF, Rooney NJ, Bradshaw JWS. Behavioural and physiological responses of dogs entering re-homing kennels. Physiol Behav. United States; 2006;89: 385–391. doi: 10.1016/j.physbeh.2006.07.012 16905163
25. Haverbeke A, Diederich C, Depiereux E, Giffroy JM. Cortisol and behavioral responses of working dogs to environmental challenges. Physiol Behav. 2008;93: 59–67. doi: 10.1016/j.physbeh.2007.07.014 17868751
26. Ott S, Soler L, Moons CPH, Kashiha MA, Bahr C, Vandermeulen J, et al. Different stressors elicit different responses in the salivary biomarkers cortisol, haptoglobin, and chromogranin A in pigs. Res Vet Sci. 2014;97: 124–128. doi: 10.1016/j.rvsc.2014.06.002 24957408
27. Ivković N, Božović Đ, Račić M, Popović-Grubač D, Davidović B. Biomarkers of Stress in Saliva / Biomarkeri stresa u pljuvački. Acta Fac Medicae Naissensis. 2015;32: 91–99. doi: 10.1515/afmnai-2015-0010
28. Kennedy B, Dillon E, Mills PJ, Ziegler MG. Catecholamines in human saliva. Life Sci. 2001;69: 87–99. doi: 10.1016/s0024-3205(01)01111-0 11411808
29. Akiyoshi H, Aoki M, Shimada T, Noda K, Kumagai D, Saleh N, et al. Measurement of plasma chromogranin A concentrations for assessment of stress responses in dogs with insulin-induced hypoglycemia. Am J Vet Res. 2005;66: 1830–1835. 16273918
30. Toda M, Morimoto K. Effect of lavender aroma on salivary endocrinological stress markers. Arch Oral Biol. 2008;53: 964–968. doi: 10.1016/j.archoralbio.2008.04.002 18635155
31. Escribano D, Soler L, Gutiérrez AM, Martínez-Subiela S, Cerón JJ. Measurement of chromogranin A in porcine saliva: validation of a time-resolved immunofluorometric assay and evaluation of its application as a marker of acute stress. Animal. 2013;7: 640–647. doi: 10.1017/S1751731112002005 23137814
32. Ninomiya S, Sato S. The assessment of the effect of presenting a companion’s face picture on social isolation stress using saliva sampling in cows. Anim Sci J. 2011;82: 787–791. doi: 10.1111/j.1740-0929.2011.00896.x 22111636
33. Kanai K, Hino M, Hori Y, Nakao R, Hoshi F, Itoh N, et al. Circadian variations in salivary chromogranin a concentrations during a 24-hour period in dogs. J Vet Sci. 2008;9: 421–423. doi: 10.4142/jvs.2008.9.4.421 19052500
34. Skandakumar S, Stodulski G, Hau J. Salivary IgA: a Possible Stress Marker In Dogs. Anim Welf. 1995;4: 339–350.
35. Hsu Y, Serpell JA. Development and validation of a questionnaire for measuring behavior and temperament traits in pet dogs. J Am Vet Med Assoc. 2003;223: 1293–1300. doi: 10.2460/javma.2003.223.1293 14621216
36. Mills DS. Medical paradigms for the study of problem behaviour: a critical review. Appl Anim Behav Sci. 2003;81: 265–277. doi: 10.1016/S0168-1591(02)00286-1
37. Duffy DL, Serpell JA. Predictive validity of a method for evaluating temperament in young guide and service dogs. Appl Anim Behav Sci. 2012;138: 99–109. doi: 10.1016/j.applanim.2012.02.011
38. Foyer P, Bjällerhag N, Wilsson E, Jensen P. Behaviour and experiences of dogs during the first year of life predict the outcome in a later temperament test. Appl Anim Behav Sci. 2014;155: 93–100. doi: 10.1016/j.applanim.2014.03.006
39. van den Berg L, Schilder MB, de Vries H, Leegwater PA, van Oost BA. Phenotyping of aggressive behavior in golden retriever dogs with a questionnaire. Behav Genet. 2006;36: 882–902. doi: 10.1007/s10519-006-9089-0 16799833
40. Sandri M, Colussi A, Perrotta MG, Stefanon B. Salivary cortisol concentration in healthy dogs is affected by size, sex, and housing context. J Vet Behav. 2015;10: 302–306.
41. Glenk LM, Kothgassner OD, Stetina BU, Palme R, Kepplinger B, Baran H. Salivary cortisol and behavior in therapy dogs during animal-assisted interventions: A pilot study. J Vet Behav. 2014;9: 98–106. doi: 10.1016/j.jveb.2014.02.005
42. Salimetrics LLC, SalivaBio. Saliva collection and handling advice, 3rd edition [Internet]. 2013 [cited 10 Feb 2014]. Available: https://www.salimetrics.com/assets/documents/Saliva_Collection_Handbook.
43. Hennessy MB T. Williams M, Miller DD, Douglas CW, Voith VL. Influence of male and female petters on plasma cortisol and behaviour: can human interaction reduce the stress of dogs in a public animal shelter? Appl Anim Behav Sci. 1998;61: 63–77. doi: 10.1016/s0168-1591(98)00179-8
44. Kobelt AJ, Hemsworth PH, Barnett JL, Butler KL. Sources of sampling variation in saliva cortisol in dogs. Res Vet Sci. 2003;75: 157–161. doi: 10.1016/s0034-5288(03)00080-8 12893165
45. Ricos C. Biologic variation and desirable specifications for QC [Internet]. [cited 9 Feb 2015]. Available: https://www.westgard.com/guest17.htm
46. Schalke E, Stichnoth J, Ott S, Jones-Baade R. Clinical signs caused by the use of electric training collars on dogs in everyday life situations. Appl Anim Behav Sci. 2007;105: 369–380. doi: 10.1016/j.applanim.2006.11.002
47. King T, Hemsworth PH, Coleman GJ. Fear of novel and startling stimuli in domestic dogs. Appl Anim Behav Sci. 2003;82: 45–64. doi: 10.1016/s0168-1591(03)00040-6
49. Hennessy MB. Using hypothalamic-pituitary-adrenal measures for assessing and reducing the stress of dogs in shelters: A review. Appl Anim Behav Sci. 2013;149: 1–12. doi: 10.1016/j.applanim.2013.09.004
50. Beerda B, Schilder MBH, Van Hooff JARAM, De Vries HW, Mol JA. Behavioural and hormonal indicators of enduring environmental stress in dogs. Anim Welf. 2000;9: 49–62.
51. Beerda B, Schilder MBH, Bernadina W, Van Hooff JARAM, De Vries HW, Mol JA. Chronic Stress in Dogs Subjected to Social and Spatial Restriction. II. Hormonal and Immunological Responses. Physiol Behav. 1999;66: 243–254. doi: 10.1016/s0031-9384(98)00290-x 10336150
52. Rosado B, García-Belenguer S, León M, Chacón G, Villegas A, Palacio J. Blood concentrations of serotonin, cortisol and dehydroepiandrosterone in aggressive dogs. Appl Anim Behav Sci. 2010;123: 124–130. doi: 10.1016/j.applanim.2010.01.009
53. Kikkawa A, Uchida Y, Suwa Y, Taguchi K. A novel method for estimating the adaptive ability of guide dogs using salivary sIgA. J Vet Med Sci. 2005;67: 707–712. doi: 10.1292/jvms.67.707 16082119
54. Bosch JA, Ring C, de Geus EJC, Veerman ECI, Nieuw Amerongen AV. Stress and secretory immunity. Int Rev Neurobiol. 2002; 52: 213–253. doi: 10.1016/S0074-7742(02)52011-0 12498106
55. Nakane H, Asami O, Yamada Y, Harada T, Matsui N, Kanno T, et al. Salivary Chromogranin A as an index of psychosomatic stress response. Biomed Res. 1998;19: 401–406. doi: 10.2220/biomedres.19.401
56. Svobodová I, Chaloupková H, Končel R, Bartoš L, Hradecká L, Jebavý L. Cortisol and Secretory Immunoglobulin A Response to Stress in German Shepherd Dogs. PLoS One. 2014;9: e90820. doi: 10.1371/journal.pone.0090820 24637917
57. Uher J. Individual behavioral phenotypes: An integrative meta-theoretical framework. Why “behavioral syndromes” are not analogs of “personality.” Dev Psychobiol. 2011;53: 521–548. doi: 10.1002/dev.20544 21432848
58. Day MJ. Immune System Development in the Dog and Cat. J Comp Pathol. 2007;137: S10–S15. doi: 10.1016/j.jcpa.2007.04.005 17560591
59. Cobb ML, Iskandarani K, Chinchilli VM, Dreschel NA. A systematic review and meta-analysis of salivary cortisol measurement in domestic canines. Domest Anim Endocrinol. 2016;57: 31–42. doi: 10.1016/j.domaniend.2016.04.003 27315597
60. Appleby D, Pluijmakers J. Separation anxiety in dogs. The function of homeostasis in its development and treatment. Vet Clin North Am Small Anim Pract. 2003;33: 321–344. doi: 10.1016/S0195-5616(02)00101-8 12701515
61. Bennett A, Hayssen V. Measuring cortisol in hair and saliva from dogs: coat color and pigment differences. Domest Anim Endocrinol. 2010/08/14. 2010;39: 171–180. doi: 10.1016/j.domaniend.2010.04.003 20705413
62. Siniscalchi M, McFarlane JR, Kauter KG, Quaranta A, Rogers LJ. Cortisol levels in hair reflect behavioural reactivity of dogs to acoustic stimuli. Res Vet Sci. 2013;94: 49–54. doi: 10.1016/j.rvsc.2012.02.017 22449333