An inbreeding model of associative overdominance during a population bottleneck

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

Stresses affect inbreeding depression in complex ways: disentangling stress-specific genetic effects from effects of initial size in plants

The magnitude of inbreeding depression (ID) varies unpredictably among environments. ID often increases in stressful environments suggesting that these expose more deleterious alleles to selection or increase their effects. More simply, ID could increase under conditions that amplify phenotypic variation (CV²), e.g., by accentuating size hierarchies among plants. These mechanisms are difficult to distinguish when stress increases both ID and phenotypic variation. We grew in- and outbred progeny of Mimulus guttatus under six abiotic stress treatments (control, waterlogging, drought, nutrient deficiency, copper addition, and clipping) with and without competition by the grass Poa palustris. ID differed greatly among stress treatments with δ varying from 7% (control) to 61% (waterlogging) but did not consistently increase with stress intensity. Poa competition increased ID under nutrient deficiency but not other stresses. Analyzing effects of initial size on performance of outbred plants suggests that under some conditions (low N, clipping) competition increased ID by amplifying initial size differences. In other cases (e.g., high ID under waterlogging), particular environments amplified the deleterious genetic effects of inbreeding suggesting differential gene expression. Interestingly, conditions that increased the phenotypic variability of inbred progeny regularly increased ID whereas variability among outbred progeny showed no relationship to ID. Our study reconciles the stress- and phenotypic variability hypotheses by demonstrating how specific conditions (rather than stress per se) act to increase ID. Analyzing CV² separately in inbred and outbred progeny while including effects of initial plant size improve our ability to predict how ID and gene expression vary across environments.

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.

Addressing Darwin's dilemma: Can pseudo-overdominance explain persistent inbreeding depression and load?

Darwin spent years investigating the effects of self-fertilization, concluding that "nature abhors perpetual self-fertilization." Given that selection purges inbred populations of strongly deleterious mutations and drift fixes mild mutations, why does inbreeding depression (ID) persist in highly inbred taxa and why do no purely selfing taxa exist? Background selection, associations and interference among loci, and drift within small inbred populations all limit selection while often increasing fixation. These mechanisms help to explain why more inbred populations in most species consistently show more fixed load. This drift load is manifest in the considerable heterosis regularly observed in between-population crosses. Such heterosis results in subsequent high ID, suggesting a mechanism by which small populations could retain variation and inbreeding load. Multiple deleterious recessive mutations linked in repulsion generate pseudo-overdominance. Many tightly linked load loci could generate a balanced segregating load high enough to sustain ID over many generations. Such pseudo-overdominance blocks (or "PODs") are more likely to occur in regions of low recombination. They should also result in clear genetic signatures including genomic hotspots of heterozygosity; distinct haplotypes supporting alleles at intermediate frequency; and high linkage disequilibrium in and around POD regions. Simulation and empirical studies tend to support these predictions. Additional simulations and comparative genomic analyses should explore POD dynamics in greater detail to resolve whether PODs exist in sufficient strength and number to account for why ID and load persist within inbred lineages.

Partial Selfing Can Reduce Genetic Loads While Maintaining Diversity During Experimental Evolution

Partial selfing, whereby self- and cross- fertilization occur in populations at intermediate frequencies, is generally thought to be evolutionarily unstable. Yet, it is found in natural populations. This could be explained if populations with partial selfing are able to reduce genetic loads and the possibility for inbreeding depression while keeping genetic diversity that may be important for future adaptation. To address this hypothesis, we compare the experimental evolution of Caenorhabditis elegans populations under partial selfing, exclusive selfing or predominant outcrossing, while they adapt to osmotically challenging conditions. We find that the ancestral genetic load, as measured by the risk of extinction upon inbreeding by selfing, is maintained as long as outcrossing is the main reproductive mode, but becomes reduced otherwise. Analysis of genome-wide single-nucleotide polymorphisms (SNPs) during experimental evolution and among the inbred lines that survived enforced inbreeding indicates that populations with predominant outcrossing or partial selfing maintained more genetic diversity than expected with neutrality or purifying selection. We discuss the conditions under which this could be explained by the presence of recessive deleterious alleles and/or overdominant loci. Taken together, our observations suggest that populations evolving under partial selfing can gain some of the benefits of eliminating unlinked deleterious recessive alleles and also the benefits of maintaining genetic diversity at partially dominant or overdominant loci that become associated due to variance of inbreeding levels.

Seascape genomics of eastern oyster (Crassostrea virginica) along the Atlantic coast of Canada

