Skip to main content
  • Research article
  • Open access
  • Published:

Context-dependent associations between heterozygosity and immune variation in a wild carnivore



A multitude of correlations between heterozygosity and fitness proxies associated with disease have been reported from wild populations, but the genetic basis of these associations is unresolved. We used a longitudinal dataset on wild Galapagos sea lions (Zalophus wollebaeki) to develop a relatively new perspective on this problem, by testing for associations between heterozygosity and immune variation across age classes and between ecological contexts.


Homozygosity by locus was negatively correlated with serum immunoglobulin G production in pups (0–3 months of age), suggesting that reduced genetic diversity has a detrimental influence on the early development of immune defence in the Galapagos sea lion. In addition, homozygosity by locus was positively correlated with total circulating leukocyte concentration in juveniles (6–24 months of age), but only in a colony subject to the anthropogenic environmental impacts of development, pollution and introduced species, which suggests that reduced genetic diversity influences mature immune system activity in circumstances of high antigen exposure.


These findings demonstrate the environmental context-dependency of the phenotypic expression of immune variation, which is implicit in the theory of ecoimmunology, but which has been rarely demonstrated in the wild. They also indicate that heterozygosity may be linked to the maintenance of heterogeneity in mammalian immune system development and response to infection, adding to the body of evidence on the nature of the mechanistic link between heterozygosity and fitness.


Associations between indices of multi-locus heterozygosity and many fitness-related traits, from birth weight [1] to song complexity [2], have been reported in wild animal populations [3]. Heterozygosity is of interest because it quantifies within individual genetic diversity, and is likely to be related to inbreeding in some way – although in exactly which way is contentious. Inbreeding has been known to have deleterious effects on fitness for a long time [4, 5], and these effects have been well documented in captive, laboratory and domesticated animals [6]. Historically, the occurrence of inbreeding depression in natural populations was disputed, as evidence was scarcer from the wild than other contexts. It is now thought that this is likely to have been due to detection difficulties [7], and that inbreeding can affect fitness in natural settings [811]. The consequences of inbreeding in the wild are of particular interest to biologists working to conserve declining populations or threatened species [9, 10, 12]. Reduced genetic diversity may interact with extrinsic stressors, such as disease, to influence population dynamics [1215]. Moreover, the extent to which the effects of reduced genetic diversity vary across ecological contexts within species has been relatively poorly explored.

An obstacle to the study of inbreeding in the wild is that there is no problem-free way to estimate levels of inbreeding in a population about which we have incomplete information. For example, the estimation of inbreeding coefficients from pedigrees is complicated by the assumption that founding individuals are outbred and unrelated [16], by the difficulty of accounting for the chance events of Mendelian segregation [17], and by the rarity of multi-generation pedigrees for natural populations [18]. An alternative to calculation from pedigrees is to use data from neutral genetic markers to either summarise genetic diversity as multi-locus heterozygosity (MLH), or to estimate inbreeding coefficients directly using relatedness algorithms, as a proxy for inbreeding value. These latter methods have the advantage of being calculable for individuals sampled from a single cohort, and are therefore feasible in many study systems, especially given recent advances in sequencing technology and the availability of resources for the study of genetics.

However, whether indices of MLH calculated using small numbers of loci represent variation across the genome well enough for inbreeding depression to be invoked as the cause of heterozygosity-fitness correlations (HFCs) is unresolved [1921]. Empirical comparisons in populations for which both pedigrees and marker data are available have shown that pedigree-derived inbreeding estimates and indices of MLH are not well correlated, especially when the mean and variance of inbreeding value are low [2226]. However, whether this suggests that MLH indices are less accurate than pedigree-derived estimates or that both are inaccurate in different ways is unknown [17].

An alternative explanation of HFCs to inbreeding depression is that of linkage between one or a few neutral markers and functional genes under balancing selection, which could give rise to the frequently observed pattern of heterosis [19, 2730]. A reappraisal noted, however, that the uneven contribution of loci to HFCs is implicit in inbreeding theory [16], as weak inbreeding is not expected to lead to detectable identity disequilibrium [16, 31]. Therefore, the effects of weak inbreeding may be present even if heterozygosity is not correlated across loci – an observation that has previously been used to rule out the presence of inbreeding – and the finding that HFCs are driven by variation at one or a few loci does not preclude weak inbreeding as their underlying cause. In addition, the frequentist statistical tests that have been used to identify single-locus contributions to HFCs suffer from problems of power [30] and non-independence [16].

To date, variation in the detection and strength of HFCs has largely been attributed to methodology, and there is consensus that increases in numbers of samples and markers will help delineate the mechanisms that give rise to HFCs [3, 19]. Another way in which HFC studies are working towards this same goal is through the inclusion of ecological heterogeneity in their study designs [3]. This approach co-opts natural variation in the expression of the consequences of inbreeding, for example, through episodic heterozygote advantage [32] or the effects of environmental stress [33, 34], to further understanding of the mechanisms that drive them in natural settings. This study takes such an approach, by examining the relationships between heterozygosity, growth, body condition and immune variation in the Galapagos sea lion.

The Galapagos sea lion (Zalophus wollebaeki) is endemic to the Galapagos archipelago and has an estimated population size of 20,000–40,000 animals [35], which is spatially and genetically structured amongst small colonies (20–500 animals) [36, 37]. It has a polygynous mating system [38, 39] and is philopatric [80]. Under these conditions we might expect the inevitable imprecision of marker-based inbreeding estimates to be minimised, as these are both traits that limit the degree to which alleles are mixed within populations, and may therefore promote the necessary variation in heterozygosity required for the statistical detection of underlying patterns [16, 19, 24].

