Heterozygosity is unrelated to adult fitness measures in a large, noninbred population of great tits (Parus major)

The extent to which heterozygosity-fitness correlations (HFCs) are expected in wild populations is an important and unresolved question in evolutionary biology, because it relates to our understanding of the genetic architecture of fitness. Here, we report a study of HFCs in a wild, noninbred population of great tits (Parus major), based on a sample comprising 281 individuals typed at 26 markers, resulting in a data set comprising over 5600 genotypes. We regressed pedigree-derived f-score and multilocus genetic diversity against eight life-history traits known to be associated with fitness in this population, including lifetime reproductive success (LRS), as well as several morphological traits under weak selection. We found no evidence for either multilocus or single-locus HFCs for any morphological or fitness trait, and further found no evidence that effect sizes were stronger for those life-history traits more closely associated with reproductive fitness. This result may, in part, be explained by the fact that we found no evidence that our set of 26 markers had any power to infer genome-wide heterozygosity in this population and that marker-derived heterozygosity was uncorrelated with pedigree-derived f-score. Overall, these results emphasize the fact that the often-reported strong HFCs detected in small, inbred populations do not reflect a general phenomenon of increasing individual reproductive fitness with increasing heterozygosity.     

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.     

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.     

Contrasting heterozygosity‐fitness correlations across life in a long‐lived seabird

Selection is a central force underlying evolutionary change and can vary in strength and direction, for example across time and space. The fitness consequences of individual genetic diversity have often been investigated by testing for multilocus heterozygosity‐fitness correlations (HFCs), but few studies have been able to assess HFCs across life stages and in both sexes. Here, we test for HFCs using a 26‐year longitudinal individual‐based data set from a large population of a long‐lived seabird (the common tern, Sterna hirundo), where 7,974 chicks and breeders of known age were genotyped at 15 microsatellite loci and sampled for life‐history traits over the complete life cycle. Heterozygosity was not correlated with fledging or post‐fledging prospecting probabilities, but was positively correlated with recruitment probability. For breeders, annual survival was not correlated with heterozygosity, but annual fledgling production was negatively correlated with heterozygosity in males and highest in intermediately heterozygous females. The contrasting HFCs among life stages and sexes indicate differential selective processes and emphasize the importance of assessing fitness consequences of traits over complete life histories.     

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.     

Estimating genome-wide heterozygosity: effects of demographic history and marker type

Heterozygosity–fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.     

Contrasting heterozygosity-fitness correlations between populations of a self-compatible shrub in a fragmented landscape

The mechanisms underlying heterozygosity-fitness correlations (HFCs) are subject of intense debates, especially about how important population features such as size or degree of isolation influence HFCs. Here, we report variation in HFCs between Large and Small populations of a self-compatible shrub (Myrtus communis) occurring within an extremely fragmented landscape. In each of the five study populations, we obtained data on both heterozygosity and fitness for 9-12 maternal families (i.e. offspring from the same mother plant). Whereas heterozygosity explained most of the variance (60-86?%) in growth rate of seedling families within Large populations, this relationship was absent within Small populations. Our results suggest that inbreeding may explain the observed HFCs within Large populations, and that different genetic processes (such as genetic drift and/or selection) could have overridden HFCs within Small populations. While it is difficult to draw general conclusions from five populations, we think our results open new research perspectives on how different genetic processes underlie variation in HFCs under different population contexts. Our study also points to a need for further attention on the complex relationships between heterozygosity in self-compatible plants and their progeny in relation to mating system variation. Finally, our results provide interesting new insights into how population genetic diversity is maintained or lost in a highly fragmented landscape.     

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.     

A quantitative review of heterozygosity–fitness correlations in animal populations

The ease of obtaining genotypic data from wild populations has renewed interest in the relationship between individual genetic diversity and fitness-related traits (heterozygosity–fitness correlations, or HFC). Here we present a comprehensive meta-analysis of HFC studies using powerful multivariate techniques which account for nonindependence of data. We compare these findings with those from univariate techniques, and test the influence of a range of factors hypothesized to influence the strength of HFCs. We found small but significantly positive effect sizes for life-history, morphological, and physiological traits; while theory predicts higher mean effect sizes for life-history traits, effect size did not differ consistently with trait type. Newly proposed measures of variation were no more powerful at detecting relationships than multilocus heterozygosity, and populations predicted to have elevated inbreeding variance did not exhibit higher mean effect sizes. Finally, we found evidence for publication bias, with studies reporting weak, nonsignificant effects being under-represented in the literature. In general, our review shows that HFC studies do not generally reveal patterns predicted by population genetic theory, and are of small effect (less than 1% of the variance in phenotypic characters explained). Future studies should use more genetic marker data and utilize sampling designs that shed more light on the biological mechanisms that may modulate the strength of association, for example by contrasting the strength of HFCs in mainland and island populations of the same species, investigating the role of environmental stress, or by considering how selection has shaped the traits under investigation.     

Heterozygosity–fitness correlations in blue tit nestlings (<Emphasis Type="Italic">Cyanistis caeruleus</Emphasis>) under contrasting rearing conditions

Understanding the relation between genetic variation and fitness remains a key question in evolutionary biology. Although heterozygosity has been reported to correlate with many fitness-related traits, the strength of the heterozygosity–fitness correlations (HFCs) is usually weak and it is still difficult to assess the generality of these associations in natural populations. It has been suggested that HFCs may become meaningful only under particular environmental conditions. Moreover, existing evidence suggests that HFCs may also differ between sexes. The aim of this study was to investigate correlations between heterozygosity in neutral markers (microsatellites) and fitness-related traits in a natural population of blue tits (Cyanistes caeruleus). Additionally, we tested whether sex and environmental conditions may influence the magnitude and direction of HFCs. We found a positive relationship between heterozygosity and body mass of 14 days post-hatching nestlings, but only among females. Our results suggest that the correlation between heterozygosity and nestling body mass observed among female offspring could be attributed to within-brood effects. We failed to find any evidence that environmental conditions as simulated by brood size manipulation affect HFCs.     

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.     

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 declining seabird population

Loss of genetic diversity is thought to lead to increased risk of extinction in endangered populations due to decreasing fitness of homozygous individuals. Here, we evaluated the presence of inbreeding depression in a long‐lived seabird, the European shag (Phalacrocorax aristotelis), after a severe decline in population size by nearly 70%. During three reproductive seasons, 85 breeders were captured and genotyped at seven microsatellite loci. Nest sites were monitored during the breeding season to estimate reproductive success as the number of chicks surviving to full‐size‐grown per nest. Captured birds were tagged with a ring with an individual code, and resighting data were collected during 7‐year period. We found a strong effect of multilocus heterozygosity on female reproductive performance, and a significant, although weaker, effect on breeder survival. However, our matrix population model suggests that this relatively small effect of genetic diversity on breeder survival may have a profound effect on fitness. This highlights the importance of integrating life history consequences in HFC studies. Importantly, heterozygosity was correlated across loci, suggesting that genomewide effects, rather than single loci, are responsible for the observed HFCs. Overall, the HFCs are a worrying symptom of genetic erosion in this declining population. Many long‐lived species are prone to extinction, and future studies should evaluate the magnitude of fitness impact of genetic deterioration on key population parameters, such as survival of breeders.     

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.     

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.     

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.     

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.     

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.     

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.