Heterozygosity-fitness Correlations Research Articles

Linking pathogen-microbiome-host interactions to explain amphibian population dynamics.

Symbiotic interactions can determine the evolutionary trajectories of host species, influencing genetic variation through selection and changes in demography. In the context of strong selective pressures such as those imposed by infectious diseases, symbionts providing defences could contribute to increase host fitness upon pathogen emergence. Here, we generated genome-wide data of an amphibian species to find evidence of evolutionary pressures driven by two skin symbionts: a batrachochytrid fungal pathogen and an antifungal bacterium. Using demographic modelling, we found evidence of decreased effective population size, probably due to pathogen infections. Additionally, we investigated host genetic associations with infection status, antifungal bacterium abundance and overall microbiome diversity using structural equation models. We uncovered relatively lower nucleotide diversity in infected frogs and potential heterozygote advantage to recruit the candidate beneficial symbiont and fight infections. Our models indicate that environmental conditions have indirect effects on symbiont abundance through both host body traits and microbiome diversity. Likewise, we uncovered a potential offsetting effect among host heterozygosity-fitness correlations, plausibly pointing to different ecological and evolutionary processes among the three species due to dynamic interactions. Our findings revealed that evolutionary pressures not only arise from the pathogen but also from the candidate beneficial symbiont, and both interactions shape the genetics of the host. Our results advance knowledge about multipartite symbiotic relationships and provide a framework to model ecological and evolutionary dynamics in wild populations. Finally, our study approach can be applied to inform conservation actions such as bioaugmentation strategies for other imperilled amphibians affected by infectious diseases.

Parapatric Genetic Lineages Persist in a Multiply Introduced Non-native Bush-Cricket

To understand colonization success of an invasive species we need to know the origin of the founders, where and when they were introduced, and how they spread from the introduction site(s) through the landscape. Admixture of different genetic lineages from multiple introductions is generally hypothesized to be beneficial to invasive species thanks to adaptive variation and heterozygosity-fitness correlations. In this study, population genetic and landscape data was gathered for Roesel’s bush-cricket,Roeseliana roeseliia small bush-cricket common in central and eastern Europe that currently is expanding its range in northern Europe. We examined how colonization history and landscape structure affect the spread of the species and its population genetic structure, as a consequence of multiple introductions. Using comprehensive information of the species ecology and dispersal, together with genetic structure inferred from samples from 29 locations in central Sweden (we employed data published by Preuss et al., 2015), we found that two parapatric founding lineages have coexisted with very little gene flow during a long time span. An isolation-by-distance pattern and a decrease of genetic diversity toward marginal areas were more pronounced in the lineage situated in forest dominated landscapes. Our findings are in strong contrast to the hypothesis that different genetic lineages will admix when introduced to the same area. The presence of the separate lineages decades after introduction and without physical barriers for gene flow shows that some mechanism prevents them from admixture. One possibility is that the lineages with different genetic setups have adapted independently to local conditions and their admixture resulted in loss of locally adapted genotypes and hybrid offspring, less viable than the respective ancestral genotypes. However, an alternative post-mating reproductive barrier and hybrid breakdown phenomenon should also be considered. Our data indicate that besides landscape characteristics, human transportation of agricultural goods may play an important role for the overall spatial genetic pattern of the species in the study area by aiding the spread of the species.

Local selection signals in the genome of blue tits emphasize regulatory and neuronal evolution.

Understanding the genomic landscape of adaptation is central to understanding microevolution in wild populations. Genomic targets of selection and the underlying genomic mechanisms of adaptation can be elucidated by genome-wide scans for past selective sweeps or by scans for direct fitness associations. We sequenced and assembled 150haplotypes of 75 blue tits (Cyanistes caeruleus) of a single Central European population by a linked-read technology. We used these genome data in combination with coalescent simulations (i) to estimate an historical effective population size of ~250,000, which recently declined to ~10,000, and (ii) to identify genome-wide distributed selective sweeps of beneficial variants probably originating from standing genetic variation (soft sweeps). The genes linked to these soft sweeps, but also those linked to hard sweeps based on new beneficial mutants, showed a significant enrichment for functions associated with gene expression and transcription regulation. This emphasizes the importance of regulatory evolution in the population's adaptive history. Soft sweeps were further enriched for genes related to axon and synapse development, indicating the significance of neuronal connectivity changes in the brain potentially linked to behavioural adaptations. A previous scan of heterozygosity-fitness correlations revealed a consistent negative effect on arrival date at the breeding site for a single microsatellite in the MDGA2gene. Here, we used the haplotype structure around this microsatellite to explain the effect as a local and direct outbreeding effect of a gene involved in synapse development.

