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.     

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.     

Population demography and heterozygosity–fitness correlations in natural guppy populations: An examination using sexually selected fitness traits

Heterozygosity–fitness correlations (HFCs) have been examined in a wide diversity of contexts, and the results are often used to infer the role of inbreeding in natural populations. Although population demography, reflected in population‐level genetic parameters such as allelic diversity or identity disequilibrium, is expected to play a role in the emergence and detectability of HFCs, direct comparisons of variation in HFCs across many populations of the same species, with different genetic histories, are rare. Here, we examined the relationship between individual microsatellite heterozygosity and a range of sexually selected traits in 660 male guppies from 22 natural populations in Trinidad. Similar to previous studies, observed HFCs were weak overall. However, variation in HFCs among populations was high for some traits (although these variances were not statistically different from zero). Population‐level genetic parameters, specifically genetic diversity levels (number of alleles, observed/expected heterozygosity) and measures of identity disequilibrium (g2 and heterozygosity–heterozygosity correlations), were not associated with variation in population‐level HFCs. This latter result indicates that these metrics do not necessarily provide a reliable predictor of HFC effect sizes across populations. Importantly, diversity and identity disequilibrium statistics were not correlated, providing empirical evidence that these metrics capture different essential characteristics of populations. A complex genetic architecture likely underpins multiple fitness traits, including those associated with male fitness, which may have reduced our ability to detect HFCs in guppy populations. Further advances in this field would benefit from additional research to determine the demographic contexts in which HFCs are most likely to occur.     

HETEROZYGOSITY-FITNESS CORRELATIONS REVEALED BY NEUTRAL AND CANDIDATE GENE MARKERS IN ROE DEER FROM A LONG-TERM STUDY

Heterozygosity-fitness correlations (HFCs) are increasingly reported but the underlying mechanisms causing HFCs are generally poorly understood. Here, we test for HFCs in roe deer ( Capreolus capreolus ) using 22 neutral microsatellites widely distributed in the genome and four microsatellites in genes that are potentially under selection. Juvenile survival was used as a proxy for individual fitness in a population that has been intensively studied for 30 years in northeastern France. For 222 juveniles, we computed two measures of genetic diversity: individual heterozygosity ( H ), and mean d ² (relatedness of parental genomes). We found a relationship between genetic diversity and fitness both for the 22 neutral markers and two candidate genes: IGF1 (Insulin-like Growth Factor I) and NRAMP (natural resistance-associated macrophage protein). Statistical evidence and the size of genetic effects on juvenile survival were comparable to those reported for early development and cohort variation, suggesting a substantial influence of genetic components on fitness in this roe deer population. For the 22 neutral microsatellites, a correlation with fitness was revealed for mean d ², but not for H , suggesting a possible outbreeding advantage. This heterosis effect could have been favored by introduction of genetically distant (Hungarian) roe deer to the population in recent times and, possibly, by the structuring of the population into distinct clans. The locus-specific correlations with fitness may be driven by growth rate advantages and resistance to diseases known to exist in the studied population. Our analyses of neutral and candidate gene markers both suggest that the observed HFCs are likely mainly due to linkage with dominant or overdominant loci that affect fitness ("local" effect) rather than to a genome-wide relationship with homozygosity due to inbreeding ("general" effect).     

HETEROZYGOSITY‐FITNESS CORRELATIONS: A TIME FOR REAPPRAISAL

Owing to the remarkable progress of molecular techniques, heterozygosity‐fitness correlations (HFCs) have become a popular tool to study the impact of inbreeding in natural populations. However, their underlying mechanisms are often hotly debated. Here we argue that these “debates” rely on verbal arguments with no basis in existing theory and inappropriate statistical testing, and that it is time to reconcile HFC with its historical and theoretical fundaments. We show that available data are quantitatively and qualitatively consistent with inbreeding‐based theory. HFC can be used to estimate the impact of inbreeding in populations, although such estimates are bound to be imprecise, especially when inbreeding is weak. Contrary to common belief, linkage disequilibrium is not an alternative to inbreeding, but rather comes with some forms of inbreeding, and is not restricted to closely linked loci. Finally, the contribution of local chromosomal effects to HFC, while predicted by inbreeding theory, is expected to be small, and has rarely if ever proven statistically significant using adequate tests. We provide guidelines to safely interpret and quantify HFCs, and present how HFCs can be used to quantify inbreeding load and unravel the structure of natural populations.     

No correlation between multi‐locus heterozygosity and fitness in the common buzzard despite heterozygote advantage for plumage colour