Interactions between environmental factors and complex life-history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean F ST = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype-environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.

A high-quality pedigree and genetic markers both reveal inbreeding depression for quality but not survival in a cooperative mammal

Inbreeding depression, the reduced fitness of offspring of closely related parents, is commonplace in both captive and wild populations and has important consequences for conservation and mating system evolution. However, because of the difficulty of collecting pedigree and life-history data from wild populations, relatively few studies have been able to compare inbreeding depression for traits at different points in the life cycle. Moreover, pedigrees give the expected proportion of the genome that is identical by descent (IBD_g ) whereas in theory with enough molecular markers realized IBD_g can be quantified directly. We therefore investigated inbreeding depression for multiple life-history traits in a wild population of banded mongooses using pedigree-based inbreeding coefficients (f_ped ) and standardized multilocus heterozygosity (sMLH) measured at 35-43 microsatellites. Within an information theoretic framework, we evaluated support for either f_ped or sMLH as inbreeding terms and used sequential regression to determine whether the residuals of sMLH on f_ped explain fitness variation above and beyond f_ped . We found no evidence of inbreeding depression for survival, either before or after nutritional independence. By contrast, inbreeding was negatively associated with two quality-related traits, yearling body mass and annual male reproductive success. Yearling body mass was associated with f_ped but not sMLH, while male annual reproductive success was best explained by both f_ped and residual sMLH. Thus, our study not only uncovers variation in the extent to which different traits show inbreeding depression, but also reveals trait-specific differences in the ability of pedigrees and molecular markers to explain fitness variation and suggests that for certain traits, genetic markers may capture variation in realized IBD_g above and beyond the pedigree expectation.

Resolving the Conflict Between Associative Overdominance and Background Selection

In small populations, genetic linkage between a polymorphic neutral locus and loci subject to selection, either against partially recessive mutations or in favor of heterozygotes, may result in an apparent selective advantage to heterozygotes at the neutral locus (associative overdominance) and a retardation of the rate of loss of variability by genetic drift at this locus. In large populations, selection against deleterious mutations has previously been shown to reduce variability at linked neutral loci (background selection). We describe analytical, numerical, and simulation studies that shed light on the conditions under which retardation vs. acceleration of loss of variability occurs at a neutral locus linked to a locus under selection. We consider a finite, randomly mating population initiated from an infinite population in equilibrium at a locus under selection. With mutation and selection, retardation occurs only when S, the product of twice the effective population size and the selection coefficient, is of order 1. With S >> 1, background selection always causes an acceleration of loss of variability. Apparent heterozygote advantage at the neutral locus is, however, always observed when mutations are partially recessive, even if there is an accelerated rate of loss of variability. With heterozygote advantage at the selected locus, loss of variability is nearly always retarded. The results shed light on experiments on the loss of variability at marker loci in laboratory populations and on the results of computer simulations of the effects of multiple selected loci on neutral variability.

Local Heterozygosity Effects on Nestling Growth and Condition in the Great Cormorant

Under inbreeding, heterozygosity at neutral genetic markers is likely to reflect genome-wide heterozygosity and, thus, is expected to correlate with fitness. There is, however, growing evidence that some of heterozygosity-fitness correlations (HFCs) can be explained by ‘local effects’, where noncoding loci are at linkage disequilibrium with functional genes. The aim of this study was to investigate correlations between heterozygosity at seven microsatellite loci and two fitness-related traits, nestling growth rate and nutritional condition, in a recently bottlenecked population of great cormorant Phalacrocorax carbo sinensis. We found that heterozygosity was positively associated with both nestling traits at the between-brood level, but the individual (within-brood) effects of heterozygosity were non-significant. We also found that only one locus per trait was primarily responsible for the significant multi-locus HFCs, suggesting a linkage disequilibrium with non-identified functional loci. The results give support for ‘local effect’ hypothesis, confirming that HFCs may not only be interpreted as evidence of inbreeding and that genetic associations between functional and selectively neutral markers could be much more common in natural populations than previously thought.

On the Consequences of Purging and Linkage on Fitness and Genetic Diversity

Using computer simulation we explore the consequences of linkage on the inbreeding load of an equilibrium population, and on the efficiency of purging and the loss of genetic diversity after a reduction in population size. We find that linkage tends to cause increased inbreeding load due to the build up of coupling groups of (partially) recessive deleterious alleles. It also induces associative overdominance at neutral sites but rarely causes increased neutral genetic diversity in equilibrium populations. After a reduction in population size, linkage can cause some delay both for the expression of the inbreeding load and the corresponding purging. However, reasonable predictions can be obtained for the evolution of fitness under inbreeding and purging by using empirical estimates of the inbreeding depression rate. Purging selection against homozygotes for deleterious alleles affects the population's pedigree. Furthermore, it can slow the loss of genetic diversity compared to that expected from the variance of gametic contributions to the breeding group and even from pedigree inbreeding. Under some conditions, this can lead to a smaller loss of genetic diversity, even below that expected from population size in the absence of selection.

