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.     

Heterozygosity-fitness correlations in zebra finches: microsatellite markers can be better than their reputation

Numerous studies have reported associations between heterozygosity in microsatellite markers and fitness-related traits (heterozygosity-fitness correlations, HFCs). However, it has often been questioned whether HFCs reflect general inbreeding depression, because a small panel of microsatellite markers does not reflect very well an individual's inbreeding coefficient (F) as calculated from a pedigree. Here, we challenge this prevailing view. Because of chance events during Mendelian segregation, an individual's realized proportion of the genome that is identical by descent (IBD) may substantially deviate from the pedigree-based expectation (i.e. F). This Mendelian noise may result in a weak correlation between F and multi-locus heterozygosity, but this does not imply that multi-locus heterozygosity is a bad estimator of realized IBD. We examined correlations between 11 fitness-related traits measured in up to 1192 captive zebra finches and three measures of inbreeding: (i) heterozygosity across 11 microsatellite markers, (ii) heterozygosity across 1359 single-nucleotide polymorphism (SNP) markers and (iii) F, based on a 5th-generation pedigree. All 11 phenotypic traits showed positive relationships with measures of heterozygosity, especially traits that are most closely related to fitness. Remarkably, the small panel of microsatellite markers produced equally strong HFCs as the large panel of SNP markers. Both marker-based approaches produced stronger correlations with phenotypes than the pedigree-based F, and this did not seem to result from the shortness of our pedigree. We argue that a small panel of microsatellites with high allelic richness may better reflect an individual's realized IBD than previously appreciated, especially in species like the zebra finch, where much of the genome is inherited in large blocks that rarely experience cross-over during meiosis.     

The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species

In nonpedigreed wild populations, inbreeding depression is often quantified through the use of heterozygosity-fitness correlations (HFCs), based on molecular estimates of relatedness. Although such correlations are typically interpreted as evidence of inbreeding depression, by assuming that the marker heterozygosity is a proxy for genome-wide heterozygosity, theory predicts that these relationships should be difficult to detect. Until now, the vast majority of empirical research in this area has been performed on generally outbred, nonbottlenecked populations, but differences in population genetic processes may limit extrapolation of results to threatened populations. Here, we present an analysis of HFCs, and their implications for the interpretation of inbreeding, in a free-ranging pedigreed population of a bottlenecked species: the endangered takahe (Porphyrio hochstetteri). Pedigree-based inbreeding depression has already been detected in this species. Using 23 microsatellite loci, we observed only weak evidence of the expected relationship between multilocus heterozygosity and fitness at individual life-history stages (such as survival to hatching and fledging), and parameter estimates were imprecise (had high error). Furthermore, our molecular data set could not accurately predict the inbreeding status of individuals (as 'inbred' or 'outbred', determined from pedigrees), nor could we show that the observed HFCs were the result of genome-wide identity disequilibrium. These results may be attributed to high variance in heterozygosity within inbreeding classes. This study is an empirical example from a free-ranging endangered species, suggesting that even relatively large numbers (>20) of microsatellites may give poor precision for estimating individual genome-wide heterozygosity. We argue that pedigree methods remain the most effective method of quantifying inbreeding in wild populations, particularly those that have gone through severe bottlenecks.     

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–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.     

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.     

Age-dependent genetic effects on a secondary sexual trait in male Alpine ibex, Capra ibex