Growth and condition have been empirically linked with survival probability and reproductive success in many taxa [40], including marine mammals [41, 42]. These are likely to be important components of fitness in the Galapagos sea lion as they are related to reproduction, but also because they are related to the ability to resist starvation, which is important in an ecosystem in which fluctuations in marine productivity are driven by unpredictable environmental variation [43, 44]. The concentration of the most common class of circulating mammalian antibody, Immunoglobulin G (IgG), has been linked to survival probability in the grey seal (Halichoerus grypus) [42], and the capacity to regulate immunity effectively may be an important determinant of Galapagos sea lion fitness given disease risk concerns [35, 36, 45]. In addition, a recent study reported a link between diversity at an immune gene locus (MHC-DRB) and fitness in the Galapagos sea lion [46]. The data presented here complement this immunogenetic approach, through the inclusion of phenotypic variation in immunity.

The many influences on immune activity in natural populations mean that phenotypic measures of immune variation are unlikely to be related to fitness in straightforward ways [47]. In this study we make use of a previously described longitudinal dataset and sample archive on immune variation in the Galapagos sea lion [48, 49] to take within individual variation, and covariation with other aspects of life history, into account. This previous work describes immune variation in the context of growth and development over the first two years of life in the Galapagos sea lion [48, 49]. Here we build on this by testing for relationships between heterozygosity and immune activity while controlling for the potentially confounding influences of growth, body condition and age class in two breeding colonies, one located in a town and subject to high levels of anthropogenic environmental impact, and the other in a relatively undisturbed habitat [48, 49].



We extracted genomic DNA from 166 Galapagos sea lion skin samples and amplified 23 polymorphic microsatellite loci previously developed for the Galapagos sea lion and other pinniped species [39, 67]. All genetic analyses were carried out at the University of Bielefeld, Germany, as described in [39, 67]. Sequencing was performed on an Applied Biosystems 3130 Sequencer (Life Technologies) and genotyping was performed in GENEMARKER (Soft Genetics, USA). We tested for deviations from Hardy-Weinberg equilibrium (HWE) using Monte Carlo simulations in the adegenet package [68] in R 2.14.1 [69], and for identity disequilibrium using the g 2 measure [31].

Main statistical analysis

Samples were collected in two Galapagos sea lion colonies: one on an uninhabited island (control colony, Santa Fe), and one located in a town (human-impacted colony, San Cristobal) [48]; and from sea lions in two ecologically distinct age-classes: pups (0–3 months old), which are restricted to land and are dependent on their mothers for nutrition; and juveniles (6–24 months old), which swim out to sea and are capable of foraging independently [70]. We sampled pups shortly after birth and at 3 months of age, between which two ages they undergo a growth spurt [44]; and juveniles at 6, 12, 18 and 24 months of age. We captured sea lions using hoop nets and briefly restrained them in a prone position without the use of chemical immobilization, and without causing harm, by following the capture protocol in [70]. All work was approved by the Zoological Society of London Ethics Committee, and carried out under Galapagos National Park permits PC-18-09, N046-2009-PNG, N101-2010-PNG and N032-2010-PNG, which covered all fieldwork, capture and sample protocols. We used two measures of immune variation: immunoglobulin G (IgG) concentration and total leukocyte concentration, as they were highlighted in previous analyses as most likely to vary meaningfully with other aspects of Galapagos sea lion life history [48, 49].

Galapagos sea lion pups undergo rapid growth and physiological development during the sampled period of their early development [44]. In order to take these changes into account, and given that pups were only sampled at two time points during this period, we calculated absolute changes in body mass (kg), body length (cm), IgG concentration (mg/ml) and total leukocyte concentration (109/l) between shortly after birth and 3 months of age. The possibility of phenotypic correlation [54] means that growth and changes in immune measures may covary, and we have shown that the direction of these associations varies between colonies in the Galapagos sea lion [49]. Therefore, in fast-growing pups, we partitioned variation in changes in each immune measure into subsets that were correlated with changes in body length and body mass in different ways using principal components analyses, carried out separately for each colony. For each colony and immune measure we fit generalised linear models with principal components that explained ≥ 5 % of the variation as response variables to homozygosity weighted by locus (HL) [52], sex and their interaction as explanatory terms, removing interactions if they were non-significant [71]. This amounted to eight statistical models fitted to pup data: two principal components, from two immune measures, in two colonies. This approach addresses the problem of the potentially confounding influence of phenotypic correlation on associations between HL and changes in immune measures in fast-growing pups, as it partitions the variation in changes in immune measures into components that are correlated with different kinds of growth, and allows for comparison between their association with homozygosity.

We chose HL as the most appropriate measure of heterozygosity for the main statistical analyses, so that they could be compared with other published results (e.g. [46]), and because the distribution of the variation in HL was amenable to modelling in a generalised linear model (GLM) framework. However, we undertook a detailed exploration of the biases inherent in different estimates of heterozygosity and inbreeding using simulation analysis to provide context for these results, and other analyses that use measures of heterozygosity more generally (Additional file 1: Supplementary Text 1.1–2, Table S2–3, Figure S1).

In comparison with pups, relatively little growth occurs in juvenile Galapagos sea lions between the ages of 6 and 24 months [44], fewer physiological changes take place, and body mass and length are more closely correlated than in younger animals [49]. In addition, we sampled juveniles at four rather than two time points. The nature of the juvenile data, therefore, allowed us to take a simpler approach to correcting for phenotypic correlation, which we did by including body mass as an explanatory variable. Separately for each colony, we fitted generalised linear mixed models (GLMMs) with each immune measure as a response variable and HL, body mass, sex and the interaction between HL and sex as explanatory terms. We included individual identity as a random effect to account for the pseudoreplication implicit in the repeated sampling of individuals. This amounted to four statistical models fitted to juvenile data, which covered two colonies and two immune measures. We compared models with and without the interaction between sex and HL using likelihood ratios tests [72]. The analysis of juvenile data was therefore analogous to that of pup data, but did not require partitioning by principal components analysis. We checked all models for signs of heteroscedasticity, heterogeneity of variance, non-normality of error and the disproportionate influence of outliers.


