Genetic variability is unrelated to growth and parasite infestation in natural populations of the European eel (Anguilla anguilla)

Positive correlations between individual genetic heterozygosity and fitness-related traits (HFCs) have been observed in organisms as diverse as plants, marine bivalves, fish or mammals. HFCs are not universal and the strength and stability of HFCs seem to be variable across species, populations and ages. We analysed the relationship between individual genetic variability and two different estimators of fitness in natural samples of European eel, growth rate (using back-calculated length-at-age 1, 2 and 3) and parasite infestation by the swimbladder nematode Anguillicola crassus . Despite using a large data set of 22 expressed sequence tags-derived microsatellite loci and a large sample size of 346 individuals, no heterozygote advantage was observed in terms of growth rate or parasite load. The lack of association was evidenced by (i) nonsignificant global HFCs, (ii) a Multivariate General Linear Model showing no effect of heterozygosity on fitness components, (iii) single-locus analysis showing a lower number of significant tests than the expected false discovery rate, (iv) sign tests showing only a significant departure from expectations at one component, and, (v) a random distribution of significant single-locus HFCs that was not consistent across fitness components or sampling sites. This contrasts with the positive association observed in farmed eels in a previous study using allozymes, which can be explained by the nature of the markers used, with the allozyme study including many loci involved in metabolic energy pathways, while the expressed sequence tags-linked microsatellites might be located in genes or in the proximity of genes uncoupled with metabolism/growth.     

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.     

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.     

Linking genetic mechanisms of heterozygosity-fitness correlations to footprints of selection at single loci

Investigations of heterozygosity-fitness correlations (HFCs) are central to the understanding how genetic diversity is maintained in natural populations. Advanced genome-wide approaches will enrich the number of functional loci to be tested. We argue that a combined analysis of the genetic mechanisms of HFCs and selection signals at single loci will allow researchers to better understand the micro-evolutionary basis of HFCs. Different dominance relationships among the alleles at the locus can lead to positive, negative or null HFCs depending on the allele frequency distribution. These scenarios differ in the temporal stability of the HFCs and in the patterns of allele frequency changes over time. Here, we describe a simple theoretical framework that links the analyses of heterozygosity-fitness associations (ecological timescale) with tests for selection signals (evolutionary timescale). Different genomic footprints of selection can be expected for the different underlying genetic mechanisms of HFCs, and this information can be independently used for the classification of HFCs. We suggest that in addition to inbreeding and single-locus overdominant effects also loci under directional selection could play a significant role in the development of heterozygosity-fitness effects in large natural populations under recent or fluctuating ecological changes.     

Heterozygosity–fitness correlations and associative overdominance: new detection method and proof of principle in the Iberian wild boar

Heterozygosity-fitness correlations (HFC) may result from a genome-wide process — inbreeding — or local effects within the genome. The majority of empirical studies reporting HFCs have attributed correlations to inbreeding depression. However, HFCs are unlikely to be caused by inbreeding depression because heterozygosity measured at a small number of neutral markers is unlikely to accurately capture a genome-wide pattern. Testing the strengths of localized effects caused by associative overdominance has proven challenging. In their current paper, Amos and Acevedo-Whitehouse present a novel test for local HFCs. Using stochastic simulations, they determine the conditions under which single-locus HFCs arise, before testing the strength of the correlation between the neutral marker and a linked gene under selection in their simulations. They used insights gained from simulation to statistically investigate the likely cause of correlations between heterozygosity and disease status using data on bovine tuberculosis infections in a wild boar population. They discover that a single microsatellite marker is an excellent predictor of tuberculosis progression in infected individuals. The results are relevant for wild boar management but, more generally, they demonstrate how single-locus HFCs could be used to identify coding loci under selection in free-living populations.     

Sex-Specific Correlations of Individual Heterozygosity,Parasite Load,and Scalation Asymmetry in a Sexually Dichromatic Lizard

Heterozygosity-fitness correlations (HFCs) provide insights into the genetic bases of individual fitness variation in natural populations. However, despite decades of study, the biological significance of HFCs is still under debate. In this study, we investigated HFCs in a large population of the sexually dimorphic lizard Takydromus viridipunctatus (Lacertidae). Because of the high prevalence of parasitism from trombiculid mites in this lizard, we expect individual fitness (i.e., survival) to decrease with increasing parasite load. Furthermore, because morphological asymmetry is likely to influence individuals'' mobility (i.e., limb asymmetry) and male biting ability during copulation (i.e., head asymmetry) in this species, we also hypothesize that individual fitness should decrease with increasing morphological asymmetry. Although we did not formally test the relationship between morphological asymmetry and fitness in this lizard, we demonstrated that survival decreased with increasing parasite load using a capture-mark-recapture data set. We used a separate sample of 140 lizards to test the correlations between individual heterozygosity (i.e., standardized mean d² and HL based on 10 microsatellite loci) and the two fitness traits (i.e., parasite load and morphological asymmetry). We also evaluated and excluded the possibility that single-locus effects produced spurious HFCs. Our results suggest male-only, negative correlations between individual heterozygosity and parasite load and between individual heterozygosity and asymmetry, suggesting sex-specific, positive HFCs. Male T. viridipunctatus with higher heterozygosity tend to have lower parasite loads (i.e., higher survival) and lower asymmetry, providing a rare example of HFC in reptiles.     

