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

Individual genetic diversity and probability of infection by avian malaria parasites in blue tits (Cyanistes caeruleus)

Understanding the importance of host genetic diversity for coping with parasites and infectious diseases is a long‐standing goal in evolutionary biology. Here, we study the association between probability of infection by avian malaria (Plasmodium relictum) and individual genetic diversity in three blue tit (Cyanistes caeruleus) populations that strongly differ in prevalence of this parasite. For this purpose, we screened avian malaria infections and genotyped 789 blue tits across 26 microsatellite markers. We used two different arrays of markers: 14 loci classified as neutral and 12 loci classified as putatively functional. We found a significant relationship between probability of infection and host genetic diversity estimated at the subset of neutral markers that was not explained by strong local effects and did not differ among the studied populations. This relationship was not linear, and probability of infection increased up to values of homozygosity by locus (HL) around 0.15, reached a plateau at values of HL from 0.15 to 0.40 and finally declined among a small proportion of highly homozygous individuals (HL > 0.4). We did not find evidence for significant identity disequilibrium, which may have resulted from a low variance of inbreeding in the study populations and/or the small power of our set of markers to detect it. A combination of subtle positive and negative local effects and/or a saturation threshold in the association between probability of infection and host genetic diversity in combination with increased resistance to parasites in highly homozygous individuals may explain the observed negative quadratic relationship. Overall, our study highlights that parasites play an important role in shaping host genetic variation and suggests that the use of large sets of neutral markers may be more appropriate for the study of heterozygosity–fitness correlations.     

Recent admixture generates heterozygosity–fitness correlations during the range expansion of an invading species

Admixture, the mixing of historically isolated gene pools, can have immediate consequences for the genetic architecture of fitness traits. Admixture may be especially important for newly colonized populations, such as during range expansion and species invasions, by generating heterozygosity that can boost fitness through heterosis. Despite widespread evidence for admixture during species invasions, few studies have examined the demographic history leading to admixture, how admixture affects the heterozygosity and fitness of invasive genotypes, and whether such fitness effects are maintained through time. We address these questions using the invasive plant Silene vulgaris, which shows evidence of admixture in both its native Europe and in North America where it has invaded. Using multilocus genotype data in conjunction with approximate Bayesian computation analysis of demographic history, we showed that admixture during the invasion of North America was independent from and much younger than admixture in the native range of Europe. We tested for fitness consequences of admixture in each range and detected a significant positive heterozygosity–fitness correlation (HFC) in North America; in contrast, no HFC was present in Europe. The lack of HFC in Europe may reflect the longer time since admixture in the native range, dissipating associations between heterozygosity at markers and fitness loci. Our results support a key short‐term role for admixture during the early stages of invasion by generating HFCs that carry populations past the threat of extinction from inbreeding and demographic stochasticity.     

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.     

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.     

Less can be more: loss of MHC functional diversity can reflect adaptation to novel conditions during fish invasions

The ability of invasive species to adapt to novel conditions depends on population size and environmental mismatch, but also on genetic variation. Away from their native range, invasive species confronted with novel selective pressures may display different levels of neutral versus functional genetic variation. However, the majority of invasion studies have only examined genetic variation at neutral markers, which may reveal little about how invaders adapt to novel environments. Salmonids are good model systems to examine adaptation to novel pressures because they have been translocated all over the world and represent major threats to freshwater biodiversity in the Southern Hemisphere, where they have become invasive. We examined patterns of genetic differentiation at seven putatively neutral (microsatellites) loci and one immune‐related major histocompatibility complex (MHC class II‐β) locus among introduced rainbow trout living in captivity (farmed) or under natural conditions (naturalized) in Chilean Patagonia. A significant positive association was found between differentiation at neutral and functional markers, highlighting the role of neutral evolutionary forces in shaping genetic variation at immune‐related genes in salmonids. However, functional (MHC) genetic diversity (but not microsatellite diversity) decreased with time spent in the wild since introduction, suggesting that there was selection against alleles associated with captive rearing of donor populations that do not provide an advantage in the wild. Thus, although high genetic diversity may initially enhance fitness in translocated populations, it does not necessarily reflect invasion success, as adaptation to novel conditions may result in rapid loss of functional MHC diversity.     

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.     

