Comparing molecular measures for detecting inbreeding depression

Correlations between heterozygosity and components of fitness have been investigated in natural populations for over 20 years. Positive correlations between a trait of interest and heterozygosity (usually measured at allozyme loci) are generally recognized as evidence of inbreeding depression. More recently, molecular markers such as microsatellites have been employed for the same purpose. A typical study might use around five to ten markers. In this paper we use a panel of 71 microsatellite loci to: (1) Compare the efficacy of heterozygosity and a related microsatellite‐specific variable, mean d², in detecting inbreeding depression; (2) Examine the statistical power of heterozygosity to detect such associations. We performed our analyses in a wild population of red deer (Cervus elaphus) in which inbreeding depression in juvenile traits had previously been detected using a panel of nine markers. We conclude that heterozygosity‐based measures outperform mean d²‐based measures, but that power to detect heterozygosity‐fitness associations is nonetheless low when ten or fewer markers are typed.     

Inbreeding coefficient and heterozygosity–fitness correlations in unhatched and hatched song sparrow nestmates

Heterozygosity–fitness correlations use molecular measures of heterozygosity as proxy estimates of individual inbreeding coefficients (f) to examine relationships between inbreeding and fitness traits. Heterozygosity–fitness correlations partly depend on the assumption that individual heterozygosity and f are strongly and negatively correlated. Although theory predicts that this relationship will be strongest when mean f and variance in f are high, few studies of heterozygosity–fitness correlations include estimates of f based on pedigrees, which allow for more thorough examinations of the relationship between f, heterozygosity and fitness in nature. We examined relationships between pedigree‐based estimates of f, multilocus heterozygosity (MLH) and the probability of survival to hatch in song sparrow nestmates. f and MLH were weakly, but significantly negatively correlated. Inbreeding coefficient predicted the probability of survival to hatch. In contrast, MLH did not predict the probability of survival to hatch nor did it account for residual variation in survival to hatch after statistically controlling for the effects of f. These results are consistent with the expectation that heterozygosity–f correlations will be weak when mean and variance in f are low. Our results also provide empirical support for recent simulation studies, which show that variation in MLH among siblings with equal f can be large and may obscure MLH–fitness relationships.     

Cell‐mediated immunity and multi‐locus heterozygosity in bluethroat nestlings

Recent evidence suggests that marker‐based heterozygosity‐fitness correlations may be driven by only one or a few markers, indicating local heterozygosity effects caused by linkage disequilibrium with functional genes. In this study, we investigated the relationship between microsatellite heterozygosity and a measure of cell‐mediated immunity (phytohaemagglutinin; PHA) in bluethroat (Luscinia s. svecica) nestlings using a full‐sibling design. We found significant positive associations between PHA response and two different indices of microsatellite heterozygosity, i.e. multi‐locus heterozygosity and mean d². However, model comparisons disclosed that both associations were more likely caused by local effects rather than general effects and that the two local effects appeared to be realized through two different genetic mechanisms. Our results indicate that both the random assortment of parental chromosomes during meiosis as well as inbreeding can drive heterozygosity‐fitness correlations.     

Genetic diversity-fitness correlation revealed by microsatellite analyses in European alpine marmots (Marmota marmota)

The relationship between individual genetic diversity and fitness-related traits are poorly understood in the wild. The availability of highly polymorphic molecular markers, such as microsatellites, has made research on this subject more feasible. We used three microsatellite-based measures of genetic diversity, individual heterozygosity H, mean d ² and mean d ² _outbreeding to test for a relationship between individual genetic diversity and important fitness trait, juvenile survival, in a population of alpine marmots (Marmota marmota), after controlling for the effects of ecological, social and physiological parameters that potentially influence juvenile survival in marmots. Analyses were conducted on 158 juveniles, and revealed a positive association between juvenile survival and genetic diversity measured by mean H. No association was found with mean d ² and with mean d ² _outbreeding. This suggests a fitness disadvantage to less heterozygous juveniles. The genetic diversity-fitness correlation (GDFC) was somewhat stronger during years with poor environmental conditions (i.e. wet summers). The stressful environmental conditions of this high mountain population might enhance inbreeding depression and make this association between genetic diversity and fitness detectable. Moreover the mating system, allowing extra pair copulation by occasional immigrants, as well as close inbreeding, favours a wide range of individual genetic diversity (mean H ranges from 0.125 to 1), which also may have facilitated the detection of the GDFC. The results further suggest that the observed GDFC is likely to be explained by the “local effect” hypothesis rather than by the “general effect” hypothesis.     