Heterozygosity-fitness correlations in a declining seabird population

Loss of genetic diversity is thought to lead to increased risk of extinction in endangered populations due to decreasing fitness of homozygous individuals. Here, we evaluated the presence of inbreeding depression in a long-lived seabird, the European shag (Phalacrocorax aristotelis), after a severe decline in population size by nearly 70%. During three reproductive seasons, 85 breeders were captured and genotyped at seven microsatellite loci. Nest sites were monitored during the breeding season to estimate reproductive success as the number of chicks surviving to full-size-grown per nest. Captured birds were tagged with a ring with an individual code, and resighting data were collected during 7-year period. We found a strong effect of multilocus heterozygosity on female reproductive performance, and a significant, although weaker, effect on breeder survival. However, our matrix population model suggests that this relatively small effect of genetic diversity on breeder survival may have a profound effect on fitness. This highlights the importance of integrating life history consequences in HFC studies. Importantly, heterozygosity was correlated across loci, suggesting that genomewide effects, rather than single loci, are responsible for the observed HFCs. Overall, the HFCs are a worrying symptom of genetic erosion in this declining population. Many long-lived species are prone to extinction, and future studies should evaluate the magnitude of fitness impact of genetic deterioration on key population parameters, such as survival of breeders.

Estimating relatedness and inbreeding using molecular markers and pedigrees: the effect of demographic history

Estimates of inbreeding and relatedness are commonly calculated using molecular markers, although the accuracy of such estimates has been questioned. As a further complication, in many situations, such estimates are required in populations with reduced genetic diversity, which is likely to affect their accuracy. We investigated the correlation between microsatellite- and pedigree-based coefficients of inbreeding and relatedness in laboratory populations of Drosophila melanogaster that had passed through bottlenecks to manipulate their genetic diversity. We also used simulations to predict expected correlations between marker- and pedigree-based estimates and to investigate the influence of linkage between loci and null alleles. Our empirical data showed lower correlations between marker- and pedigree-based estimates in our control (nonbottleneck) population than were predicted by our simulations or those found in similar studies. Correlations were weaker in bottleneck populations, confirming that extreme reductions in diversity can compromise the ability of molecular estimates to detect recent inbreeding events. However, this result was highly dependent on the strength of the bottleneck and we did not observe or predict any reduction in correlations in our population that went through a relatively severe bottleneck of N = 10 for one generation. Our results are therefore encouraging, as molecular estimates appeared robust to quite severe reductions in genetic diversity. It should also be remembered that pedigree-based estimates may not capture realized identity-by-decent and that marker-based estimates may actually be more useful in certain situations.

Estimating genome-wide heterozygosity: effects of demographic history and marker type

Heterozygosity-fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.

The genetic basis and experimental evolution of inbreeding depression in Caenorhabditis elegans

Determining the genetic basis of inbreeding depression is important for understanding the role of selection in the evolution of mixed breeding systems. Here, we investigate how androdioecy (a breeding system characterized by partial selfing and outcrossing) and dioecy (characterized by obligatory outcrossing) influence the experimental evolution of inbreeding depression in Caenorhabditis elegans. We derived inbred lines from ancestral and evolved populations and found that the dioecious lineages underwent more extinction than androdioecious lineages. For both breeding systems, however, there was selection during inbreeding because the diversity patterns of 337 single-nucleotide polymorphisms (SNPs) among surviving inbred lines deviated from neutral expectations. In parallel, we also followed the evolution of embryo to adult viability, which revealed similar starting levels of inbreeding depression in both breeding systems, but also outbreeding depression. Under androdioecy, diversity at a neutral subset of 134 SNPs correlated well with the viability trajectories, showing that the population genetic structure imposed by partial selfing affected the opportunity for different forms of selection. Our findings suggest that the interplay between the disruptions of coevolved sets of loci by outcrossing, the efficient purging of deleterious recessive alleles with selfing and overdominant selection with outcrossing can help explain mixed breeding systems.

Postcopulatory selection for dissimilar gametes maintains heterozygosity in the endangered North Atlantic right whale