Of 23 amplified loci, one (ZcwG06) [50] showed significant deviation from Hardy-Weinberg equilibrium (Additional file 1: Table S1) and was omitted from all analyses. All individuals were genotyped twice, samples that produced allele mismatches were excluded from the analyses (n = 2), and the data are provided in Additional file 2. Homozoygosity by locus (HL) did not vary significantly with sex (GLM: n = 52, estimate = −0.50, standard error, SE = 0.42, p = 0.23), colony (GLM: n = 52, estimate = −0.52, SE = 1.57, p = 0.74) or estimated birthdate (GLM: n = 52, estimate = 0.012, SE = 0.045, p = 0.78) in pups; or with sex (GLM: n = 73, estimate = 0.03, SE = 0.02, p = 0.137) or colony (GLM: n = 73, estimate = −0.03, SE = 0.02, p = 0.09) in juveniles. There was no evidence of identity disequilibrium (g 2 = 0.0033, standard deviation, SD = 0.0027, p < 0.888, 1000 iterations).

The mean change in pup body mass from shortly after birth until 3 months of age was 6.01 kg (SD = 1.81 kg) and mean change in pup body length over the same period was 14.59 cm (SD = 5.19 cm). Change in pup body length was greater in the human-impacted colony than in the control colony (GLM: n = 40, estimate = −3.93, SE = 1.54, p = 0.014), but there was no significant difference between sexes (GLM: n = 40, estimate = 3.14, SE = 1.59, p = 0.056). Change in body mass was greater in male pups than female pups (GLM: n = 40, estimate = 1.39, SE = 0.54, p = 0.013), but there was no significant difference between colonies (GLM: n = 40, estimate = −0.99, SE = 0.56, p = 0.083).

For IgG, the first component produced by principal component analyses for each colony was positively correlated with changes in IgG concentration, changes in body mass, and changes in body length (0.93, 0.33, 0.09 and 0.92, 0.33, 0.15 in the human-impacted and control colonies respectively; Table 1); and the second component was positively correlated with changes in IgG concentration but negatively correlated with changes in body mass and length (0.34, −0.83, −0.43 and 0.36, −0.91, −0.17 in the human-impacted and control colonies respectively; Table 1). This demonstrates the separation of variation in IgG into that which is positively correlated with growth through phenotypic correlation, represented by the first principal components, and that which is not, represented by the second principal components. For total leukocyte concentration, the principal component analyses produced a similar pattern to IgG, and for both IgG and total leukocyte concentration, the differences between colonies were minimal (Table 1; Additional file 1: Figure S2).

Table 1 Principal component analysis results for IgG (total immunoglobulin G concentration, mg/ml) in the human-impacted

In pups homozoygosity by locus (HL) was negatively associated with the second principal component derived from IgG in the control colony (GLM: estimate = 16.75, SE = 7.58, p = 0.039; Tables 2 and 3; Fig. 1a). This principal component was most strongly correlated with changes in body length (−0.91), next with changes in IgG (0.36) and lastly with changes in mass (−0.17; Table 1). The pattern of variation observed with HL suggests that individuals with high levels of heterozygosity grew less but produced more IgG. Sex was included in the statistical models to account for its potentially confounding influence but was never significant as part of an interaction with HL. In juveniles the interaction between HL and sex was not retained in any model, and HL was retained in a single model, in which it was positively correlated with total leukocyte concentration in the human-impacted colony (GLMM: estimate = 5.98, SE = 2.21; likelihood ratio test, LRT: χ2 = 7.29, p = 0.006; Tables 2 and 3 and Fig. 1b).

Table 2 Statistical associations between homozygosity weighted by locus (HL) and principal components derived from immune and growth variables and in pups (linear models), and immune variables in juveniles (linear mixed effect models)
Table 3 Statistical associations between homozygosity weighted by locus (HL) and principal components derived from immune and growth variables and immune variables in juveniles (linear mixed effect models)
Fig. 1
figure 1

a The relationship between homozygosity weighted by locus (HL) and the second principal component derived from pup data on changes in body length, mass and IgG concentration in the control colony (GLM: estimate = 16.75, SE = 7.58, p = 0.039; Tables 2 and 3); (b) the relationship between homozygosity weighted by locus (HL) and total leukocyte concentration (109/l) in juveniles in the human-impacted colony (GLMM: estimate = 5.98, SE = 2.21; LRT: χ2 = 7.29, p = 0.006; Tables 2 and 3)


This study has two main findings: first, two different kinds of immune variation were associated with heterozygosity in a wild mammal; second, the occurrence and nature of these associations varied between ecologically distinct circumstances. The latter allows for the possibility that human impact may influence the relationship between genotype and immune phenotype in this species [53]. These results – as is necessarily the case with studies of species such as the Galapagos sea lion in the wild – are correlative, and preclude the inference of direct causation. Also, given that the colony on San Cristobal (human-impacted) is unique – as discussed in detail in [48, 49] – it was not possible to test for interactions between heterozygosity and human-impact on immune variation across a larger sample of colonies. Nevertheless, the results reported here answer a call to incorporate more ecologically relevant variation into the exploration of the impacts of variation in heterozygosity in natural populations [3]. In addition, given recent interest in the relationship between variation at immune gene loci and fitness in wild organisms, including in the Galapagos sea lion [46], the quantification of physiological variation in immunity may prove a valuable tool for understanding the complex mechanisms that give rise to correlations between such immunogenetic variation and fitness.

Immunoglobulin G (IgG) and heterozygosity

Physiological measures of immune variation may be condition-dependent [54], and variation in IgG concentration is correlated with growth in Galapagos sea lion pups [49]. We partitioned variation in growth and IgG production into principal components to take this co-linearity into account. The majority of Galapagos sea lion pups grow [44] and produce IgG [48] during their first three months of life. Given the lack of evidence for infection or disease in the sampled Galapagos sea lion pups (Additional file 1: Supplementary text 1.3–1.4), and a field experiment that showed IgG is not passed from mother to pup during the sampled stage of development in the Galapagos sea lion [48], the observed IgG production is likely to be driven by post-natal antigen exposure, which stimulates the build-up of protective baseline populations of lymphocytes and antibodies [42, 55].