Correlations between heterozygosity and fitness are frequently found but rarely well understood. Fitness can be affected by single loci of large effect which correlate with neutral markers via linkage disequilibrium, or as a result of variation in genome‐wide heterozygosity following inbreeding. We explored these alternatives in the common buzzard, a raptor species in which three colour morphs differ in their lifetime reproductive success. Using 18 polymorphic microsatellite loci, we evaluated potential genetic differences among the morphs which may lead to subpopulation structuring and tested for correlations between three fitness‐related traits and heterozygosity, both genome wide and at each locus separately. Despite their assortative mating pattern, the buzzard morphs were found to be genetically undifferentiated. Multilocus heterozygosity was only found to be correlated with a single fitness‐related trait, infection with the blood parasite, Leucocytozoon buteonis, and this was via interactions with vole abundance and age. One locus also showed a significant relationship with blood parasite infection and ectoparasite infestation. The vicinity of this locus contains two genes, one of which is potentially implicated in the immune system of birds. We conclude that genome‐wide heterozygosity is unlikely to be a major determinant of parasite burden and body condition in the polymorphic common buzzard.     

Growth rate correlates to individual heterozygosity in the European eel, Anguilla anguilla L

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme over-dominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.     

Negative effects of individual heterozygosity on reproductive success in a wild bird population

The evolutionary consequences of individual genetic diversity are frequently studied by assessing heterozygosity–fitness correlations (HFCs). The prevalence of positive and negative HFCs and the predominance of general versus local effects in wild populations are far from understood, partly because comprehensive studies testing for both inbreeding and outbreeding depression are lacking. We studied a genetically diverse population of blue tits in southern Germany using a genome‐wide set of 87 microsatellites to investigate the relationship between proxies of reproductive success and measures of multilocus and single‐locus individual heterozygosity (MLH and SLH). We used complimentary measures of MLH and partitioned markers into functional categories according to their position in the blue tit genome. HFCs based on MLH were consistently negative for functional loci, whereas correlations were rather inconsistent for loci found in nonfunctional areas of the genome. Clutch size was the only reproductive variable showing a general effect. We found evidence for local effects for three measures of reproductive success: arrival date at the breeding site, the probability of breeding at the study site and male reproductive success. For these, we observed consistent, and relatively strong, negative effects at one functional locus. Remarkably, this marker had a similar effect in another blue tit population from Austria (~400 km to the east). We suggest that a genetic local effect on timing of arrival might be responsible for most negative HFCs detected, with carry‐over effects on other reproductive traits. This effect could reflect individual differences in the distance between overwintering areas and breeding sites.     

Heterozygosity–fitness correlations and associative overdominance: new detection method and proof of principle in the Iberian wild boar

Heterozygosity-fitness correlations (HFC) may result from a genome-wide process — inbreeding — or local effects within the genome. The majority of empirical studies reporting HFCs have attributed correlations to inbreeding depression. However, HFCs are unlikely to be caused by inbreeding depression because heterozygosity measured at a small number of neutral markers is unlikely to accurately capture a genome-wide pattern. Testing the strengths of localized effects caused by associative overdominance has proven challenging. In their current paper, Amos and Acevedo-Whitehouse present a novel test for local HFCs. Using stochastic simulations, they determine the conditions under which single-locus HFCs arise, before testing the strength of the correlation between the neutral marker and a linked gene under selection in their simulations. They used insights gained from simulation to statistically investigate the likely cause of correlations between heterozygosity and disease status using data on bovine tuberculosis infections in a wild boar population. They discover that a single microsatellite marker is an excellent predictor of tuberculosis progression in infected individuals. The results are relevant for wild boar management but, more generally, they demonstrate how single-locus HFCs could be used to identify coding loci under selection in free-living populations.     

Heterozygosity–fitness correlations and their relevance to studies on inbreeding depression in threatened species

The majority of reported multilocus heterozygosity–fitness correlations (HFCs) are from large, outbred populations, and their relevance to studies on inbreeding depression in threatened populations is often stressed. The results of such HFC studies conducted on outbred populations may be of limited application to threatened population management, however, as bottlenecked populations exhibit increased incidence of inbreeding, increased linkage disequilibrium, reduced genetic diversity and possible effects of historical inbreeding such as purging. These differences may affect both our ability to detect inbreeding depression in threatened species, and our interpretation of the underlying mechanisms for observed heterozygosity–fitness relationships. The study of HFCs in outbred populations is of interest in itself, but the results may not translate directly to threatened populations that have undergone severe bottlenecks.     

DOES LINKAGE DISEQUILIBRIUM GENERATE HETEROZYGOSITY‐FITNESS CORRELATIONS IN GREAT REED WARBLERS?

