The deer ked (Lipoptena cervi), a hematophagous ectoparasite of cervids, is currently spreading in Scandinavia, and the moose (Alces alces) is its main host. However, little is known about the impact of deer keds on moose. We analyzed the hair cortisol concentration (HCC) from 262 moose harvested in the fall in relation to age class, sex, body mass (BM), and deer ked infestation intensity, and BM in relation to age class, sex, and infestation intensity. We found that HCC decreased with increasing deer ked intensity at low ked intensities, but for the higher levels of ked intensities, there was a positive relationship between HCC and ked intensity. The HCC was higher in males than in females and lower in yearlings than in calves and adults. Our failure to find any association between BM and deer ked intensity suggested a negligible impact of deer ked infestation on moose foraging and metabolism at the level of infestation observed early in the infestation, but did not exclude an effect later in winter. Our findings suggested that moose generally tolerated moderate parasitism by keds. However, the increase in HCC at higher ked intensities suggested that the tolerance strategy could be disrupted with further increases in intensities and consequently may negatively affect animal health and welfare.

The deer ked (Lipoptena cervi) is a hematophagous louse fly (Hippoboscidae) found on moose (Alces alces), red deer (Cervus elaphus), and other species of deer in Europe, Asia, and North America (Bequaert 1942). In many areas, the swarming keds constitute a major obstacle for human outdoor activities, but the parasite has been regarded as harmless for cervids (Allan 2001; Paakkonen et al. 2012). However, Madslien et al. (2011) reported an outbreak of severe alopecia in moose associated with massive infestation of deer keds in Norway, and Kynkäänniemi et al. (2014) stated that reindeer (Rangifer tarandus tarandus) showed signs of distress and pruritus after infestation with keds.

Hair analysis is becoming increasingly popular as a method to examine glucocorticoid (GC) and sex steroid hormone concentrations in wild mammals (Koren et al. 2019). Glucocorticoid concentrations are regarded as physiological indices of stress and the relative condition or health of individuals and populations (Bonier et al. 2009).

We aimed to 1) examine the relationship between infestation intensity of deer keds and hair cortisol concentration (HCC) in moose and 2) test whether there was a negative association between infestation intensity of deer keds and moose body mass (BM). Our hypothesis was that increasing intensities of deer keds results in increased HCC and decreased BM.

Skin samples were collected from 262 moose: 36 calves (16 males/20 females), 92 yearlings (56 males/36 females), and 134 adults (67 males/67 females) in Hedmark and Akershus counties, in southeastern Norway, during the first week of the hunting season (5–12 October 2010). Hunters recorded sex and age class (calf, yearling, adult) of the moose, measured the carcass weight (about 50% of full BM; Wallin et al. 1996) and collected a standardized 20×20-cm skin sample from the neck area (Madslien et al. 2012).

In the laboratory, guard hairs for HCC analysis were collected by shaving a skin area of 2×2 cm. The hair samples were stored in paper envelopes at room temperature until analysis. Finally, deer ked infestation intensities were calculated from all skin samples (Madslien et al. 2012). The cortisol analysis was performed at the Toxicology Centre, University of Saskatchewan (Saskatoon, Saskatchewan, Canada). Only guard hairs with the follicles removed were used to determine HCC. Surface contamination was removed by washing hairs with methanol (five 3-min washes). Between 100 and 200 mg of washed and dried hair was ground to a fine powder, using a ball mill, and weighed. Extraction of cortisol from hair and cortisol concentration (pg/mg powdered hair) were determined as described by Macbeth et al. (2010).

Generalized additive modeling with gamma distribution and inverse link was run to analyze the effect of ked intensity on HCC. Ked intensity was transformed (natural log) to normalize and stabilize the variance. Age class, sex, and BM were then added to the baseline model (i.e., forward selection), using the small-sample–corrected Akaike information criterion (Burnham and Anderson 2003). We used generalized linear regression to analyze the effect of age class, sex, and ked intensity on moose BM. Statistical analyses were performed using the R statistical software (R version 2.14.1; R Development Core Team 2019) and the generalized additive modeling was run in the R package mgcv.