Although small populations are expected to lose genetic diversity through genetic drift and inbreeding, a number of mechanisms exist that could minimize this genetic decline. Examples include mate choice for unrelated mates and fertilization patterns biased toward genetically dissimilar gametes. Both processes have been widely documented, but the long-term implications have received little attention. Here, we combined over 25 years of field data with high-resolution genetic data to assess the long-term impacts of biased fertilization patterns in the endangered North Atlantic right whale. Offspring have higher levels of microsatellite heterozygosity than expected from this gene pool (effect size = 0.326, P < 0.011). This pattern is not due to precopulatory mate choice for genetically dissimilar mates (P < 0.600), but instead results from postcopulatory selection for gametes that are genetically dissimilar (effect size = 0.37, P < 0.003). The long-term implication is that heterozygosity has slowly increased in calves born throughout the study period, as opposed to the slight decline that was expected. Therefore, this mechanism represents a natural means through which small populations can mitigate the loss of genetic diversity over time.

Population Genomics of sub-saharan Drosophila melanogaster: African diversity and non-African admixture

Drosophila melanogaster has played a pivotal role in the development of modern population genetics. However, many basic questions regarding the demographic and adaptive history of this species remain unresolved. We report the genome sequencing of 139 wild-derived strains of D. melanogaster, representing 22 population samples from the sub-Saharan ancestral range of this species, along with one European population. Most genomes were sequenced above 25X depth from haploid embryos. Results indicated a pervasive influence of non-African admixture in many African populations, motivating the development and application of a novel admixture detection method. Admixture proportions varied among populations, with greater admixture in urban locations. Admixture levels also varied across the genome, with localized peaks and valleys suggestive of a non-neutral introgression process. Genomes from the same location differed starkly in ancestry, suggesting that isolation mechanisms may exist within African populations. After removing putatively admixed genomic segments, the greatest genetic diversity was observed in southern Africa (e.g. Zambia), while diversity in other populations was largely consistent with a geographic expansion from this potentially ancestral region. The European population showed different levels of diversity reduction on each chromosome arm, and some African populations displayed chromosome arm-specific diversity reductions. Inversions in the European sample were associated with strong elevations in diversity across chromosome arms. Genomic scans were conducted to identify loci that may represent targets of positive selection within an African population, between African populations, and between European and African populations. A disproportionate number of candidate selective sweep regions were located near genes with varied roles in gene regulation. Outliers for Europe-Africa F(ST) were found to be enriched in genomic regions of locally elevated cosmopolitan admixture, possibly reflecting a role for some of these loci in driving the introgression of non-African alleles into African populations.

Within- and among-population impact of genetic erosion on adult fitness-related traits in the European tree frog Hyla arborea

Assessing in wild populations how fitness is impacted by inbreeding and genetic drift is a major goal for conservation biology. An approach to measure the detrimental effects of inbreeding on fitness is to estimate correlations between molecular variation and phenotypic performances within and among populations. Our study investigated the effect of individual multilocus heterozygosity on body size, body condition and reproductive investment of males (that is, chorus attendance) and females (that is, clutch mass and egg size) in both small fragmented and large non-fragmented populations of European tree frog (Hyla arborea). Because adult size and/or condition and reproductive investment are usually related, genetic erosion may have detrimental effects directly on reproductive investment, and also on individual body size and condition that in turn may affect reproductive investment. We confirmed that the reproductive investment was highly size-dependent for both sexes. Larger females invested more in offspring production, and larger males attended the chorus in the pond more often. Our results did not provide evidence for a decline in body size, condition and reproductive effort with decreased multilocus heterozygosity both within and among populations. We showed that the lack of heterozygosity-fitness correlations within populations probably resulted from low inbreeding levels (inferior to ca. 20% full-sib mating rate), even in the small fragmented populations. The detrimental effects of fixation load were either low in adults or hidden by environmental variation among populations. These findings will be useful to design specific management actions to improve population persistence.

The effect and relative importance of neutral genetic diversity for predicting parasitism varies across parasite taxa

