Selection maintains MHC diversity through a natural population bottleneck

A perceived consequence of a population bottleneck is the erosion of genetic diversity and concomitant reduction in individual fitness and evolutionary potential. Although reduced genetic variation associated with demographic perturbation has been amply demonstrated for neutral molecular markers, the effective management of genetic resources in natural populations is hindered by a lack of understanding of how adaptive genetic variation will respond to population fluctuations, given these are affected by selection as well as drift. Here, we demonstrate that selection counters drift to maintain polymorphism at a major histocompatibility complex (MHC) locus through a population bottleneck in an inbred island population of water voles. Before and after the bottleneck, MHC allele frequencies were close to balancing selection equilibrium but became skewed by drift when the population size was critically low. MHC heterozygosity generally conformed to Hardy-Weinberg expectations except in one generation during the population recovery where there was a significant excess of heterozygous genotypes, which simulations ascribed to strong differential MHC-dependent survival. Low allelic diversity and highly skewed frequency distributions at microsatellite loci indicated potent genetic drift due to a strong founder affect and/or previous population bottlenecks. This study is a real-time examination of the predictions of fundamental evolutionary theory in low genetic diversity situations. The findings highlight that conservation efforts to maintain the genetic health and evolutionary potential of natural populations should consider the genetic basis for fitness-related traits, and how such adaptive genetic diversity will vary in response to both the demographic fluctuations and the effects of selection.     

Haplotype structure in Ashkenazi Jewish BRCA1 and BRCA2 mutation carriers

Three founder mutations in BRCA1 and BRCA2 contribute to the risk of hereditary breast and ovarian cancer in Ashkenazi Jews (AJ). They are observed at increased frequency in the AJ compared to other BRCA mutations in Caucasian non-Jews (CNJ). Several authors have proposed that elevated allele frequencies in the surrounding genomic regions reflect adaptive or balancing selection. Such proposals predict long-range linkage disequilibrium (LD) resulting from a selective sweep, although genetic drift in a founder population may also act to create long-distance LD. To date, few studies have used the tools of statistical genomics to examine the likelihood of long-range LD at a deleterious locus in a population that faced a genetic bottleneck. We studied the genotypes of hundreds of women from a large international consortium of BRCA1 and BRCA2 mutation carriers and found that AJ women exhibited long-range haplotypes compared to CNJ women. More than 50% of the AJ chromosomes with the BRCA1 185delAG mutation share an identical 2.1?Mb haplotype and nearly 16% of AJ chromosomes carrying the BRCA2 6174delT mutation share a 1.4?Mb haplotype. Simulations based on the best inference of Ashkenazi population demography indicate that long-range haplotypes are expected in the context of a genome-wide survey. Our results are consistent with the hypothesis that a local bottleneck effect from population size constriction events could by chance have resulted in the large haplotype blocks observed at high frequency in the BRCA1 and BRCA2 regions of Ashkenazi Jews.     

Redistribution of gene frequency and changes of genetic variation following a bottleneck in population size

Although the distribution of frequencies of genes influencing quantitative traits is important to our understanding of their genetic basis and their evolution, direct information from laboratory experiments is very limited. In theory, different models of selection and mutation generate different predictions of frequency distributions. When a large population at mutation-selection balance passes through a rapid bottleneck in size, the frequency distribution of genes is dramatically altered, causing changes in observable quantities such as the mean and variance of quantitative traits. We investigate the gene frequency distribution of a population at mutation-selection balance under a joint-effect model of real stabilizing and pleiotropic selection and its redistribution and thus changes of the genetic properties of metric and fitness traits after the population passes a rapid bottleneck and expands in size. If all genes that affect the trait are neutral with respect to fitness, the additive genetic variance (VA) is always reduced by a bottleneck in population size, regardless of their degree of dominance. For genes that have been under selection, VA increases following a bottleneck if they are (partially) recessive, while the dominance variance increases substantially for any degree of dominance. With typical estimates of mutation parameters, the joint-effect model can explain data from laboratory experiments on the effect of bottlenecking on fitness and morphological traits, providing further support for it as a plausible mechanism for maintenance of quantitative genetic variation.     

A population-genetic test of founder effects and implications for Ashkenazi Jewish diseases

A founder effect can account for the presence of an allele at an unusually high frequency in an isolated population if the allele is selectively neutral and if all copies are identical by descent with a copy that either was carried by a founder individual or arose by mutation later. Here, a statistical test of both aspects of the founder-effect hypothesis is developed. The test is performed by a modified version of a program that implements the Slatkin-Bertorelle test of neutrality. The test is applied to several disease-associated alleles found predominantly in Ashkenazi Jews. Despite considerable uncertainty about the demographic history of Ashkenazi Jews and their ancestors, available genetic data are consistent with a founder effect resulting from a severe bottleneck in population size between a.d. 1100 and a.d. 1400 and an earlier bottleneck in a.d. 75, at the beginning of the Jewish Diaspora. The relatively high frequency of alleles causing four different lysosomal storage disorders, including Tay-Sachs disease and Gaucher disease, can be accounted for if the disease-associated alleles are recessive in their effects on reproductive fitness.     