Heterozygosity at a single locus explains a large proportion of variation in two fitness‐related traits in great tits: a general or a local effect?

In natural populations, mating between relatives can have important fitness consequences due to the negative effects of reduced heterozygosity. Parental level of inbreeding or heterozygosity has been also found to influence the performance of offspring, via direct and indirect parental effects that are independent of the progeny own level of genetic diversity. In this study, we first analysed the effects of parental heterozygosity and relatedness (i.e. an estimate of offspring genetic diversity) on four traits related to offspring viability in great tits (Parus major) using 15 microsatellite markers. Second, we tested whether significant heterozygosity–fitness correlations (HFCs) were due to ‘local’ (i.e. linkage to genes influencing fitness) and/or ‘general’ (genome‐wide heterozygosity) effects. We found a significant negative relationship between parental genetic relatedness and hatching success, and maternal heterozygosity was positively associated with offspring body size. The characteristics of the studied populations (recent admixture, polygynous matings) together with the fact that we found evidence for identity disequilibrium across our set of neutral markers suggest that HFCs may have resulted from genome‐wide inbreeding depression. However, one locus (Ase18) had disproportionately large effects on the observed HFCs: heterozygosity at this locus had significant positive effects on hatching success and offspring size. It suggests that this marker may lie near to a functional locus under selection (i.e. a local effect) or, alternatively, heterozygosity at this locus might be correlated to heterozygosity across the genome due to the extensive ID found in our populations (i.e. a general effect). Collectively, our results lend support to both the general and local effect hypotheses and reinforce the view that HFCs lie on a continuum from inbreeding depression to those strictly due to linkage between marker loci and genes under selection.     

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

HETEROZYGOSITY-FITNESS CORRELATIONS REVEALED BY NEUTRAL AND CANDIDATE GENE MARKERS IN ROE DEER FROM A LONG-TERM STUDY

Heterozygosity-fitness correlations (HFCs) are increasingly reported but the underlying mechanisms causing HFCs are generally poorly understood. Here, we test for HFCs in roe deer ( Capreolus capreolus ) using 22 neutral microsatellites widely distributed in the genome and four microsatellites in genes that are potentially under selection. Juvenile survival was used as a proxy for individual fitness in a population that has been intensively studied for 30 years in northeastern France. For 222 juveniles, we computed two measures of genetic diversity: individual heterozygosity ( H ), and mean d ² (relatedness of parental genomes). We found a relationship between genetic diversity and fitness both for the 22 neutral markers and two candidate genes: IGF1 (Insulin-like Growth Factor I) and NRAMP (natural resistance-associated macrophage protein). Statistical evidence and the size of genetic effects on juvenile survival were comparable to those reported for early development and cohort variation, suggesting a substantial influence of genetic components on fitness in this roe deer population. For the 22 neutral microsatellites, a correlation with fitness was revealed for mean d ², but not for H , suggesting a possible outbreeding advantage. This heterosis effect could have been favored by introduction of genetically distant (Hungarian) roe deer to the population in recent times and, possibly, by the structuring of the population into distinct clans. The locus-specific correlations with fitness may be driven by growth rate advantages and resistance to diseases known to exist in the studied population. Our analyses of neutral and candidate gene markers both suggest that the observed HFCs are likely mainly due to linkage with dominant or overdominant loci that affect fitness ("local" effect) rather than to a genome-wide relationship with homozygosity due to inbreeding ("general" effect).     

Heterozygosity-fitness correlations in the gilthead sea bream Sparus aurata using microsatellite loci from unknown and gene-rich genomic locations

Heterozygosity-fitness correlations (HFC) were assessed for a sample of a gilthead sea bream Sparus aurata population. Two hundred and seventy-one fish were genotyped at 22 known and novel microsatellite loci, from which correlations between the multilocus heterozygosity index (I(MLH) ) and various fitness traits (fork length, mass and specific growth rates) were calculated. Significant global HFCs were found in this sample (0·02 ≤r(2) ≤ 0·08). In addition, all the significant correlations found in this work were negative, indicating that heterozygotes had lower fitness than their homozygote counterparts. Marker location could not explain the observed HFCs. Evidence of inbreeding, outbreeding or population and family structuring was not found in this work. The presence of undetected general effects that may lead to the appearance of HFCs, however, cannot be ruled out. These results seem to be best explained by the occurrence of local effects (due to linkage) or even by possible direct locus advantages.     

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.     