Bioindicator snake shows genomic signatures of natural and anthropogenic barriers to gene flow.

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant 'islands'. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which 'islands' of habitat are most in need of management and protection.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population.

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.

Heterozygosity-Fitness Correlations in a Continental Island Population of Thorn-Tailed Rayadito.

Heterozygosity-fitness correlations (HFCs) have been used to monitor the effects of inbreeding in threatened populations. HFCs can also be useful to investigate the potential effects of inbreeding in isolated relict populations of long-term persistence and to better understand the role of inbreeding and outbreeding as drivers of changes in genetic diversity. We studied a continental island population of thorn-tailed rayadito (Aphrastura spinicauda) inhabiting the relict forest of Fray Jorge National Park, north-central Chile. This population has experienced a long-term, gradual process of isolation since the end of the Tertiary. Using 10 years of field data in combination with molecular techniques, we tested for HFCs to assess the importance of inbreeding depression. If inbreeding depression is important, we predict a positive relationship between individual heterozygosity and fitness-related traits. We genotyped 183 individuals at 12 polymorphic microsatellite loci and used 7 measures of reproductive success and estimates of apparent survival to calculate HFCs. We found weak to moderate statistical support (P-values between 0.05 and 0.01) for a linear effect of female multi-locus heterozygosity (MLH) on clutch size and nonlinear effects on laying date and fledging success. While more heterozygous females laid smaller clutches, nonlinear effects indicated that females with intermediate values of MLH started laying earlier and had higher fledging success. We found no evidence for effects of MLH on annual fecundity or on apparent survival. Our results along with the long-term demographic stability of the study population contradict the hypothesis that inbreeding depression occurs in this population.

Negative effects of individual heterozygosity on reproductive success in a wild bird population.

The evolutionary consequences of individual genetic diversity are frequently studied by assessing heterozygosity-fitness correlations (HFCs). The prevalence of positive and negative HFCs and the predominance of general versus local effects in wild populations are far from understood, partly because comprehensive studies testing for both inbreeding and outbreeding depression are lacking. We studied a genetically diverse population of blue tits in southern Germany using a genome-wide set of 87 microsatellites to investigate the relationship between proxies of reproductive success and measures of multilocus and single-locus individual heterozygosity (MLH and SLH). We used complimentary measures of MLH and partitioned markers into functional categories according to their position in the blue tit genome. HFCs based on MLH were consistently negative for functional loci, whereas correlations were rather inconsistent for loci found in nonfunctional areas of the genome. Clutch size was the only reproductive variable showing a general effect. We found evidence for local effects for three measures of reproductive success: arrival date at the breeding site, the probability of breeding at the study site and male reproductive success. For these, we observed consistent, and relatively strong, negative effects at one functional locus. Remarkably, this marker had a similar effect in another blue tit population from Austria (~400km to the east). We suggest that a genetic local effect on timing of arrival might be responsible for most negative HFCs detected, with carry-over effects on other reproductive traits. This effect could reflect individual differences in the distance between overwintering areas and breeding sites.

Investigating inbreeding in a free-ranging, captive population of an Australian marsupial