A general population genetic framework for antagonistic selection that accounts for demography and recurrent mutation

Antagonistic selection--where alleles at a locus have opposing effects on male and female fitness ("sexual antagonism") or between components of fitness ("antagonistic pleiotropy")--might play an important role in maintaining population genetic variation and in driving phylogenetic and genomic patterns of sexual dimorphism and life-history evolution. While prior theory has thoroughly characterized the conditions necessary for antagonistic balancing selection to operate, we currently know little about the evolutionary interactions between antagonistic selection, recurrent mutation, and genetic drift, which should collectively shape empirical patterns of genetic variation. To fill this void, we developed and analyzed a series of population genetic models that simultaneously incorporate these processes. Our models identify two general properties of antagonistically selected loci. First, antagonistic selection inflates heterozygosity and fitness variance across a broad parameter range--a result that applies to alleles maintained by balancing selection and by recurrent mutation. Second, effective population size and genetic drift profoundly affect the statistical frequency distributions of antagonistically selected alleles. The "efficacy" of antagonistic selection (i.e., its tendency to dominate over genetic drift) is extremely weak relative to classical models, such as directional selection and overdominance. Alleles meeting traditional criteria for strong selection (N(e)s > 1, where N(e) is the effective population size, and s is a selection coefficient for a given sex or fitness component) may nevertheless evolve as if neutral. The effects of mutation and demography may generate population differences in overall levels of antagonistic fitness variation, as well as molecular population genetic signatures of balancing selection.     

Comparing the effects of genetic drift and fluctuating selection on genotype frequency changes in the scarlet tiger moth

One of the recurring arguments in evolutionary biology is whether evolution occurs principally through natural selection or through neutral processes such as genetic drift. A 60-year-long time series of changes in the genotype frequency of a colour polymorphism of the scarlet tiger moth, Callimorpha dominula, was used to compare the relative effects of genetic drift and variable natural selection. The analysis showed that most of the variation in frequency was the result of genetic drift. In addition, although selection was acting, mean fitness barely increased. This supports the 'Red Queen's hypothesis' that long-term improvements in fitness may not occur, because populations have to keep pace with changes in the environment.     

Bad coat, ripped genes: cryptic selection on coat colour varies with ontogeny in Soay sheep

A fundamental goal for evolutionary biologists is to connect the processes of natural selection and genetic drift with the maintenance of genetic variation in the wild. There are now scores of examples in which mapping phenotypes to genotypes has identified the molecular basis of traits in natural populations, but documenting the fitness consequences of these genes has proven more elusive. A long-term study of several thousand unmanaged sheep on the remote Scottish island of Hirta in the St. Kilda archipelago provides a unique opportunity to make these connections between genotype, phenotype and fitness. Since 1985, there has been a steady decline in the frequency of sheep with a rare uniform ('self'-) pigmentation pattern as opposed to the more common wild-type pattern of dark upper body and pale belly. In this issue of Molecular Ecology, Gratten et al. (2012) link variation at the agouti signalling protein (ASIP) gene with 25 years of pedigree information and individual fitness data to investigate the contributions of natural selection and genetic drift to this decline. Consistent with the decrease in frequency of self-type sheep, Gratten et al. demonstrate that recessive genetic variants at ASIP are associated with reduced lifetime fitness. But surprisingly, these variants show the opposite trend to the self-type phenotype by increasing throughout the study period. This paradoxical result occurs because heterozygotes that harbour the majority of self-type alleles have increased in frequency. This pattern would not have been detectable if only phenotype had been monitored because the pigmentation pattern of heterozygotes is indistinguishable from that of wild-type homozygotes. The study provides an excellent example of the importance of measuring selection at both the phenotypic and genetic level, and demonstrates how long-term pedigrees can be used to link alleles to phenotypes and ultimately to fitness.     

MEASUREMENTS OF NATURAL SELECTION ON FLORAL TRAITS IN WILD RADISH (RAPHANUS RAPHANISTRUM). II. SELECTION THROUGH LIFETIME MALE AND TOTAL FITNESS

It has often been suggested that selection on floral traits in hermaphroditic plants should occur primarily through differences in male fitness. However, measurements of selection on floral traits through differences in lifetime male fitness have been lacking. We measured selection on a variety of wild radish floral traits using lifetime male fitness measures derived from genetic paternity analysis. These male fitness estimates were then combined with estimates of lifetime female fitness of the same plants to produce measurements of selection based on lifetime total fitness. Contrary to the prediction above, there was no strong evidence for selection on floral morphology through male fitness differences in any of the three years of the study, but there was strong selection for increased flower size through female fitness differences in one year. The main determinant of both male and female fitness in all years was flower number; this lead to moderately positive correlations between male and female fitness in all three years.     