Heterozygosity–fitness correlations and their relevance to studies on inbreeding depression in threatened species

The majority of reported multilocus heterozygosity–fitness correlations (HFCs) are from large, outbred populations, and their relevance to studies on inbreeding depression in threatened populations is often stressed. The results of such HFC studies conducted on outbred populations may be of limited application to threatened population management, however, as bottlenecked populations exhibit increased incidence of inbreeding, increased linkage disequilibrium, reduced genetic diversity and possible effects of historical inbreeding such as purging. These differences may affect both our ability to detect inbreeding depression in threatened species, and our interpretation of the underlying mechanisms for observed heterozygosity–fitness relationships. The study of HFCs in outbred populations is of interest in itself, but the results may not translate directly to threatened populations that have undergone severe bottlenecks.     

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

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

The effects of multilocus balancing selection on neutral variability

We studied the effect of multilocus balancing selection on neutral nucleotide variability at linked sites by simulating a model where diallelic polymorphisms are maintained at an arbitrary number of selected loci by means of symmetric overdominance. Different combinations of alleles define different genetic backgrounds that subdivide the population and strongly affect variability. Several multilocus fitness regimes with different degrees of epistasis and gametic disequilibrium are allowed. Analytical results based on a multilocus extension of the structured coalescent predict that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. Our simulation results show that although variability increases with the number of genetic backgrounds that are maintained in the population, it is reduced by random fluctuations in the frequencies of those backgrounds and does not reach high levels even in very large populations. We also show that previous results on balancing selection in single-locus systems do not extend to the multilocus scenario in a straightforward way. Different patterns of linkage disequilibrium and of the frequency spectrum of neutral mutations are expected under different degrees of epistasis. Interestingly, the power to detect balancing selection using deviations from a neutral distribution of allele frequencies seems to be diminished under the fitness regime that leads to the largest increase of variability over the neutral case. This and other results are discussed in the light of data from the Mhc.     

Growth rate correlates to individual heterozygosity in the European eel, Anguilla anguilla L

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme over-dominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.     

HETEROZYGOSITY‐FITNESS CORRELATIONS: A TIME FOR REAPPRAISAL

Owing to the remarkable progress of molecular techniques, heterozygosity‐fitness correlations (HFCs) have become a popular tool to study the impact of inbreeding in natural populations. However, their underlying mechanisms are often hotly debated. Here we argue that these “debates” rely on verbal arguments with no basis in existing theory and inappropriate statistical testing, and that it is time to reconcile HFC with its historical and theoretical fundaments. We show that available data are quantitatively and qualitatively consistent with inbreeding‐based theory. HFC can be used to estimate the impact of inbreeding in populations, although such estimates are bound to be imprecise, especially when inbreeding is weak. Contrary to common belief, linkage disequilibrium is not an alternative to inbreeding, but rather comes with some forms of inbreeding, and is not restricted to closely linked loci. Finally, the contribution of local chromosomal effects to HFC, while predicted by inbreeding theory, is expected to be small, and has rarely if ever proven statistically significant using adequate tests. We provide guidelines to safely interpret and quantify HFCs, and present how HFCs can be used to quantify inbreeding load and unravel the structure of natural populations.     

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

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

Does heterozygosity estimate inbreeding in real populations?

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

THE EVOLUTION OF STRONG REPRODUCTIVE ISOLATION

Felsenstein distinguished two ways by which selection can directly strengthen isolation. First, a modifier that strengthens prezygotic isolation can be favored everywhere. This fits with the traditional view of reinforcement as an adaptation to reduce deleterious hybridization by strengthening assortative mating. Second, selection can favor association between different incompatibilities, despite recombination. We generalize this "two allele" model to follow associations among any number of incompatibilities, which may include both assortment and hybrid inviability. Our key argument is that this process, of coupling between incompatibilities, may be quite different from the usual view of reinforcement: strong isolation can evolve through the coupling of any kind of incompatibility, whether prezygotic or postzygotic. Single locus incompatibilities become coupled because associations between them increase the variance in compatibility, which in turn increases mean fitness if there is positive epistasis. Multiple incompatibilities, each maintained by epistasis, can become coupled in the same way. In contrast, a single-locus incompatibility can become coupled with loci that reduce the viability of haploid hybrids because this reduces harmful recombination. We obtain simple approximations for the limits of tight linkage, and strong assortment, and show how assortment alleles can invade through associations with other components of reproductive isolation.     

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

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