Correlations between heterozygosity and reproductive success in the blue tit (Cyanistes caeruleus): an analysis of inbreeding and single locus effects

To understand the mechanisms behind heterozygosity-fitness correlations (HFC), it is necessary to employ large numbers of markers with known function and independently estimate the variation in inbreeding in the population. Here we genotyped 794 blue tits with 79 microsatellites that were distributed across 25 chromosomes and that were classified either as "functional" (N= 58) or "neutral" (N= 21). We found a positive effect of individual heterozygosity at multiple loci on clutch size, on the number of eggs sired by males, and on the number of recruits produced by males and females. We documented the occurrence of some consanguineous matings and found evidence for a particular type of population structure that can contribute to the occurrence of inbreeding. As the set of "neutral" loci provided more power to detect HFC and identity disequilibrium, we argue that "neutral" markers are better predictors of the effects of inbreeding. The number of significant effects at single loci did not exceed the expected number of false positives and no strong effects were associated with heterozygosity at "functional" markers. Thus, the HFC found here cannot be attributed to strong effects of the loci under study.     

Population genomic evidence for adaptive differentiation in the Baltic Sea herring

Detecting and estimating the degree of genetic differentiation among populations of highly mobile marine fish having pelagic larval stages is challenging because their effective population sizes can be large, and thus, little genetic drift and differentiation is expected in neutral genomic sites. However, genomic sites subject to directional selection stemming from variation in local environmental conditions can still show substantial genetic differentiation, yet these signatures can be hard to detect with low‐throughput approaches. Using a pooled RAD‐seq approach, we investigated genomewide patterns of genetic variability and differentiation within and among 20 populations of Atlantic herring in the Baltic Sea (and adjacent Atlantic sites), where previous low‐throughput studies and/or studies based on few populations have found limited evidence for genetic differentiation. Stringent quality control was applied in the filtering of 1 791 254 SNPs, resulting in a final data set of 68 182 polymorphic loci. Clear differentiation was identified between Atlantic and Baltic populations in many genomic sites, while differentiation within the Baltic Sea area was weaker and geographically less structured. However, outlier analyses – whether including all populations or only those within the Baltic Sea – uncovered hundreds of directionally selected loci in which variability was associated with either salinity, temperature or both. Hence, our results support the view that although the degree of genetic differentiation among Baltic Sea herring populations is low, there are many genomic regions showing elevated divergence, apparently as a response to temperature‐ and salinity‐related natural selection. As such, the results add to the increasing evidence of local adaptation in highly mobile marine organisms, and those in the young Baltic Sea in particular.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Adaptive divergence along environmental gradients in a climate‐change‐sensitive mammal

In the face of predicted climate change, a broader understanding of biotic responses to varying environments has become increasingly important within the context of biodiversity conservation. Local adaptation is one potential option, yet remarkably few studies have harnessed genomic tools to evaluate the efficacy of this response within natural populations. Here, we show evidence of selection driving divergence of a climate‐change‐sensitive mammal, the American pika (Ochotona princeps), distributed along elevation gradients at its northern range margin in the Coast Mountains of British Columbia (BC), Canada. We employed amplified‐fragment‐length‐polymorphism‐based genomic scans to conduct genomewide searches for candidate loci among populations inhabiting varying environments from sea level to 1500 m. Using several independent approaches to outlier locus detection, we identified 68 candidate loci putatively under selection (out of a total 1509 screened), 15 of which displayed significant associations with environmental variables including annual precipitation and maximum summer temperature. These candidate loci may represent important targets for predicting pika responses to climate change and informing novel approaches to wildlife conservation in a changing world.     

Multiple sexual ornaments signal heterozygosity in male blue tits