Heterozygote advantage in Tay-Sachs carriers?

Chi-square analyses of new data as well as data previously reported by Myrianthopoulos have shown that grandparents of Tay-Sachs carriers die from proportionally the same causes as grandparents of noncarriers. It is unlikely that there is any advantage to being a Tay-Sachs carrier insofar as resistance to tuberculosis is concerned. Our results are further evidence to support Fraikor's claim that the high carrier frequency of the allele in Ashkenazi Jews is probably caused by a combination of founder effect, genetic drift, and differential immigration patterns.     

Evolution of Fitness in Experimental Populations of Vesicular Stomatitis Virus

The evolution of fitness in experimental clonal populations of vesicular stomatitis virus (VSV) has been compared under different genetic (fitness of initial clone) and demographic (population dynamics) regimes. In spite of the high genetic heterogeneity among replicates within experiments, there is a clear effect of population dynamics on the evolution of fitness. Those populations that went through strong periodic bottlenecks showed a decreased fitness in competition experiments with wild type. Conversely, mutant populations that were transferred under the dynamics of continuous population expansions increased their fitness when compared with the same wild type. The magnitude of the observed effect depended on the fitness of the original viral clone. Thus, high fitness clones showed a larger reduction in fitness than low fitness clones under dynamics with included periodic bottleneck. In contrast, the gain in fitness was larger the lower the initial fitness of the viral clone. The quantitative genetic analysis of the trait ``fitness' in the resulting populations shows that genetic variation for the trait is positively correlated with the magnitude of the change in the same trait. The results are interpreted in terms of the operation of MULLER's ratchet and genetic drift as opposed to the appearance of beneficial mutations.     

Genetic drift during the spread phase of a biological invasion

Recent theoretical and experimental models have revealed the role played by evolution during species spread, and in particular have questioned the influence of genetic drift at range edges. By investigating the spread of an aquatic invader in patchy habitats, we quantified genetic drift and explored its consequences for genetic diversity and fitness. We examined the interplay of gene flow and genetic drift in 36 populations of the red swamp crayfish, Procambarus clarkii, in a relatively recently invaded wetland area (30 years, Brière, northwest France). Despite the small spatial scale of our study (15 km²), populations were highly structured according to the strong barrier of land surfaces and revealed a clear pattern of colonization through watercourses. Isolated populations exhibited small effective sizes and low dispersal rates that depended on water connectivity, suggesting that genetic drift dominated in the evolution of allele frequencies in these populations. We also observed a significant decrease in the genetic diversity of isolated populations over only a 2‐year period, but failed to demonstrate an associated fitness cost using fluctuating asymmetry. This study documents the possible strong influence of genetic drift during the spread of a species, and such findings provide critical insights into the current context of profound rearrangements in species distributions due to global change.     

Simulating natural selection in landscape genetics