Secondary sexual traits, such as horns in ungulates, may be good indicators of genetic quality because they are costly to develop. Genetic effects on such traits may be revealed by examining correlations between multilocus heterozygosity (MLH) and trait value. Correlations between MLH and fitness traits, termed heterozygosity-fitness correlations (HFC), may reflect inbreeding depression or associative overdominance of neutral microsatellite loci with loci directly affecting fitness traits. We investigated HFCs for horn growth, body mass and faecal counts of nematode eggs in wild Alpine ibex (Capra ibex). We also tested if individual inbreeding coefficients (f') estimated from microsatellite data were more strongly correlated with fitness traits than MLH. MLH was more strongly associated with trait variation than f'. We found HFC for horn growth but not for body mass or faecal counts of nematode eggs. The effect of MLH on horn growth was age-specific. The slope of the correlation between MLH and yearly horn growth changed from negative to positive as males aged, in accordance with the mutation accumulation theory of the evolution of senescence. Our results suggest that the horns of ibex males are an honest signal of genetic quality.     

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.     

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.     

Does heterozygosity estimate inbreeding in real populations?

Many recent studies report that individual heterozygosity at a handful of apparently neutral microsatellite markers is correlated with key components of fitness, with most studies invoking inbreeding depression as the likely underlying mechanism. The implicit assumption is that an individual's inbreeding coefficient can be estimated reliably using only 10 or so markers, but the validity of this assumption is unclear. Consequently, we have used individual-based simulations to examine the conditions under which heterozygosity and inbreeding are likely to be correlated. Our results indicate that the parameter space in which this occurs is surprisingly narrow, requiring that inbreeding events are both frequent and severe, for example, through selfing, strong population structure and/or high levels of polygyny. Even then, the correlations are strong only when large numbers of loci (~200) can be deployed to estimate heterozygosity. With the handful of markers used in most studies, correlations only become likely under the most extreme scenario we looked at, namely 20 demes of 20 individuals coupled with strong polygyny. This finding is supported by the observation that heterozygosity is only weakly correlated among markers within an individual, even in a dataset comprising 400 markers typed in diverse human populations, some of which favour consanguineous marriages. If heterozygosity and inbreeding coefficient are generally uncorrelated, then heterozygosity-fitness correlations probably have little to do with inbreeding depression. Instead, one would need to invoke chance linkage between the markers used and one or more gene(s) experiencing balancing selection. Unfortunately, both explanations sit somewhat uncomfortably with current understanding. If inbreeding is the dominant mechanism, then our simulations indicate that consanguineous mating would have to be vastly more common than is predicted for most realistic populations. Conversely, if heterosis provides the answer, there need to be many more polymorphisms with major fitness effects and higher levels of linkage disequilibrium than are generally assumed.     

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.     

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.     

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.     

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.     

TESTING THE INFLUENCE OF FAMILY STRUCTURE AND OUTBREEDING DEPRESSION ON HETEROZYGOSITY‐FITNESS CORRELATIONS IN SMALL POPULATIONS

Theory predicts that positive heterozygosity‐fitness correlations (HFCs) arise as a consequence of inbreeding, which is often assumed to have a strong impact in small, fragmented populations. Yet according to empirical data, HFC in such populations seem highly variable and unpredictable. We here discuss two overlooked phenomena that may contribute to this variation. First, in a small population, each generation may consist of a few families. This generates random correlations between particular alleles and fitness (AFCs, allele‐fitness correlations) and results in too liberal tests for HFC. Second, in some contexts, small populations receiving immigrants may be more impacted by outbreeding depression than by inbreeding depression, resulting in negative rather than positive HFC. We investigated these processes through a case study in tadpole cohorts of Pelodytes punctatus living in small ponds. We provide evidence for a strong family structure and significant AFC in this system, as well as an example of negative HFC. By simulations, we show that this negative HFC cannot be a spurious effect of family structure, and therefore reflects outbreeding depression in the studied population. Our example suggests that a detailed examination of AFC and HFC patterns can provide valuable insights into the internal genetic structure and sources of fitness variation in small populations.     

Heterozygosity and survival in blue tits (Cyanistes caeruleus): contrasting effects of presumably functional and neutral loci

The relationship between genetic diversity and fitness has important implications in evolutionary and conservation biology. This relationship has been widely investigated at the individual level in studies of heterozygosity-fitness correlations (HFC). General effects caused by inbreeding and/or local effects at single loci have been used as explanations of HFC, but the debate about the causes of HFC in open, natural populations is still ongoing. Study designs that control for variation in the inbreeding level of the individuals, and knowledge on the function and location of the markers used to measure heterozygosity, are fundamental to understand the causes of HFC. Here we investigated correlations between individual heterozygosity and estimates of survival at different life-history stages in an open population of blue tits (Cyanistes caeruleus). For survival at the embryo, nestling and fledgling stage, we used a full-sibling approach, i.e. we controlled for the level of inbreeding. We genotyped 1496 individuals with 79 microsatellites mapped across 25 chromosomes in the zebra finch (Taeniopygia guttata) that were classified either as potentially functional (58 loci) or as neutral (21 loci). We found different effects of standardized multilocus heterozygosity (SH): SH(functional) had a negative effect on the probability of hatching and local recruitment of females, whereas SH(neutral) had a positive effect on adult survival. The negative effects of functional loci are better explained by local effects, whereas the positive effects of neutral markers could reflect inbreeding effects in the population. Our results highlight the importance of considering the characteristics of the markers used in HFC studies and confirm the mixed effects of heterozygosity in different contexts (e.g. sex and life-history stage).     

SELECTION ON OVERDOMINANT GENES MAINTAINS HETEROZYGOSITY ALONG MULTIPLE CHROMOSOMES IN A CLONAL LINEAGE OF HONEY BEE

Correlations between fitness and genome‐wide heterozygosity (heterozygosity‐fitness correlations, HFCs) have been reported across a wide range of taxa. The genetic basis of these correlations is controversial: do they arise from genome‐wide inbreeding (“general effects”) or the “local effects” of overdominant loci acting in linkage disequilibrium with neutral loci? In an asexual thelytokous lineage of the Cape honey bee (Apis mellifera capensis), the effects of inbreeding have been homogenized across the population, making this an ideal system in which to detect overdominant loci, and to make inferences about the importance of overdominance on HFCs in general. Here we investigate the pattern of zygosity along two chromosomes in 42 workers from the clonal Cape honey bee population. On chromosome III (which contains the sex‐locus, a gene that is homozygous‐lethal) and chromosome IV we show that the pattern of zygosity is characterized by loss of heterozygosity in short regions followed by the telomeric restoration of heterozygosity. We infer that at least four selectively overdominant genes maintain heterozygosity on chromosome III and three on chromosome IV via local effects acting on neutral markers in linkage disequilibrium. We conclude that heterozygote advantage and local effects may be more common and evolutionarily significant than is generally appreciated.     

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.     

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).