Heterozygosity-fitness correlations (HFCs) at noncoding genetic markers are commonly assumed to reflect fitness effects of heterozygosity at genomewide distributed genes in partially inbred populations. However, in populations with much linkage disequilibrium (LD), HFCs may arise also as a consequence of selection on fitness loci in the local chromosomal vicinity of the markers. Recent data suggest that relatively high levels of LD may prevail in many ecological situations. Consequently, LD may be an important factor, together with partial inbreeding, in causing HFCs in natural populations. In the present study, we evaluate whether LD can generate HFCs in a small and newly founded population of great reed warblers (Acrocephalus arundinaceus). For this purpose dyads of full siblings of which only one individual survived to adult age (i.e., returned to breed at the study area) were scored at 19 microsatellite loci, and at a gene region of hypothesized importance for survival, the major histocompatibility complex (MHC). By examining siblings, we controlled for variation in the inbreeding coefficient and thus excluded genome-wide fitness effects in our analyses. We found that recruited individuals had significantly higher multilocus heterozygosity (MLH), and mean d2 (a microsatellite-specific variable), than their nonrecruited siblings. There was a tendency for the survivors to have a more diverse MHC than the nonsurvivors. Single-locus analyses showed that the strength of the genotype-survival association was especially pronounced at four microsatellite loci. By using genotype data from the entire breeding population, we detected significant LD between five of 162 pairs of microsatellite loci after accounting for multiple tests. Our present finding of a significant within-family multilocus heterozygosity-survival association in a nonequilibrium population supports the view that LD generates HFCs in natural populations.     

On the correlation between heterozygosity and fitness in natural populations

Three primary hypotheses currently prevail for correlations between heterozygosity at a set of molecular markers and fitness in natural populations. First, multilocus heterozygosity-fitness correlations might result from selection acting directly on the scored loci, such as at particular allozyme loci. Second, significant levels of linkage disequilibrium, as in recently bottlenecked-and-expanded populations, might cause associations between the markers and fitness loci in the local chromosomal vicinity. Third, in partially inbred populations, heterozygosity at the markers might reflect variation in the inbreeding coefficient and might associate with fitness as a result of effects of homozygosity at genome-wide distributed loci. Despite years of research, the relative importance of these hypotheses remains unclear. The screening of heterozygosity at polymorphic DNA markers offers an opportunity to resolve this issue, and relevant empirical studies have now emerged. We provide an account of the recent progress on the subject, and give suggestions on how to distinguish between the three hypotheses in future studies.     

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.     

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.     

Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee

Inbreeding (the mating between closely related individuals) often has detrimental effects that are associated with loss of heterozygosity at overdominant loci, and the expression of deleterious recessive alleles. However, determining which loci are detrimental when homozygous, and the extent of their phenotypic effects, remains poorly understood. Here, we utilize a unique inbred population of clonal (thelytokous) honey bees, Apis mellifera capensis, to determine which loci reduce individual fitness when homozygous. This asexual population arose from a single worker ancestor approximately 20 years ago and has persisted for at least 100 generations. Thelytokous parthenogenesis results in a 1/3 of loss of heterozygosity with each generation. Yet, this population retains heterozygosity throughout its genome due to selection against homozygotes. Deep sequencing of one bee from each of the three known sub‐lineages of the population revealed that 3,766 of 10,884 genes (34%) have retained heterozygosity across all sub‐lineages, suggesting that these genes have heterozygote advantage. The maintenance of heterozygosity in the same genes and genomic regions in all three sub‐lineages suggests that nearly every chromosome carries genes that show sufficient heterozygote advantage to be selectively detrimental when homozygous.     

Negative heterozygosity-fitness correlations observed with microsatellites located in functional areas of the genome

Genome-wide heterozygosity inferred from neutral markers such as microsatellites is often expected to (i) reflect individual inbreeding and (ii) covary positively with fitness, generating positive heterozygosity-fitness correlations (HFCs). The often forgotten other end of the inbreeding-outbreeding continuum is outbreeding depression: past a certain degree of heterozygosity, heterozygotes tend to have lower fitness than homozygotes. Outbreeding depression arises from the breakup of co-adapted gene complexes and/or the introgression of nonlocally adapted genes. Provided that a correlation in heterozygosity exists across loci, outbreeding depression will be reflected in negative HFCs. In this issue, Olano-Marin et al. (2011a) describe negative heterozygosity-fitness correlations (HFCs) in blue tits Cyanistes caeruleus (Fig. 1), whereby heterozygosity has a significant, negative effect on female hatching success and recruitment. This study, together with a similar study by the same authors published in Evolution (Olano-Marin et al. 2011b), forms an original contribution in two respects. First, in the same population, positive and negative HFCs were recorded, revealing both inbreeding and outbreeding depression depending on the trait studied (whereby both processes were reliant on unknown, and possibly different, sets of coding loci). Second, a large number of microsatellite markers were split into two functional groups: microsatellite markers were either designed using zebra finch expressed sequence tags (ESTs) or derived using traditional cloning methods and presumed to be neutral. Contrasting large classes of loci and their varying levels of polymorphism, rather than looking for one locus that would stand out among tens of randomly selected markers, pave the way for a more elegant and powerful approach to explore how HFCs vary across traits and among regions of the genome. [Figure: see text].     

Heterozygosity-fitness correlations in the gilthead sea bream Sparus aurata using microsatellite loci from unknown and gene-rich genomic locations

Heterozygosity-fitness correlations (HFC) were assessed for a sample of a gilthead sea bream Sparus aurata population. Two hundred and seventy-one fish were genotyped at 22 known and novel microsatellite loci, from which correlations between the multilocus heterozygosity index (I(MLH) ) and various fitness traits (fork length, mass and specific growth rates) were calculated. Significant global HFCs were found in this sample (0·02 ≤r(2) ≤ 0·08). In addition, all the significant correlations found in this work were negative, indicating that heterozygotes had lower fitness than their homozygote counterparts. Marker location could not explain the observed HFCs. Evidence of inbreeding, outbreeding or population and family structuring was not found in this work. The presence of undetected general effects that may lead to the appearance of HFCs, however, cannot be ruled out. These results seem to be best explained by the occurrence of local effects (due to linkage) or even by possible direct locus advantages.     

Local heterozygosity-fitness correlations with global positive effects on fitness in threespine stickleback

The complex interactions between genetic diversity and evolution have important implications in many biological areas including conservation, speciation, and mate choice. A common way to study these interactions is to look at heterozygosity-fitness correlations (HFCs). Until recently, HFCs based on noncoding markers were believed to result primarily from global inbreeding effects. However, accumulating theoretical and empirical evidence shows that HFCs may often result from genes being linked to the markers used (local effect). Moreover, local effect HFCs could differ from global inbreeding effects in their direction and occurrence. Consequently, the investigation of the structure and consequences of local HFCs is emerging as a new important goal in evolutionary biology. In this study of a wild threespine stickleback (Gasterosteus aculeatus) population, we first tested the presence of significant positive or negative local effects of heterozygosity at 30 microsatellites loci on five fitness components: survival, mating success, territoriality, length, and body condition. Then, we evaluated the direction and shape of total impact of local HFCs, and estimated the magnitude of the impacts on fitness using regression coefficients and selection differentials. We found that multilocus heterozygosity was not a reliable estimator of individual inbreeding coefficient, which supported the relevance of single-locus based analyses. Highly significant and temporally stable local HFCs were observed. These were mainly positive, but negative effects of heterozygosity were also found. Strong and opposite effects of heterozygosity are probably present in many populations, but may be blurred in HFC analyses looking for global effects only. In this population, both negative and positive HFCs are apparently driving mate preference by females, which is likely to contribute to the maintenance of both additive and nonadditive genetic variance.     

Molecular and pedigree measures of relatedness provide similar estimates of inbreeding depression in a bottlenecked population

Individual‐based estimates of the degree of inbreeding or parental relatedness from pedigrees provide a critical starting point for studies of inbreeding depression, but in practice wild pedigrees are difficult to obtain. Because inbreeding increases the proportion of genomewide loci that are identical by descent, inbreeding variation within populations has the potential to generate observable correlations between heterozygosity measured using molecular markers and a variety of fitness related traits. Termed heterozygosity‐fitness correlations (HFCs), these correlations have been observed in a wide variety of taxa. The difficulty of obtaining wild pedigree data, however, means that empirical investigations of how pedigree inbreeding influences HFCs are rare. Here, we assess evidence for inbreeding depression in three life‐history traits (hatching and fledging success and juvenile survival) in an isolated population of Stewart Island robins using both pedigree‐ and molecular‐derived measures of relatedness. We found results from the two measures were highly correlated and supported evidence for significant but weak inbreeding depression. However, standardized effect sizes for inbreeding depression based on the pedigree‐based kin coefficients (k) were greater and had smaller standard errors than those based on molecular genetic measures of relatedness (RI), particularly for hatching and fledging success. Nevertheless, the results presented here support the use of molecular‐based measures of relatedness in bottlenecked populations when information regarding inbreeding depression is desired but pedigree data on relatedness are unavailable.     

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.