Conservation breeding programs, such as fenced reserves, conserve species by removing populations from key threatening processes and providing supplementary resources to support species preservation. This form of management often has the eventual aim to reintroduce individuals back to their native habitat. A key issue with captive management is the small breeding pool, due to generally small population sizes and limited migration, which may lead to inbreeding and consequent inbreeding depression. We investigated inbreeding and its potential fitness consequences in a large, free-ranging, captive population of an Australian marsupial, the long-nosed potoroo (Potorous tridactylus), a species that has been recently impacted by severe bushfire. We used reduced-representation sequencing to genotype single nucleotide polymorphisms and reconstruct a pedigree for the population. We detected low levels of inbreeding in our partial pedigree that incorporated 50.7% (34/67) of sampled individuals. Inbreeding depression was investigated using heterozygosity–fitness correlations of six fitness measurements, but we found no evidence of inbreeding depression in the population despite strong molecular evidence for inter-individual variation in inbreeding coefficient. Our findings suggest that the potoroo population at Tidbinbilla Nature Reserve (TNR) maintains a high level of diversity. As the population is descendant of animals from divergent sources, the TNR population also helps conservationists to learn more about potential consequences of admixture, and plan recovery strategies for the species. For TNR potoroos, we suggest that current management protocols have been effective at maintaining diversity to date. More generally, our results show how modern molecular techniques can efficiently characterise the genetic profile of a free-ranging, captive population to help inform management guidelines for wild recovery.

PSVIII-22 Microsatellite-based heterozygosity fitness correlations in reindeer

Abstract Inbred individuals often show reduced fitness due to increased number of loci identical by descent as the loss of heterozygosity may lead to the unmasking of deleterious recessive alleles. Inbreeding level can be measured either based on pedigree or genetic markers. The latter can be used for heterozygosity fitness correlation (HFC) studies in populations which pedigrees are unknown. Here, we estimated individual heterozygosity level in reindeer based on ten STR and investigated heterozygosity fitness correlations. A total of 64 individuals were divided into two groups: males (M, n = 16) and females (F, n = 48). The influence of inbreeding on height at withers (HaW), chest depth (CD), chest width (CW), chest girth (CG), wrist girth (WG), body length(BL), head length (HL), body weight (BW), hip width (HW) and loin width (LW) was examined. R package “inbreedR” was used for statistical analysis. Identity disequilibrium point estimate - g2 was calculated to estimate the correlations between inbreeding, multilocus heterozygosity (MLH) and fitness. For M group the measure of g2 was 0.03390±0.02564, p (g2 &amp;gt;0) = 0.017, and for F group - 0.03973±0.02377, p (g2 &amp;gt; 0) = 0.006. Expected r2 between inbreeding level and heterozygosity was 0.49795 and 0.40812 for M and F, respectively. The highest HFCs in M group were detected for CG - 0.35, WG – 0.26, BW – 0.18 and HW – 0.14 i.e. traits that are characteristics for body composition. The highest HFCs in F group were for LW – 0.22, CW – 0.16 and HaW – 0.15 i.e. traits related to reproductivity. For all the other traits HFCs were close to zero. Our study demonstrated that STR can be used for calculations of MLH and inbreeding correlations as well as for detecting traits that are influenced by inbreeding in reindeer. The study was supported by Russian Science Foundation (Pr. no.16-16-10068).

Heterozygosity fitness correlations and generation interval of the Norway lobster in the Aegean Sea, eastern Mediterranean

BackgroundComprehensively detailed information on population dynamics for benthic species is crucial since potential admixture of individuals could shift the genetic subdivision and age structure during a full breeding period. The apparent genetic impact of the potential recruitment strategy of Norway lobster Nephrops norvegicus is still under research. For this reason the present study was focused on genetic variation of the species over a given continuous year period in a semi-enclosed gulf of the Aegean Sea.ResultsAnalyses revealed that the relative smaller size class in females and the apparent faster growth of males may represent a key-role differential strategy for the two sexes, whereas females tend to mature slower. Heterozygosity fitness correlations (HFCs) showed substantially significant associations suggesting that inbreeding depression for females and outbreeding depression for males are the proximate fitness mechanisms, respectively.ConclusionsNephrops norvegicus uniformal genetic composition (background of high gene flow), could be attributed to potential population recolonization, due to a hypothesized passive larval movement from deeper waters, which may suggest that some offspring of local residents and potential male non-breeders from other regions admixture randomly.

Heterozygosity-Fitness Correlations Reveal Inbreeding Depression in Neonatal Body Size in a Critically Endangered Rock Iguana.