MICROSATELLITE LOCI REVEAL SEX-DEPENDENT RESPONSES TO INBREEDING AND OUTBREEDING IN RED DEER CALVES

Mean d² is a recently devised microsatellite-based measure that is hypothesised to allow the detection of inbreeding depression and heterosis in free-living populations. Two studies that have investigated the measure have both demonstrated an association between mean d² and traits related to fitness. Here we present an association between mean d¹ and an important component of fitness, first-year overwinter survival, in a population of red deer on the Isle of Rum, Scotland. The association between survival and mean d² differed between males and females. As predicted, outbred female calves (high mean d²) survived better than those that were inbred (low mean d²). However, the association was in the opposite direction in male calves. We suggest that this difference is due to different early growth strategies between the sexes. The association between mean d² and survival was not significantly influenced by any single locus. Decomposition of mean d² into a recent inbreeding component and an outbreeding component showed that it was the degree of outbreeding that influenced survival in males and both the degree of outbreeding and recent inbreeding that influenced survival in females. Our analyses suggest that mean d² is an easy-to-calculate measure of inbreeding and degree of outbreeding that can reveal interesting interactions between genetics and ecology.     

Recent admixture generates heterozygosity–fitness correlations during the range expansion of an invading species

Admixture, the mixing of historically isolated gene pools, can have immediate consequences for the genetic architecture of fitness traits. Admixture may be especially important for newly colonized populations, such as during range expansion and species invasions, by generating heterozygosity that can boost fitness through heterosis. Despite widespread evidence for admixture during species invasions, few studies have examined the demographic history leading to admixture, how admixture affects the heterozygosity and fitness of invasive genotypes, and whether such fitness effects are maintained through time. We address these questions using the invasive plant Silene vulgaris, which shows evidence of admixture in both its native Europe and in North America where it has invaded. Using multilocus genotype data in conjunction with approximate Bayesian computation analysis of demographic history, we showed that admixture during the invasion of North America was independent from and much younger than admixture in the native range of Europe. We tested for fitness consequences of admixture in each range and detected a significant positive heterozygosity–fitness correlation (HFC) in North America; in contrast, no HFC was present in Europe. The lack of HFC in Europe may reflect the longer time since admixture in the native range, dissipating associations between heterozygosity at markers and fitness loci. Our results support a key short‐term role for admixture during the early stages of invasion by generating HFCs that carry populations past the threat of extinction from inbreeding and demographic stochasticity.     

Adult survival and microsatellite diversity in possums: effects of major histocompatibility complex-linked microsatellite diversity but not multilocus inbreeding estimators

Adult survival is perhaps the fitness parameter most important to population growth in long-lived species. Intrinsic and extrinsic covariates of survival are therefore likely to be important drivers of population dynamics. We used long-term mark-recapture data to identify genetic, individual and environmental covariates of local survival in a natural population of mountain brushtail possums (Trichosurus cunninghami). Rainfall and intra-individual diversity at microsatellite DNA markers were associated with increased local survival of adults and juveniles. We contrasted the performance of several microsatellite heterozygosity measures, including internal relatedness (IR), homozygosity by loci (HL) and the mean multilocus estimate of the squared difference in microsatellite allele sizes within an individual (mean d ²). However, the strongest effect on survival was not associated with multilocus microsatellite diversity (which would indicate a genome-wide inbreeding effect), but a subset of two loci. This included a major histocompatibility complex (MHC)-linked marker and a putatively neutral microsatellite locus. For both loci, diversity measures incorporating allele size information had stronger associations with survival than measures based on heterozygosity, whether or not allele frequency information was included (such as IR). Increased survival was apparent among heterozygotes at the MHC-linked locus, but the benefits of heterozygosity to survival were reduced in heterozygotes with larger differences in allele size. The effect of heterozygosity on fitness-related traits was supported by data on endoparasites in a subset of the individuals studied in this population. There was no apparent density dependence in survival, nor an effect of sex, age or immigrant status. Our findings suggest that in the apparent absence of inbreeding, variation at specific loci can generate strong associations between fitness and diversity at linked markers.     

Heterozygosity–fitness correlations in a wild mammal population: accounting for parental and environmental effects

