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

On the correlation between heterozygosity and fitness in natural populations

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

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

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.     

Are extra-pair young better than within-pair young? A comparison of survival and dominance in alpine marmot

1. In socially monogamous species, females may seek extra-pair copulation to gain genetic benefits. In order to test this 'genetic quality' hypothesis, one must compare the performance of extra-pair young (EPY) and within-pair young (WPY). Such tests, however, are scarce and results published so far are inconclusive. 2. Here, we test the 'genetic quality' hypothesis using multistate capture-recapture models to compare age-specific survival and access to dominance between EPY and WPY in the alpine marmot Marmota marmota, a socially monogamous mammal showing extra-pair paternities. 3. When compared with WPY, survival of EPY was higher by 15%, 10% and 30%, for juveniles, yearlings and 2-year-old individuals, respectively. Survival at older ages did not differ. 4. Survival corresponded to true survival for yearlings and juveniles as dispersal does not occur before 2 years of age in marmots. For older individuals, survival estimates included a mixture of survival and dispersal. The 30% increase of the 2-year-old EPY survival might reflect delayed dispersal rather than high survival of EPY as compared with WPY. 5. WPY and EPY had the same probability (0.28) to access dominance at 2 years of age, but EPY were more successful at older ages than WPY (0.46 vs. 0.10). 6. Both survival and reproductive performance were higher in EPY than in WPY. The fitness advantages of adopting such a mixed mating tactic are thus likely to be high for marmot females. We suggest that obtaining genetic benefits is the main evolutionary force driving extra-pair paternity in alpine marmots.     

Parasites as mediators of heterozygosity-fitness correlations in the Great Tit (Parus major)

Positive correlations between heterozygosity and fitness traits are frequently observed, and it has been hypothesized, but rarely tested experimentally, that parasites play a key role in mediating the heterozygosity-fitness association. We evaluated this hypothesis in a wild great tit (Parus major) population by testing the prediction that the heterozygosity-fitness association would appear in broods experimentally infested with a common ectoparasite, but not in parasite-free broods. We simultaneously assessed the effects of parental and offspring heterozygosity on nestling growth and found that body mass of nestlings close to independence, which is a strong predictor of post-fledging survival, increased significantly with nestling levels of heterozygosity in experimentally infested nests, but not in parasite-free nests. Heterozygosity level of the fathers also showed a significant positive correlation with offspring body mass under an experimental parasite load, whereas there was no correlation with the mothers' level of heterozygosity. Thus, our results indicate a key role for parasites as mediators of the heterozygosity-fitness correlations.     

Major histocompatibility complex heterozygosity enhances reproductive success

We investigated how heterozygosity at the major histocompatibility complex (MHC) affects fitness in wild-derived (F2) house mice (Mus musculus musculus). To compare and control for potential confounding effects from close inbreeding and genome-wide heterozygosity, we used mice that were systematically outbred. We assessed how heterozygosity at MHC and background loci (using 15 microsatellite markers on 11 different chromosomes) affects individual survival and reproductive success (RS) in large, semi-natural population enclosures. We found that overall heterozygosity significantly increased RS, and this correlation was entirely explained by heterozygosity at two MHC loci. Moreover, we found that the effects of MHC heterozygosity depend on the level of background heterozygosity, and the benefits of maximal MHC heterozygosity show a curvilinear effect with increasing background heterozygosity. The enhanced RS from MHC heterozygosity was not because of increased survival, and although MHC heterozygosity was correlated with body mass, body mass did not correlate with RS when heterozygosity is controlled. Breeders were more MHC heterozygous than nonbreeders for both sexes, indicating that MHC heterozygosity enhanced fecundity, mating success or both. Our results show that (i) MHC heterozygosity enhances fitness among wild, outbred as well as congenic laboratory mice; (ii) heterozygosity-fitness correlations can potentially be explained by a few loci, such as MHC; (iii) MHC heterozygosity can increase fitness, even without affecting survival, by increasing mating and RS; and (iv) MHC effects depend on background genes, and maximal MHC heterozygosity is most beneficial at intermediate or optimal levels of background heterozygosity.     

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.     

A quantitative model of the relationship between phenotypic variance and heterozygosity at marker loci under partial selfing.

Negative relationships between allozyme heterozygosity and morphological variance have often been observed and interpreted as evidence for increased developmental stability in heterozygotes. However, inbreeding can also generate such relationships by decreasing heterozygosity at neutral loci and redistributing genetic variance at the same time. I here provide a quantitative genetic model of this process by analogy with heterozygosity-fitness relationships. Inbreeding generates negative heterozygosity-variance relationships irrespective of the genetic architecture of the trait. This holds for fitness traits as well as neutral traits, the effect being stronger for fitness traits under directional dominance or overdominance. The order of magnitude of heterozygosity-variance regressions is compatible with empirical data even with very low inbreeding. Although developmental stability effects cannot be excluded, inbreeding is a parsimonious explanation that should be seriously considered to explain correlations between heterozygosity and both mean and variance of phenotypes in natural populations.     

Multilocus heterozygosity and inbreeding in the Siberian jay