Inbreeding depression, though challenging to identify in nature, may play an important role in regulating the dynamics of small and isolated populations. Conversely, greater expression of genetic load can enhance opportunities for natural selection. Conditional expression concentrates these opportunities for selection and may lead to failure of detection. This study investigates the possibility for age-dependent expression of inbreeding depression in a critically endangered population of rock iguanas, Cyclura nubila caymanensis. We employ heterozygote-fitness correlations to examine the contributions of individual genetic factors to body size, a fitness-related trait. Nonsignificant reductions in homozygosity (up to 7%) were detected between neonates and individuals surviving past their first year, which may reflect natural absorption of inbreeding effects by this small, fecund population. The majority of variation in neonate body size was attributed to maternal or environmental effects (i.e., clutch identity and incubation length); however, heterozygosity across 22 microsatellite loci also contributed significantly and positively to model predictions. Conversely, effects of heterozygosity on fitness were not detectable when adults were examined, suggesting that inbreeding depression in body size may be age dependent in this taxon. Overall, these findings emphasize the importance of taking holistic, cross-generational approaches to genetic monitoring of endangered populations.

Effective population size and heterozygosity-fitness correlations in a population of the Mediterranean lagoon ecotype of long-snouted seahorse Hippocampus guttulatus

The management of endangered species is complicated in the marine environment owing to difficulties to directly access, track and monitor in situ. Population genetics provide a genuine alternative to estimate population size and inbreeding using non-lethal procedures. The long-snouted seahorse, Hippocampus guttulatus, is facing multiple threats such as human disturbance or by-catch, and has been listed in the red list of IUCN. One large population is found in the Thau lagoon, in the south of France. A recent study has shown this population belongs to a genetic lineage only found in Mediterranean lagoons that can be considered as an Evolutionarily Significant Unit (ESU) and should be managed with dedicated conservation strategies. In the present study, we used genetic analysis of temporal samples to estimate the effective population size of the Thau population and correlations between individual multilocus heterozygosity and fitness traits to investigate the possible expression of inbreeding depression in the wild. Non-invasive sampling of 172 seahorses for which profiles were pictured and biometric data recorded were genotyped using 291 informative SNPs. Genetic diversity remained stable over a 7-year time interval. In addition, very low levels of close relatedness and inbreeding were observed, with only a single pair of related individuals in 2008 and two inbreds in 2013. We did not detect departure from identity equilibrium. The effective population size was estimated to be Ne = 2742 (~ 40 reproductive seahorses per km2), larger than previously thought. No correlation was observed between heterozygosity and fluctuating asymmetry or other morphometric traits, suggesting a population with low variance in inbreeding. Together these results suggest this population does not meet conventional genetic criteria of an endangered population, as the population seems sufficiently large to avoid inbreeding and its detrimental effects. This study paves the way for the genetic monitoring of this recently discovered ESU of a species with patrimonial and conservation concerns.

Heterozygosity at neutral and immune loci is not associated with neonatal mortality due to microbial infection in Antarctic fur seals.

Numerous studies have reported correlations between the heterozygosity of genetic markers and fitness. These heterozygosity–fitness correlations (HFCs) play a central role in evolutionary and conservation biology, yet their mechanistic basis remains open to debate. For example, fitness associations have been widely reported at both neutral and functional loci, yet few studies have directly compared the two, making it difficult to gauge the relative contributions of genome‐wide inbreeding and specific functional genes to fitness. Here, we compared the effects of neutral and immune gene heterozygosity on death from bacterial infection in Antarctic fur seal (Arctocephalus gazella) pups. We specifically developed a panel of 13 microsatellites from expressed immune genes and genotyped these together with 48 neutral loci in 234 individuals, comprising 39 pups that were classified at necropsy as having most likely died of bacterial infection together with a five times larger matched sample of healthy surviving pups. Identity disequilibrium quantified from the neutral markers was positive and significant, indicative of variance in inbreeding within the study population. However, multilocus heterozygosity did not differ significantly between healthy and infected pups at either class of marker, and little evidence was found for fitness associations at individual loci. These results support a previous study of Antarctic fur seals that found no effects of heterozygosity at nine neutral microsatellites on neonatal survival and thereby help to refine our understanding of how HFCs vary across the life cycle. Given that nonsignificant HFCs are underreported in the literature, we also hope that our study will contribute toward a more balanced understanding of the wider importance of this phenomenon.

Sibling competition in Blue Tits (Cyanistes caeruleus) is stronger in homozygous broods

Abstract Heterozygosity affects mate selection and can modulate interactions among family members and their fitness-related decisions. We studied whether nestling heterozygosity affected parent–offspring interactions and sib–sib competition in the Blue Tit (Cyanistes caeruleus) while controlling for the degree of relatedness among nestlings. Demanding environmental conditions might make the detection of heterozygosity-fitness correlations easier. Thus, we also investigated whether the decision rules of family members according to offspring heterozygosity were affected by brood size, as a proxy of the strength of sibling conflict. We found that chick individual heterozygosity was positively although weakly associated with individual body mass. Mean brood heterozygosity did not predict fledging success, but broods that fledged more chicks showed a higher number of less common alleles. Interestingly, fathers, but not mothers, favored heterozygous broods with many nestlings, that is, heterozygous broods with higher potential for sibling conflict. Moreover, the lower the mean brood heterozygosity the stronger the begging intensity when parents were absent, regardless of brood size. Finally, the degree of relatedness among nestlings was not associated with any behavioral parameter, supporting a more prominent role for heterozygosity in shaping intra-family interactions. Our findings suggest that offspring heterozygosity determines sex-specific rules of parental care and that genetic diversity is associated with lower sibling competition.

Contrasting heterozygosity-fitness correlations across life in a long-lived seabird.

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

Genetic signatures of small effective population sizes and demographic declines in an endangered rattlesnake, Sistrurus catenatus.

Endangered species that exist in small isolated populations are at elevated risk of losing adaptive variation due to genetic drift. Analyses that estimate short‐term effective population sizes, characterize historical demographic processes, and project the trajectory of genetic variation into the future are useful for predicting how levels of genetic diversity may change. Here, we use data from two independent types of genetic markers (single nucleotide polymorphisms [SNPs] and microsatellites) to evaluate genetic diversity in 17 populations spanning the geographic range of the endangered eastern massasauga rattlesnake (Sistrurus catenatus). First, we use SNP data to confirm previous reports that these populations exhibit high levels of genetic structure (overall Fst = 0.25). Second, we show that most populations have contemporary Ne estimates <50. Heterozygosity–fitness correlations in these populations provided no evidence for a genetic cost to living in small populations, though these tests may lack power. Third, model‐based demographic analyses of individual populations indicate that all have experienced declines, with the onset of many of these declines occurring over timescales consistent with anthropogenic impacts (<200 years). Finally, forward simulations of the expected loss of variation in relatively large (Ne = 50) and small (Ne = 10) populations indicate they will lose a substantial amount of their current standing neutral variation (63% and 99%, respectively) over the next 100 years. Our results argue that drift has a significant and increasing impact on levels of genetic variation in isolated populations of this snake, and efforts to assess and mitigate associated impacts on adaptive variation should be components of the management of this endangered reptile.

Heterozygosity–behavior and heterozygosity–fitness correlations in a salamander with limited dispersal

AbstractInbreeding and genetic drift can decrease genetic heterozygosity, and this low heterozygosity can depress fitness, resulting in heterozygosity–fitness correlations (HFCs). HFCs are typically small in magnitude, a result often attributed to power of the analyses. Animal behaviors often affect fitness and are often heritable to some degree. We hypothesized that heterozygosity influences behavior, which, in turn, potentially influences fitness. Specifically, in red‐backed salamanders, Plethodon cinereus, which have limited dispersal and the potential for inbreeding, we tested whether heterozygosity, as estimated from six microsatellite loci, affected home range size, juvenile growth, and survival. We found that salamanders with higher heterozygosity had larger home ranges and grew faster, which is indicative of reproductive success. However, we found no effects of heterozygosity on survival. We conclude that, because activity in P. cinereus is tightly linked to food uptake and mass gain, heterozygosity influences growth via effects on foraging behavior. Future research should investigate how the relationship between heterozygosity, behavior, and fitness may be affected or mediated by endocrine or immune systems.