Mean ked intensity (mean number of keds/cm2) in skin samples and HCC varied among sex and age classes (Table 1). The HCC decreased with increasing infestation intensity at low ked intensities, but for the higher levels of ked intensity, there was a positive relationship between HCC and infestation intensity (Fig. 1A). Yearlings had lesser HCC than adults (P=0.003) and males had greater HCC than females (P<0.001; Table 2). When analyzing moose with deer ked intensities greater or equal to the mean, we found a slight increase in HCC with increasing ked intensities (Table 3 and Fig. 2). Variation in BM was best explained by age class and sex (Fig. 1B and Table 4).

Table 1

Summary statistics for deer ked (Lipoptena cervi) intensity (keds/cm2), hair cortisol concentration (HCC; pg/mg), and body mass (kg) of moose (Alces alces; n=262) infested with deer keds in southeastern Norway.

Summary statistics for deer ked (Lipoptena cervi) intensity (keds/cm2), hair cortisol concentration (HCC; pg/mg), and body mass (kg) of moose (Alces alces; n=262) infested with deer keds in southeastern Norway.
Summary statistics for deer ked (Lipoptena cervi) intensity (keds/cm2), hair cortisol concentration (HCC; pg/mg), and body mass (kg) of moose (Alces alces; n=262) infested with deer keds in southeastern Norway.
Figure 1

Prediction of hair cortisol concentration (HCC; pg/mg) shown as a spline function of deer ked (Lipoptena cervi) intensity (keds/cm2) (A) and body mass (kg) (B) in 262 Norwegian moose (Alces alces). Prediction lines together with 95% confidence envelopes are shown for males (blue) and females (green) at the level of (A) overall mean of deer ked intensity and (B) mean body mass of adult females and adult males, respectively. Raw data of deer ked intensity and corresponding HCC are imposed as points of various symbols for the three age classes.

Figure 1

Prediction of hair cortisol concentration (HCC; pg/mg) shown as a spline function of deer ked (Lipoptena cervi) intensity (keds/cm2) (A) and body mass (kg) (B) in 262 Norwegian moose (Alces alces). Prediction lines together with 95% confidence envelopes are shown for males (blue) and females (green) at the level of (A) overall mean of deer ked intensity and (B) mean body mass of adult females and adult males, respectively. Raw data of deer ked intensity and corresponding HCC are imposed as points of various symbols for the three age classes.

Close modal
Table 2

Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) in Norwegian moose (Alces alces; n=262), using gamma distribution with inverse link. Because of the inverse link of the gamma distribution, a significant negative estimate refers to a positive association with HCC. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is adult females at mean ln deer ked intensity (Intensity). The nonlinear relationship between deer ked intensity and HCC is modeled as a smooth spline (s), as shown in Figure 1A, and reported by the estimated degrees of freedom and the significance of the smooth term.

Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) in Norwegian moose (Alces alces; n=262), using gamma distribution with inverse link. Because of the inverse link of the gamma distribution, a significant negative estimate refers to a positive association with HCC. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is adult females at mean ln deer ked intensity (Intensity). The nonlinear relationship between deer ked intensity and HCC is modeled as a smooth spline (s), as shown in Figure 1A, and reported by the estimated degrees of freedom and the significance of the smooth term.
Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) in Norwegian moose (Alces alces; n=262), using gamma distribution with inverse link. Because of the inverse link of the gamma distribution, a significant negative estimate refers to a positive association with HCC. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is adult females at mean ln deer ked intensity (Intensity). The nonlinear relationship between deer ked intensity and HCC is modeled as a smooth spline (s), as shown in Figure 1A, and reported by the estimated degrees of freedom and the significance of the smooth term.
Table 3

Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) for the subset of moose (Alces alces; n=93) with deer ked intensity higher or equal to mean intensity. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1.

Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) for the subset of moose (Alces alces; n=93) with deer ked intensity higher or equal to mean intensity. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1.
Hair cortisol concentration (HCC; pg/mg) as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) for the subset of moose (Alces alces; n=93) with deer ked intensity higher or equal to mean intensity. The column ΔAICc refers to the change in the Akaike information criterion corrected for small sample sizes when the corresponding variable is excluded from the model. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1.
Figure 2

