Soft Shoulders Ahead: Spurious Signatures of Soft and Partial Selective Sweeps Result from Linked Hard Sweeps

Characterizing the nature of the adaptive process at the genetic level is a central goal for population genetics. In particular, we know little about the sources of adaptive substitution or about the number of adaptive variants currently segregating in nature. Historically, population geneticists have focused attention on the hard-sweep model of adaptation in which a de novo beneficial mutation arises and rapidly fixes in a population. Recently more attention has been given to soft-sweep models, in which alleles that were previously neutral, or nearly so, drift until such a time as the environment shifts and their selection coefficient changes to become beneficial. It remains an active and difficult problem, however, to tease apart the telltale signatures of hard vs. soft sweeps in genomic polymorphism data. Through extensive simulations of hard- and soft-sweep models, here we show that indeed the two might not be separable through the use of simple summary statistics. In particular, it seems that recombination in regions linked to, but distant from, sites of hard sweeps can create patterns of polymorphism that closely mirror what is expected to be found near soft sweeps. We find that a very similar situation arises when using haplotype-based statistics that are aimed at detecting partial or ongoing selective sweeps, such that it is difficult to distinguish the shoulder of a hard sweep from the center of a partial sweep. While knowing the location of the selected site mitigates this problem slightly, we show that stochasticity in signatures of natural selection will frequently cause the signal to reach its zenith far from this site and that this effect is more severe for soft sweeps; thus inferences of the target as well as the mode of positive selection may be inaccurate. In addition, both the time since a sweep ends and biologically realistic levels of allelic gene conversion lead to errors in the classification and identification of selective sweeps. This general problem of “soft shoulders” underscores the difficulty in differentiating soft and partial sweeps from hard-sweep scenarios in molecular population genomics data. The soft-shoulder effect also implies that the more common hard sweeps have been in recent evolutionary history, the more prevalent spurious signatures of soft or partial sweeps may appear in some genome-wide scans.     

Inferring the effects of demography and selection on Drosophila melanogaster populations from a chromosome-wide scan of DNA variation

Identifying regions of the Drosophila melanogaster genome that have been recent targets of positive Darwinian selection will provide evidence for adaptations that have helped this species to colonize temperate habitats. We have begun a search for such genomic regions by analyzing multiple loci (about 250) dispersed across the X chromosome in a putatively ancestral population from East Africa and a derived European population. For both populations we found evidence for past changes in population size. We estimated that a major bottleneck associated with the colonization of Europe occurred about 3,500-16,000 years ago. We also found that while this bottleneck can account for most of the reduction in variation observed in the European sample, there is a deficit of polymorphism in some genomic regions that cannot be explained by demography alone.     

Detecting strong positive selection in the genome

New statistical tests have been developed in the past decade that enable us to infer evidence of recent strong positive selection from genome-wide data on single-nucleotide polymorphism and to localize the targets of selection in the genome. Based on these tests, past demographic events that led to distortions of the site-frequency spectrum of variation can be distinguished from selection, in particular if linkage disequilibrium is taken into account. These methods have been successfully applied to species from which complete sequence information and polymorphism data are available, including Drosophila melanogaster, humans, and several plant species. To make full use of the available data, however, the tests that were primarily designed for panmictic populations need to be extended to spatially structured populations.     

The impact of epistatic selection on the genomic traces of selection

The rapid accumulation of genomic data has led to an explosion of studies searching for signals of past selection left within DNA sequences. Yet the majority of theoretical studies investigating the traces of selection have assumed a simple form of selection, without interactions among selectively fixed sites. Fitness interactions—‘epistasis’—are commonplace, however, and take on a myriad of forms ( Whitlock et al. 1995 ; Segrèet al. 2005 ; Phillips 2008 ). It is thus important to determine how such epistasis would influence selective sweeps. On p. 5018 of this issue, Takahasi (2009) explores the effect of epistasis on genetic variation neighbouring two sites that interact in determining fitness, finding that such epistasis has a dramatic impact on the genetic variability in regions surrounding the interacting sites.     

Limited Evidence for Classic Selective Sweeps in African Populations

While hundreds of loci have been identified as reflecting strong-positive selection in human populations, connections between candidate loci and specific selective pressures often remain obscure. This study investigates broader patterns of selection in African populations, which are underrepresented despite their potential to offer key insights into human adaptation. We scan for hard selective sweeps using several haplotype and allele-frequency statistics with a data set of nearly 500,000 genome-wide single-nucleotide polymorphisms in 12 highly diverged African populations that span a range of environments and subsistence strategies. We find that positive selection does not appear to be a strong determinant of allele-frequency differentiation among these African populations. Haplotype statistics do identify putatively selected regions that are shared across African populations. However, as assessed by extensive simulations, patterns of haplotype sharing between African populations follow neutral expectations and suggest that tails of the empirical distributions contain false-positive signals. After highlighting several genomic regions where positive selection can be inferred with higher confidence, we use a novel method to identify biological functions enriched among populations’ empirical tail genomic windows, such as immune response in agricultural groups. In general, however, it seems that current methods for selection scans are poorly suited to populations that, like the African populations in this study, are affected by ascertainment bias and have low levels of linkage disequilibrium, possibly old selective sweeps, and potentially reduced phasing accuracy. Additionally, population history can confound the interpretation of selection statistics, suggesting that greater care is needed in attributing broad genetic patterns to human adaptation.     