Linking landscape effects to key evolutionary processes through individual organism movement and natural selection is essential to provide a foundation for evolutionary landscape genetics. Of particular importance is determining how spatially-explicit, individual-based models differ from classic population genetics and evolutionary ecology models based on ideal panmictic populations in an allopatric setting in their predictions of population structure and frequency of fixation of adaptive alleles. We explore initial applications of a spatially-explicit, individual-based evolutionary landscape genetics program that incorporates all factors--mutation, gene flow, genetic drift and selection--that affect the frequency of an allele in a population. We incorporate natural selection by imposing differential survival rates defined by local relative fitness values on a landscape. Selection coefficients thus can vary not only for genotypes, but also in space as functions of local environmental variability. This simulator enables coupling of gene flow (governed by resistance surfaces), with natural selection (governed by selection surfaces). We validate the individual-based simulations under Wright-Fisher assumptions. We show that under isolation-by-distance processes, there are deviations in the rate of change and equilibrium values of allele frequency. The program provides a valuable tool (cdpop v1.0; http://cel.dbs.umt.edu/software/CDPOP/) for the study of evolutionary landscape genetics that allows explicit evaluation of the interactions between gene flow and selection in complex landscapes.     

Selection and microevolution of coat pattern are cryptic in a wild population of sheep

Understanding the maintenance of genetic variation in natural populations is a core aim of evolutionary genetics. Insight can be gained by quantifying selection at the level of the genotype, as opposed to the phenotype. Here, we show that in a natural population of Soay sheep which is polymorphic for coat pattern, recessive genetic variants at the causal gene, agouti signalling protein (ASIP) are associated with reduced lifetime fitness. This was due primarily to a reduction in juvenile survival of uniformly coloured (self-type) sheep, which are homozygous recessive, and occurs despite significantly higher reproductive success in surviving self-type adults. Consistent with their relatively low fitness, we show that the frequency of self-type individuals has declined from 1985 to 2008. Remarkably though, the frequency of the underlying self-allele has increased, because the frequency of heterozygous individuals (who harbour the majority of all self-alleles) has increased. Indeed, the ratio of observed/expected heterozygous individuals has increased during the study, such that there is now a significant excess of heterozygotyes. By employing gene-dropping simulations, we show that microevolutionary trends in the frequency and excess of ASIP heterozygotes are too pronounced to be caused by genetic drift. Studying this polymorphism at the level of phenotype rather than underlying genotype would have failed to detect cryptic fitness differences. We would also have been unable to rule out genetic drift as an evolutionary force driving genetic change. This highlights the importance of resolving the underlying genetic basis of phenotypic variation in explaining evolutionary dynamics.     

The effects of genetic drift in experimental evolution

Experimental evolution is characterized by exponential or logistic growth and periodic population bottlenecks. The fate of a rare beneficial mutation in these systems is influenced both by the bottleneck effect and by genetic drift. This paper explores the effects of incorporating genetic drift into models of fixation probability in populations with periodic bottlenecks. To model the inherent stochasticity during the growth phase in these populations, we assume a Poisson distribution of offspring. An analytical solution is developed to calculate the fixation probability and a computer simulation is used to verify the results. We find that the fixation rate of a favourable mutant is significantly lower when genetic drift is considered; fixation probability is reduced by over 25% for realistic experimental protocols. Our method is valid for both weak and strong selection; since very large selection coefficients have been reported in the experimental literature, this is an important improvement over previous results.     

The impact of random frequency-dependent mutations on the average population fitness

ABSTRACT: BACKGROUND: In addition to selection, the process of evolution is accompanied by stochastic effects, such as changing environmental conditions, genetic drift and mutations. Commonly it is believed that without genetic drift, advantageous mutations quickly fixate in a halpoid population due to strong selection and lead to a continuous increase of the average fitness. This conclusion is based on the assumption of constant fitness. However, for frequency dependent fitness, where the fitness of an individual depends on the interactions with other individuals in the population, this does not hold. RESULTS: We propose a mathematical model that allows to understand the consequences of random frequency dependent mutations on the dynamics of an infinite large population. The frequencies of different types change according to the replicator equations and the fitness of a mutant is random and frequency dependent. To capture the interactions of different types, we employ a payoff matrix of variable size and thus are able to accommodate an arbitrary number of mutations. We assume that at most one mutant type arises at a time. The payoff entries to describe the mutant type are random variables obeying a probability distribution which is related to the fitness of the parent type. CONCLUSIONS: We show that a random mutant can decrease the average fitness under frequency dependent selection, based on analytical results for two types, and on simulations for n types. Interestingly, in the case of at most two types the probabilities to increase or decrease the average fitness are independent of the concrete probability density function. Instead, they only depend on the probability that the payoff entries of the mutant are larger than the payoff entries of the parent type.     

Natural selection and the origin and maintenance of standard genetic marker systems

Natural selection has always been assumed to be the major force of evolution, but its presence has been difficult to demonstrate. A review of the evidence for selective differences among genotypes for most human genetic polymorphisms indicates there is little of a direct nature. Indirect theoretical evidence, however, seems to support a major role for natural selection, and it does not seem to support the hypothesis that most amino acid substitutions within the human species are neutral. Among small isolates, most of the gene frequency differences are most likely due to genetic drift or the founder effect, and the principal counterbalancing force is gene flow or migration. But genetic differences among the major human subdivisions do not seem to be due to the same interacting forces. One reason for the inability to detect selection has been an oversimplified view of its operation, which assigns genotypes a constant fitness in every generation. Many recent theoretical developments of more complicated kinds of selection may lead to a resolution of the problem and suggest better interpretations of the enormous amount of data on human genetic variation that is rapidly accumulating.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Rapid loss of MHC class II variation in a bottlenecked population is explained by drift and loss of copy number variation

Population bottlenecks may reduce genetic variation and potentially increase the risk of extinction. Here, we present the first study to use historic samples to analyse loss of variation at the major histocompatibility complex (MHC), which plays a central role in vertebrate disease resistance. Balancing selection acts on the MHC and could moderate the loss of variation expected from drift; however, in a Wisconsin population of greater prairie-chickens (Tympanuchus cupido), the number of MHC class II B alleles per individual declined by 44% following a population bottleneck, compared to a loss of only 8% at microsatellites. Simulations indicate that drift likely reduced MHC variation at the population level, as well as within individuals by reducing the number of gene copies per individual or by fixing the same alleles across multiple loci. These multiple effects of genetic drift on MHC variation could have important implications for immunity and fitness.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.