Prediction of hair cortisol concentration (HCC; pg/mg) shown as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) for Norwegian moose (Alces alces) with above mean infestation intensity (prediction lines together with 95% confidence envelopes). Raw data of deer ked intensity and corresponding HCC are imposed as points of various symbols for the three age classes.

Figure 2

Prediction of hair cortisol concentration (HCC; pg/mg) shown as a function of deer ked (Lipoptena cervi) intensity (keds/cm2) for Norwegian moose (Alces alces) with above mean infestation intensity (prediction lines together with 95% confidence envelopes). Raw data of deer ked intensity and corresponding HCC are imposed as points of various symbols for the three age classes.

Close modal
Table 4

Body mass (kg) in Norwegian moose (Alces alces; n=262), modelled as a spline of ln deer ked (Lipoptena cervi) intensity (keds/cm2), using gamma distribution with a log link. When accounting for sex and age class, there was no significant association with body mass and deer ked intensity. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is female adults at population average of ln deer ked intensity. The smooth term (spline function [s]) of deer ked intensity (see Fig. 1B) is reported by the estimated degrees of freedom and the significance.

Body mass (kg) in Norwegian moose (Alces alces; n=262), modelled as a spline of ln deer ked (Lipoptena cervi) intensity (keds/cm2), using gamma distribution with a log link. When accounting for sex and age class, there was no significant association with body mass and deer ked intensity. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is female adults at population average of ln deer ked intensity. The smooth term (spline function [s]) of deer ked intensity (see Fig. 1B) is reported by the estimated degrees of freedom and the significance.
Body mass (kg) in Norwegian moose (Alces alces; n=262), modelled as a spline of ln deer ked (Lipoptena cervi) intensity (keds/cm2), using gamma distribution with a log link. When accounting for sex and age class, there was no significant association with body mass and deer ked intensity. Deer ked intensity is ln transformed, centered around mean, and scaled to variance 1. The reference level (Intercept) is female adults at population average of ln deer ked intensity. The smooth term (spline function [s]) of deer ked intensity (see Fig. 1B) is reported by the estimated degrees of freedom and the significance.

Varying results have been found in studies evaluating the relationship between parasite intensity and GCs (Goldstein et al. 2005; Carlsson et al. 2016). We found that males have greater HCC than females and that HCC was positively associated with BM larger than 150 kg (i.e., adults; Fig. 1B). Similarly, Di Francesco et al. (2017) found significantly greater HCCs in muskoxen males than females. Body mass was found to be negatively associated with HCC for BM less than 125 kg (i.e., yearlings; Fig. 1B). This association agrees with Mislan et al. (2016), who reported polar bears in poorer body condition had greater HCC.

According to the cort-fitness hypothesis, high baseline GC concentrations indicate an individual or population in poorer condition and with lower relative fitness (Bonier et al. 2009). However, the baseline or stress-induced GC concentrations in plasma from laboratory and wild animals did not change in a consistent manner in response to chronic stress (Bonier et al. 2009; Dickens and Romero 2013) also reported inconsistent relationships between plasma GCs and fitness. Pawluski et al. (2017) suggested that the reported variability of GC concentrations under chronic stress may be due to the species studied, type of stressor or welfare measured, duration of stress, and techniques used to analyze GC concentrations.

Most moose in the present study had low to moderate deer ked intensities; only a few moose had high ked intensities. Paakkonen et al. (2012) found only minor effects of high infestation levels of deer ked on a range of systemic physiological measures, including plasma cortisol concentration (PCC). This contrasts with our findings. An explanation might be that the HCC measured in the current study is a more sensitive measurement of chronic stress than is PCC. In addition, HCC may remain less affected by variations in GC secretion in response to acute stress caused by hunting, which has the potential to mask small changes in GC caused by chronic stress compared to PCC.

Hair may accumulate GC hormones over weeks to months, and although HCC is thought to be insensitive to the impact of acute stress, some evidence suggests that capture method may influence HCC (Cattet et al. 2014). Koren et al. (2019) found correlation between hair and serum cortisol in samples from moose. However, several studies have suggested that local cortisol production may contribute to GC integration in hair, because of a parallel, but peripheral, stress axis within the hair follicles (Ito et al. 2005; Keckeis et al. 2012). This does not imply that steroids measured in hair reflect local production only, but that this process may be independently influenced by both the central and peripheral stress axes (Cattet et al. 2014). Salaberger et al. (2016) found that mechanical irritation significantly increased HCC in sheep. Skin irritation caused by deer keds may have induced local GC production and metabolism in hair follicle cells regulated by locally expressed HPA mediators (Stubsjøen et al. 2015).