Understanding factors that determine heterogeneity in levels of parasitism across individuals is a major challenge in disease ecology. It is known that genetic makeup plays an important role in infection likelihood, but the mechanism remains unclear as does its relative importance when compared to other factors. We analyzed relationships between genetic diversity and macroparasites in outbred, free-ranging populations of raccoons (Procyon lotor). We measured heterozygosity at 14 microsatellite loci and modeled the effects of both multi-locus and single-locus heterozygosity on parasitism using an information theoretic approach and including non-genetic factors that are known to influence the likelihood of parasitism. The association of genetic diversity and parasitism, as well as the relative importance of genetic diversity, differed by parasitic group. Endoparasite species richness was better predicted by a model that included genetic diversity, with the more heterozygous hosts harboring fewer endoparasite species. Genetic diversity was also important in predicting abundance of replete ticks (Dermacentor variabilis). This association fit a curvilinear trend, with hosts that had either high or low levels of heterozygosity harboring fewer parasites than those with intermediate levels. In contrast, genetic diversity was not important in predicting abundance of non-replete ticks and lice (Trichodectes octomaculatus). No strong single-locus effects were observed for either endoparasites or replete ticks. Our results suggest that in outbred populations multi-locus diversity might be important for coping with parasitism. The differences in the relationships between heterozygosity and parasitism for the different parasites suggest that the role of genetic diversity varies with parasite-mediated selective pressures.

The opportunity for balancing selection in experimental populations of Caenorhabditis elegans

The role of balancing selection in maintaining diversity during the evolution of sexual populations to novel environments is poorly understood. To address this issue, we studied the impact of two mating systems, androdioecy and dioecy, on genotype distributions during the experimental evolution of Caenorhabditis elegans. We analyzed the temporal trajectories of 334 single nucleotide polymorphisms, covering 1/3 of the genome, and found extensive allele frequency changes and little loss of heterozygosities after 100 generations. As modeled with numerical simulations, SNP differentiation was consistent with genetic drift and average fitness effects of 2%, assuming that selection acted independently at each locus. Remarkably, inbreeding by self-fertilization was of little consequence to SNP differentiation. Modeling selection on deleterious recessive alleles suggests that the initial evolutionary dynamics can be explained by associative overdominance, but not the later stages because much lower heterozygosities would be maintained during experimental evolution. By contrast, models with selection on true overdominant loci can explain the heterozygote excess observed at all periods, particularly when negative epistasis or independent fitness effects were considered. Overall, these findings indicate that selection at single loci, including purging of recessive alleles, underlies most of the genetic differentiation accomplished during the experiment. Nonetheless, they also imply that maintenance of genetic diversity may in large part be due to balancing selection at multiple loci.

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

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

Heterozygosity-fitness correlations among wild populations of European tree frog (Hyla arborea) detect fixation load

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.

Inbreeding depression accumulation across life-history stages of the endangered Takahe

Studies evaluating the impact of inbreeding depression on population viability of threatened species tend to focus on the effects of inbreeding at a single life-history stage (e.g., juvenile survival). We examined the effects of inbreeding across the full life-history continuum, from survival up to adulthood, to subsequent reproductive success, and to the recruitment of second-generation offspring, in wild Takahe ( Porphyrio hochstetteri ) by analyzing pedigree and fitness data collected over 21 breeding seasons. Although the effect size of inbreeding at individual life-history stages was small, inbreeding depression accumulated across multiple life-history stages and ultimately reduced long-term fitness (i.e., successful recruitment of second-generation offspring). The estimated total lethal equivalents (2B) summed across all life-history stages were substantial (16.05, 95% CI 0.08-90.8) and equivalent to an 88% reduction in recruitment of second-generation offspring for closely related pairs (e.g., sib-sib pairings) relative to unrelated pairs (according to the pedigree). A history of small population size in the Takahe could have contributed to partial purging of the genetic load and the low level of inbreeding depression detected at each single life-history stage. Nevertheless, our results indicate that such "purged" populations can still exhibit substantial inbreeding depression, especially when small but negative fitness effects accumulate across the species' life history. Because inbreeding depression can ultimately affect population viability of small, isolated populations, our results illustrate the importance of measuring the effects of inbreeding across the full life-history continuum.

Inbreeding coefficient and heterozygosity-fitness correlations in unhatched and hatched song sparrow nestmates

Heterozygosity-fitness correlations use molecular measures of heterozygosity as proxy estimates of individual inbreeding coefficients (f) to examine relationships between inbreeding and fitness traits. Heterozygosity-fitness correlations partly depend on the assumption that individual heterozygosity and f are strongly and negatively correlated. Although theory predicts that this relationship will be strongest when mean f and variance in f are high, few studies of heterozygosity-fitness correlations include estimates of f based on pedigrees, which allow for more thorough examinations of the relationship between f, heterozygosity and fitness in nature. We examined relationships between pedigree-based estimates of f, multilocus heterozygosity (MLH) and the probability of survival to hatch in song sparrow nestmates. f and MLH were weakly, but significantly negatively correlated. Inbreeding coefficient predicted the probability of survival to hatch. In contrast, MLH did not predict the probability of survival to hatch nor did it account for residual variation in survival to hatch after statistically controlling for the effects of f. These results are consistent with the expectation that heterozygosity-f correlations will be weak when mean and variance in f are low. Our results also provide empirical support for recent simulation studies, which show that variation in MLH among siblings with equal f can be large and may obscure MLH-fitness relationships.

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.