Age-dependent, negative heterozygosity-fitness correlations and local effects in an endangered Caribbean reptile, Iguana delicatissima.

Inbreeding depression can have alarming impacts on threatened species with small population sizes. Assessing inbreeding has therefore become an important focus of conservation research. In this study, heterozygosity–fitness correlations (HFCs) were measured by genotyping 7 loci in 83 adult and 184 hatchling Lesser Antillean Iguanas, Iguana delicatissima, at a communal nesting site in Dominica to assess the role of inbreeding depression on hatchling fitness and recruitment to the adult population in this endangered species. We found insignificant correlations between multilocus heterozygosity and multiple fitness proxies in hatchlings and adults. Further, multilocus heterozygosity did not differ significantly between hatchlings and adults, which suggests that the survivorship of homozygous hatchlings does not differ markedly from that of their heterozygous counterparts. However, genotypes at two individual loci were correlated with hatching date, a finding consistent with the linkage between specific marker loci and segregating deleterious recessive alleles. These results provide only modest evidence that inbreeding depression influences the population dynamics of I. delicatissima on Dominica.

Evidence of opposing fitness effects of parental heterozygosity and relatedness in a critically endangered marine turtle?

How individual genetic variability relates to fitness is important in understanding evolution and the processes affecting populations of conservation concern. Heterozygosity-fitness correlations (HFCs) have been widely used to study this link in wild populations, where key parameters that affect both variability and fitness, such as inbreeding, can be difficult to measure. We used estimates of parental heterozygosity and genetic similarity ('relatedness') derived from 32 microsatellite markers to explore the relationship between genetic variability and fitness in a population of the critically endangered hawksbill turtle, Eretmochelys imbricata. We found no effect of maternal MLH (multilocus heterozygosity) on clutch size or egg success rate, and no single-locus effects. However, we found effects of paternal MLH and parental relatedness on egg success rate that interacted in a way that may result in both positive and negative effects of genetic variability. Multicollinearity in these tests was within safe limits, and null simulations suggested that the effect was not an artefact of using paternal genotypes reconstructed from large samples of offspring. Our results could imply a tension between inbreeding and outbreeding depression in this system, which is biologically feasible in turtles: female-biased natal philopatry may elevate inbreeding risk and local adaptation, and both processes may be disrupted by male-biased dispersal. Although this conclusion should be treated with caution due to a lack of significant identity disequilibrium, our study shows the importance of considering both positive and negative effects when assessing how variation in genetic variability affects fitness in wild systems.

Population demography and heterozygosity-fitness correlations in natural guppy populations: An examination using sexually selected fitness traits.

Heterozygosity-fitness correlations (HFCs) have been examined in a wide diversity of contexts, and the results are often used to infer the role of inbreeding in natural populations. Although population demography, reflected in population-level genetic parameters such as allelic diversity or identity disequilibrium, is expected to play a role in the emergence and detectability of HFCs, direct comparisons of variation in HFCs across many populations of the same species, with different genetic histories, are rare. Here, we examined the relationship between individual microsatellite heterozygosity and a range of sexually selected traits in 660 male guppies from 22 natural populations in Trinidad. Similar to previous studies, observed HFCs were weak overall. However, variation in HFCs among populations was high for some traits (although these variances were not statistically different from zero). Population-level genetic parameters, specifically genetic diversity levels (number of alleles, observed/expected heterozygosity) and measures of identity disequilibrium (g2 and heterozygosity-heterozygosity correlations), were not associated with variation in population-level HFCs. This latter result indicates that these metrics do not necessarily provide a reliable predictor of HFC effect sizes across populations. Importantly, diversity and identity disequilibrium statistics were not correlated, providing empirical evidence that these metrics capture different essential characteristics of populations. A complex genetic architecture likely underpins multiple fitness traits, including those associated with male fitness, which may have reduced our ability to detect HFCs in guppy populations. Further advances in this field would benefit from additional research to determine the demographic contexts in which HFCs are most likely to occur.