Homozygosity weighted by locus (HL) was not associated with the first principal component derived from changes in IgG concentration, body mass and body length. This suggests that heterozygosity may not be straightforwardly related to fitness, as it was not correlated with the axis of variation that most likely describes phenotypic correlation between traits (Tables 2 and 3) [54]. The second principal component was negatively correlated with HL, which, given the loadings (Tables 2 and 3), suggests that relatively heterozygous individuals grew less but produced more IgG, which is congruent with findings from other systems [14, 56]. Such a pattern could arise through variation in resource allocation or acquisition between sea lion pups with different levels of heterozygosity [57]. However, this does not explain why relatively homozygous pups appeared to grow more. The answer to this may lie in the relative magnitude of the length and mass loadings on the two components: there is a small difference between these in the first component (0.33 and 0.15 for length and mass respectively; Tables 2 and 3), as we would expect under phenotypic correlation; but length loads more heavily on the second component than mass (−0.91 and −0.17 for length and mass respectively; Tables 2 and 3). This suggests that the first component better represents overall growth that is more likely to be positively correlated with fitness, while the second represents elongation, which corresponds to a decrease in body condition [40].

Pups from the human-impacted colony varied less in terms of heterozygosity than pups from the control colony (Additional file 1: Figure S1a-b), so the power to detect an effect may have been lower in the former, which could explain why the above pattern was only observed in the control colony. An alternative explanation is that pups (or mothers) in the human-impacted colony did not vary resource allocation according to heterozygosity in the same way that those in the control colony did. Pups in the human-impacted colony produced more IgG than pups from the control colony during the sampled period of development, which may be due to higher post-natal antigen exposure in the human-impacted colony [48]. If this is the case, environmental stimulation of IgG production in the human-impacted colony may have over-ridden genetic influence on IgG production, which could obscure the statistical signal of genetic influence on IgG production in the human-impacted colony.

Leukocytes and heterozygosity

Total leukocyte concentration in mammals is likely to vary over shorter timescales than IgG concentration, and to be more prone to fluctuation, as leukocytes are less cumulative in the blood than antibodies [58]. In young mammals, though, (including sea lions) total leukocyte concentration changes with age, as cell populations develop during post-natal growth [48, 59]. Therefore, in pups we used the same principal component analysis approach as we did with IgG. However, heterozygosity was not statistically associated with any of the principal components derived from total leukocyte concentration in pups (Table 1).

The positive correlation between homozygosity and total leukocyte concentration in juveniles in the human-impacted colony suggests that low within individual genetic diversity may be associated with an increase in the number of circulating leukocytes, under certain ecological circumstances. Total leukocyte concentration in juveniles is more likely to represent mature immune system activity than early-developmental leukocyte production, given that the latter happens during only a short period following birth in pinnipeds [60, 61]. Low heterozygosity has been associated with increased susceptibility to parasitism [1, 13, 62], so may have compromised aspects of innate immunity, e.g. [14, 56], in these juvenile sea lions, which could have led to relatively high infection rates, e.g. [63]. This may have been observed only in the human-impacted colony as this is where infection risk is likely to be relatively high, given the threat of pathogen exposure [36, 45, 52]; an effect that may be compounded by relatively high stress and pollution levels in the human-impacted colony [48].

Heterozygosity and inbreeding

A general assumption of heterozygosity-fitness-correlation (HFC) studies is that phenotypic variation is linearly related to fitness and therefore under directional selection [3]. Immune variables such as IgG concentration are unlikely to be so, due to the damage caused by immunopathology and other trade-offs involving immunity [57]; moreover, different kinds of antibody responses have been shown to be under different kinds of selection: in blue tits (Parus caeruleus) “…primary responsiveness to diphtheria was subject to stabilizing selection, whereas secondary responsiveness to tetanus was subject to positive directional selection” [47]. Correlations between such traits and heterozygosity are nonetheless useful, as when, and under what ecological circumstances, they are detectable can offer insight into their mechanisms [3]. Immune traits are unlikely to be influenced by as many genes as complex composite traits such as survival and lifetime reproductive success (to which they contribute), and therefore may be more amenable to a candidate gene approach to the investigation of HFC mechanisms, e.g. [64]. However, this also implies that immune traits are less likely to be influenced by genome-wide inbreeding [16], which raises the question of whether heterozygosity is acting as a proxy for inbreeding in the Galapagos sea lion.

As g 2 was not significantly different from zero (g 2 = 0.0033, SD = 0.0027, p < 0.888, 1000 iterations), identity disequilibrium of our marker set was low [31]. This suggests an absence of inbreeding (but see [21]) in the Galapagos sea lion, which agrees with other assessments using different proxies for inbreeding in the Galapagos sea lion [37, 46]. If the correlations between heterozygosity and immune variation described here are not due to inbreeding, it is possible they are the result of the other widely suggested explanation for HFCs: linkage between one or a few presumed-neutral markers to functional genes under balancing selection [19]. A variety of statistical approaches have been proposed to identify heterogeneities in the contributions of loci to HFCs [19, 27, 28, 65], but problems of power and non-independence mean their results are difficult to assess [16, 30]. In addition, unbalanced contributions of one or a few loci to HFCs are implicit in inbreeding theory, as weak inbreeding is not expected to lead to detectable identity disequilibrium [16, 31]. Therefore, an absence of identity disequilibrium (or other equivalent measure of correlation of heterozygosity across loci) does not preclude the presence of weak inbreeding [16].