Efficiency of selection, as measured by single nucleotide polymorphism variation, is dependent on inbreeding rate in Drosophila melanogaster

It is often hypothesized that slow inbreeding causes less inbreeding depression than fast inbreeding at the same absolute level of inbreeding. Possible explanations for this phenomenon include the more efficient purging of deleterious alleles and more efficient selection for heterozygote individuals during slow, when compared with fast, inbreeding. We studied the impact of inbreeding rate on the loss of heterozygosity and on morphological traits in Drosophila melanogaster. We analysed five noninbred control lines, 10 fast inbred lines and 10 slow inbred lines; the inbred lines all had an expected inbreeding coefficient of approximately 0.25. Forty single nucleotide polymorphisms in DNA coding regions were genotyped, and we measured the size and shape of wings and counted the number of sternopleural bristles on the genotyped individuals. We found a significantly higher level of genetic variation in the slow inbred lines than in the fast inbred lines. This higher genetic variation was resulting from a large contribution from a few loci and a smaller effect from several loci. We attributed the increased heterozygosity in the slow inbred lines to the favouring of heterozygous individuals over homozygous individuals by natural selection, either by associative over-dominance or balancing selection, or a combination of both. Furthermore, we found a significant polynomial correlation between genetic variance and wing size and shape in the fast inbred lines. This was caused by a greater number of homozygous individuals among the fast inbred lines with small, narrow wings, which indicated inbreeding depression. Our results demonstrated that the same amount of inbreeding can have different effects on genetic variance depending on the inbreeding rate, with slow inbreeding leading to higher genetic variance than fast inbreeding. These results increase our understanding of the genetic basis of the common observation that slow inbred lines express less inbreeding depression than fast inbred lines. In addition, this has more general implications for the importance of selection in maintaining genetic variation.

Association between SNP heterozygosity and quantitative traits in the Framingham Heart Study

Associations between multilocus heterozygosity and fitness traits, also termed heterozygosity and fitness correlations (HFCs), have been reported in numerous organisms. These studies, in general, indicate a positive relationship between heterozygosity and fitness traits. We studied the association between genome-wide heterozygosity at 706 non-synonymous and synonymous SNPs and 19 quantitative traits, including morphological, biochemical and fitness traits in the Framingham Heart Study. Statistically significant association was found between heterozygosity and systolic and diastolic blood pressures as well as left ventricular diameter and wall thickness. These results suggest that heterozygosity may be associated with traits, such as blood pressure that closely track environmental variations. Balancing selection may be operating in the maintenance of heterozygosity and the major components of blood pressure and hypertension. Genome wide SNP heterozygosity may be used to understand the phenomenon of dominance as well as the evolutionary basis of many quantitative traits in humans.

Single Nucleotide polymorphisms and their relationship to codon usage bias in the Pacific oyster Crassostrea gigas

DNA sequence polymorphism and codon usage bias were investigated in a set of 41 nuclear loci in the Pacific oyster Crassostrea gigas. Our results revealed a very high level of DNA polymorphism in oysters, in the order of magnitude of the highest levels reported in animals to date. A total of 290 single nucleotide polymorphisms (SNPs) were detected, 76 of which being localised in exons and 214 in non-coding regions. Average density of SNPs was estimated to be one SNP every 60 bp in coding regions and one every 40 bp in non-coding regions. Non-synonymous substitutions contributed substantially to the polymorphism observed in coding regions. The non-synonymous to silent diversity ratio was 0.16 on average, which is fairly higher to the ratio reported in other invertebrate species recognised to display large population sizes. Therefore, purifying selection does not appear to be as strong as it could have been expected for a species with a large effective population size. The level of non-synonymous diversity varied greatly from one gene to another, in accordance with varying selective constraints. We examined codon usage bias and its relationship with DNA polymorphism. The table of optimal codons was deduced from the analysis of an EST dataset, using EST counts as a rough assessment of gene expression. As recently observed in some other taxa, we found a strong and significant negative relationship between codon bias and non-synonymous diversity suggesting correlated selective constraints on synonymous and non-synonymous substitutions. Codon bias as measured by the frequency of optimal codons for expression might therefore provide a useful indicator of the level of constraint upon proteins in the oyster genome.