HFCs (heterozygosity–fitness correlations) measure the direct relationship between an individual's genetic diversity and fitness. The effects of parental heterozygosity and the environment on HFCs are currently under‐researched. We investigated these in a high‐density U.K. population of European badgers (Meles meles), using a multimodel capture–mark–recapture framework and 35 microsatellite loci. We detected interannual variation in first‐year, but not adult, survival probability. Adult females had higher annual survival probabilities than adult males. Cubs with more heterozygous fathers had higher first‐year survival, but only in wetter summers; there was no relationship with individual or maternal heterozygosity. Moist soil conditions enhance badger food supply (earthworms), improving survival. In dryer years, higher indiscriminate mortality rates appear to mask differential heterozygosity‐related survival effects. This paternal interaction was significant in the most supported model; however, the model‐averaged estimate had a relative importance of 0.50 and overlapped zero slightly. First‐year survival probabilities were not correlated with the inbreeding coefficient (f); however, small sample sizes limited the power to detect inbreeding depression. Correlations between individual heterozygosity and inbreeding were weak, in line with published meta‐analyses showing that HFCs tend to be weak. We found support for general rather than local heterozygosity effects on first‐year survival probability, and g2 indicated that our markers had power to detect inbreeding. We emphasize the importance of assessing how environmental stressors can influence the magnitude and direction of HFCs and of considering how parental genetic diversity can affect fitness‐related traits, which could play an important role in the evolution of mate choice.     

Evidence of opposing fitness effects of parental heterozygosity and relatedness in a critically endangered marine turtle?

How individual genetic variability relates to fitness is important in understanding evolution and the processes affecting populations of conservation concern. Heterozygosity–fitness correlations (HFCs) have been widely used to study this link in wild populations, where key parameters that affect both variability and fitness, such as inbreeding, can be difficult to measure. We used estimates of parental heterozygosity and genetic similarity (‘relatedness’) derived from 32 microsatellite markers to explore the relationship between genetic variability and fitness in a population of the critically endangered hawksbill turtle, Eretmochelys imbricata. We found no effect of maternal MLH (multilocus heterozygosity) on clutch size or egg success rate, and no single‐locus effects. However, we found effects of paternal MLH and parental relatedness on egg success rate that interacted in a way that may result in both positive and negative effects of genetic variability. Multicollinearity in these tests was within safe limits, and null simulations suggested that the effect was not an artefact of using paternal genotypes reconstructed from large samples of offspring. Our results could imply a tension between inbreeding and outbreeding depression in this system, which is biologically feasible in turtles: female‐biased natal philopatry may elevate inbreeding risk and local adaptation, and both processes may be disrupted by male‐biased dispersal. Although this conclusion should be treated with caution due to a lack of significant identity disequilibrium, our study shows the importance of considering both positive and negative effects when assessing how variation in genetic variability affects fitness in wild systems.     

Assessment of identity disequilibrium and its relation to empirical heterozygosity fitness correlations: a meta‐analysis

Heterozygosity fitness correlations (HFCs) have frequently been used to detect inbreeding depression, under the assumption that genome‐wide heterozygosity is a good proxy for inbreeding. However, meta‐analyses of the association between fitness measures and individual heterozygosity have shown that often either no correlations are observed or the effect sizes are small. One of the reasons for this may be the absence of variance in inbreeding, a requisite for generating general‐effect HFCs. Recent work has highlighted identity disequilibrium (ID) as a measure that may capture variance in the level of inbreeding within a population; however, no thorough assessment of ID in natural populations has been conducted. In this meta‐analysis, we assess the magnitude of ID (as measured by the g₂ statistic) from 50 previously published HFC studies and its relationship to the observed effect sizes of those studies. We then assess how much power the studies had to detect general‐effect HFCs, and the number of markers that would have been needed to generate a high expected correlation (r² = 0.9) between observed heterozygosity and inbreeding. Across the majority of studies, g₂ values were not significantly different than zero. Despite this, we found that the magnitude of g₂ was associated with the average effect sizes observed in a population, even when point estimates were nonsignificant. These low values of g₂ translated into low expected correlations between heterozygosity and inbreeding and suggest that many more markers than typically used are needed to robustly detect HFCs.     

Using heterozygosity–fitness correlations to study inbreeding depression in an isolated population of white‐tailed deer founded by few individuals