Although the links between heterozygosity and inbreeding are unresolved [20], the context-dependent expression of patterns of covariance between heterozygosity and other traits can yield useful and interesting information. We offer these results as examples of covariation between aspects of genotype and phenotype in a natural population, which demonstrate the context-dependency of immune variation expression predicted by ecoimmunology, and which may provide new perspective on heterozygosity-fitness correlations, by allowing for the formation of increasingly detailed hypotheses on the mechanisms that drive them. Irrespective of their causes and mechanisms, and in light of recent findings in the Antarctic fur seal (Arctocephalus gazelle) [12], our results may also have implications for the management of declining populations, as they allow for the possibility that genetic effects interact with the two major threats to the Galapagos sea lion: climate-driven changes in food availability [35, 44, 66] through effects on growth and condition; and disease threat [36], through effects on immunity; both of which could act to compound the risk of decline and extinction.

Availability of supporting data

The genetic data used in the analyses are available in Additional file 1.



Multi-locus heterozygosity


Heterozygosity-fitness correlations


Immunoglobulin G


Hardy-Weinberg equilibrium


Homozygosity weighted by locus


Generalised linear model


Generalised linear mixed model


Standard error


Standard deviation


Likelihood ratio test


  1. Coltman DW, Bowen WD, Wright JM. Birth weight and neonatal survival of harbour seal pups are positively correlated with genetic variation measured by microsatellites. Proc R Soc B Biol Sci. 1998;265:803–9. doi:10.1098/rspb.1998.0363.

    Article  CAS  Google Scholar 

  2. Marshall RC, Buchanan KL, Catchpole CK. Sexual selection and individual genetic diversity in a songbird. Proc R Soc B Biol Sci. 2003;270:S248–50. doi:10.1098/rsbl.2003.0081.

    Article  Google Scholar 

  3. Chapman JR, Nakagawa S, Coltman DW, Slate J, Sheldon BC. A quantitative review of heterozygosity-fitness correlations in animal populations. Mol Ecol. 2009;18:2746–65. doi:10.1111/j.1365-294X.2009.04247.x.

    Article  CAS  PubMed  Google Scholar 

  4. Darwin CR. The effects of cross and self fertilization in the vegetable kingdom. London, UK: J. Murray & Co.; 1876.

    Google Scholar 

  5. Wright S. Coefficients of inbreeding and relationship. Am Nat. 1922;56:330–8.

    Article  Google Scholar 

  6. Frankham R. Inbreeding and extinction: a threshold effect. Conserv Biol. 1995;9(4):792–9. doi:10.1046/j.1523-1739.1995.09040792.x.

    Article  Google Scholar 

  7. Frankham R, Ballou J, Briscoe D. An introduction to conservation genetics. Cambridge, UK: Cambridge University Press. 2002

  8. Crnokrak P, Roff DA. Inbreeding depression in the wild. Heredity. 1999;83:260–70. doi:10.1038/sj.hdy.6885530.

    Article  PubMed  Google Scholar 

  9. Keller LF, Waller DM. Inbreeding effects in wild populations. Trends Ecol Evol. 2002;17(5):230–41.

    Article  Google Scholar 

  10. Brekke P, Bennett PM, Wang J, Pettorelli N, Ewen JG. Sensitive males: inbreeding depression in an endangered bird. Proc R Soc B Biol Sci. 2010;277(1700):3677–84. doi:10.1098/rspb.2010.1144.

    Article  Google Scholar 

  11. Walling CA, Nussey DH, Morris A, Clutton-Brock TH, Kruuk LEB, Pemberton JM. Inbreeding depression in red deer calves. BMC Evol Biol. 2011;11:318–31. doi:10.1186/1471-2148-11-318.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Forcada J, Hoffman JI. Climate change selects for heterozygosity in a declining fur seal population. Nature. 2014;511(7510):462–5. doi:10.1038/nature13542.

    Article  CAS  PubMed  Google Scholar 

  13. Acevedo-Whitehouse K, Gulland FMD, Greig D, Amos W. Disease susceptibility in California sea lions. Nature. 2003;422:35. doi:10.1038/422035a.

    Article  CAS  PubMed  Google Scholar 

  14. Reid JM, Arcese P, Keller LF, Elliott KH, Sampson L, Hasselquist D. Inbreeding effects on immune response in free-living song sparrows (Melospiza melodia). Proc R Soc B Biol Sci. 2007;274:697–706. 10.1098/rspb.2006.0092.

    Article  Google Scholar 

  15. Acevedo-Whitehouse K, Duffus ALJ. Effects of environmental change on wildlife health. Philos Trans R Soc B. 2009;364:3429–38. doi:10.1098/rstb.2009.0128.

    Article  Google Scholar 

  16. Szulkin M, Bierne N, David P. Heterozygosity-fitness correlations: a time for reappraisal. Evolution. 2010;64(5):1202–17. doi:10.1111/j.1558-5646.2010.00966.x.

    PubMed  Google Scholar 

  17. Forstmeier W, Schielzeth H, Mueller JC, Ellegren H, Kempenaers B. Heterozygosity-fitness correlations in zebra finches: microsatellite markers can be better than their reputation. Mol Ecol. 2012;21:3237–49. doi:10.1111/j.1365-294X.2012.05593.x.

    Article  PubMed  Google Scholar 

  18. Frentiu FD, Clegg SM, Chittock J, Burke T, Blows MW, Owens IPF. Pedigree-free animal models: the relatedness matrix reloaded. Proc R Soc B Biol Sci. 2008;275:639–47. 10.1098/rspb.2007.1032.

    Article  Google Scholar 

  19. Balloux F, Amos W, Coulson T. Does heterozygosity estimate inbreeding in real populations? Mol Ecol. 2004;13:3021–31. doi:10.1111/j.1365-294X.2004.02318.x.

    Article  CAS  PubMed  Google Scholar 

  20. Pemberton JM. Wild pedigrees: the way forward. Proc R Soc B Biol Sci. 2008;275:613–21. doi:10.1098/rspb.2007.1531.

    Article  CAS  Google Scholar 

  21. Miller JM, Coltman DW. Assessment of identity disequilibrium and its relation to empirical heterozygosity fitness correlations: a meta-analysis. Molecular Ecology, in press. 2014. doi:10.1111/mec.12707

  22. Slate J, David P, Dodds KG, Veenvliet BA, Glass BC, Broad TE, et al. Understanding the relationship between the inbreeding coefficient and multilocus heterozygosity: theoretical expectations and empirical data. Heredity. 2004;93:255–65. doi:10.1038/sj.hdy.6800485.

    Article  CAS  PubMed  Google Scholar 

  23. Overall A, Byrne K, Pilkington J, Pemberton J. Heterozygosity, inbreeding and neonatal traits in Soay sheep on St Kilda. Mol Ecol. 2005;14:3383–93. doi:10.1111/j.1365-294X.2005.02682.x.

    Article  CAS  PubMed  Google Scholar 

  24. Grueber CE, Waters JM, Jamieson IG. The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species. Mol Ecol. 2011;20:67–79. doi:10.1111/j.1365-294X.2010.04930.x.

    Article  PubMed  Google Scholar 

  25. Wetzel DP, Stewart IRK, Westneat DF. Heterozygosity predicts clutch and egg size but not plasticity in a house sparrow population with no evidence of inbreeding. Mol Ecol. 2012;21:406–20. doi:10.1111/j.1365-294X.2011.05380.x.

    Article  CAS  PubMed  Google Scholar 

  26. Hoffman JI, Simpson F, David P, Rijks JM, Kuiken T, Thorne MAS, et al. High-throughput sequencing reveals inbreeding depression in a natural population. Proceedings of the National Academy of Sciences, online February 28th 2014. 2014. doi:10.1073/pnas.1318945111.

    Google Scholar 

  27. Acevedo-Whitehouse K, Spraker TR, Lyons E, Melin SR, Gulland F, Delong RL, et al. Contrasting effects of heterozygosity on survival and hookworm resistance in California sea lion pups. Mol Ecol. 2006;15(7):1973–82. doi:10.1111/j.1365-294X.2006.02903.x.

    Article  CAS  PubMed  Google Scholar 

  28. Acevedo-Whitehouse K, Petetti L, Duignan P, Castinel A. Hookworm infection, anaemia and genetic variability of the New Zealand sea lion. Proc R Soc B Biol Sci. 2009;276:3523–9. doi:10.1098/rspb.2009.1001.

    Article  CAS  Google Scholar 

  29. Hoffman JI, Forcada J, Amos W. Exploring the mechanisms underlying a heterozygosity-fitness correlation for canine size in the Antarctic fur seal Arctocephalus gazella. J Hered. 2010;101(5):539–52. doi:10.1093/jhered/esq046.

    Article  CAS  PubMed  Google Scholar 

  30. Hoffman JI, Hanson N, Forcada J, Trathan PN, Amos W. Getting long in the tooth: a strong positive correlation between canine size and heterozygosity in Antarctic fur seals Arctocephalus gazella. J Hered. 2010;101(5):527–38. doi:10.1093/jhered/esq045.

    Article  CAS  PubMed  Google Scholar 

  31. David P, Pujol B, Viard F, Castella V, Goudet J. Reliable selfing rate estimates from imperfect population genetic data. Mol Ecol. 2007;16:2474–87. doi:10.1111/j.1365-294X.2007.03330.x.

    Article  CAS  PubMed  Google Scholar 

  32. Samollow P, Soulé M. A case of stress related heterozygote superiority in nature. Evolution. 1983;37(3):646–9.

    Article  Google Scholar 

  33. Hoffman A, Hercus M. Environmental stress as an evolutionary force. Bioscience. 2000;50(3):217–26. doi:[0217:ESAAEF]2.3.CO;2.

    Article  Google Scholar 

  34. Enders LS, Nunney L. Seasonal stress drives predictable changes in inbreeding depression in field-tested captive populations of Drosophila melanogaster. Proceedings of The Royal Society B: Biological Sciences. 2012. doi: 10.1098/rspb.2012.1018

  35. Aurioles D, Trillmich F. Zalophus wollebaeki. IUCN Red List of Threatened Species. Version 2011.1. 2008.

    Google Scholar 

  36. Alava JJ, Salazar S. Status and conservation of Otariids in Ecuador and the Galapagos Islands. In: Trites AW, Atkinson SK, De Master DP, Fritz LW, Gelatt TS, Rea LD, Wynne KM, editors. Sea lions of the world. Alaska Sea Grant College Program: Fairbanks; 2006. p. 495–520.

    Chapter  Google Scholar 

  37. Wolf JBW, Harrod C, Brunner S, Salazar S, Trillmich F, Tautz D. Tracing early stages of species differentiation: ecological, morphological and genetic divergence of Galapagos sea lion populations. BMC Evol Biol. 2008;8:14. doi:10.1186/1471-2148-8-150.

    Article  Google Scholar 

  38. Wolf JBW, Kauermann G, Trillmich F. Males in the shade: habitat use and sexual segregation in the Galápagos sea lion (Zalophus californianus wollebaeki). Behav Ecol Sociobiol. 2005;59:293–302. doi:10.1007/s00265-005-0042-7.

    Article  Google Scholar 

  39. Pörschmann U, Trillmich F, Mueller B, Wolf JBW. Male reproductive success and its behavioural correlates in a polygynous mammal, the Galápagos sea lion (Zalophus wollebaeki). Mol Ecol. 2010;19(12):2574–86. doi:10.1111/j.1365-294X.2010.04665.x.

    PubMed  Google Scholar 

  40. Schulte-Hostedde AI, Zinner B, Millar JS, Hickling GJ. Restitution of mass-size residuals: validating body condition indices. Ecology. 2005;86(1):155–63. doi:

    Article  Google Scholar 

  41. Hall AJ, McConnell BJ, Barker RJ. Factors affecting first-year survival in grey seals and their implications for life history strategy. J Anim Ecol. 2001;70:138–49. doi:10.1111/j.1365-2656.2001.00468.x.

    Article  Google Scholar 

  42. Hall AJ, McConnell BJ, Barker RJ. The effect of total immunoglobulin levels, mass and condition on the first-year survival of Grey Seal pups. Funct Ecol. 2002;16(4):462–74. doi:10.1046/j.1365-2435.2002.00649.x.

    Article  Google Scholar 

  43. Trillmich F, Limberger D. Drastic effects of El Nino on Galapagos pinnipeds. Oecologia. 1985;67(1):19–22. doi:10.1007/BF00378445.

    Article  Google Scholar 

  44. Mueller B, Pörschmann U, Wolf JBW, Trillmich F. Growth under uncertainty: the influence of marine variability on early development of Galapagos sea lions. Mar Mamm Sci. 2011;27(2):350–65. doi:10.1111/j.1748-7692.2010.00404.x.

    Article  Google Scholar 

  45. Cunningham AA, Daszak P, Rodriguez J. Pathogen pollution: defining a parasitological threat to biodiversity conservation. J Parasitol. 2003;89:S78–83.

    Google Scholar 

  46. Lenz TL, Mueller B, Trillmich F, Wolf JBW. Divergent allele advantage at MHC-DRB through direct and maternal genotypic effects and its consequences for allele pool composition and mating. Proc R Soc B. 2013;280(1762):20130714. doi:10.1098/rspb.2013.0714.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Råberg L, Stjernman M, Hasselquist D. Immune responsiveness in adult blue tits: heritability and effects of nutritional status during ontogeny. Oecologia. 2003;136(3):360–4. doi:10.1007/s00442-003-1287-3.

    Article  PubMed  Google Scholar 

  48. Brock PM, Hall AJ, Goodman SJ, Cruz M, Acevedo-Whitehouse K. Applying the tools of ecological immunology to conservation: a test case in the Galapagos sea lion. Anim Conserv. 2013;16(1):19–31. doi:10.1111/j.1469-1795.2012.00567.x.

    Article  Google Scholar 

  49. Brock PM, Hall AJ, Goodman SJ, Cruz M, Acevedo-Whitehouse K. Immune Activity, Body Condition and Human-Associated Environmental Impacts in a Wild Marine Mammal. PLoS ONE. 2013;8(6), e67132. doi:10.1371/journal.pone.0067132.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Wolf JBW, Tautz D, Trillmich F. Galápagos and Californian sea lions are separate species: genetic analysis of the genus Zalophus and its implications for conservation management. Front Zool. 2007;4:20. doi:10.1186/1742-9994-4-20.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Wang J. Triadic IBD coefficients and applications to estimating pairwise relatedness. Genet Res. 2007;89(3):135–53. doi:10.1017/S0016672307008798.

    Article  CAS  PubMed  Google Scholar 

  52. Aparicio JM, Ortego J, Cordero PJ. What should we weigh to estimate heterozygosity, alleles or loci? Mol Ecol. 2006;15:4659–65. doi:10.1111/j.1365-294X.2006.03111.x.

    Article  CAS  PubMed  Google Scholar 

  53. Pedersen AB, Babayan SA. Wild immunology. Mol Ecol. 2011;20:872–80. doi:10.1111/j.1365-294X.2010.04938.x.

    Article  CAS  PubMed  Google Scholar 

  54. McDade TW. Life history theory and the immune system: steps toward a human ecological immunology. Am J Phys Anthropol. 2003;46:100–25. doi:10.1002/ajpa.10398.

    Article  Google Scholar 

  55. Freitas A, Viale A, Sundblad A, Heusser C, Coutinho A. Normal serum immunoglobulins participate in the selection of peripheral B-cell repertoires. Proc Natl Acad Sci U S A. 1991;88:5640–4. doi:10.1073/pnas.88.13.5640.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Whiteman NK, Matson KD, Bollmer JL, Parker PG. Disease ecology in the Galápagos Hawk (Buteo galapagoensis): host genetic diversity, parasite load and natural antibodies. Proc R Soc Lond B Biol Sci. 2006;273(1588):797–804. doi:10.1098/rspb.2005.3396.

    Article  Google Scholar 

  57. Schmid-Hempel P. Ecological immunology. In: Evolutionary Parasitology. Oxford: Oxford University Press; 2011. p. 98–140.

    Google Scholar 

  58. Tizard IR. Veterinary Immunology (8th Edition). St Louis (MO), USA: Saunders Elsevier. 2009

  59. Keogh MJ, Maniscalco JM, Atkinson S. Steller sea lion (Eumetopias jubatus) pups undergo a decrease in circulating white blood cells and the ability of T cells to proliferate during early post-natal development. Vet Immunol Immunopathol. 2010;137(3–4):298–304. doi:10.1016/j.vetimm.2010.05.008.

    Article  CAS  PubMed  Google Scholar 

  60. Day MJ. Immune system development in the dog and cat. J Comp Pathol. 2007;137:S10–5. doi:10.1016/j.jcpa.2007.04.005.

    Article  CAS  PubMed  Google Scholar 

  61. Bossart GB, Reidarson TH, Dierauf LA, Duffield DA. Clinical pathology of marine mammals. In: Dierauf LA, Gulland MD, editors. Marine Mammal Medicine. 2nd ed. Boca Raton, FL, USA: CRC Press; 2001. p. 383–430.

    Google Scholar 

  62. Townsend AK, Clark AB, McGowan KJ, Miller AD, Buckles EL. Condition, innate immunity and disease mortality of inbred crows. Proc R Soc Lond B Biol Sci. 2010;277:2875–83. doi:10.1098/rspb.2010.0480.

    Article  Google Scholar 

  63. Coltman D, Pilkington J, Smith J, Pemberton J. Parasite-mediated selection against inbred Soay sheep in a free-living island population. Evolution. 1999;53(4):1259–67.

    Article  Google Scholar 

  64. McCarthy AJ, Shaw M, Jepson PD, Brasseur SMJM, Reijnders PJH, Goodman SJ. Variation in European harbour seal immune response genes and susceptibility to phocine distemper virus (PDV). Infect Genet Evol. 2011;11(7):1616–23. doi:10.1016/j.meegid.2011.06.002.

    Article  CAS  PubMed  Google Scholar 

  65. Amos W, Acevedo-Whitehouse K. A new test for genotype-fitness associations reveals a single microsatellite allele that strongly predicts the nature of tuberculosis infections in wild boar. Mol Ecol Resour. 2009;9(4):1102–11. doi:10.1111/j.1755-0998.2009.02560.x.

    Article  CAS  PubMed  Google Scholar 

  66. Salazar S, Denkinger J. Possible effects of climate change on the populations of Galapagos pinnipeds. Galapagos Research (Noticias Galapagos). 2010;67:45–9.

    Google Scholar 

  67. Wolf JBW, Tautz D, Caccone A, Steinfartz S. Development of new microsatellite loci and evaluation of loci from other pinniped species for the Galapagos sea lion (Zalophus californianus wollebaeki). Conserv Genet. 2006;7:461–5. doi:10.1007/s10592-005-9045-1.

    Article  Google Scholar 

  68. Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5. doi:10.1093/bioinformatics/btn129.

    Article  CAS  PubMed  Google Scholar 

  69. Development Core Team R. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2011.

    Google Scholar 

  70. Jeglinski JWE, Werner C, Robinson PW, Costa DP, Trillmich F. Age, body mass and environmental variation shape the foraging ontogeny of Galapagos sea lions. Mar Ecol Prog Ser. 2012;453:279–96.

    Article  Google Scholar 

  71. Crawley M J. The R book. Chichester, UK: Wiley. 2007

  72. Zuur A, Ieno E, Walker N, Saveliev A, Smith G. Mixed effects models and extensions in ecology with R. New York: NY: Springer; 2009.

    Book  Google Scholar 

  73. Wang J. COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients. Mol Ecol Resour. 2011;11(1):141–5. doi:10.1111/j.1755-0998.2010.02885.x.

    Article  PubMed  Google Scholar 

  74. Ritland K. Estimators for pairwise relatedness and individual inbreeding coefficients. Genet Res. 1996;67:175–85. doi:

  75. Lynch M, Ritland K. Estimation of pairwise relatedness with molecular markers. Genetics. 1999;152:1753–66.

    PubMed Central  CAS  PubMed  Google Scholar 

  76. Milligan BG. Maximum-likelihood estimation of relatedness. Genetics. 2003;163(3):1153–67.

    PubMed Central  PubMed  Google Scholar 

  77. Amos W, Wilmer JW, Fullard K, Burg TM, Croxall JP, Bloch D, et al. The influence of parental relatedness on reproductive success. Proc R Soc B Biol Sci. 2001;268(1480):2021–7. doi:10.1098/rspb.2001.1751.

    Article  CAS  Google Scholar 

  78. Alho JS, Välimäki K, Merilä J. Rhh: an R extension for estimating multilocus heterozygosity and heterozygosity-heterozygosity correlation. Mol Ecol Resour. 2010;10(4):720–2. doi:10.1111/j.1755-0998.2010.02830.x.

    Article  PubMed  Google Scholar 

  79. Dailey MD. Parasitic diseases. In: Dierauf LA, Gulland MD, editors. Marine Mammal Medicine. 2nd ed. Boca Raton, FL, USA: CRC Press; 2001. p. 383–430.

    Google Scholar 

  80. Wolf JBW, Trillmich F. Beyond habitat requirements: individual fine-scale site fidelity in a colony of the Galapagos sea lion (Zalophus wollebaeki) creates conditions for social structuring. Oecologia. 2007;152:553–67. doi:10.1007/s00442-007-0665-7.

    Article  PubMed  Google Scholar 

