Heterozygosity–behaviour correlations in nine‐spined stickleback (Pungitius pungitius) populations: contrasting effects at random and functional loci

The study of heterozygosity‐fitness correlations (HFCs) has a long history in the fields of ecology and evolutionary biology but remains controversial. Recently, it has been shown that the genetic distance of markers from functional loci can be an important factor to be considered in addition to marker numbers and variability. In this study, we investigated the correlation between individual heterozygosity and behaviour (aggression, boldness and feeding activity) in nine‐spined stickleback (Pungitius pungitius) individuals originating from four populations in two contrasting environments. Offspring of full‐sib families raised in a common garden setting were assessed for behaviour and genotyped using 84 microsatellite markers that were either located within or near behaviourally or physiologically important genes (termed ‘functional’) or were randomly selected. No associations were detected with any behavioural trait in any population or over all populations when genetic variability was measured using all 84 markers combined. However, when the markers were separated into three functional categories (behavioural, physiological and random), several significant associations were observed both with functional markers and random markers in one of the four populations assessed. Interestingly, contrasting correlations with behaviour were observed when using physiological gene (negative) and random (positive) markers. Upon dividing the physiological gene markers into further subcategories based on their specific physiological functions, a strong relationship between the heterozygosity of markers linked to osmoregulation‐related genes, and behaviour was revealed in the brackish water population. Our results indicate that both local (physiological) and general (neutral) effects are important in shaping behaviour and that heterozygosity–behaviour correlations are population dependent.     

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.     

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.     

Heterozygosity correlates with body size,nest site quality and productivity in a colonial waterbird,the whiskered tern (Chlidonias hybrida,Aves: Sternidae)

Positive effects of individual heterozygosity on naturally selected traits have been reported in wild populations of many animal taxa. The aim of this study was to test whether heterozygosity predicts the quality of acquired nest sites and productivity in a colonially breeding waterbird, the whiskered tern (Chlidonias hybrida). For this purpose, 40 adult terns from a small, recently established population in Central Poland were typed at eight microsatellite loci. We demonstrate that individual heterozygosity is positively related to hatching success. We hypothesize that this association could be mediated by direct effects of heterozygosity on the competitive abilities of individuals. We found that more heterozygous terns tended to breed in better protected central parts of the colony, suggesting that they had capabilities of outcompeting less heterozygous individuals and relegating them to the less attractive peripheries of the colony. It was also demonstrated that the link between heterozygosity and individual abilities to acquire more attractive nest site could be mediated by the larger size of heterozygous individuals. Although no correlations between heterozygosity and different components of condition were found, there was a positive association between female heterozygosity and both clutch size and egg size. We suggest that demonstrated heterozygosity‐fitness correlations could be primarily caused by inbreeding depression in the studied whiskered tern population.     

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.     

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.     

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.     

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.     

Specific alleles at immune genes,rather than genome‐wide heterozygosity,are related to immunity and survival in the critically endangered Attwater's prairie‐chicken

The negative effects of inbreeding on fitness are serious concerns for populations of endangered species. Reduced fitness has been associated with lower genome‐wide heterozygosity and immune gene diversity in the wild; however, it is rare that both types of genetic measures are included in the same study. Thus, it is often unclear whether the variation in fitness is due to the general effects of inbreeding, immunity‐related genes or both. Here, we tested whether genome‐wide heterozygosity (20 990 SNPs) and diversity at nine immune genes were better predictors of two measures of fitness (immune response and survival) in the endangered Attwater's prairie‐chicken (Tympanuchus cupido attwateri). We found that postrelease survival of captive‐bred birds was related to alleles of the innate (Toll‐like receptors, TLRs) and adaptive (major histocompatibility complex, MHC) immune systems, but not to genome‐wide heterozygosity. Likewise, we found that the immune response at the time of release was related to TLR and MHC alleles, and not to genome‐wide heterozygosity. Overall, this study demonstrates that immune genes may serve as important genetic markers when monitoring fitness in inbred populations and that in some populations specific functional genes may be better predictors of fitness than genome‐wide heterozygosity.     

Heterozygosity-fitness correlations in a bottlenecked island species: a case study on the Seychelles warbler

