Inbreeding effects on immune response in free-living song sparrows (Melospiza melodia)

The consequences of inbreeding for host immunity to parasitic infection have broad implications for the evolutionary and dynamical impacts of parasites on populations where inbreeding occurs. To rigorously assess the magnitude and the prevalence of inbreeding effects on immunity, multiple components of host immune response should be related to inbreeding coefficient (f) in free-living individuals. We used a pedigreed, free-living population of song sparrows (Melospiza melodia) to test whether individual responses to widely used experimental immune challenges varied consistently with f. The patagial swelling response to phytohaemagglutinin declined markedly with f in both females and males in both 2002 and 2003, although overall inbreeding depression was greater in males. The primary antibody response to tetanus toxoid declined with f in females but not in males in both 2004 and 2005. Primary antibody responses to diphtheria toxoid were low but tended to decline with f in 2004. Overall inbreeding depression did not solely reflect particularly strong immune responses in outbred offspring of immigrant–native pairings or weak responses in highly inbred individuals. These data indicate substantial and apparently sex-specific inbreeding effects on immune response, implying that inbred hosts may be relatively susceptible to parasitic infection to differing degrees in males and females.


INTRODUCTION
The relationship between a host's inbreeding level and its immunity to parasitic infection is of broad interest in evolutionary and population ecology. A decline in immunity with inbreeding would indicate a dominance genetic component to parasite resistance, bearing on the evolutionary consequences of host-parasite interactions and the maintenance of genetic variation in natural populations (Potts & Wakeland 1990;Coltman et al. 1999;Westerdahl et al. 2005). Such declines would also suggest that parasites could mediate inbreeding depression in reproduction and survival, thereby threatening the persistence of small, inbred populations (Altizer et al. 2003;Tompkins et al. 2006;Whiteman et al. 2006). Quantifying the magnitude and the prevalence of inbreeding depression in immunity is therefore integral to predicting the genetic, evolutionary and population dynamic consequences of host-parasite interactions (Coltman et al. 1999;Reid et al. 2003;Spielman et al. 2004). Furthermore, in general, the evolutionary and dynamical impacts of parasites on host populations depend on whether the fitness cost of parasite exposure is constant or variable, and on the host sex, age or stage classes that are most affected (Murdoch et al. 1997;Coltman et al. 1999). Comprehensive studies of inbreeding effects on host immunity should therefore quantify whether effects act consistently or are exhibited only occasionally or in specific categories of population members (Giese & Hedrick 2003;Cô té et al. 2005).
One common approach to measuring inbreeding effects on immunity is to relate host inbreeding level to the degree of infection or mortality after natural or experimental exposure to specific natural parasites (e.g. Cassinello et al. 2001;Haag et al. 2003;Spielman et al. 2004). A complementary approach is to relate inbreeding level to the host's response to a novel, experimental immune challenge (Giese & Hedrick 2003;Reid et al. 2003;Hawley et al. 2005). This approach has the advantage of eliminating confounding variation in a host's prior exposure to the focal parasite, thereby allowing inbreeding effects on immunity per se to be clearly distinguished from inbreeding effects on exposure (Norris & Evans 2000;Staszewski & Boulinier 2004). Such novel immune challenges are also ecologically and evolutionarily pertinent, since newly emerging parasites may exert particularly severe selection on naive host populations (Daszak et al. 2000;Altizer et al. 2003). Studies using this approach should apply multiple challenges targeting different components of host immunity, since different components may not show correlated responses (Norris & Evans 2000;Haag et al. 2003;Adamo 2004). Furthermore, since inbreeding effects on phenotypic traits are often less severe in captivity than in the wild (Hedrick & Kalinowski 2000;Joron & Brakefield 2003), studies should ideally be carried out in free-living individuals; the evolutionary and dynamical impacts of parasites on natural populations might otherwise be inaccurately inferred. However, despite their broad relevance, there remain few data describing relationships between host inbreeding per se and multiple components of immune response in freeliving individuals, or the extent to which inbreeding effects vary among seasons or categories of population members.
We used a free-living, pedigreed population of song sparrows, Melospiza melodia, inhabiting Mandarte Island, Canada, to test whether individual responses to two widely used experimental immune challenges varied with an individual's coefficient of inbreeding ( f ). Specifically, we related individual variation in the patagial swelling response to phytohaemagglutinin (PHA), and in the humoral antibody response to diphtheria-tetanus vaccine (DTV), to individual f. We repeated each experimental challenge in two different years on different samples of males and females of varying ages, and thereby investigate whether inbreeding effects on immunity varied between years or with host age or sex.

