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.     

ON THE LOW HERITABILITY OF LIFE-HISTORY TRAITS

Life-history traits such as longevity and fecundity often show low heritability. This is usually interpreted in terms of Fisher's fundamental theorem to mean that populations are near evolutionary equilibrium and genetic variance in total fitness is low. We develop the causal relationship between metric traits and life-history traits to show that a life-history trait is expected to have a low heritability whether or not the population is at equilibrium. This is because it is subject to all the environmental variation in the metric traits that affect it plus additional environmental variation. There is no simple prediction regarding levels of additive genetic variance in life-history traits, which may be high at equilibrium. Several other patterns in the inheritance of life-history traits are readily predicted from the causal model. These include the strength of genetic correlations between life-history traits, levels of nonadditive genetic variance, and the inevitability of genotype-environment interaction.     

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.     

Microsatellite measures of inbreeding: a meta-analysis

Abstract Meta-analyses of published and unpublished correlations between phenotypic variation and two measures of genetic variation at microsatellite loci, multilocus heterozygosity (MLH) and mean d², revealed that the strength of these associations are generally weak (mean r < 0.10). Effects on life-history trait variation were significantly greater than zero for both measures over all reported effect sizes ( r = 0. 0856 and 0.0479 for MLH and mean d ², respectively), whereas effects on morphometric traits were not ( r = 0.0052 and r = 0.0038), which is consistent with the prediction that life-history traits exhibit greater inbreeding depression than morphometric traits. Effect sizes reported using mean d ² were smaller and more variable than those reported using MLH, suggesting that MLH may be a better metric for capturing inbreeding depression most of the time. However, analyses of paired effect sizes reported using both measures from the same data did not differ significantly. Several lines of evidence suggest that published effects sizes are upwardly biased. First, effect sizes from published studies were significantly higher than those reported in unpublished studies. Second, fail-safe numbers for reported effect sizes were generally quite low, with the exception of correlations between MLH and life-history traits. Finally, the slope of the regression of effect size on sample size was negative for most sets of traits. Taken together, these results suggest that studies designed to detect inbreeding depression on a life-history trait using microsatellites will need to sample in excess of 600 individuals to detect an average effect size ( r = 0.10) with reasonable statistical power (0.80). Very few published studies have used samples sizes approaching this value.     

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.     

Genetic correlations do not constrain the evolution of sexual dimorphism in the moss Ceratodon purpureus

The trajectory of phenotypic evolution is constrained in the short term by genetic correlations among traits. However, the extent to which genetic correlations impose a lasting constraint is generally unknown. Here, I examine the genetic architecture of life-history variation in male and female gametophytes from two populations of the moss Ceratodon purpureus, focusing on genetic correlations within and between the sexes. A significant negative correlation between allocation to vegetative and reproductive tissue was evident in males of both populations, but not females. All traits showed between-sex correlations of significantly less than one, indicating additive genetic variance for sexual dimorphism. The degree of dimorphism for traits was significantly negatively associated with the strength of the between-sex correlation. The structure of genetic correlations among life-history traits was more divergent between the two populations in females than in males. Collectively, these results suggest that genetic correlations do not impose a lasting constraint on the evolution of life-history variation in the species.     

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.     

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

Genetic analysis of life-history constraint and evolution in a wild ungulate population

Trade-offs among life-history traits are central to evolutionary theory. In quantitative genetic terms, trade-offs may be manifested as negative genetic covariances relative to the direction of selection on phenotypic traits. Although the expression and selection of ecologically important phenotypic variation are fundamentally multivariate phenomena, the in situ quantification of genetic covariances is challenging. Even for life-history traits, where well-developed theory exists with which to relate phenotypic variation to fitness variation, little evidence exists from in situ studies that negative genetic covariances are an important aspect of the genetic architecture of life-history traits. In fact, the majority of reported estimates of genetic covariances among life-history traits are positive. Here we apply theory of the genetics and selection of life histories in organisms with complex life cycles to provide a framework for quantifying the contribution of multivariate genetically based relationships among traits to evolutionary constraint. We use a Bayesian framework to link pedigree-based inference of the genetic basis of variation in life-history traits to evolutionary demography theory regarding how life histories are selected. Our results suggest that genetic covariances may be acting to constrain the evolution of female life-history traits in a wild population of red deer Cervus elaphus: genetic covariances are estimated to reduce the rate of adaptation by about 40%, relative to predicted evolutionary change in the absence of genetic covariances. Furthermore, multivariate phenotypic (rather than genetic) relationships among female life-history traits do not reveal this constraint.     

Evaluating the role of inbreeding depression in heterozygosity‐fitness correlations: how useful are tests for identity disequilibrium?