Higher individual genetic quality has been hypothesized to be associated with the expression of conspicuous ornaments. However, the relationship between multicomponent sexual signals and heterozygosity is poorly understood. In this study, we examined whether different ornaments, including song (repertoire size and bout length) and plumage coloration (yellow breast and blue crown), reflect individual genetic diversity in male blue tits (Aves: Cyanistes caeruleus). We estimated genetic diversity using 26 microsatellite markers that were classified as putatively functional (12 loci) and neutral (14 loci). We found that yellow breast carotenoid chroma, blue crown brightness, bout length and body condition were positively associated with heterozygosity at functional loci, but not with genetic diversity estimated at all typed loci or the subset of neutral markers. The lack of strong single‐locus effects and the presence of identity disequilibrium in our population suggest that the observed heterozygosity‐phenotype associations are driven by loci widely distributed across the genome. The predominant role of putatively functional loci evidences that the expression of secondary sexual characters is more tightly reflected by heterozygosity at genomic regions containing coding genes that are being actively expressed, a fact that may make ornamental traits more reliable indicators of the genetic quality of individuals. Overall, this study shows that multiple secondary sexual characters reflect male genetic diversity and lends support to the good‐genes‐as‐heterozygosity hypothesis. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 362–375.     

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.     

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host–pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1′665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen‐mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host–pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.     

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

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

Genetic divergence and isolation by thermal environment in geothermal populations of an aquatic invertebrate

Temperature is one of the most influential forces of natural selection impacting all biological levels. In the face of increasing global temperatures, studies over small geographic scales allowing investigations on the effects of gene flow are of great value for understanding thermal adaptation. Here, we investigated genetic population structure in the freshwater gastropod Radix balthica originating from contrasting thermal habitats in three areas of geothermal activity in Iceland. Snails from 32 sites were genotyped at 208 AFLP loci. Five AFLPs were identified as putatively under divergent selection in Lake Mývatn, a geothermal lake with an almost 20 °C difference in mean temperature across a distance of a few kilometres. In four of these loci, variation across all study populations was correlated with temperature. We found significant population structure in neutral markers both within and between the areas. Cluster analysis using neutral markers classified the sites mainly by geography, whereas analyses using markers under selection differentiated the sites based on temperature. Isolation by distance was stronger in the neutral than in the outlier loci. Pairwise differences based on outlier F_ST were significantly correlated with temperature at different spatial scales, even after correcting for geographic distance or neutral pairwise F_ST differences. In general, genetic variation decreased with increasing environmental temperature, possibly suggesting that natural selection had reduced the genetic diversity in the warm origin sites. Our results emphasize the influence of environmental temperature on the genetic structure of populations and suggest local thermal adaptation in these geothermal habitats.     

Demographic responses to weather fluctuations are context dependent in a long‐lived amphibian

Weather fluctuations have been demonstrated to affect demographic traits in many species. In long‐lived organisms, their impact on adult survival might be buffered by the evolution of traits that reduce variation in interannual adult survival. For example, skipping breeding is an effective behavioral mechanism that may limit yearly variation in adult survival when harsh weather conditions occur; however, this in turn would likely lead to strong variation in recruitment. Yet, only a few studies to date have examined the impact of weather variation on survival, recruitment and breeding probability simultaneously in different populations of the same species. To fill this gap, we studied the impact of spring temperatures and spring rainfall on survival, on reproductive skipping behavior and on recruitment in five populations of a long‐lived amphibian, the yellow‐bellied toad (Bombina variegata). Based on capture–recapture data, our findings demonstrate that survival depends on interactions between age, population and weather variation. Varying weather conditions in the spring result in strong variation in the survival of immature toads, whereas they have little effect on adult toads. Breeding probability depends on both the individual's previous reproductive status and on the weather conditions during the current breeding season, leading to high interannual variation in recruitment. Crucially, we found that the impact of weather variation on demographic traits is largely context dependent and may thus differ sharply between populations. Our results suggest that studies predicting the impact of climate change on population dynamics should be taken with caution when the relationship between climate and demographic traits is established using only one population or few populations. We therefore highly recommend further research that includes surveys replicated in a substantial number of populations to account for context‐dependent variation in demographic processes.     

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.     

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

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

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.