Lipid remodeling in wild and selectively bred hard clams at low temperatures in relation to genetic and physiological parameters

A temperature decrease usually induces an ordering effect in membrane phospholipids, which can lead to membrane dysfunction. Poikilotherms inhabiting eurythermal environments typically counteract this temperature effect by remodeling membrane lipids as stipulated in the homeoviscous adaptation theory (HVA). Hard clams, Mercenaria mercenaria, can suffer high overwintering mortalities in the Gulf of St Lawrence, Canada. The selectively bred M. mercenaria var. notata can have higher overwintering mortalities than the wild species, thus suggesting that the two varieties have different degrees of adaptation to low temperatures. The objective of this study was to investigate the changes in lipid composition of soft tissues in wild and selected hard clams in relation to their metabolic and genetic characteristics. Clams were placed at the northern limit of their distribution from August 2003 to May 2004; they were exposed to a gradual temperature decrease and then maintained at <0 degrees C for 3.5 months. This study is the first to report a major remodeling of lipids in this species as predicted by HVA; this remodeling involved a sequential response of the phospholipid to sterol ratio as well as in levels of 22:6n-3 and non-methylene interrupted dienoic fatty acids. Hard clams showed an increase in 20:5n-3 as temperature decreased, but this was not maintained during overwintering, which suggests that 20:5n-3 may have been used for eicosanoid biosynthesis as a stress response to environmental conditions. Selectively bred hard clams were characterized by a higher metabolic demand and a deviation from Hardy-Weinberg equilibrium at several genetic loci due to a deficit in heterozygote frequency compared with wild clams, which is believed to impose additional stress and render these animals more vulnerable to overwintering mortality. Finally, an intriguing finding is that the lower metabolic requirements of wild animals coincide with a lower unsaturation index of their lipids, as predicted by Hulbert's theory of membranes as pacemakers of metabolism.

Reliable selfing rate estimates from imperfect population genetic data

Genotypic frequencies at codominant marker loci in population samples convey information on mating systems. A classical way to extract this information is to measure heterozygote deficiencies (FIS) and obtain the selfing rate s from FIS = s/(2 - s), assuming inbreeding equilibrium. A major drawback is that heterozygote deficiencies are often present without selfing, owing largely to technical artefacts such as null alleles or partial dominance. We show here that, in the absence of gametic disequilibrium, the multilocus structure can be used to derive estimates of s independent of FIS and free of technical biases. Their statistical power and precision are comparable to those of FIS, although they are sensitive to certain types of gametic disequilibria, a bias shared with progeny-array methods but not FIS. We analyse four real data sets spanning a range of mating systems. In two examples, we obtain s = 0 despite positive FIS, strongly suggesting that the latter are artefactual. In the remaining examples, all estimates are consistent. All the computations have been implemented in a open-access and user-friendly software called rmes (robust multilocus estimate of selfing) available at http://ftp.cefe.cnrs.fr, and can be used on any multilocus data. Being able to extract the reliable information from imperfect data, our method opens the way to make use of the ever-growing number of published population genetic studies, in addition to the more demanding progeny-array approaches, to investigate selfing rates.

Genotypic selection in Daphnia populations consisting of inbred sibships

The genetic basis of fitness reduction associated with inbreeding is still poorly understood. Here we use associations between allozyme genotypes and fitness to investigate the genetic basis of inbreeding depression in experimental outdoor populations of the water flea, Daphnia magna. In Daphnia, a phase of clonal reproduction follows hatching from sexually produced resting eggs, and changes in genotype frequencies during the clonal phase can be used to estimate fitness. Our experiment resembles natural colonization of ponds in that single clones colonize an empty pool, expand asexually and produce sexual offspring by selfing (sisters mate with their clonal brothers). These offspring diapause and form populations consisting of selfed sibships in the following spring. In 12 of 13 experimental populations, genotypes of selfed hatchlings after diapause conformed to Mendelian expectations. During the subsequent ca. 10 asexual generations, however, genotype frequencies changed significantly at 19 of 27 single loci studied within populations, mostly in favour of heterozygotes, with heterozygosity at multiple loci affecting the change in genotype frequency multiplicatively. Because variance in heterozygosity among siblings at a given marker reflects only heterozygosity in the chromosomal region around this marker, our results suggest that selection at fitness-associated loci in the chromosomal regions near the markers were responsible for these changes. The genotype frequency changes were more consistent with selection acting on linked loci than on the allozymes themselves. Taken together, the evidence for abundant selection in the chromosomal regions of the markers and the fact that changes in genotype frequencies became apparent only after several generations of clonal selection, point to a genetic load consisting of many alleles of small or intermediate effects, which is consistent with the strong genetic differentiation and repeated genetic bottlenecks in the metapopulation from which the animals for this study were obtained.