Download references


We would like to thank the editors and reviewers that have helped improve this article, and P Martinez, J Torres, A Velez, V Cedeño, W Tapia, J Murillo, E Espinoza, R Quinnell, J Jeglinski, J Hoffman, F Trillmich, T Meise, E Hippauf, P Brekke & J Wang for their help. All work was carried out under Galapagos National Park permits PC-18-09, N°046-2009-PNG, N°101-2010-PNG and N°032-2010-PNG, and DEFRA permits POAO/2008/925 and POAO/2010/136. PMB was funded by NERC, UFAW, SMRU and IOZ.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Patrick M. Brock or Karina Acevedo-Whitehouse.

Additional information

Competing interests

The authors declare no competing interests.

Authors’ contributions

KAW and PMB conceived and designed the study. PMB and MC collected the data in the field, and AJH and MC contributed reagents and materials. PMB analysed the data and wrote the paper, in discussion with KAW, SJG and AJH. All authors read and contributed to the final manuscript.

Additional files

Additional file 1: Table S1.

Tests for deviations from Hardy-Weinberg equilibrium. Supplementary text 1.1: Simulation methods. Supplementary text 1.2: Simulation results. Tables S2. & S3. Correlations with simulated inbreeding coefficients. Figure S1. Triadic likelihood estimator vs. homozygosity weighted by locus. Figure S2. Principal component analysis results. Supplementary text 1.3: Parasites methods. Supplementary text 1.4: Parasite results. (DOCX 95 kb)

Additional file 2:

Microsatellite data. (CSV 30 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Brock, P.M., Goodman, S.J., Hall, A.J. et al. Context-dependent associations between heterozygosity and immune variation in a wild carnivore. BMC Evol Biol 15, 242 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: