INBREEDING INFLUENCES WITHIN‐BROOD HETEROZYGOSITY‐FITNESS CORRELATIONS (HFCS) IN AN ISOLATED PASSERINE POPULATION

Molecular estimates of inbreeding may be made using genetic markers such as microsatellites, however the interpretation of resulting heterozygosity‐fitness correlations (HFCs) with respect to inbreeding depression is not straightforward. We investigated the relationship between pedigree‐determined inbreeding coefficients (f) and HFCs in a closely monitored, reintroduced population of Stewart Island robins (Petroica australis rakiura) on Ulva Island, New Zealand. Using a full sibling design, we focused on differences in juvenile survival associated specifically with individual sibling variation in standardized multilocus heterozygosity (SH) when expected f was identical. We found that within broods, siblings with higher SH at microsatellite loci experienced a higher probability of juvenile survival. This effect, however, was detected primarily within broods that experienced inbreeding or when inbreeding had occurred in their pedigree histories (i.e., at the parents’ level). Thus we show, for the first time in a wild population, that the strength of an HFC is partially dependent on the presence of inbreeding events in the recent pedigree history. Our results illustrate the importance of realized effects of inbreeding on genetic variation and fitness and the value of full‐sibling designs for the study of HFCs in the context of small, inbred populations.     

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.     

A high‐quality pedigree and genetic markers both reveal inbreeding depression for quality but not survival in a cooperative mammal

Inbreeding depression, the reduced fitness of offspring of closely related parents, is commonplace in both captive and wild populations and has important consequences for conservation and mating system evolution. However, because of the difficulty of collecting pedigree and life‐history data from wild populations, relatively few studies have been able to compare inbreeding depression for traits at different points in the life cycle. Moreover, pedigrees give the expected proportion of the genome that is identical by descent (IBD_g) whereas in theory with enough molecular markers realized IBD_g can be quantified directly. We therefore investigated inbreeding depression for multiple life‐history traits in a wild population of banded mongooses using pedigree‐based inbreeding coefficients (f_ped) and standardized multilocus heterozygosity (sMLH) measured at 35–43 microsatellites. Within an information theoretic framework, we evaluated support for either f_ped or sMLH as inbreeding terms and used sequential regression to determine whether the residuals of sMLH on f_ped explain fitness variation above and beyond f_ped. We found no evidence of inbreeding depression for survival, either before or after nutritional independence. By contrast, inbreeding was negatively associated with two quality‐related traits, yearling body mass and annual male reproductive success. Yearling body mass was associated with f_ped but not sMLH, while male annual reproductive success was best explained by both f_ped and residual sMLH. Thus, our study not only uncovers variation in the extent to which different traits show inbreeding depression, but also reveals trait‐specific differences in the ability of pedigrees and molecular markers to explain fitness variation and suggests that for certain traits, genetic markers may capture variation in realized IBD_g above and beyond the pedigree expectation.     

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.     

Analysis of the relationship between allozyme heterozygosity and fitness in the rare Gentiana pneumonanthe L.

Especially for rare species occurring in small populations, which are prone to loss of genetic variation and inbreeding, detailed knowledge of the relationship between heterozygosity and fitness is generally lacking. After reporting on allozyme variation and fitness in relation to population size in the rare plant Gentiana pneumonanthe, we present a more detailed analysis of the association between heterozygosity and individual fitness. The aim of this study was to test whether increased fitness of more heterozygous individuals is explained best by the ‘inbreeding’ hypothesis or by the ‘overdominance’ hypothesis. Individual fitness was measured during 8 months of growth in the greenhouse as the performance for six life-history parameters. PCA reduced these parameters to four main Fitness Components. Individual heterozygosity was scored for seven polymorphic allozyme loci. For some of these loci (e.g. Aat3, Pgm1 and 6Pgdh2) heterozygotes showed a significantly higher relative fitness than homozygotes. To test the inbreeding model, regression analyses were performed between each Fitness Component and the number of heterozygous loci per individual. Multiple regressions with the adaptive distance of five loci as independent variables were used to test the overdominance model. Only the inbreeding model was a statistically significant explanation for the relationship between heterozygosity and fitness in G. pneumonanthe. The number of heterozygous loci was significantly negatively correlated with the coefficients of variation of three of the six initially measured fitness parameters. This suggests a lower developmental stability among more homozygous plants and may explain the higher phenotypic variation in small populations of the species observed earlier. The importance of the results for conservation biology is discussed.     

Age‐dependent,negative heterozygosity–fitness correlations and local effects in an endangered Caribbean reptile,Iguana delicatissima

Inbreeding depression can have alarming impacts on threatened species with small population sizes. Assessing inbreeding has therefore become an important focus of conservation research. In this study, heterozygosity–fitness correlations (HFCs) were measured by genotyping 7 loci in 83 adult and 184 hatchling Lesser Antillean Iguanas, Iguana delicatissima, at a communal nesting site in Dominica to assess the role of inbreeding depression on hatchling fitness and recruitment to the adult population in this endangered species. We found insignificant correlations between multilocus heterozygosity and multiple fitness proxies in hatchlings and adults. Further, multilocus heterozygosity did not differ significantly between hatchlings and adults, which suggests that the survivorship of homozygous hatchlings does not differ markedly from that of their heterozygous counterparts. However, genotypes at two individual loci were correlated with hatching date, a finding consistent with the linkage between specific marker loci and segregating deleterious recessive alleles. These results provide only modest evidence that inbreeding depression influences the population dynamics of I. delicatissima on Dominica.     

Association between host's genetic diversity and parasite burden in damselflies

Recent research indicates that low genetic variation in individuals can increase susceptibility to parasite infection, yet evidence from natural invertebrate populations remains scarce. Here, we studied the relationship between genetic heterozygosity, measured as AFLP‐based inbreeding coefficient f_AFLP, and gregarine parasite burden from eleven damselfly, Calopteryx splendens, populations. We found that in the studied populations, 5–92% of males were parasitized by endoparasitic gregarines (Apicomplexa: Actinocephalidae). Number of parasites ranged from none to 47 parasites per male, and parasites were highly aggregated in a few hosts. Mean individual f_AFLP did not differ between populations. Moreover, we found a positive association between individual's inbreeding coefficient and parasite burden. In other words, the more homozygous the individual, the more parasites it harbours. Thus, parasites are likely to pose strong selection pressure against inbreeding and homozygosity. Our results support the heterozygosity‐fitness correlation hypothesis, which suggests the importance of heterozygosity for an individual's pathogen resistance.     

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.     

Effects of inbreeding on fitness‐related traits in a small isolated moose population

Inbreeding can affect fitness‐related traits at different life history stages and may interact with environmental variation to induce even larger effects. We used genetic parentage assignment based on 22 microsatellite loci to determine a 25 year long pedigree for a newly established island population of moose with 20–40 reproducing individuals annually. We used the pedigree to calculate individual inbreeding coefficients and examined for effects of individual inbreeding (f) and heterozygosity on fitness‐related traits. We found negative effects of f on birth date, calf body mass and twinning rate. The relationship between f and calf body mass and twinning rate were found to be separate but weaker after accounting for birth date. We found no support for an inbreeding effect on the age‐specific lifetime reproductive success of females. The influence of f on birth date was related to climatic conditions during the spring prior to birth, indicating that calves with a low f were born earlier after a cold spring than calves with high f. In years with a warm spring, calf f did not affect birth date. The results suggest that severe inbreeding in moose has both indirect effects on fitness through delayed birth and lower juvenile body mass, as well as separate direct effects, as there still was a significant relationship between f and twinning rate after accounting for birth date and body mass as calf. Consequently, severe inbreeding as found in the study population may have consequences for population growth and extinction risk.     

The use (or misuse) of microsatellite allelic distances in the context of inbreeding and conservation genetics

In line with inbreeding theory, genetic diversity at a set of molecular markers may explain variation in fitness‐associated traits in partially inbred populations, and such associations will appear as ‘genotype–fitness correlations’. An individual genetic diversity index specifically used for microsatellites is ‘mean d²’, i.e. the mean squared distance between alleles. The original hypothesis for mean d²–fitness correlations assumes that mean d² captures fitness effects at both ends of the inbreeding–outbreeding spectrum. This hypothesis received strong criticism from work showing that even a plain diversity estimate such as multi‐locus heterozygosity (MLH) outperforms mean d² as a predictor of the inbreeding coefficient and fitness in most realistic situations. Despite this critique, the mean d²‐approach is still used frequently in ecological and evolutionary research, producing results suggesting that mean d² sometimes provides a stronger prediction of fitness than does MLH. In light of the critique, such results are unexpected, but potential explanations for them may exist (at least hypothetically), including scenarios based on close linkage and recent admixture. Nevertheless, a major caveat is that it is very difficult to predict a priori if mean d² will improve the genotype–fitness correlation, which in turn makes objective interpretations difficult. Mean d²–fitness associations are potentially interesting, but the fact that we cannot easily understand them is problematic and should be thoroughly addressed in each study. Therefore, instead of hastily reached interpretations of mean d²–fitness correlations, conclusions need support from complementary analyses, e.g. verifying admixture of genetically structured populations.     

Differential allocation in a lekking bird: females lay larger eggs and are more likely to have male chicks when they mate with less related males

The differential allocation hypothesis predicts increased investment in offspring when females mate with high-quality males. Few studies have tested whether investment varies with mate relatedness, despite evidence that non-additive gene action influences mate and offspring genetic quality. We tested whether female lekking lance-tailed manakins (Chiroxiphia lanceolata) adjust offspring sex and egg volume in response to mate attractiveness (annual reproductive success, ARS), heterozygosity and relatedness. Across 968 offspring, the probability of being male decreased with increasing parental relatedness but not father ARS or heterozygosity. This correlation tended to diminish with increasing lay-date. Across 162 offspring, egg volume correlated negatively with parental relatedness and varied with lay-date, but was unrelated to father ARS or heterozygosity. Offspring sex and egg size were unrelated to maternal age. Comparisons of maternal half-siblings in broods with no mortality produced similar results, indicating differential allocation rather than covariation between female quality and relatedness or sex-specific inbreeding depression in survival. As males suffer greater inbreeding depression, overproducing females after mating with related males may reduce fitness costs of inbreeding in a system with no inbreeding avoidance, while biasing the sex of outbred offspring towards males may maximize fitness via increased mating success of outbred sons.     

Multilocus heterozygosity and inbreeding depression in an insular house sparrow metapopulation

Inbreeding causes reduction of genetic variability that may have severe fitness consequences. In spite of its potentially huge impact on viability and evolutionary processes especially in small populations, quantitative demonstrations of genetic and demographic effects of inbreeding in natural populations are few. Here, we examine the relationship between individual inbreeding coefficients (F) and individual standardized multilocus heterozygosity (H) in an insular metapopulation of house sparrows (Passer domesticus) in northern Norway in order to evaluate whether H is a good predictor for F. We then relate variation in fitness (i.e. the probability of surviving from fledging to recruitment) to F and H, which enables us to examine whether inbreeding depression is associated with a reduction in genetic variability. The average level of inbreeding in the house sparrow metapopulation was high, and there was large inter-individual variation in F. As expected, standardized multilocus heterozygosity decreased with the level of inbreeding. The probability of recruitment was significantly negatively related to F, and, accordingly, increased with H. However, H explained no significant additional variation in recruitment rate than was explained by F. This suggests that H is a good predictor for F in this metapopulation, and that an increase in F is likely to be associated with a general increase in the level of homozygosity on loci across the genome, which has severe fitness consequences.     

Heterozygosity-fitness correlations in a migratory bird: an analysis of inbreeding and single-locus effects

Studies in a multitude of taxa have described a correlation between heterozygosity and fitness and usually conclude that this is evidence for inbreeding depression. Here, we have used multilocus heterozygosity (MLH) estimates from 15 microsatellite markers to show evidence of heterozygosity-fitness correlations (HFCs) in a long-distance migratory bird, the light-bellied Brent goose. We found significant, positive heterozygosity-heterozygosity correlations between random subsets of the markers we employed, and no evidence that a model containing all loci as individual predictors in a multiple regression explained significantly more variation than a model with MLH as a single predictor. Collectively, these results lend support to the hypothesis that the HFCs we have observed are a function of inbreeding depression. However, we do find that fitness correlations are only detectable in years where population-level productivity is high enough for the reproductive asymmetry between high and low heterozygosity individuals to become apparent. We suggest that lack of evidence of heterozygosity-fitness correlations in animal systems may be because heterozygosity is a poor proxy measure of inbreeding, especially when employing low numbers of markers, but alternatively because the asymmetries between individuals of different heterozygosities may only be apparent when environmental effects on fitness are less pronounced.     

Direct and indirect causal effects of heterozygosity on fitness-related traits in Alpine ibex

Heterozygosity–fitness correlations (HFCs) are a useful tool to investigate the effects of inbreeding in wild populations, but are not informative in distinguishing between direct and indirect effects of heterozygosity on fitness-related traits. We tested HFCs in male Alpine ibex (Capra ibex) in a free-ranging population (which suffered a severe bottleneck at the end of the eighteenth century) and used confirmatory path analysis to disentangle the causal relationships between heterozygosity and fitness-related traits. We tested HFCs in 149 male individuals born between 1985 and 2009. We found that standardized multi-locus heterozygosity (MLH), calculated from 37 microsatellite loci, was related to body mass and horn growth, which are known to be important fitness-related traits, and to faecal egg counts (FECs) of nematode eggs, a proxy of parasite resistance. Then, using confirmatory path analysis, we were able to show that the effect of MLH on horn growth was not direct but mediated by body mass and FEC. HFCs do not necessarily imply direct genetic effects on fitness-related traits, which instead can be mediated by other traits in complex and unexpected ways.     

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.     

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.     

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.