MATERIAL AND METHODS
(a) Study population Mandarte Island, approximately 6 ha in size, lies 25 km northeast of Victoria, British Columbia, Canada. Its resident song sparrow population, which numbered 16-27 breeding pairs during 2002-2005, has been studied intensively since 1975 (Smith et al. 2006). During this long-term study, all song sparrows fledged on Mandarte have been individually colourringed, almost always before leaving their natal territory. All immigrants to the breeding population (1.1 yr K1 on average) have been colour-ringed soon after settling. All population members are therefore individually identifiable (Smith et al. 2006). Detailed observations of breeding behaviour have allowed a complete social pedigree to be constructed, covering all sparrows hatched since 1981 (Keller 1998). We used standard algorithms to estimate each individual's f directly from the pedigree (Falconer & Mackay 1996;Keller 1998).
The value of f reflects the probability that two homologous alleles will be identical by descent and estimates an individual's genome-wide homozygosity relative to the baseline population (Falconer & Mackay 1996). While immigrants to Mandarte are themselves of unknown f, they are genetically distinguishable from Mandarte's existing population at the time of arrival (assignment test based on nine microsatellite loci, p!0.02 for all immigrants during 1993-1996, see also Keller et al. 2001). Offspring of immigrant-native pairings can therefore be defined as outbred ( fZ0.000, Marr et al. 2002). Extra-pair fertilizations (epfs) occur on Mandarte (approx. 25% of offspring hatched during 1993-1996, O'Connor et al. 2006 and in song sparrows more widely (e.g. 11% of offspring, Major & Barber 2004). These epfs introduce error into the social pedigree and therefore into estimates of f. However, there is no evidence that epfs occur systematically with respect to relatedness in song sparrows on Mandarte. Briefly, across 283 broods hatched during 1993-1996, epfs were not more frequent in females that were more closely related to their social mates ( pZ0.52,O'Connor et al. in preparation). Furthermore, females were no more or less closely related to their extra-pair mate than to their social mate (nZ100 triads, pZ0.65, O'Connor et al. in preparation). Epfs are therefore likely to introduce error but not substantial bias into estimates of f. Such error is expected to cause inbreeding depression in phenotypic traits to be underestimated Kruuk et al. 2002;Marr et al. in press). Since f reflects the relatedness between an individual's parents rather than an individual's relatedness to its own offspring, epfs will not necessarily introduce more error into estimates of f in males than in females. Such sex-biased error in f, and consequently in the estimated magnitude of inbreeding depression in phenotypic traits, could only arise given a biased sex ratio in extrapair young. On Mandarte, the sex ratio of extra-pair young identified during 1993-1996 did not differ from 1 : 1 (65 males, 54 females, c 2 1 Z0.4, pZ0.53). Furthermore, sex ratio did not differ between within-pair and extra-pair young ( pZ0.61, Generalized linear mixed model controlling for year and brood). Therefore, although paternity of the song sparrows exposed to recent experimental immune challenges has not yet been fully verified, data from previous years do not lead us to expect sex-biased error in estimates of f or inbreeding depression.
(b) PHA response The patagial (wing-web) swelling response to subcutaneous injection of PHA is a widely used measure of avian immune responsiveness (Goto et al. 1978;Tella et al. 2002;Martin et al. 2006). PHA response is often interpreted as a measure of 'cell-mediated immunity', but may in fact reflect multiple immunological processes manifested as the infiltration and proliferation of multiple cell types over varying time periods . PHA response is thought to trade-off against other life-history components and to vary with individual condition and 'quality' (Tella et al. 2002;Martin et al. 2006), is often positively correlated with survival (Møller & Saino 2003) and has been related to MHC genotype (Taylor et al. 1987;Bonneaud et al. 2005). Therefore, while the underlying immunological processes may be more complex than often recognized, PHA response remains a useful measure of immunological 'condition' in the context of immunoecology.
During February 2002, September 2002 and September 2003, we measured PHA response in song sparrows on Mandarte. Methodology followed Reid et al. (2003). Briefly, sparrows were mist-netted after 11.00 and wing length, tarsus length and mass were recorded. Patagium thicknesses were measured three times using a modified dial calliper (Mitutoyo, Japan). Sparrows were injected with 30 ml 2 mg ml K1 PHA (L9132, Sigma, St Louis, Missouri) in phosphate buffer solution (PBS, 9.7 gl K1 D5773, Sigma, St Louis, Missouri) in the right patagium and 30 ml PBS in the left patagium and roosted overnight in individual enclosures with ad libitum food and water. Left and right patagial thicknesses were remeasured approximately 18 h after injection and sparrows were released. PHA response was estimated as the difference in increase in thickness between right and left patagia over the experimental period. Patagium thickness measurements were highly repeatable within individuals (rO0.94, p!0.0001, Reid et al. 2003).
(c) Diphtheria-tetanus response The specific humoral antibody response to foreign antigens constitutes one major component of avian acquired immunity (Roitt et al. 1998). During September 2004 and September 2005, we measured humoral immunity as a song sparrow's primary antibody response to tetanus and diphtheria toxoids (Hasselquist et al. 1999(Hasselquist et al. , 2001Owen-Ashley et al. 2004). Sparrows were mist-netted, biometrics were recorded as above, and approximately 100 ml blood was collected by brachial venipuncture. Individuals were vaccinated with 70 ml human DTV(2 Lf diphtheria toxoid, 5 Lf tetanus toxoid adsorbed in aluminium phosphate, Aventis Pasteur Ltd, Toronto, Canada) in the pectoral muscle and released. Primary antibody responses peak 9-15 days after vaccination in song sparrows and other passerines (Hasselquist et al. 1999;Owen-Ashley et al. 2004). Correspondingly, we attempted to recapture sparrows 10-12 days after vaccination. Recaptured sparrows were reweighed, blood sampled as before and released. Since we could not control the exact timing of recapture of free-flying individuals and weather conditions impeded mist-netting on some days, we blood sampled all individuals recaptured 8-14 days after vaccination and subsequently controlled statistically for intersample period (see §2d ). Blood samples were placed immediately on ice and centrifuged for 4 min at 3000rpm. within 5 h. Plasma was separated off, stored buried in ice and frozen at K208C within 72 h. Enzyme-linked immunosorbent assays (ELISAs) were subsequently used to quantify tetanus and diphtheria antibody titres in pre-and postvaccination plasma samples. Protocols followed those previously developed for song sparrows and other passerines (Hasselquist et al. 1999(Hasselquist et al. , 2001Owen-Ashley et al. 2004). ELISA plates held an individual's pre-and post-vaccination plasma samples in duplicate. Titres were standardized for any among-plate variation in conditions or reagents within each year by reference to serially diluted standard samples that were included on each plate. Each individual's primary antibody response to tetanus and diphtheria toxoids was estimated as the difference between post-and pre-vaccination standardized antibody titres (Hasselquist et al. 1999;Owen-Ashley et al. 2004). Owing to within-year standardization of titres and slight between-year differences in ELISA assays, between-year differences in mean antibody response may not solely reflect true biological variation. In naive hosts, pre-vaccination antibody titres are expected to be low but non-zero due to natural antibodies that can bind antigens without previous exposure (Lee et al. 2006).

(d) Analysis
We used general linear models (GLMs) to test whether PHA, tetanus or diphtheria responses varied with individual f within each year in which each immune challenge was applied, or across all data combined. We additionally modelled effects of individual age, body condition, maternal f and paternal f as covariates and sex and year as fixed factors. Sexes were determined by observing adult breeding behaviour or PCR amplification of sex-linked CHD1 genes (Smith et al. 2006;Heinrich et al. in preparation). Body condition was estimated as the residual of mass on the cube of the first principal component of wing and tarsus length. Age was calculated to the nearest month from ringing data. We additionally modelled challenge date (within season) and inter-sample period as linear and quadratic covariates. Challenge date was eliminated from all models (all pO0.35). Some individuals vaccinated with DTV in 2005 were offspring of mothers that had been vaccinated in 2004. Since humoral immunity shows inter-generational effects of maternal vaccination in song sparrows (Reid et al. 2006), we included maternal vaccination history as a binary fixed factor ('vaccinated' or 'unvaccinated') in analyses of antibody responses in 2005. Since samples included sets of siblings, we initially modelled 'family' as a random factor. However, since family effects were not significant, models reduced to final GLMs.
A small number of individuals experienced the same immune challenge twice in consecutive years. Since repeat immune responses may not be independent (Roitt et al. 1998), we restricted analyses to each individual's first exposure to each challenge. We did not apply both PHA and DTV to the same individuals in the same year because mounting multiple simultaneous immune responses can reduce survival (Hansson et al. 2004). Ten individuals that experienced DTV during 2004had experienced PHA during 2002. Diphtheria and tetanus responses did not differ between individuals that had and had not experienced PHA ( pO0.6).
Since the Mandarte song sparrow pedigree is relatively deep, the only individuals that are considered completely outbred ( fZ0.000) are the offspring of immigrant-native pairings. Such F1 offspring of between-population crosses may show enhanced phenotypes due to heterosis (Falconer & Mackay 1996;Marr et al. 2002). Our samples also included some highly inbred offspring of pairings between first-or second-order relatives. Such pairings may involve a phenotypically non-random subset of parents (Reid et al. in preparation) and form relative outliers with respect to f. Therefore, to test whether overall declines in immune response with inbreeding were solely driven by heterosis in outbred offspring or by particularly weak responses in highly inbred individuals, we investigated whether final models or parameter estimates were altered by excluding outbred individuals, or individuals with f R0.125.
Dependent variables were log transformed to reduce deviations from normality. Estimated slopes can consequently be directly interpreted as the inbreeding load or number of lethal equivalents (the slope of a log (trait) on f regression, Falconer & Mackay 1996). B f , B m and B are the estimated inbreeding loads for females, males and all individuals (assuming no sex!f interaction) and are presented with 95% confidence limits. h 2 is the partial effect size for f. Independent variables were not tightly correlated (all r!0.4). Residuals were approximately normally distributed and were not correlated or heteroscedastic with respect to predicted values. Immune response data were collected blind to individual f and other characteristics. Analyses were run in Pedigree Viewer (http://www-personal.une.edu.au/ bkinghor/pedigree.htm), SPSS (v. 14.0) and R (v. 2.2.1). All tests were two tailed. Variables were retained in models if p%0.10. Two-and three-way interactions involving individual f were tested and eliminated except where stated. Sample sizes vary among models because body condition was not calculated for two individuals in wing moult during September 2002 and f was unknown for immigrant parents. Analyses of PHA response in 2002 differ from those reported previously (Reid et al. 2003) in that data from adult and juvenile sparrows are combined. Fieldwork was approved by the University of British Columbia Animal Care Committee.
Inbreeding and immunity J. M. Reid et al. 699 3. RESULTS Table 1 (table 3). However, the main effect of sex and the sex!f interaction were significant in both 2004 and 2005 and across all data; females mounted slightly higher mean tetanus responses than males, and tetanus response declined with f in females but not in males (table 3, figure 2). Tetanus response also varied with inter-sample period and maternal vaccination history and tended to vary with paternal f in 2005. Sparrows whose mothers had been vaccinated in 2004 or whose fathers were relatively inbred showed stronger antibody responses (although the latter effect was weak, table 3, see also Reid et al. 2006). The final model remained quantitatively similar after excluding nine outbred individuals (B f ZK7.4 (95% CI K12.8 to K2.0), B m ZK1.5 (95% CI K5.0 to 2.0)) and qualitatively similar after excluding five highly inbred individuals (although confidence intervals around estimated inbreeding loads increased substantially in this case, B f ZK11.2 (95% CI K20.5 to K1.9), B m Z2.3 (95% CI K5.6 to 10.3)). Diphtheria response tended to decline with increasing f (figure 3) and increased then declined with age (age: FZ18.2, p!0.001; age 2 : FZ8.3, pZ0.007). Females mounted marginally higher mean diphtheria responses than males (FZ3.2, pZ0.082), but the sex!f interaction was not significant ( pZ0.30). Other covariates were eliminated (all pO0.1; final model: F 4,37 Z7.5, p!0.001, R 2 Z0.41). The final model remained quantitatively similar after excluding a single outbred individual. After excluding four highly inbred individuals, the estimated inbreeding load increased to K4.4, although this effect no longer approached statistical significance due to increased confidence limits (95% CL K13.4 to 4.7, pZ0.33).

DISCUSSION
In free-living song sparrows, PHA, tetanus and diphtheria responses all varied to some extent with individual f. Estimated inbreeding loads for PHA and tetanus responses did not vary markedly among years or sparrows of different ages, but apparently differed between females and males. PHA response declined with f in both sexes, but declined more markedly in males than in females in 2003 and across all data. Tetanus response declined with f in females but not in males in 2004 and 2005 and across all data. We had sufficient statistical power to detect the inbreeding load for tetanus response estimated in females had it also occurred in males (observed rz0.56, nZ40, powerz0.97, given that overall variance in f was similar in both sexes). There is no clear expectation that apparent female-specific inbreeding depression should have arisen because measurement error in tetanus response or f was consistently greater in males (e.g. due to paternity error, see §2), and collinearity of explanatory variables did not differ between the sexes. Therefore, while it remains possible that the repeatable sex-specific inbreeding depression observed in tetanus response simply reflects sampling variance or stochastic male-biased error in f, our Table 1. Descriptive statistics of song sparrows challenged with phytohaemagglutinin (PHA) and diphtheria-tetanus vaccine (DTV) during 2002-2005. (N t , N juv  Although the sex!f interaction no longer explained a significant proportion of variation in PHA response after excluding inbred individuals, and confidence intervals around parameter estimates were in some cases increased, overall patterns of declining immune responses with increasing f remained clear. Apparent inbreeding depression in immune response was therefore not solely attributable to heterosis in offspring of immigrants-native pairings, or to particularly weak immunity in the few highly inbred individuals that were tested. Inbreeding depression is thought to reflect increased genetic homozygosity and consequently increased expression of deleterious recessive alleles and reduced expression of overdominance (Charlesworth & Charlesworth 1999). Inbreeding depression in immune response might therefore reflect direct effects of increased homozygosity at key response loci. However, since immune responses can be energetically demanding and depend on an individual's physiological and hormonal state (Hasselquist et al. 1999;Ots et al. 2001;Owen-Ashley et al. 2004), inbreeding depression may also reflect more general detrimental effects of inbreeding on metabolic or endocrinological pathways. PHA response declined with body condition in song sparrows, indicating condition dependence in this component of immunity (see also Alonso-Alvarez & Tella 2001;. Tetanus and diphtheria responses did not vary with body condition (see also . However, the apparent sex-specific inbreeding depression observed in tetanus response may suggest an indirect physiological mechanism. Direct phenotypic consequences of inbreeding should not differ between females and males unless key loci are sex linked (Falconer & Mackay 1996), and major MHC loci that mediate major aspects of humoral (and cellmediated) immunity are autosomal in birds (Edwards et al. 1999, although it is possible that sex-linked loci are also involved). In contrast, sex-specific variation in immune response resulting from sex-specific variation in Table 2    life-history allocation and physiology is widely predicted and documented (Zuk & McKean 1998;Hasselquist et al. 1999;Stoehr & Kokko 2006). Predicted variation, however, primarily relates to mean responses and patterns of seasonal variation. The apparent sex-specific inbreeding depression observed in immune response is therefore unexpected. It is not clear why inbreeding effects should differ between males and females or why the direction of sex-specific effects should differ among components of immunity. Since investment in immune responsiveness may trade-off against investment in other life-history components (Tella et al. 2002;Adamo 2004), one final possibility is that relatively inbred individuals strategically invest less in immunity. Such reduced allocation could be viewed in a life-history context as an adaptive response to increased intrinsic mortality rates in inbred individuals, and could conceivably be modulated in response to direct physiological or endocrinological effects of inbreeding. Such strategic investment patterns may provide one possible explanation for sex-specific allocation to immunity with respect to f. Irrespective of the underlying mechanisms, our results broadly concur with an increasing body of evidence that parasite resistance commonly declines with reduced heterozygosity and/or genetic diversity in vertebrates (e.g. Acevedo-Whitehouse et al. 2005;MacDougall-Shackleton et al. 2005;Westerdahl et al. 2005;Whiteman et al. 2006, although see Giese & Hedrick 2003Cô té et al. 2005). However, most studies, particularly those on free-living individuals, have related parasite resistance or immunity to heterozygosity or genetic diversity measured across a few microsatellite or MHC loci rather than individual f. In many cases, it is therefore uncertain whether observed effects reflect genome-wide heterozygosity or inbreeding per se, as opposed to heterozygosity at specific focal or linked loci (Balloux et al. 2004;Slate et al. 2004). Furthermore, few studies have related any measure of inbreeding, heterozygosity or genetic diversity to an individual's response to any novel, non-specific immune challenge, and no general pattern has emerged. For example, one generation of sib-sib mating did not reduce encapsulation response in bumble-bees (Bombus terrestris, Gerloff et al. 2003). PHA response but not the antibody response to sheep red blood cells (SRBC response) increased with multilocus heterozygosity in captive house finches (Carpodacus mexicanus, Hawley et al. 2005). Neither PHA nor SRBC response varied with MHC diversity in captive house sparrows (Passer domesticus, Bonneaud et al. 2005). In contrast, PHA, tetanus and diphtheria responses all declined to some extent with increasing f in song sparrows. Further studies are required to characterize the patterns and causes of this variation and identify the circumstances under which inbreeding and/or genetic diversity is most likely to influence immune response.
It is likely to be simplistic to interpret an individual's response to any single immune challenge as a measure of general parasite resistance (Norris & Evans 2000;Adamo 2004). Since only ten song sparrows experienced both PHA and DTV and these individuals experienced these challenges in different years, we have little power to assess whether PHA, diphtheria and tetanus responses were correlated within individuals. Notwithstanding these limitations, individual PHA and tetanus responses tended to be negatively correlated after controlling for year, inter-sample period and maternal vaccination history (r s ZK0.55, nZ10, pZ0.09). However, our repeated challenge experiments suggested that PHA, tetanus and diphtheria responses, which to some degree reflect the main components of avian immunity (Norris & Evans 2000;Martin et al. 2006), show broadly similar declines with increasing f in free-living song sparrows. Inbreeding    loads were estimated as 2.1-7.0, 7.2-7.8 and 2.8 lethal equivalents per gamete for PHA, tetanus (in females) and diphtheria responses, respectively. Individual f explained a remarkably high proportion of variation, particularly in PHA response. Furthermore, given the likely paternity error in the Mandarte pedigree, these figures may underestimate the true magnitude of inbreeding depression Kruuk et al. 2002;Marr et al. in press). We cannot quantitatively compare the inbreeding loads observed in PHA, diphtheria and tetanus responses, since these different components of immunity were measured in different individuals in different years. Among-component variation in inbreeding depression cannot therefore be distinguished from among-year variation, such as might result from among-year variation in environmental severity (e.g. Keller et al. 2002). However, our data indicate that the magnitude of inbreeding depression can be of the same order in multiple components of immunity, at least in females, and equals or exceeds that commonly observed in major fitness components, including in song sparrows on Mandarte (Keller 1998;Kruuk et al. 2002).
The substantial inbreeding depression observed in immune response seems likely to be biologically significant. Variation in PHA response has been linked to substantial changes in life-history allocation , and an average effect size of rZ0.43 has been reported for the relationship between immune response and survival in birds (Møller & Saino 2003; see also . Given this effect size and the inbreeding depression in immunity observed in song sparrows, even moderate inbreeding ( fz0.1) would equate to a substantial reduction in adult survival. Therefore, interpreted directly, our data indicate that even moderate inbreeding may impair immunity sufficiently to reduce fitness in free-living birds (assuming exposure to a relevant parasite and notwithstanding that fitness may not be maximized by maximizing immune responsiveness, Adamo 2004;Viney et al. 2005). These data support the general suggestion that host-parasite interactions can be influenced by host genotype and specifically highlight the likely importance of dominance genetic effects such that the parasite resistance declines with host genome-wide heterozygosity. Resulting parasite-mediated selection against inbred hosts may help maintain genetic diversity in natural populations (Potts & Wakeland 1990;Coltman et al. 1999), influence the evolution and maintenance of inter-sexual selection (Reid et al. 2005) and threaten the persistence of small populations (Altizer et al. 2003;Whiteman et al. 2006). The apparent sex-specific inbreeding depression observed in PHA response and particularly tetanus response is unexpected, and implies that further compound consequences of host-parasite interactions might arise. For example, given sex-specific inbreeding depression in immunity, parasites might cause biased sex ratios in small populations of inbred hosts, thereby reducing effective population size and further increasing inbreeding, genetic drift and stochastic extinction risk. Further studies are required to investigate the prevalence, magnitude and causes of sex-specific inbreeding depression in parasite resistance and explore the demographic and evolutionary implications.