We used capture-mark-recapture models to investigate the effects of both individual and parental heterozygosity, measured at microsatellite loci on the survival of Seychelles warblers (Acrocephalus sechellensis), an endemic island species which went through a severe population bottleneck in the middle of the last century. We found that an individual's survival was not correlated with multilocus heterozygosity, or with heterozygosity at any specific locus. However, maternal, but not paternal, multilocus heterozygosity was positively associated with offspring survival, but only in years with low survival probabilities. A nestling cross-fostering experiment showed that this was a direct maternal effect as there was an effect of the genetic mother's, but not of the social mother's, heterozygosity. Heterozygosity-fitness correlations at microsatellite markers were generally assumed to reflect genome-wide effects. Although this might be true in partially inbred populations, such correlations may also arise as a result of local effects with specific markers being closely linked to genes which determine fitness. However, heterozygosity at the individual microsatellite loci was not correlated and therefore does not seem to reflect genome-wide heterozygosity. This suggests that even in a small bottlenecked population, heterozygosity-fitness correlations may not be caused by genome-wide effects. Support for the local effects hypothesis was also equivocal; although three specific loci were associated with offspring survival, including all single-locus heterozygosities as independent predictors for the variation in survival was not supported by the data. Furthermore, in contrast to the local effects hypothesis, the loci which contributed most to the heterozygosity-survival relationship were not more polymorphic than the other loci. This study highlights the difficulties in distinguishing between the two hypotheses.     

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.     

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.     

The strength of the association between heterozygosity and probability of interannual local recruitment increases with environmental harshness in blue tits

The extent of inbreeding depression and the magnitude of heterozygosity–fitness correlations (HFC) have been suggested to depend on the environmental context in which they are assayed, but little evidence is available for wild populations. We combine extensive molecular and capture–mark–recapture data from a blue tit (Cyanistes caeruleus) population to (1) analyze the relationship between heterozygosity and probability of interannual adult local recruitment and (2) test whether environmental stress imposed by physiologically suboptimal temperatures and rainfall influence the magnitude of HFC. To address these questions, we used two different arrays of microsatellite markers: 14 loci classified as neutral and 12 loci classified as putatively functional. We found significant relationships between heterozygosity and probability of interannual local recruitment that were most likely explained by variation in genomewide heterozygosity. The strength of the association between heterozygosity and probability of interannual local recruitment was positively associated with annual accumulated precipitation. Annual mean heterozygosity increased over time, which may have resulted from an overall positive selection on heterozygosity over the course of the study period. Finally, neutral and putatively functional loci showed similar trends, but the former had stronger effect sizes and seemed to better reflect genomewide heterozygosity. Overall, our results show that HFC can be context dependent, emphasizing the need to consider the role of environmental heterogeneity as a key factor when exploring the consequences of individual genetic diversity on fitness in natural populations.     

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.     

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

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

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

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

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

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

Microsatellite primers for fungus‐growing ants

We isolated five polymorphic microsatellite loci from a library of two thousand recombinant clones of two fungus‐growing ant species, Cyphomyrmex longiscapus and Trachymyrmex cf. zeteki. Amplification and heterozygosity were tested in five species of higher attine ants using both the newly developed primers and earlier published primers that were developed for fungus‐growing ants. A total of 20 variable microsatellite loci, developed for six different species of fungus‐growing ants, are now available for studying the population genetics and colony kin‐structure of these ants.     

Effect of microsatellite selection on individual and population genetic inferences: an empirical study using cross‐specific and species‐specific amplifications

Although whole‐genome sequencing is becoming more accessible and feasible for nonmodel organisms, microsatellites have remained the markers of choice for various population and conservation genetic studies. However, the criteria for choosing microsatellites are still controversial due to ascertainment bias that may be introduced into the genetic inference. An empirical study of red deer (Cervus elaphus) populations, in which cross‐specific and species‐specific microsatellites developed through pyrosequencing of enriched libraries, was performed for this study. Two different strategies were used to select the species‐specific panels: randomly vs. highly polymorphic markers. The results suggest that reliable and accurate estimations of genetic diversity can be obtained using random microsatellites distributed throughout the genome. In addition, the results reinforce previous evidence that selecting the most polymorphic markers leads to an ascertainment bias in estimates of genetic diversity, when compared with randomly selected microsatellites. Analyses of population differentiation and clustering seem less influenced by the approach of microsatellite selection, whereas assigning individuals to populations might be affected by a random selection of a small number of microsatellites. Individual multilocus heterozygosity measures produced various discordant results, which in turn had impacts on the heterozygosity‐fitness correlation test. Finally, we argue that picking the appropriate microsatellite set should primarily take into account the ecological and evolutionary questions studied. Selecting the most polymorphic markers will generally overestimate genetic diversity parameters, leading to misinterpretations of the real genetic diversity, which is particularly important in managed and threatened populations.