Marker-based investigation of inbreeding depression in the endangered species Brassica insularis

Various methods have been proposed to estimate inbreeding depression and to assess its consequences for natural populations. As an alternative to controlled crosses, the use of molecular markers has allowed direct investigation of inbreeding depression in natural populations, but usually suffers from low statistical power. Here, we investigated the effect of inbreeding depression on survival in two populations of the rare species Brassica insularis, using both controlled crosses and a marker-based approach. We compare the respective merits of the two approaches for studying inbreeding depression. We also use information from the molecular markers to dissect in detail patterns of inbreeding depression in this species. A posteriori, we find that combining the approaches was not necessary to obtain simple point estimates of inbreeding depression. However, using molecular markers may give insight into the genetic basis of inbreeding depression, such as the occurrence of epistatic interactions among deleterious alleles or purging.

Isolation by Distance in the Eastern Oyster, Crassostrea virginica, in Chesapeake Bay

Intensive efforts are underway to restore depleted stocks of Crassostrea virginica in Chesapeake Bay. However, the extent of gene flow among local populations, an important force mediating the success of these endeavors, is poorly understood. Spatial and temporal population structures were examined in C. virginica from Chesapeake Bay using eight microsatellite loci. Deficits in heterozygosity relative to Hardy-Weinberg expectations were seen at all loci and were best explained by null alleles. Permutation tests indicated that heterozygote deficiency reduced power in tests of differentiation. Nonetheless, genotypic exact tests demonstrated significant levels of geographic differentiation overall, and a subtle pattern of isolation by distance (IBD) was observed. Comparisons between age classes failed to show differences in genotype frequencies, allelic richness, gene diversity, or differentiation as measured by F(ST), contrary to predictions made by the sweepstakes hypothesis. The IBD pattern could reflect an evolutionary equilibrium established because local gene flow predominates, or be influenced in either direction by recent anthropogenic activities. An evolutionary interpretation appears justified as more parsimonious, implying that local efforts to restore oyster populations will have local demographic payoffs, perhaps at the scale of tributaries or regional subestuaries within Chesapeake Bay.

Heterozygosity, inbreeding and neonatal traits in Soay sheep on St Kilda

We investigated whether birth weight and neonatal survival, a period within which 24% of all mortalities occur, were correlated with levels of inbreeding in St Kilda Soay sheep, using pedigree inbreeding coefficients and four marker-based estimators of inbreeding. None of the inbreeding estimators, either of the offspring, or of their mothers, explained significant variation in a lamb's birth weight or probability of surviving the neonatal period, suggesting low inbreeding depression for these traits. We evaluated the correlation between the marker-based measures of inbreeding and inbreeding coefficients obtained from the Soay pedigree, where paternal links were inferred using the same panel of microsatellite markers. Even when using a relatively complete portion of the pedigree, in which all individuals had known maternal and paternal grandparents, the correlation was found to be weak (r = -0.207, where mean f = 0.0168). These results add support to the recent prediction that when the mean and variance in inbreeding are low in a population, heterozygosity-fitness correlations can be very weak or even undetectable. The pursuit of more detailed pedigrees offers the best prospect for identifying inbreeding depression within this study population.

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.

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.

Understanding the relationship between the inbreeding coefficient and multilocus heterozygosity: theoretical expectations and empirical data

Geneticists have been interested in inbreeding and inbreeding depression since the time of Darwin. Two alternative approaches that can be used to measure how inbred an individual is involve the use of pedigree records to estimate inbreeding coefficients or molecular markers to measure multilocus heterozygosity. However, the relationship between inbreeding coefficient and heterozygosity has only rarely been investigated. In this paper, a framework to predict the relationship between the two variables is presented. In addition, microsatellite genotypes at 138 loci spanning all 26 autosomes of the sheep genome were used to investigate the relationship between inbreeding coefficient and multilocus heterozygosity. Multilocus heterozygosity was only weakly correlated with inbreeding coefficient, and heterozygosity was not positively correlated between markers more often than expected by chance. Inbreeding coefficient, but not multilocus heterozygosity, detected evidence of inbreeding depression for morphological traits. The relevance of these findings to the causes of heterozygosity--fitness correlations is discussed and predictions for other wild and captive populations are presented.

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.