A heterozygosity–fitness correlations (HFCs) may reflect inbreeding depression, but the extent to which they do so is debated. HFCs are particularly likely to occur after demographic disturbances such as population bottleneck or admixture. We here study HFC in an introduced and isolated ungulate population of white‐tailed deer Odocoileus virginianus in Finland founded in 1934 by four individuals. A total of 422 ≥ 1‐year‐old white‐tailed deer were collected in the 2012 hunting season in southern Finland and genotyped for 14 microsatellite loci. We find significant identity disequilibrium as estimated by g₂. Heterozygosity was positively associated with size‐ and age‐corrected body mass, but not with jaw size or (in males) antler score. Because of the relatively high identity disequilibrium, heterozygosity of the marker panel explained 51% of variation in inbreeding. Inbreeding explained approximately 4% of the variation in body mass and is thus a minor, although significant source of variation in body mass in this population. The study of HFC is attractive for game‐ and conservation‐oriented wildlife management because it presents an affordable and readily used approach for genetic monitoring that allowing identification of fitness costs associated with genetic substructuring in what may seem like a homogeneous population.     

Heterozygosity at a single locus explains a large proportion of variation in two fitness‐related traits in great tits: a general or a local effect?

In natural populations, mating between relatives can have important fitness consequences due to the negative effects of reduced heterozygosity. Parental level of inbreeding or heterozygosity has been also found to influence the performance of offspring, via direct and indirect parental effects that are independent of the progeny own level of genetic diversity. In this study, we first analysed the effects of parental heterozygosity and relatedness (i.e. an estimate of offspring genetic diversity) on four traits related to offspring viability in great tits (Parus major) using 15 microsatellite markers. Second, we tested whether significant heterozygosity–fitness correlations (HFCs) were due to ‘local’ (i.e. linkage to genes influencing fitness) and/or ‘general’ (genome‐wide heterozygosity) effects. We found a significant negative relationship between parental genetic relatedness and hatching success, and maternal heterozygosity was positively associated with offspring body size. The characteristics of the studied populations (recent admixture, polygynous matings) together with the fact that we found evidence for identity disequilibrium across our set of neutral markers suggest that HFCs may have resulted from genome‐wide inbreeding depression. However, one locus (Ase18) had disproportionately large effects on the observed HFCs: heterozygosity at this locus had significant positive effects on hatching success and offspring size. It suggests that this marker may lie near to a functional locus under selection (i.e. a local effect) or, alternatively, heterozygosity at this locus might be correlated to heterozygosity across the genome due to the extensive ID found in our populations (i.e. a general effect). Collectively, our results lend support to both the general and local effect hypotheses and reinforce the view that HFCs lie on a continuum from inbreeding depression to those strictly due to linkage between marker loci and genes under selection.     

Accounting for cryptic population substructure enhances detection of inbreeding depression with genomic inbreeding coefficients: an example from a critically endangered marsupial

Characterizing inbreeding depression in wildlife populations can be critical to their conservation. Coefficients of individual inbreeding can be estimated from genome‐wide marker data. The degree to which sensitivity of inbreeding coefficients to population genetic substructure alters estimates of inbreeding depression in wild populations is not well understood. Using generalized linear models, we tested the power of two frequently used inbreeding coefficients that are calculated from genome‐wide SNP markers, F_H and F^^III, to predict four fitness traits estimated over two decades in an isolated population of the critically endangered Leadbeater's possum. F_H estimates inbreeding as excess observed homozygotes relative to equilibrium expectations, whereas F^^III quantifies allelic similarity between the gametes that formed an individual, and upweights rare homozygotes. We estimated F_H and F^^III from 1,575 genome‐wide SNP loci in individuals with fitness trait data (N = 179–237 per trait), and computed revised coefficients, F_H^{by group} and F^^{IIIby group}, adjusted for population genetic substructure by calculating them separately within two different genetic groups of individuals identified in the population. Using F_H or F^^III in the models, inbreeding depression was detected for survival to sexual maturity, longevity and whether individuals bred during their lifetime. F^^{IIIby group} (but not F_H^{by group}) additionally revealed significant inbreeding depression for lifetime reproductive output (total offspring assigned to each individual). Estimates of numbers of lethal equivalents indicated substantial inbreeding load, but differing between inbreeding estimators. Inbreeding depression, declining population size, and low and declining genetic diversity suggest that genetic rescue may assist in preventing extinction of this unique Leadbeater's possum population.     

Genome‐wide regulatory deterioration impedes adaptive responses to stress in inbred populations of Drosophila melanogaster*