Natural selection at linked sites in humans

Theoretical and empirical work indicates that patterns of neutral polymorphism can be affected by linked, selected mutations. Under background selection, deleterious mutations removed from a population by purifying selection cause a reduction in linked neutral diversity. Under genetic hitchhiking, the rise in frequency and fixation of beneficial mutations also reduces the level of linked neutral polymorphism. Here we review the evidence that levels of neutral polymorphism in humans are affected by selection at linked sites. We then discuss four approaches for distinguishing between background selection and genetic hitchhiking based on (i) the relationship between polymorphism level and recombination rate for neutral loci with high mutation rates, (ii) relative levels of variation on the X chromosome and the autosomes, (iii) the frequency distribution of neutral polymorphisms, and (iv) population-specific patterns of genetic variation. Although the evidence for selection at linked sites in humans is clear, current methods and data do not allow us to clearly assess the relative importance of background selection and genetic hitchhiking in humans. These results contrast with those obtained for Drosophila, where the signals of positive selection are stronger.     

Selective sweeps in multilocus models of quantitative traits

We study the trajectory of an allele that affects a polygenic trait selected toward a phenotypic optimum. Furthermore, conditioning on this trajectory we analyze the effect of the selected mutation on linked neutral variation. We examine the well-characterized two-locus two-allele model but we also provide results for diallelic models with up to eight loci. First, when the optimum phenotype is that of the double heterozygote in a two-locus model, and there is no dominance or epistasis of effects on the trait, the trajectories of selected mutations rarely reach fixation; instead, a polymorphic equilibrium at both loci is approached. Whether a polymorphic equilibrium is reached (rather than fixation at both loci) depends on the intensity of selection and the relative distances to the optimum of the homozygotes at each locus. Furthermore, if both loci have similar effects on the trait, fixation of an allele at a given locus is less likely when it starts at low frequency and the other locus is polymorphic (with alleles at intermediate frequencies). Weaker selection increases the probability of fixation of the studied allele, as the polymorphic equilibrium is less stable in this case. When we do not require the double heterozygote to be at the optimum we find that the polymorphic equilibrium is more difficult to reach, and fixation becomes more likely. Second, increasing the number of loci decreases the probability of fixation, because adaptation to the optimum is possible by various combinations of alleles. Summaries of the genealogy (height, total length, and imbalance) and of sequence polymorphism (number of polymorphisms, frequency spectrum, and haplotype structure) next to a selected locus depend on the frequency that the selected mutation approaches at equilibrium. We conclude that multilocus response to selection may in some cases prevent selective sweeps from being completed, as described in previous studies, but that conditions causing this to happen strongly depend on the genetic architecture of the trait, and that fixation of selected mutations is likely in many instances.     

Sexual selection is stabilizing selection in pupfish [Cyprinodon pecosensis)

One of the predictions of the 'good genes' model of sexual selection is that reproductively successful males with well-developed indicator traits should show smaller variances for non-indicator traits, that are not directly associated with mating success, when compared to non-breeding males and females. Thus sexual selection should reinforce stabilizing natural selection in reducing the variance in quantitative traits. This prediction is tested by analysing variation in eight morphological traits of breeding males, non-breeding males, and females of pupfish (Cyprinodon pecosensis). Breeding males tended to be less variable than non-breeding males for all principal component factors, and for all morphological traits except for depth, although these differences were statistically significant only for PC2, and PC5 and for pelvic fin length, number of pelvic fin rays and number of preopercular and preorbital pores. Similarly, breeding males tended to be less variable than females for all principal component factors and for all morphological traits except for number of preopercular pores. These differences were statistically significant for PC2, and for depth, pelvic fin length, number of preorbital pores and pectoral fin rays. The overall pattern of reduced variability in independent traits of breeding males revealed by principal component analysis is very consistent and highly significant (P<10⁵). These results support the prediction of the 'good genes' model and show that reproductively active males are subject to more severe stabilizing selection for several quantitative traits than non-breeding males and females. Thus sexual selection, through male-male competition, female choice, or an interaction of both selective processes, results in stabilizing selection on quantitative morphological traits.     

Reduced variation around drug-resistant dhfr alleles in African Plasmodium falciparum