The lack of association between deer ked intensity and BM suggests that keds have minimal impact on moose foraging and metabolism early in the infestation period. However, we cannot rule out the possibility that deer keds may significantly impact moose during late winter when food is limited, the weather is harsh, and the parasites have exploited the host for a longer period of time.

Our results suggest that moderately intense deer ked infestations do not cause major chronic stress or loss of BM during the early infestation period of the deer ked. However, the increase in HCC at higher ked intensities even during this early parasitic phase of the infestation may be indicative of an increasingly negative response later in the infestation, affecting animal welfare negatively.

This work was supported by the National Health Surveillance Program for Cervids and Muskox (HOP) and the Norwegian Environment Agency. We appreciate the superb sampling performed by voluntary moose hunters. The authors thank Marthe Opland, Norwegian Veterinary Institute, and Bryan Sarauer, University of Saskatchewan, for the excellent technical assistance in the laboratory.

Allan
SA.
2001
.
Biting flies (class Insecta: order Diptera).
In
:
Parasitic disease of wild mammals
,
Samuel
WM
,
Pybus
MJ
,
Kocan
AA
, editors.
Iowa State University Press
,
Ames, Iowa
, pp.
18
45
.
Bequaert
J.
1942
.
A monograph of the Melophaginae or ked flies of sheep, goats, deer and antelopes (Diptera, Hippoboscidae).
Entomol Am
22
:
1
220
.
Bonier
F
,
Martin
PR
,
Moore
IT
,
Wingfield
JC
.
2009
.
Do baseline glucocorticoids predict fitness?
Trends Ecol Evol
24
:
634
642
.
Burnham
KP
,
Anderson
DR
.
2003
.
Model selection and multimodel inference: A practical information-theoretic approach.
2nd Ed.
Springer
,
New York, New York
,
488
pp.
Carlsson
AM
,
Mastromonaco
G
,
Vandervalk
E
,
Kutz
S
.
2016
.
Parasites, stress and reindeer: Infection with abomasal nematodes is not associated with elevated glucocorticoid levels in hair or faeces.
Conserv Physiol
4
:
cow058
.
Cattet
M
,
Macbeth
BJ
,
Janz
DM
,
Zedrosser
A
,
Swenson
JE
,
Dumond
M
,
Stenhouse
GB
.
2014
.
Quantifying long-term stress in brown bears with the hair cortisol concentration: A biomarker that may be confounded by rapid changes in response to capture and handling.
Conserv Physiol
2
:
cou026
.
Di Francesco
J
,
Navarro-Gonzalez
N
,
Wynne-Edwards
K
,
Peacock
S
,
Leclerc
LM
,
Tomaselli
M
,
Davison
T
,
Carlsson
A
,
Kutz
S
.
2017
.
Qiviut cortisol in muskoxen as a potential tool for informing conservation strategies.
Conserv Physiol
5
:
cox052
.
Dickens
MJ
,
Romero
LM
.
2013
.
A consensus endocrine profile for chronically stressed wild animals does not exist.
Gen Comp Endocrinol
191
:
177
189
.
Goldstein
EJ
,
Millspaugh
JJ
,
Washburn
BE
,
Brundige
GC
,
Raedeke
KJ
.
2005
.
Relationships among fecal lungworm loads, fecal glucocorticoid metabolites, and lamb recruitment in free-ranging Rocky Mountain bighorn sheep.
J Wildl Dis
41
:
416
425
.
Ito
N
,
Ito
T
,
Kromminga
A
,
Bettermann
A
,
Takigawa
M
,
Kees
F
,
Straub
RH
,
Paus
R
.
2005
.
Human hair follicles display a functional equivalent of the hypothalamic-pituitary-adrenal axis and synthesize cortisol.
FASEB J
19
:
1332
1334
.
Keckeis
K
,
Lepschy
M
,
Schöpper
H
,
Moser
L
,
Troxler
J
,
Palme
R
.
2012
.
Hair cortisol: A parameter of chronic stress? Insights from a radiometabolism study in guinea pigs.
J Comp Physiol B
182
:
985
996
.
Koren
L
,
Bryan
H
,
Matas
D
,
Tinman
S
,
Fahlman
Å
,
Whiteside
D
,
Smits
J
,
Wynne-Edwards
K
.
2019
.
Towards the validation of endogenous steroid testing in wildlife hair.
J Appl Ecol
56
:
547
561
.
Kynkäänniemi
SM
,
Kettu
M
,
Kortet
R
,
Häarköonen
L
,
Kaitala
A
,
Paakkonen
T
,
Mustonen
AM
,
Nieminen
P
,
Härkönen
S
,
Ylönen
H
, et al.
2014
.
Acute impacts of the deer ked (Lipoptena cervi) infestation on reindeer (Rangifer tarandus tarandus) behaviour.
Parasitol Res
113
:
1489
1497
.
Macbeth
BJ
,
Cattet
MRL
,
Stenhouse
GB
,
Gibeau
ML
,
Janz
DM
.
2010
.
Hair cortisol concentration as a noninvasive measure of long-term stress in free-ranging grizzly bears (Ursus arctos): Considerations with implications for other wildlife.
Can J Zool
88
:
935
949
.
Madslien
K
,
Ytrehus
B
,
Vikøren
T
,
Malmsten
J
,
Isaksen
K
,
Hygen
HO
,
Solberg
EJ
.
2011
.
Hair-loss epizootic in moose (Alces alces) associated with massive deer ked (Lipoptena cervi) infestation.
J Wildl Dis
47
:
893
906
.
Madslien
K
,
Ytrehus
B
,
Viljugrein
H
,
Solberg
EJ
,
Bråten
KR
,
Mysterud
A
.
2012
.
Factors affecting deer ked (Lipoptena cervi) prevalence and infestation intensity in moose (Alces alces) in Norway.
Parasit Vectors
5
:
251
.
Mislan
P
,
Derocher
AE
, St
Louis
VL
,
Richardson
E
,
Lunn
NJ
,
Janz
DM
.
2016
.
Assessing stress in Western Hudson Bay polar bears using hair cortisol concentration as a biomarker.
Ecol Indic
71
:
47
54
.
Paakkonen
T
,
Mustonen
AM
,
Käkelä
R
,
Laaksonen
S
,
Solismaa
M
,
Aho
J
,
Puukka
K
,
Nieminen
P
.
2012
.
The effects of an abundant ectoparasite, the deer ked (Lipoptena cervi), on the health of moose (Alces alces) in Finland.
Parasitol Res
111
:
1223
1232
.
Pawluski
J
,
Jego
P
,
Henry
S
,
Bruchet
A
,
Palme
R
,
Coste
C
,
Hausberger
M
.
2017
.
Low plasma cortisol and fecal cortisol metabolite measures as indicators of compromised welfare in domestic horses (Equus caballus).
PLoS One
12
:
e0182257
.
R Development Core Team
.
2019
.
R: A language and environment for statistical computing.
R Foundation for Statistical Computing
,
Vienna, Austria
.
http://www.R-project.org. Accessed September 2019
.
Salaberger
T
,
Millard
M
,
Makarem
SE
,
Möstl
E
,
Grünberger
V
,
Krametter-Frötscher
R
,
Wittek
T
,
Palme
R
.
2016
.
Influence of external factors on hair cortisol concentrations.
Gen Comp Endocrinol
233
:
73
78
.
Stubsjøen
SM
,
Bohlin
J
,
Dahl
E
,
Knappe-Poindecker
M
,
Fjeldaas
T
,
Lepschy
M
,
Palme
R
,
Langbein
J
,
Ropstad
E
.
2015
.
Assessment of chronic stress in sheep (part I): The use of cortisol and cortisone in hair as non-invasive biological markers.
Small Rumin Res
132
:
25
31
.
Wallin
K
,
Cederlund
G
,
Pehrson
A
.
1996
.
Predicting body mass from chest circumference in moose Alces alces.
Wildl Biol
2
:
53
58
.