Inbreeding depression is often intensified under environmental stress (i.e., inbreeding–stress interaction). Although the fitness consequences of this phenomenon are well‐described, underlying mechanisms such as an increased expression of deleterious alleles under stress, or a lower capacity for adaptive responses to stress with inbreeding, have rarely been investigated. We investigated a fitness component (egg‐to‐adult viability) and gene‐expression patterns using RNA‐seq analyses in noninbred control lines and in inbred lines of Drosophila melanogaster exposed to benign temperature or heat stress. We find little support for an increase in the cumulative expression of deleterious alleles under stress. Instead, inbred individuals had a reduced ability to induce an adaptive gene regulatory stress response compared to controls. The decrease in egg‐to‐adult viability due to stress was most pronounced in the lines with the largest deviation in the adaptive stress response (R² = 0.48). Thus, we find strong evidence for a lower capacity of inbred individuals to respond by gene regulation to stress and that this is the main driver of inbreeding‐stress interactions. In comparison, the altered gene expression due to inbreeding at benign temperature showed no correlation with fitness and was pronounced in genomic regions experiencing the highest increase in homozygosity.     

META-ANALYSES OF THE ASSOCIATION BETWEEN MULTILOCUS HETEROZYGOSITY AND FITNESS

Meta-analyses of published correlation coefficients between multilocus heterozygosity (MLH) and two fitness surrogates, growth rate and fluctuating asymmetry, suggested that the strength of these correlations are generally weak. A variety of plants and animals was included in the meta-analyses. A statistically homogeneous group of MLH–growth rate correlation coefficients that included both plants and animals yielded a common correlation of r_z = 0.133. A common correlation of r_z = –0.170 was estimated for correlations between MLH and fluctuating asymmetry in three species of salmonid fishes. These results suggest that selection, including overdominance, has at most a weak effect at allozyme loci and cast some doubt on the widely held notion that heterozygosity and individual fitness are strongly correlated.     

On the genetic parameter determining the efficiency of purging: an estimate for Drosophila egg‐to‐pupae viability

The consequences of inbreeding on fitness can be crucial in evolutionary and conservation grounds and depend upon the efficiency of purging against deleterious recessive alleles. Recently, analytical expressions have been derived to predict the evolution of mean fitness, taking into account both inbreeding and purging, which depend on an ‘effective purging coefficient (d_e)’. Here, we explore the validity of that predictive approach and assay the strength of purging by estimating d_e for egg‐to‐pupae viability (EPV) after a drastic reduction in population size in a recently captured base population of Drosophila melanogaster. For this purpose, we first obtained estimates of the inbreeding depression rate (δ) for EPV in the base population, and we found that about 40% was due to segregating recessive lethals. Then, two sets of lines were founded from this base population and were maintained with different effective size throughout the rest of the experiment (N = 6; N = 12), their mean EPV being assayed at different generations. Due to purging, the reductions in mean EPV experienced by these lines were considerably smaller than the corresponding neutral predictions. For the 60% of δ attributable to nonlethal deleterious alleles, our results suggest an effective purging coefficient d_e > 0.02. Similarly, we obtain that d_e > 0.09 is required to roughly account for purging against the pooled inbreeding depression from lethal and nonlethal deleterious alleles. This implies that purging should be efficient for population sizes of the order of a few tens and larger, but might be inefficient against nonlethal deleterious alleles in smaller populations.     

Genomewide ancestry and divergence patterns from low‐coverage sequencing data reveal a complex history of admixture in wild baboons

Naturally occurring admixture has now been documented in every major primate lineage, suggesting its key role in primate evolutionary history. Active primate hybrid zones can provide valuable insight into this process. Here, we investigate the history of admixture in one of the best‐studied natural primate hybrid zones, between yellow baboons (Papio cynocephalus) and anubis baboons (Papio anubis) in the Amboseli ecosystem of Kenya. We generated a new genome assembly for yellow baboon and low‐coverage genomewide resequencing data from yellow baboons, anubis baboons and known hybrids (n = 44). Using a novel composite likelihood method for estimating local ancestry from low‐coverage data, we found high levels of genetic diversity and genetic differentiation between the parent taxa, and excellent agreement between genome‐scale ancestry estimates and a priori pedigree, life history and morphology‐based estimates (r² = 0.899). However, even putatively unadmixed Amboseli yellow individuals carried a substantial proportion of anubis ancestry, presumably due to historical admixture. Further, the distribution of shared vs. fixed differences between a putatively unadmixed Amboseli yellow baboon and an unadmixed anubis baboon, both sequenced at high coverage, is inconsistent with simple isolation–migration or equilibrium migration models. Our findings suggest a complex process of intermittent contact that has occurred multiple times in baboon evolutionary history, despite no obvious fitness costs to hybrids or major geographic or behavioural barriers. In combination with the extensive phenotypic data available for baboon hybrids, our results provide valuable context for understanding the history of admixture in primates, including in our own lineage.     