We have measured microsatellite diversity at 26 markers around the dhfr gene in pyrimethamine-sensitive and -resistant parasites collected in southeast Africa. Through direct comparison with diversity on sensitive chromosomes we have found significant loss of diversity across a region of 70 kb around the most highly resistant allele which is evidence of a selective sweep attributable to selection through widespread use of pyrimethamine (in combination with sulfadoxine) as treatment for malaria. Retrospective analysis through four years of direct and continuous selection from use of sulfadoxine-pyrimethamine as first-line malaria treatment on a Plasmodium falciparum population in KwaZulu Natal, South Africa, has revealed how recombination significantly narrowed the margins of the selective sweep over time. A deterministic model incorporating selection coefficients measured during the same interval indicates that the transition was toward a state of recombination-selection equilibrium. We compared loss of diversity around the same resistance allele in two populations at either extreme of the range of entomological inoculation rates (EIRs), namely, under one infective bite per year in Mpumalanga, South Africa, and more than one per day in southern Tanzania. EIRs determine effective recombination rates and are expected to profoundly influence the dimensions of the selective sweep. Surprisingly, the dimensions were broadly consistent across both populations. We conclude that despite different recombination rates and contrasting drug selection histories in neighboring countries, the region-wide movement of resistant parasites has played a key role in the establishment of resistance in these populations and the dimensions of the selective sweep are dominated by the influence of high initial starting frequencies.     

Identification of selection signatures in livestock species

The identification of regions that have undergone selection is one of the principal goals of theoretical and applied evolutionary genetics. Such studies can also provide information about the evolutionary processes involved in shaping genomes, as well as physical and functional information about genes/genomic regions. Domestication followed by breed formation and selection schemes has allowed the formation of very diverse livestock breeds adapted to a wide variety of environments and with special characteristics. The advances in genomics in the last five years have enabled the development of several methods to detect selection signatures and have resulted in the publication of a considerable number of studies involving livestock species. The aims of this review are to describe the principal effects of natural/artificial selection on livestock genomes, to present the main methods used to detect selection signatures and to discuss some recent results in this area. This review should be useful also to research scientists working with wild animals/non-domesticated species and plant biologists working with breeding and evolutionary biology.     

Low levels of nucleotide diversity in mammalian Y chromosomes

Sex chromosomes provide a useful context for the study of the relative importance of evolutionary forces affecting genetic diversity. The human Y chromosome shows levels of nucleotide diversity 20% that of autosomes, which is significantly less than expected when differences in effective population size and sex-specific mutation rates are taken into account. To study the generality of low levels of Y chromosome variability in mammalian genomes, we investigated nucleotide diversity in intron sequences of X (1.1-3.0 kb) and Y (0.7-3.5 kb) chromosome genes of five mammals: lynx, wolf, reindeer, cattle, and field vole. For all species, nucleotide diversity was found to be lower on Y than on X, with no segregating site observed in Y-linked sequences of lynx, reindeer, and cattle. For X chromosome sequences, nucleotide diversity was in the range of 1.6 x 10(-4) (lynx) to 8.0 x 10(-4) (field vole). When differences in effective population size and the extent of the male mutation bias were taken into account, all five species showed evidence of reduced levels of Y chromosome variability. Reduced levels of Y chromosome variability have also been observed in Drosophila and in plants, as well as in the female-specific W chromosome of birds. Among the different factors proposed to explain low levels of genetic variability in the sex-limited chromosome (Y/W), we note that selection is the only factor that is broadly applicable irrespective of mode of reproduction and whether there is male or female heterogamety.     

A pooling approach to detect signatures of selective sweeps in genome scans using microsatellites

We have evaluated a pooling approach that can reduce the number of polymerase chain reactions in a screen for selective sweeps by more than an order of magnitude. We show that the complex peak pattern that results from pooling of all samples from a given population is a faithful reflection of the composite pattern of the individual alleles, although with an under‐representation of the larger alleles. Candidate loci for selective sweeps can be identified by visual inspection of the pool patterns. We have also implemented a software tool, which can find suitable microsatellite loci in the vicinity of annotated genes.     

大白猪和通城猪全基因组选择性清扫分析

长期的人工选择使猪的生产性能得到显著提高,与选择相关的基因组区域也随之发生特定遗传变异表征(选择信号)。不同类型品种所受到选择强度不一,选择信号亦不相同,选择性清扫分析已逐渐成为选择信号的主要检测手段。文章基于商用型大白猪(n=45)和地方猪品种通城猪(n=45)的猪60K SNP芯片分型数据,借助遗传分化系数Fst法进行选择信号检测分析。利用gPLINK软件设定质控标准,共计34 304个SNPs被筛选出用于统计分析。使用Genepop软件包计算两个猪品种之间的遗传分化参数Fst,所得Fst平均值为0.3209。选取Fst>0.7036(即占总Fst值数目的 1%),共计344个SNPs被选择出来。SNP位置注释显示这些位点涉及到79个候选基因(Sus scrofa Build 9)。利用在线软件Ingenuity Pathway Analysis对候选基因的生物学通路进行网络分析,发现它们多与生长繁殖及免疫应答有关,如NCOA6、ERBB4、RUNX2和APOB等基因。研究结果为进行猪产肉、抗病等性状候选基因和致因突变深入挖掘提供了有益参考。