Because of its common negative association with fitness, inbreeding is a major concern in conservation biology. Traditionally it has been measured as individual inbreeding coefficient calculated from the pedigree, but recently multilocus heterozygosity estimates have become commonly used as proxies. However, theoretical and simulation studies have cast doubt on the validity of these surrogates especially when they are based on only a few molecular markers. Yet, empirical studies reporting the correlation between multilocus heterozygosity and inbreeding coefficient are rare. We studied this relationship in a wild Siberian jay (Perisoreus infaustus) population subject to a long-term field study over 30 years. The correlations between inbreeding coefficient and the employed heterozygosity measures—standardized heterozygosity and internal relatedness—based on 21 microsatellite loci were weak. These results together with results from theoretical and simulation studies caution against use of multilocus heterozygosity estimates to study inbreeding in natural populations.     

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

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

Heterozygosity-fitness correlations in 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.     

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.     

A new likelihood estimator and its comparison with moment estimators of individual genome-wide diversity

The inbreeding coefficient of an individual, F, is one of the central parameters in population genetics theory. It has found important applications in evolutionary biology, conservation and ecology, such as the study of inbreeding depression. In the absence of detailed and reliable pedigree records, researchers have developed quite a few estimators to estimate F or the genome-wide homozygosity from genetic marker data. The statistical properties and comparative performances of these metrics are rarely known, however, which impedes an informed choice of the most appropriate one in practical applications. In this investigation, I propose a new likelihood F estimator that makes efficient use of marker information and takes into account of allelic dropouts, null alleles and prior knowledge of inbreeding. I compare the likelihood estimator with three moment estimators of F and three metrics of genomic homozygosity (or heterozygosity) by analysing both simulated and empirical datasets. It is shown that the likelihood estimator invariably outperforms the other estimators and metrics across all datasets analysed. For a typical dataset in heterozygosity-fitness correlation studies involving 10-20 microsatellites and 50 individuals, the correlation between the likelihood estimator and F (the simulated true inbreeding coefficient) is about 8 ～ 35% higher than that between the moment estimators and F. A frequently applied metric, multilocus heterozygosity (MLH), and an F estimator based on the consideration of the proportion of alleles in homozygous conditions, [F R'), are shown to have particularly poor performances. The low correlation between MLH and fitness traits, which is widely observed in numerous empirical studies, might be partially caused by the adoption of this inefficient estimator of genomic inbreeding.     

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.     

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.     

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 predicts clutch and egg size but not plasticity in a house sparrow population with no evidence of inbreeding

We investigated the link between heterozygosity and the reaction norm attributes of reproductive performance in female house sparrows (Passer domesticus). We collected data on clutch size, egg size, hatching success and nestling survival in 2816 nesting attempts made by 791 marked individuals over a 16-year period. Pedigree analysis revealed no evidence of inbreeding. Neither parent-offspring regression nor an animal model revealed significant heritability in clutch or egg size. We selected 42 females that laid at least seven clutches at our study site and used a survey of 21 autosomal microsatellite loci to estimate heterozygosity for each female. We controlled for phenotypic plasticity and found that both clutch and egg size showed significant positive correlations with heterozygosity. We found no evidence that heterozygosity influenced the slope of individual reaction norms. Further analysis suggested that clutch size was affected by heterozygosity across the genome, but egg size had more complex relationships, with evidence favouring the influence of multiple loci. Given the apparent lack of inbreeding and large population size, our results suggest associative overdominance as the likely mechanism for the impact of heterozygosity, but also created a puzzle about the process producing associations between neutral markers and the genes affecting clutch size or egg size. One possible explanation is a long-term residual effect of the historical bottleneck that occurred when house sparrows were introduced into North America. The existence of heterozygosity-fitness correlations in a population with considerable phenotypic plasticity and little inbreeding implies that the effects of heterozygosity may be more significant than previously thought.     

Heterozygosity-fitness correlations of conserved microsatellite markers in Kentish plovers Charadrius alexandrinus

Heterozygosity-fitness correlations (HFCs) are frequently used to examine the relationship between genetic diversity and fitness. Most studies have reported positive HFCs, although there is a strong bias towards investigating HFCs in genetically impoverished populations. We investigated HFCs in a large genetically diverse breeding population of Kentish plovers Charadrius alexandrinus in southern Turkey. This small shorebird exhibits highly variable mating and care systems, and it is becoming an ecological model species to understand breeding system evolution. Using 11 conserved and six anonymous microsatellite markers, we tested whether and how heterozygosity was associated with chick survival, tarsus and body mass growth controlling for nongenetic effects (chick sex, hatching date, length of biparental care and site quality) that influence survival and growth. There was no genome-wide effect of heterozygosity on fitness, and we did not find any significant effects of heterozygosity on growth rates. However, two of the 11 conserved markers displayed an association with offspring survival: one marker showed a positive HFC, whereas the other marker showed a negative HFC. Heterozygosity at three further conserved loci showed significant interaction with nongenetic variables. In contrast, heterozygosity based on anonymous microsatellite loci was not associated with fitness or growth. Markers that were correlated with chick survival were not more likely to be located in exons or introns than other markers that lacked this association.