Microsatellite measures of inbreeding: a meta-analysis

Abstract Meta-analyses of published and unpublished correlations between phenotypic variation and two measures of genetic variation at microsatellite loci, multilocus heterozygosity (MLH) and mean d², revealed that the strength of these associations are generally weak (mean r < 0.10). Effects on life-history trait variation were significantly greater than zero for both measures over all reported effect sizes ( r = 0. 0856 and 0.0479 for MLH and mean d ², respectively), whereas effects on morphometric traits were not ( r = 0.0052 and r = 0.0038), which is consistent with the prediction that life-history traits exhibit greater inbreeding depression than morphometric traits. Effect sizes reported using mean d ² were smaller and more variable than those reported using MLH, suggesting that MLH may be a better metric for capturing inbreeding depression most of the time. However, analyses of paired effect sizes reported using both measures from the same data did not differ significantly. Several lines of evidence suggest that published effects sizes are upwardly biased. First, effect sizes from published studies were significantly higher than those reported in unpublished studies. Second, fail-safe numbers for reported effect sizes were generally quite low, with the exception of correlations between MLH and life-history traits. Finally, the slope of the regression of effect size on sample size was negative for most sets of traits. Taken together, these results suggest that studies designed to detect inbreeding depression on a life-history trait using microsatellites will need to sample in excess of 600 individuals to detect an average effect size ( r = 0.10) with reasonable statistical power (0.80). Very few published studies have used samples sizes approaching this value.     

Evaluating the role of inbreeding depression in heterozygosity‐fitness correlations: how useful are tests for identity disequilibrium?

Heterozygosity‐fitness correlations (HFCs) have been observed for several decades, but their causes are often elusive. Tests for identity disequilibrium (ID, correlated heterozygosity between loci) are commonly used to determine if inbreeding depression is a possible cause of HFCs. We used computer simulations to determine how often ID is detected when HFCs are caused by inbreeding depression. We also used ID in conjunction with HFCs to estimate the proportion of variation (r²) in fitness explained by the individual inbreeding coefficient (F). ID was not detected in a large proportion of populations with statistically significant HFCs (sample size = 120 individuals) unless the variance of F was high (σ²(F) ≥ 0.005) or many loci were used (100 microsatellites or 1000 SNPs). For example, with 25 microsatellites, ID was not detected in 49% of populations when HFCs were caused by six lethal equivalents and σ²(F) was typical of vertebrate populations (σ²(F) ≈ 0.002). Estimates of r² between survival and F based on ID and HFCs were imprecise unless ID was strong and highly statistically significant (P ≈ 0.01). These results suggest that failing to detect ID in HFC studies should not be taken as evidence that inbreeding depression is absent. The number of markers necessary to simultaneously detect HFC and ID depends strongly on σ²(F). Thus the mating system and demography of populations, which influence σ²(F), should be considered when designing HFC studies. ID should be used in conjunction with HFCs to estimate the correlation between fitness and F, because HFCs alone reveal little about the strength of inbreeding depression.     

Heterozygosity at neutral and immune loci is not associated with neonatal mortality due to microbial infection in Antarctic fur seals

Numerous studies have reported correlations between the heterozygosity of genetic markers and fitness. These heterozygosity–fitness correlations (HFCs) play a central role in evolutionary and conservation biology, yet their mechanistic basis remains open to debate. For example, fitness associations have been widely reported at both neutral and functional loci, yet few studies have directly compared the two, making it difficult to gauge the relative contributions of genome‐wide inbreeding and specific functional genes to fitness. Here, we compared the effects of neutral and immune gene heterozygosity on death from bacterial infection in Antarctic fur seal (Arctocephalus gazella) pups. We specifically developed a panel of 13 microsatellites from expressed immune genes and genotyped these together with 48 neutral loci in 234 individuals, comprising 39 pups that were classified at necropsy as having most likely died of bacterial infection together with a five times larger matched sample of healthy surviving pups. Identity disequilibrium quantified from the neutral markers was positive and significant, indicative of variance in inbreeding within the study population. However, multilocus heterozygosity did not differ significantly between healthy and infected pups at either class of marker, and little evidence was found for fitness associations at individual loci. These results support a previous study of Antarctic fur seals that found no effects of heterozygosity at nine neutral microsatellites on neonatal survival and thereby help to refine our understanding of how HFCs vary across the life cycle. Given that nonsignificant HFCs are underreported in the literature, we also hope that our study will contribute toward a more balanced understanding of the wider importance of this phenomenon.