Heterozygosity‐fitness correlations (HFCs) have been observed for several decades, but their causes are often elusive. Tests for identity disequilibrium (ID, correlated heterozygosity between loci) are commonly used to determine if inbreeding depression is a possible cause of HFCs. We used computer simulations to determine how often ID is detected when HFCs are caused by inbreeding depression. We also used ID in conjunction with HFCs to estimate the proportion of variation (r²) in fitness explained by the individual inbreeding coefficient (F). ID was not detected in a large proportion of populations with statistically significant HFCs (sample size = 120 individuals) unless the variance of F was high (σ²(F) ≥ 0.005) or many loci were used (100 microsatellites or 1000 SNPs). For example, with 25 microsatellites, ID was not detected in 49% of populations when HFCs were caused by six lethal equivalents and σ²(F) was typical of vertebrate populations (σ²(F) ≈ 0.002). Estimates of r² between survival and F based on ID and HFCs were imprecise unless ID was strong and highly statistically significant (P ≈ 0.01). These results suggest that failing to detect ID in HFC studies should not be taken as evidence that inbreeding depression is absent. The number of markers necessary to simultaneously detect HFC and ID depends strongly on σ²(F). Thus the mating system and demography of populations, which influence σ²(F), should be considered when designing HFC studies. ID should be used in conjunction with HFCs to estimate the correlation between fitness and F, because HFCs alone reveal little about the strength of inbreeding depression.     

QUANTITATIVE GENETICS OF SURVIVAL AND GROWTH IN IMPATIENS CAPENSIS

When variation in life-history characters is caused by many genes of small effect, then quantitative-genetic parameters may quantify constraints on rate and direction of microevolutionary change. I estimated heritabilities and genetic correlations for 16 life-history and morphological characters in two populations of Impatiens capensis, a partially self-pollinating herbaceous annual. The Madison population had little or no additive genetic variance for any of these characters, while the Milwaukee population had significant narrowsense heritabilities and genetic correlations for several traits, including adult size, which is highly correlated with fitness. All genetic correlations among fitness components were positive, hence there is no evidence for antagonistic pleiotropy among these traits. Dissimilarity of heritabilities in the two populations supports theoretical predictions that long-term changes in genetic variance-covariance patterns may occur when population sizes are small and selection is strong, as may occur in many plant species.     

The change in quantitative genetic variation with inbreeding

Inbreeding is known to reduce heterozygosity of neutral genetic markers, but its impact on quantitative genetic variation is debated. Theory predicts a linear decline in additive genetic variance (V(A)) with increasing inbreeding coefficient (F) when loci underlying the trait act additively, but a nonlinear hump-shaped relationship when dominance and epistasis are important. Predictions for heritability (h2) are similar, although the exact shape depends on the value of h2 in the absence of inbreeding. We located 22 published studies in which the level of genetic variation in experimentally inbred populations (measured by V(A) or h2) was compared with that in outbred control populations. For life-history traits, the data strongly supported a nonlinear change in genetic variation with increasing F. V(A) and h2 were, respectively, 244% and 50% higher at F = 0.4 than in outbred populations, and dominance plus epistatic variance together exceeded additive variance by a factor of four. For nonfitness traits the decline was linear and estimates of nonadditive variance were small. These results confirm that population bottlenecks frequently increase V(A) in some traits, and imply that life-history traits are underlain by substantial dominance or epistasis. However, the importance of drift-induced genetic variation in conservation or evolutionary biology is questionable, in part because inbreeding depression usually accompanies inbreeding.     

THE EVOLUTIONARY GENETICS OF MALE LIFE-HISTORY CHARACTERS IN DROSOPHILA MELANOGASTER

Alternative models of the maintenance of genetic variability, theories of life-history evolution, and theories of sexual selection and mate choice can be tested by measuring additive and nonadditive genetic variances of components of fitness. A quantitative genetic breeding design was used to produce estimates of genetic variances for male life-history traits in Drosophila melanogaster. Additive genetic covariances and correlations between traits were also estimated. Flies from a large, outbred, laboratory population were assayed for age-specific competitive mating ability, age-specific survivorship, body mass, and fertility. Variance-component analysis then allowed the decomposition of phenotypic variation into components associated with additive genetic, nonadditive genetic, and environmental variability. A comparison of dominance and additive components of genetic variation provides little support for an important role for balancing selection in maintaining genetic variance in this suite of traits. The results provide support for the mutation-accumulation theory, but not the antagonistic-pleiotropy theory of senescence. No evidence is found for the positive genetic correlations between mating success and offspring quality or quantity that are predicted by “good genes” models of sexual selection. Additive genetic coefficients of variation for life-history characters are larger than those for body weight. Finally, this set of male life-history characters exhibits a very low correspondence between estimates of genetic and phenotypic correlations.     

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.     

Evidence for Fine-Scale Natal Homing Among Island Beach Spawning Sockeye Salmon, Oncorhynchus nerka

Salmonid fishes aggregate for breeding at spatially defined, suitable habitats. These aggregations may evolve into discrete populations when precise natal homing leads to reproductive isolation, and local regimes of selection lead to adaptation. Population structure is often defined by persistent differences in selectively neutral genetic markers and in mean values of morphological and life-history traits between locations. This approach is limited by the spatial scale at which traits diverge; low levels of reproductively successful straying, combined with similar selective pressures on life-history traits resulting from similar habitat features and environmental conditions, can significantly reduce the power of these discriminatory methods. We compared data on three life-history traits and polymorphism of DNA microsatellites for evidence of population subdivision among sockeye salmon spawning on spatially discrete but physically similar beaches on islands in Iliamna Lake, Alaska. We found small but significant differences in average body length, body depth and age composition between sites as well as significant interactions between site and year. These interactions, reflecting random variation in growth or recruitment among sites, are a powerful tool for discriminating populations with similar mean trait values. These results suggest fine-scale homing to natal sites, but the microsatellite data revealed no evidence of restricted gene flow among sites. There seems to be enough straying among the populations to prevent differentiation at neutral traits but enough homing for them to be functionally distinct.     

Evolution of pleiotropic alleles for maturation and size as a consequence of predation

Understanding life-history evolution requires knowledge about genetic interactions, physiological mechanisms and the nature of selection. For platyfish, Xiphophorus maculatus, extensive information is available about genetic and physiological mechanisms influencing life-history traits. In particular, alleles at the pituitary locus have large and antagonistic effects on age and size at sexual maturation. To examine how predation affects the evolution of these antagonistic traits, I examined pituitary allele evolution in experimental populations differing in predation risk. A smaller size, earlier maturation allele increased in frequency when predators were absent, while a larger size, later maturation allele increased in frequency when predators were present. Thus, predation favours alleles for larger size, at the cost of later maturation and reproduction. These findings are interesting for several reasons. First, predation is often predicted to favour early reproduction at smaller sizes. Second, few studies have shown how selection acts on alleles that affect age and size at sexual maturation. Finally, many studies assume that trade-offs between these life-history traits result from antagonistic pleiotropy, with alleles that positively affect one trait negatively affecting another, yet this is rarely known. This study unequivocally demonstrates that genetically based trade-offs affect life-history evolution in platyfish.     

The effect of drought stress on heterozygosity–fitness correlations in pedunculate oak (Quercus robur)

Background and Aims

The interaction between forest fragmentation and predicted climate change may pose a serious threat to tree populations. In small and spatially isolated forest fragments, increased homozygosity may directly affect individual tree fitness through the expression of deleterious alleles. Climate change-induced drought stress may exacerbate these detrimental genetic consequences of forest fragmentation, as the fitness response to low levels of individual heterozygosity is generally thought to be stronger under environmental stress than under optimal conditions.

Methods

To test this hypothesis, a greenhouse experiment was performed in which various transpiration and growth traits of 6-month-old seedlings of Quercus robur differing in multilocus heterozygosity (MLH) were recorded for 3 months under a well-watered and a drought stress treatment. Heterozygosity–fitness correlations (HFC) were examined by correlating the recorded traits of individual seedlings to their MLH and by studying their response to drought stress.

Key Results

Weak, but significant, effects of MLH on several fitness traits were obtained, which were stronger for transpiration variables than for the recorded growth traits. High atmospheric stress (measured as vapour pressure deficit) influenced the strength of the HFCs of the transpiration variables, whereas only a limited effect of the irrigation treatment on the HFCs was observed.

Conclusions

Under ongoing climate change, increased atmospheric stress in the future may strengthen the negative fitness responses of trees to low MLH. This indicates the necessity to maximize individual multilocus heterozygosity in forest tree breeding programmes.     

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.     

Insular shifts and trade-offs in life-history traits in pond frogs in the Zhoushan Archipelago, China

Island and mainland populations of animal species often differ strikingly in life-history traits such as clutch size, egg size, total reproductive effort and body size. However, despite widespread recognition of insular shifts in these life-history traits in birds, mammals and reptiles, there have been no reports of such life-history shifts in amphibians. Furthermore, most studies have focused on one specific life-history trait without explicit consideration of coordinated evolution among these intimately linked life-history traits, and thus the relationships among these traits are poorly studied. Here we provide the first evidence of insular shifts and trade-offs in a coordinated suite of life-history traits for an amphibian species, the pond frog Rana nigromaculata . Life-history data were collected from eight islands in the Zhoushan Archipelago and neighboring mainland China. We found consistent, significant shifts in all life-history traits between mainland and island populations. Island populations had smaller clutch sizes, larger egg sizes, larger female body size and invested less in total reproductive effort than mainland populations. Significant negative relationships were found between egg size and clutch size and between egg size and total reproductive effort among frog populations after controlling for the effects of body size. Therefore, decreased reproductive effort and clutch size, larger egg size and body size in pond frogs on islands were selected through trade-offs as an overall life-history strategy. Our findings contribute to the formation of a broad, repeatable ecological generality for insular shifts in life-history traits across a range of terrestrial vertebrate taxa.