The effects of linkage on comparative estimators of selection
© Chan et al.; licensee BioMed Central Ltd. 2013
Received: 16 April 2013
Accepted: 29 October 2013
Published: 7 November 2013
A major goal of molecular evolution is to determine how natural selection has shaped the evolution of a gene. One approach taken by methods such as K A /K S and the McDonald-Kreitman (MK) test is to compare the frequency of non-synonymous and synonymous changes. These methods, however, rely on the assumption that a change in frequency of one mutation will not affect changes in frequency of other mutations.
We demonstrate that linkage between sites can bias measures of selection based on synonymous and non-synonymous changes. Using forward simulation of a Wright-Fisher process, we show that hitch-hiking of deleterious mutations with advantageous mutations can lead to overestimation of the number of adaptive substitutions, while background selection and clonal interference can distort the site frequency spectrum to obscure the signal for positive selection. We present three diagnostics for detecting these effects of linked selection and apply them to the human influenza (H3N2) hemagglutinin gene.
Various forms of linked selection have characteristic effects on MK-type statistics. The extent of background selection, hitch-hiking and clonal interference can be evaluated using the diagnostic statistics presented here. The diagnostics can also be used to determine how well we expect the MK statistics to perform and whether one form of the statistic may be preferable to another.
Understanding the mechanisms by which natural selection shapes the evolution of genes is one of the major aims of molecular evolution. One commonly used approach for the detection of positive selection in protein-coding sequences is based on comparing the frequency of non-synonymous or amino-acid (A) changes to the frequency of synonymous (S) changes . For simplicity, synonymous nucleotide changes that do not affect the protein are generally assumed to be neutral. In the absence of selection and accounting for the genetic code, we expect both types of changes to be equally probable so that the rate of non-synonymous substitutions per site (K A ) is equal to the rate of synonymous substitutions per site (K S ); a ratio of K A /K S >1 indicates positive selection favouring a change in the protein . However, this test is heavily conservative as proteins are generally under negative selection against amino acid changes that may affect protein function. Positive selection at a small number of sites may be masked by negative selection removing non-synonymous changes in the rest of the protein .
The McDonald-Kreitman (MK) test  attempts to account for the presence of negatively selected sites by comparing K A /K S to f, the proportion of nearly neutral sites in the sequence . If selection is strong, deleterious and beneficial mutations are expected to make little contribution to polymorphism; deleterious mutations are removed by selection and beneficial mutations reach fixation rapidly. Polymorphic sites are expected to consist largely of neutral variation, and the ratio of the number of neutral non-synonymous polymorphic sites (P A ) to the number of synonymous polymorphic sites (P S ) can be used as an estimator of f. In the MK test, positive selection is inferred when K A /K S >P A /P S . Following similar reasoning, K A /K S measured in a related sample can be used as a measure of selective constraint so that an increase in the K A /K S ratio implies positive selection [7, 8].
With the increasing availability of sequence data, various modifications of K A /K S methods have been developed to quantify the prevalence , strength [9, 10] and dynamics of positive selection [11, 12]. These methods rely on the assumption that sites segregate independently; that is, the change in frequency at one site will not affect the change in frequency at another site. In a large population with a high mutation rate, however, multiple mutations co-occur in the population and the change in frequency of one mutation also depends on selection acting at linked sites. Depending on the type of selection, linkage can have different effects; background selection, hitch-hiking and clonal interference can both increase or decrease fixation probability or polymorphism frequency relative to expected levels, which we describe below.
Background selection is the reduction in genetic variability caused by linkage to negatively selected sites . The effect of background selection on the probability of fixation is qualitatively similar to a reduction in effective population size [13–16], which implies a higher than expected value of K A /K S under negative selection and a lower than expected value of K A /K S under positive selection relative to expectations under independently segregating sites . Background selection also reduces the number of neutral polymorphic sites , and can result in a non-monotonic site-frequency spectrum, similar to the effect of continual adaptation [18, 19]. Linkage between sites introduces dependencies in the site frequency spectrum, increasing the covariance even if the mean is unchanged . Recent work with the structured coalescent  in a model of only negative selection, provides analytical expressions for the number of both neutral and deleterious mutations showing that the effective population size varies, both going back in time, and between individuals in different fitness classes.
When both positive and negative selection operate on a locus, the dynamics of linked neutral and deleterious mutations will also be affected by hitch-hiking . Birky and Walsh  showed that hitch-hiking does not affect the fixation probability at neutral sites but increases the fixation probability at negatively selected sites, which implies that K A /K S values are elevated relative to expectation under independently segregating sites. For the MK statistic, the effect of hitch-hiking depends on its effect on polymorphism relative to its effect on divergence. The effect of hitch-hiking on neutral polymorphism has been described by Braverman et al. , but has not been characterised on a selected background. Previous findings [23–25] were largely based on coalescent simulations which allow only a small number of sites to be under selection and model the trajectory of beneficial mutations deterministically. Forward simulation studies [14, 26–28] which begin with a number of positively selected sites and evolve towards mutation-selection equilibrium show that linkage affects a number of frequency-based statistics including Tajima’s D and heterozygosity.
Clonal interference (interactions between positively selected mutations) has also been predicted to reduce the fixation probability of beneficial mutations and promote the fixation of deleterious mutations; this was demonstrated in several experimental systems [29, 30]. More recently, theoretical models assuming continual adaptation with a high supply of beneficial mutations have been used to obtain analytical expressions characterising genetic diversity. These models predict a non-monotonic site frequency spectrum with a large number of both low and high-frequency mutations [18, 19, 31]. This is equivalent to large number of lineages coalescing simultaneously and is often described as multiple-mergers [18, 19, 31].
Here, we examine the joint effects of background selection, hitch-hiking and clonal interference on the K A /K S and MK statistic. Based on theoretical studies [18, 19, 21, 31], we expect different forms of distortion in the site-frequency spectrum due to these effects. Previous simulation studies [14, 26, 27, 32] have often considered these effects together, but here we distinguish between them by allowing both the strength of selection and the level of interference to vary. We do this using forward simulations with finite sites, allowing positive selection to occur at different times. Finally, we propose three diagnostic statistics to indicate the degree to which (a) hitch-hiking of deleterious mutations (b) background selection and (c) clonal interference affect a sample of protein-coding sequences.
The effect of background selection
We begin by examining the effect of negative selection and linkage without positive selection in a protein-coding region of 500 codons evolving under a Wright-Fisher process. Negative selection is described by the distribution of fitness effects (DFE) of non-synonymous changes, which are specific to each codon site. The DFE is modelled using a gamma distribution where a large value of the shape parameter β corresponds to a higher proportion of weakly deleterious mutations.
The combined effect of background selection, clonal interference and hitch-hiking
We show results for low levels of background selection (small u) in the left column and results for high levels of background selection in the right column. The (unscaled) synonymous site frequency spectrum is shown in the top row. When the effect of background selection is small, the synonymous site frequency spectrum is close to the expectation under independently segregating sites (θ/i; black dashed lines). Background selection (bold grey lines) reduces the level of synonymous variation, particularly at medium frequencies, leading to a non-monotonic distribution, but the effect is not as severe as clonal interference. Linked positive selection further reduces polymorphism levels; a slow rate of sweeps with strong selection (orange lines) primarily affects high-frequency mutations while a high supply of weak positive selection (green lines) results in smaller levels of reduction at both low and high frequencies. When both the supply rate and the strength of positive selection is strong (pink lines), the synonymous site frequency spectrum approaches θ/i2(black dotted line), which is the leading behaviour predicted for continual adaptation .
To examine how linkage affects selected mutations, we compare the ratio of the averaged frequency spectra for non-synonymous (A) and synonymous (S) sites (Figure 3, bottom row). The A/S ratio in the absence of positive selection is indicated by the bold grey line, whereas the A/S ratio for deleterious sites linked to positively selected sites is shown by coloured solid lines. The discrepancy between the grey and coloured lines reflects the effect of hitch-hiking; there is a slight increase in the A/S ratio at high-frequencies due to hitch-hiking. Note that the actual number of deleterious polymorphisms is reduced relative to simulations with no positive selection (Additional file 1: Figure S1) but the number of synonymous polymorphisms is reduced by a relatively greater proportion.
Comparing the A/S ratio with (dashed coloured lines) and without (solid coloured lines) beneficial mutations, it can be seen that beneficial mutations can segregate at all frequencies when the supply rate is high (green and pink lines), but mutations segregating at high frequencies tend to include more beneficial mutations. Comparison of the two panels in the bottom row also shows that higher levels of background selection increase the effect of both hitch-hiking (solid coloured lines) and clonal interference (dashed coloured lines), as distortions in the site-frequency spectrum tend to occur over a wider range of frequencies. Similar results are seen for larger values of β with more pronounced reductions of synonymous polymorphism due to background selection, and changes in the A/S ratio due to hitch-hiking and clonal interference are spread across a broader frequency range (Additional file 1: Figure S1).
The different panels in Figure 4 correspond to different combinations of positive and negative selection: low levels of background selection (due to strong negative selection) and infrequent positive sweeps (top left), low levels of background selection and frequent positive sweeps (bottom left), high levels of background selection and infrequent positive sweeps (top right) and high levels of background selection with frequent positive sweeps (bottom right). When the effect of background selection is large (top right), both and tend to underestimate the true number of adaptive substitutions. For low levels of background selection or frequent positive sweeps, the effect of hitch-hiking (controlled by β) and the amount of clonal interference (using the observed number of beneficial mutations as a proxy) has a consistent effect on the MK statistics. For small values of β so that low levels of hitch-hiking occur, has smaller bias than . However, for high levels of hitch-hiking is less biased, particularly when clonal interference is low. Results for different values of s b were qualitatively similar but with larger relative error for weaker positive selection.
The reason for these biases is intuitively clear from the site frequency spectrum. differs from only in that it does not use polymorphism data at high frequency. Therefore, is more robust against clonal interference (Figure 4, bottom row) as beneficial mutations are more likely to segregate at high frequencies. However, when weakly deleterious effects are prevalent (Figure 4, solid pink lines), is upwardly biased as it does not account for the relaxation of selective constraint due to positive selection. This is confirmed by the similar values obtained for and , suggesting that removal of high and low frequency polymorphisms in the context of linked selection has a similar effect to that expected under independently segregating sites, namely the removal of both positively and negatively selected mutations. The correction of Bhatt et al.  does not perform well when there are high levels of background selection as distortions in the site frequency spectrum are spread across a wider range of frequencies than without background selection.
Diagnostics for linkage effects
In the previous section, we showed that much of the bias in the comparative estimators can be explained in terms of background selection, hitch-hiking and clonal interference. In order to detect these effects using samples of protein-coding sequences, we construct and examine three diagnostic statistics.
If weak deleterious effects are rare, then we expect that most deleterious mutations are immediately removed from the population. In this case, most polymorphisms would be selectively neutral and we would expect that the ratio of non-synonymous to synonymous polymorphisms, at any frequency range, is simply determined by the mutational bias. The difference of the two ratios in D1is therefore expected to equal zero in the absence of weak deleterious effects and large values are indicative of a high frequency of weak deleterious mutations, which results in susceptibility to hitch-hiking.
Deleterious mutations are not expected to persist to medium frequencies, so polymorphisms at medium and high frequencies can be assumed to be neutral or beneficial. Under assumptions of neutrality and independently segregating sites, the expected number of polymorphic sites that occur at frequency i is given by E(p A (i))=θ A /i, where θ A =2u N L c/(c+1), giving an expectation of D2=0. Values of D2<0 can, therefore, indicate strong linkage effects due to an excess of beneficial or deleterious mutations.
To evaluate whether D1, D2and D3differ from zero, we use a non-parametric bootstrap, recalculating statistics after resampling with replacement from the original sequence sample. The scaling factor for mutation bias c, which is omitted from D1, is automatically accounted for by this method. Confidence intervals for D1were constructed from the bootstraps using the 2.5 to 97.5 percentiles. As D2is slightly biased, confidence intervals for D2and D3were constructed using the BCA method provided in R .
Application of diagnostics to human influenza A (H3N2)
It has long been known that linkage influences polymorphism frequencies and fixation probabilities which can adversely affect methods that assume independent segregation of sites. The MK statistic, which compares the number of polymorphic sites rather than using only frequency information, is generally considered to be more robust to linkage effects than frequency-based statistics [20, 27]. In this study, we show that the MK statistic can be affected, depending on the mode of linked selection and its characteristic effect of the site-frequency spectrum.
When background selection has a large effect, distortions in the site-frequency spectrum result in a downward bias in the MK statistics. However, when the effect of negative selection is low compared to positive selection, the performance of the MK statistics as a quantitative estimator of the number of adaptive substitutions is determined by the extent of hitch-hiking and clonal interference. When negative selection is strong so that levels of hitch-hiking are low, tends to perform better. Specifically, it is more robust against distortions of the site frequency spectrum at high frequencies caused by background selection or clonal interference. However, is more robust to hitch-hiking which occurs when weak negative selection is pervasive.
Our results are consistent with that of a recent study by Messer and Petrov  showing that
tends to underestimate the number of adaptive substitutions. We find this primarily occurs when background selection has a large effect and positive selection is infrequent. However, when positive selection is strong, hitch-hiking can also lead to overestimation as suggested in some empirical studies . When interactions between deleterious polymorphisms dominate the dynamics of the populations, the asymptotic correction proposed by Messer and Petrov  can be used to correct for underestimation due to low and medium frequency deleterious polymorphisms. This method corrects for deleterious mutations, as the relative abundance of deleterious mutations is reduced at higher frequencies, but cannot be applied for beneficial mutations which have increased relative abundance at higher frequencies.
Here, we show that, when background selection is relatively weak, choosing the appropriate form of the MK statistic can reduce estimation bias. Messer and Petrov’s  results apply for organisms with large genomes and many weakly deleterious mutations, but when genomes are small and selective effects are broadly distributed, as is the case viral populations [38, 39], the considerations raised in this study apply.
Understanding the effects of linked selection also affects our interpretation of these estimators. The number of adaptive substitutions cannot be directly related to either the strength of selection or the supply of beneficial mutations, but it is a combination of both of these factors. For example, Strelkowa and Lässig  and Koelle et al.  raised alternative hypotheses concerning whether periodic positive sweeps in human influenza were due to a limiting supply of beneficial mutations, or by a high supply rate with competition between beneficial mutations limiting the fixation rate.
The selective regime is important, then, for both the application and interpretation of these estimators. We present three statistics for evaluating the effects of linked selection. D1signals the presence of weak deleterious mutations that can potentially cause hitch-hiking and is based on identifying an excess of non-synonymous low-frequency polymorphisms. More sophisticated methods to characterise the DFE have been used [9, 41], but these methods rely on a number of assumptions and have given conflicting results. Consequently, it is useful to have a simple diagnostic that flags when hitch-hiking might be an issue. We have not attempted to use standard site-frequency based indicators of hitch-hiking (e.g ) which test for an excess of low and high frequency polymorphisms. As demonstrated by Kim , the excess of high frequency polymorphisms may not occur in recurrent sweeps. In addition, this effect can be complicated by clonal interference. If no comparative information is available, the excess of low frequency polymorphisms cannot be distinguished from a model of population growth .
Our second diagnostic, D2detects an excess of high-frequency non-synonymous polymorphisms signalling strong linkage effects on the population, either due to a large number of deleterious mutations causing background selection, or a large number of beneficial mutations causing clonal interference. In both cases, MK statistics are biased and estimators assuming independently segregating sites must be interpreted with care. We can distinguish between the effects of background selection and clonal interference by using a third statistic, D3. The diagnostic for clonal interference presented here follow a similar reasoning to the method used by Strelkowa and Lässig  in testing for an excess of high frequency non-synonymous polymorphisms. Our method has more general applicability as it accounts for the effect of deleterious mutations, and can be used even when it is not known which region of the sequence is positively selected.
We have applied these diagnostics to the dataset used by Strelkowa and Lässig . Our results provide further support for their conclusion that clonal interference occurs in human influenza A. The authors also raised the question of whether strong and frequent positive selection would promote the fixation of deleterious mutations, increasing the brittleness of the protein. The values obtained for D1, however, suggest that strong negative selection is predominant, so that hitch-hiking in the HA1 region is rare; this is in agreement with Shih et al. , who showed that few non-synonymous substitutions occurred outside antigenic epitopes. It is also consistent with a phylogenetic study by Nielsen and Yang  that estimated the DFE shape parameter β in that region as 0.373, indicative of low sensitivity to hitch-hiking in our model. The combination of clonal interference and robustness against hitch-hiking suggests that the modification used by Bhatt et al.  is appropriate for application to the HA1 region.
In this study, we have not considered the effect of population size changes, which are known to affect site-frequency based statistics. However, we expect D1and D3to be relatively robust, as they are based on comparisons between the non-synonymous and synonymous site frequency spectra . Population expansions, which are expected to have the strongest effect on low-frequency mutations , should have minimal effect on D2and D3. Population bottlenecks, however, can remove medium frequency polymorphisms, potentially elevating the magnitude of both D2and D3and giving false positives for clonal interference. We have also not examined the effect of selection against synonymous mutations. We expect, however, that D1and D3should not be strongly affected as long as selection against synonymous mutations is weaker than against non-synonymous mutations. D2does not use information from the synonymous site frequency spectra and should not be affected by selection against synonymous mutations, but negative values of D2may also result from background selection at synonymous sites. These effects could be considered in more detail in future simulation studies.
We have shown that linked selection is responsible for biases in the MK statistics, causing underestimation when there are high levels of interference between selected mutations, and overestimation for strong non-interfering sweeps. The statistics presented in this study can be applied to samples of protein-coding sequences to evaluate the influence of linked selection for the parameter range studied here. Values of D1that are significantly greater than zero signal susceptibility to hitch-hiking while values of D2significantly smaller than zero indicate the occurrence of multiple mergers. Multiple mergers due to clonal interference can be distinguished from background selection when D3>0.
Based on our simulations, when D2<0, D3>0 and D1≈0, we recommend using a statistic such as , which excludes both low- and high-frequency polymorphisms. On the the hand, when high values of D1(signalling hitch-hiking) are obtained, it is more appropriate to use which uses medium and high-frequency polymorphisms, accounting for change in selective constraint due to hitchhiking. In cases where D2≪0 and D3<=0, both and are expected to perform poorly.
Simulation of sequence evolution under linkage
We simulate the evolution of a population, represented as a sequence of length L=500 codons (nucleotide triplet). Each codon site is associated with a selection coefficient, s d , which is drawn from the distribution of fitness effects (DFE; seeDistribution of deleterious effects, below). The DFE affects both the extent of background selection and hitch-hiking. To model a well-adapted population, each simulation is initialised so that all non-synonymous changes from the ancestral sequence are negatively selected, reducing fitness by a factor of 1−s d . All synonymous changes are neutral. Throughout the simulation, positive selection is introduced at a specified number of sites at fixed times. After the introduction of positive selection, an individual carrying a non-synonymous change from the ancestral sequence at the positively selected site undergoes a change of fitness by a factor of 1+s b . The timing of the introduction of positive selection and the strength of selection (see Positive selection, below) control the extent of clonal interference. The extent of hitch-hiking is determined by the interaction between the DFE and positive selection.
Each simulation is initialised with a haploid population of N=104monomorphic individuals. The mutation process follows a Kimura two-parameter model , with the transition-transversion ratio fixed at κ=3. Ancestral sequences are generated randomly assuming that the base frequency of all 61 non-stop codons are equal, and all 27 one-step mutations at a codon are allowed. For κ>1, the mutation probabilities are not equal. Individuals carrying stop-codons have fitness set to zero.
In each generation, the total number of mutations introduced into the population follows a Poisson distribution with mean uNL, where the mutation rate per site per generation is u=10−6or u=10−5and occurs uniformly across all sites and all sequences. We assume non-overlapping generations and individuals reproduce by multinomial sampling with probability proportional to their fitness, as in a Wright-Fisher process.
Distribution of deleterious effects
where β is the shape parameter and is the mean selective coefficient. We consider shape parameters of β=0.25,0.5,1,2, which is similar to the range used by Charlesworth and Eyre-Walker . Estimated values in the literature range from 0.23  to 3.22 . The mean strength of selection was set at , each of which in combination with the respective β value above gives ω0≈0.1 in the presence of linkage for u=10−6.
The shape parameter β controls the proportion of weakly deleterious mutations, and therefore the extent of hitch-hiking, and in combination with u, the amount of background selection. For small values of β, the distribution of selection coefficients is broadly distributed with a larger proportion of both nearly neutral and strongly deleterious mutations; large values of β give a more strongly peaked DFE centred at nearly neutral to weakly deleterious values. Background selection is primarily mediated by the deleterious mutations that are sufficiently weakly selected that they are able to persist to appreciable frequencies but accumulate to increase the extinction probability of linked neutral and beneficial mutations. This range of selective coefficients is given approximately by 0.5<U d /s d <5 , where U d is the genomic mutation rate at selected sites. Equating U d with the genomic mutation rate gives a range of 6.7×10−5<s d <6.7×10−4for u=10−6, but U d is generally smaller than U for finite values of β. For β=0.25, less than 5% of sites lie within this range so that strong negative selection dominates and most deleterious mutations are rapidly removed from the population. For u=10−5, all mutations with 6.7×10−4<s d <6.7×10−3contribute to background selection, which covers the range around 1/N, so that much higher levels of background can be observed. Similarly, the extent of hitch-hiking is controlled by the proportion of sites with weak deleterious effects relative to the strength of positive selection, with the specific range varying according to the strength and prevalence of positive selection.
To examine the effect of linked positive selection, we introduce positive selection at a small number of codon sites in the sequence. Unlike negatively selected sites that individually have small effects but cumulatively can have a strong effect due to the large number of negatively selected sites, positive selection is expected to be rare, but a single site can have a strong effect. Thus we model all positively selected sites to have the same fixed selective effect s b .
At regular time intervals, we randomly choose a site and change the selective coefficient to s b to generate recurrent sweeps. This models a scenario of continuous positive selection, with beneficial mutations arising at different times. By varying the interval between each introduction of positive selection, we can model full selective sweeps that occur successively  or interfering sweeps . Note that unlike coalescent simulations [43, 49], we control the rate at which beneficial mutations are introduced rather than the sweep rate. The selective sweep may occur considerably later than the time at which positive selection is introduced because genetic drift, background selection and hitch-hiking can affect the time required for beneficial mutations to reach establishment.
For a low supply rate of beneficial mutations, we expect beneficial mutations to fix primarily in successive sweeps with rare occurrences of clonal interference, whereas clonal interference will occur with high probability when the supply rate of beneficial mutations is high. The expected time for a beneficial mutation to become established in the population is given by t est =1/(u NL b s b ) ; after establishment, the beneficial mutation behaves almost deterministically, increasing rapidly in frequency and is expected to fix in generations . For population size N=104and u=10−6, a single beneficial mutation of strength s b =0.01 is expected to have establishment and fixation times of t est ≈2857 and t fix ≈460 generations. To obtain a high supply rate of beneficial mutations, we introduce positive selection at high frequency, specifically at one site in every τ=1000 generations, which is faster than the rate of establishment. For a low supply rate of beneficial mutations, we set τ=10000 generations, so that establishment and fixation of one beneficial mutation is likely to occur before a second positively selected site is introduced. Note that varying the timing of positive selection controls the supply rate of beneficial mutations (generally parameterised as U b N=uL b N) indirectly. After positive selection is introduced at a site, L b is increased by one; however, L b is also decreased when a beneficial mutation reaches fixation.
Simulations under independently segregating sites
To compare sequence statistics obtained under complete linkage with those obtained under the assumption of independently segregating sites, we simulate the number of polymorphic and divergent sites according to the Poisson Random Field (PRF) model . The PRF model assumes a Wright-Fisher population at equilibrium with an infinite number of sites so that all new mutations occur on distinct sites. Under these assumptions, Sawyer and Hartl  showed that number of sites carrying a derived mutation follows a Poisson random field, with expectations that are functions of the mutation and selection parameters. We use the PRF as it is the basis of a number of inference methods [6, 9, 10, 41], and therefore provides a better reference than a finite-site model with independently segregating sites.
In the PRF framework , the number of derived sites can be simulated as independent Poisson variables. We can then use the number of divergent and polymorphic sites to calculate sequence statistics , and as described in the main text. In the following section, we give the equations used to calculate the mean number of divergent and polymorphic sites.
denotes the Hurwitz zeta function which is provided in the GNU scientific library . When L b >0 sites are positively selected, we generate the number of divergent non-synonymous sites over the deleterious portion of the sequence using Equation (7) and the number of divergent beneficial sites is generated from a truncated Poisson distribution with mean u A L b t ω(s b ,N), capped at L b . This allows comparison with the finite sites model, which explicitly does not allow recurrent positive selection at a single site.
denotes the incomplete beta function.
This is obtained by taking the ratio between fixation probabilities of a selected and a neutral mutation . In the case where positive selection is not operating, the value of ω summed across the entire sequence is equal to the proportion of effectively neutral sites, denoted f.
The MK statistic does not require a scaling factor c, as it is given in units of the number of non-synonymous substitutions. The offset (+1) term in the denominator means that this estimator is defined in all cases. Smith and Eyre-Walker  found that the offset does not introduce noticeable bias.
The ratio in Equation (21) is an estimator of f, under the assumption that all segregating polymorphisms are selectively neutral. This assumption is valid when selection is strong so that selected mutations immediately reach fixation or extinction, but not when weak selection is frequent. This problem is further compounded in the context of linked selection as linkage has the effect of weakening the effective strength of selection so that both deleterious and beneficial mutations can potentially segregate for longer prior to extinction or fixation. Here, we examine two modifications of the MK statistic.
In fact, is not a true estimator as ω0is a fixed value (treating f as known) rather than a measurable quantity. Here, ω0is obtained using the median value of based on simulations with linkage and the same values of β and but no positive selection (ω0=0.09,0.09,0.11,0.12 for u=10−6and ω0=0.10,0.13,0.23,0.33 for u=10−5). We used simulations rather than the theoretical expectation of f to account for background selection. In practice, ω0cannot be estimated from divergence information unless there is a period where it is known positive selection has not occurred. However, we use , as it provides a comparison showing how and differ in their estimation of f.
This work was supported by a Discovery Grant DP110100465 from the Australian Research Council and by the Australian Postgraduate Award to CHSC. We also thank two anonymous reviewers for their constructive comments.
- Nielsen R: Molecular signatures of natural selection. Annu Rev Genet. 2005, 39: 197-218. 10.1146/annurev.genet.39.073003.112420.PubMedView ArticleGoogle Scholar
- Hughes AL, Nei M: Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988, 335 (6186): 167-170. 10.1038/335167a0.PubMedView ArticleGoogle Scholar
- Hughes AL: Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level. Heredity. 2007, 99 (4): 364-373. 10.1038/sj.hdy.6801031.PubMedView ArticleGoogle Scholar
- McDonald JH, Kreitman M: Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991, 351 (6328): 652-654. 10.1038/351652a0.PubMedView ArticleGoogle Scholar
- Eyre-Walker A, Keightley PD: The distribution of fitness effects of new mutations. Nat Rev Genet. 2007, 8 (8): 610-618.PubMedView ArticleGoogle Scholar
- Smith NGC, Eyre-Walker A: Adaptive protein evolution in Drosophila. Nature. 2002, 415 (6875): 1022-1024. 10.1038/4151022a.PubMedView ArticleGoogle Scholar
- Toll-Riera M, Laurie S, Albà M M: Lineage-specific variation in intensity of natural selection in mammals. Mol Biol Evol. 2011, 28: 383-398. 10.1093/molbev/msq206.PubMedView ArticleGoogle Scholar
- Czelusniak J, Goodman M, Hewett-Emmett D, Weiss ML, Venta PJ, Tashian RE: Phylogenetic origins and adaptive evolution of avian and mammalian haemoglobin genes. Nature. 1982, 298 (5871): 297-300. 10.1038/298297a0.PubMedView ArticleGoogle Scholar
- Nielsen R, Yang Z: Estimating the distribution of selection coefficients from phylogenetic data with applications to mitochondrial and viral DNA. Mol Biol Evol. 2003, 20 (8): 1231-1239. 10.1093/molbev/msg147.PubMedView ArticleGoogle Scholar
- Sawyer SA, Hartl DL: Population genetics of polymorphism and divergence. Genetics. 1992, 132 (4): 1161-1176.PubMed CentralPubMedGoogle Scholar
- Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15 (5): 568-573. 10.1093/oxfordjournals.molbev.a025957.PubMedView ArticleGoogle Scholar
- Lemey P, Kosakovsky Pond SL, Drummond AJ, Pybus OG, Shapiro B, Barroso H, Taveira N, Rambaut A: Synonymous substitution rates predict HIV disease progression as a result of underlying replication dynamics. PLoS Comput Biol. 2007, 3 (2): e29-10.1371/journal.pcbi.0030029.PubMed CentralPubMedView ArticleGoogle Scholar
- Charlesworth B, Morgan MT, Charlesworth D: The effect of deleterious mutations on neutral molecular variation. Genetics. 1993, 134: 1289-1303.PubMed CentralPubMedGoogle Scholar
- Birky CW, Walsh JB: Effects of linkage on rates of molecular evolution. Proc Natl Acad Sci U S A. 1988, 85 (17): 6414-6418. 10.1073/pnas.85.17.6414.PubMed CentralPubMedView ArticleGoogle Scholar
- Peck JR: A ruby in the rubbish: beneficial mutations, deleterious mutations and the evolution of sex. Genetics. 1994, 137 (2): 597-606.PubMed CentralPubMedGoogle Scholar
- Charlesworth B: The effect of background selection against deleterious mutations on weakly selected, linked variants. Gent Res. 1994, 63 (03): 213-227. 10.1017/S0016672300032365.View ArticleGoogle Scholar
- Charlesworth D, Charlesworth B, Morgan MT: The pattern of neutral molecular variation under the background selection model. Genetics. 1995, 141 (4): 1619-PubMed CentralPubMedGoogle Scholar
- Neher RA, Hallatschek O: Genealogies of rapidly adapting populations. Proc Natl Acad Sci U S A. 2013, 110 (2): 437-442. 10.1073/pnas.1213113110.PubMed CentralPubMedView ArticleGoogle Scholar
- Neher RA, Shraiman BI: Genetic draft and quasi-neutrality in large facultatively sexual populations. Genetics. 2011, 188 (4): 975-996. 10.1534/genetics.111.128876.PubMed CentralPubMedView ArticleGoogle Scholar
- Bustamante CD, Wakeley J, Sawyer S, Hartl DL: Directional selection and the site-frequency spectrum. Genetics. 2001, 159 (4): 1779-1788.PubMed CentralPubMedGoogle Scholar
- Walczak AM, Nicolaisen LE, Plotkin JB, Desai MM: The structure of genealogies in the presence of purifying selection: a fitness-class coalescent. Genetics. 2012, 190 (2): 753-779. 10.1534/genetics.111.134544.PubMed CentralPubMedView ArticleGoogle Scholar
- Maynard-Smith J, Haigh J, et al: The hitch-hiking effect of a favourable gene. Genet Res. 1974, 23: 23-35. 10.1017/S0016672300014634.View ArticleGoogle Scholar
- Braverman JM, Hudson RR, Kaplan NL, Langley CH, Stephan W: The hitchhiking effect on the site frequency spectrum of DNA polymorphisms. Genetics. 1995, 140 (2): 783-796.PubMed CentralPubMedGoogle Scholar
- Kim Y, Stephan W: Joint effects of genetic hitchhiking and background selection on neutral variation. Genetics. 2000, 155 (3): 1415-1427.PubMed CentralPubMedGoogle Scholar
- Kim Y, Wiehe T: Simulation of DNA sequence evolution under models of recent directional selection. Brief Bioinform. 2009, 10: 84-96.PubMed CentralPubMedView ArticleGoogle Scholar
- Li W: Models of nearly neutral mutations with particular implications for nonrandom usage of synonymous codons. J Mol Evol. 1987, 24 (4): 337-345. 10.1007/BF02134132.PubMedView ArticleGoogle Scholar
- Comeron JM, Kreitman M: Population, evolutionary and genomic consequences of interference selection. Genetics. 2002, 161: 389-410.PubMed CentralPubMedGoogle Scholar
- Zeng K, Charlesworth B: The effects of demography and linkage on the estimation of selection and mutation parameters. Genetics. 2010, 186 (4): 1411-1424. 10.1534/genetics.110.122150.PubMed CentralPubMedView ArticleGoogle Scholar
- Miralles R, Gerrish PJ, Moya A, Elena SF: Clonal interference and the evolution of RNA viruses. Science. 1999, 285 (5434): 1745-1747. 10.1126/science.285.5434.1745.PubMedView ArticleGoogle Scholar
- Rozen DE, de Visser JG, Gerrish PJ: Fitness effects of fixed beneficial mutations in microbial mopulations. Curr Biol. 2002, 12 (12): 1040-1045. 10.1016/S0960-9822(02)00896-5.PubMedView ArticleGoogle Scholar
- Desai MM, Walczak AM, Fisher DS: Genetic diversity and the structure of genealogies in rapidly adapting populations. Genetics. 2013, 193 (2): 565-585. 10.1534/genetics.112.147157.PubMed CentralPubMedView ArticleGoogle Scholar
- Messer PW, Petrov DA: Frequent adaptation and the McDonald-Kreitman test. Proc Natl Acad Sci U S A. 2013, 110 (21): 8615-8620. 10.1073/pnas.1220835110.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhatt S, Holmes EC, Pybus OG: The genomic rate of molecular adaptation of the human influenza A virus. Mol Biol Evol. 2443, 28 (9):Google Scholar
- Charlesworth J, Eyre-Walker A: The McDonald-Kreitman test and slightly deleterious mutations. Mol Biol Evol. 2008, 25 (6): 1007-1015. 10.1093/molbev/msn005.PubMedView ArticleGoogle Scholar
- Efron B: Better bootstrap confidence intervals. J Am Stat Assoc. 1987, 82 (397): 171-185. 10.1080/01621459.1987.10478410.View ArticleGoogle Scholar
- Strelkowa N, Lässig M: Clonal interference in the evolution of influenza. Genetics. 2012, 192 (2): 671-682. 10.1534/genetics.112.143396.PubMed CentralPubMedView ArticleGoogle Scholar
- Fay JC: Weighing the evidence for adaptation at the molecular level. Trends Genet. 2011, 27 (9): 343-349. 10.1016/j.tig.2011.06.003.PubMed CentralPubMedView ArticleGoogle Scholar
- Seger J, Smith WA, Perry JJ, Hunn J, Kaliszewska ZA, Sala LL, Pozzi L, Rowntree VJ, Adler FR: Gene genealogies strongly distorted by weakly interfering mutations in constant environments. Genetics. 2010, 184 (2): 529-545. 10.1534/genetics.109.103556.PubMed CentralPubMedView ArticleGoogle Scholar
- Sanjuán R, Moya A, Elena SF: The distribution of fitness effects caused by single-nucleotide substitutions in an RNA virus. Proc Natl Acad Sci U S A. 2004, 101 (22): 8396-8401. 10.1073/pnas.0400146101.PubMed CentralPubMedView ArticleGoogle Scholar
- Koelle K, Cobey S, Grenfell B, Pascual M: Epochal evolution shapes the phylodynamics of interpandemic influenza A (H3N2) in humans. Science. 2006, 314 (5807): 1898-1903. 10.1126/science.1132745.PubMedView ArticleGoogle Scholar
- Piganeau G, Eyre-Walker A: Estimating the distribution of fitness effects from DNA sequence data: Implications for the molecular clock. Proc Natl Acad Sci U S A. 2003, 100 (18): 10335-10340. 10.1073/pnas.1833064100.PubMed CentralPubMedView ArticleGoogle Scholar
- Fay JC, Wu CI: Hitchhiking under positive Darwinian selection. Genetics. 2000, 155 (3): 1405-1413.PubMed CentralPubMedGoogle Scholar
- Kim Y: Allele frequency distribution under recurrent selective sweeps. Genetics. 2006, 172 (3): 1967-1978.PubMed CentralPubMedView ArticleGoogle Scholar
- Hahn MW, Rausher MD, Cunningham CW: Distinguishing between selection and population expansion in an experimental lineage of bacteriophage T7. Genetics. 2002, 161: 11-20.PubMed CentralPubMedGoogle Scholar
- Shih ACC, Hsiao TC, Ho MS, Li WH: Simultaneous amino acid substitutions at antigenic sites drive influenza A hemagglutinin evolution. Proc Natl Acad Sci U S A. 2007, 104 (15): 6283-6288. 10.1073/pnas.0701396104.PubMed CentralPubMedView ArticleGoogle Scholar
- Maruyama T, Fuerst PA: Population bottlenecks and nonequilibrium models in population genetics. II. Number of alleles in a small population that was formed by a recent bottleneck. Genetics. 1985, 111 (3): 675-689.PubMed CentralPubMedGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16 (2): 111-120. 10.1007/BF01731581.PubMedView ArticleGoogle Scholar
- Eyre-Walker A, Woolfit M, Phelps T: The distribution of fitness effects of new deleterious amino acid mutations in humans. Genetics. 2006, 173 (2): 891-900. 10.1534/genetics.106.057570.PubMed CentralPubMedView ArticleGoogle Scholar
- Coop G, Ralph P: Patterns of neutral diversity under general models of selective sweeps. Genetics. 2012, 192: 205-224. 10.1534/genetics.112.141861.PubMed CentralPubMedView ArticleGoogle Scholar
- Desai MM, Fisher DS: Beneficial mutation-selection balance and the effect of linkage on positive selection. Genetics. 2007, 176 (3): 1759-1798. 10.1534/genetics.106.067678.PubMed CentralPubMedView ArticleGoogle Scholar
- Galassi M, Davies J, Theiler J, Gough B, Jungman G, Booth M, Rossi F: GNU Scientific Library Reference Manual. Bristol: Network Theory Limited. 2003Google Scholar
- Wright S: The distribution of gene frequencies under irreversible mutation. Proc Natl Acad Sci U S A. 1938, 24 (7): 253-10.1073/pnas.24.7.253.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang Z: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Kimura M: On the probability of fixation of mutant genes in a population. Genetics. 1962, 47 (6): 713-719.PubMed CentralPubMedGoogle Scholar
- Fay JC, Wyckoff GJ, Wu CI: Positive and negative selection on the human genome. Genetics. 2001, 158 (3): 1227-1234.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.