Hitchhiking Under Positive Darwinian Selection

Corresponding author: Chung-I Wu, Department of Ecology and Evolution, University of Chicago, 1101 E. 57th St., Chicago, IL 60637., ciwu{at}uchicago.edu (E-mail)

Communicating editor: W. STEPHAN

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

Positive selection can be inferred from its effect on linked neutral variation. In the restrictive case when there is no recombination, all linked variation is removed. If recombination is present but rare, both deterministic and stochastic models of positive selection show that linked variation hitchhikes to either low or high frequencies. While the frequency distribution of variation can be influenced by a number of evolutionary processes, an excess of derived variants at high frequency is a unique pattern produced by hitchhiking (derived refers to the nonancestral state as determined from an outgroup). We adopt a statistic, H, to measure an excess of high compared to intermediate frequency variants. Only a few high-frequency variants are needed to detect hitchhiking since not many are expected under neutrality. This is of particular utility in regions of low recombination where there is not much variation and in regions of normal or high recombination, where the hitchhiking effect can be limited to a small (<1 kb) region. Application of the H test to published surveys of Drosophila variation reveals an excess of high frequency variants that are likely to have been influenced by positive selection.

THE extent to which positive Darwinian selection shapes molecular evolution has been a dominant issue in the last three decades (KIMURA 1983 ; NEI 1987 ). Under neutrality, the amount of DNA variation within a species is proportional to the amount of divergence between species (KIMURA 1983 ). In contrast, the spread of a positively selected mutation through a population removes linked neutral variation without affecting neutral divergence (MAYNARD SMITH and HAIGH 1974 ). Surveys of DNA variation have revealed a correlation between levels of variation and rates of recombination in plants (STEPHAN and LANGLEY 1998 ), flies (STEPHAN and LANGLEY 1989 ; BEGUN and AQUADRO 1992 ; MARTIN-CAMPOS et al. 1992 ; STEPHAN and MITCHELL 1992 ), mice (NACHMAN 1997 ), and humans (NACHMAN et al. 1998 ). Reduced levels of variation are also found within small regions of some genes (MCDONALD 1998 ). However, there are a number of alternative explanations for reduced levels of variation. In particular, "background selection," or the elimination of deleterious mutations segregating in a population, predicts low levels of variation in regions of low recombination (CHARLESWORTH et al. 1993 ). Distinguishing positive selection from these alternative explanations has become a significant empirical and theoretical pursuit.

The frequency spectrum of variation provides a means of detecting positive selection independent of levels of variation. In the absence of recombination, hitchhiking eliminates all linked variation and a population in recovery is characterized by an excess of new mutations at low frequency. In the presence of recombination, hitchhiking is incomplete since not all variation is removed and its direct effect on the frequency spectrum is not known. Unless the hitchhiking effect is very strong, even intragenic variation will experience an incomplete rather than complete hitchhiking effect. For example, a selection coefficient of 10^-3 and a recombination rate of 10^-8/bp (about the average in Drosophila melanogaster; ASHBURNER 1989 ) would still leave 10% of the heterozygosity 500 bp away from a site under selection (STEPHAN et al. 1992 ). In a region where recombination is reduced 10-fold, such as the tip of the X chromosome (AGUADE et al. 1989 ), the corresponding distance would be 5 kb.

The incomplete removal of variation may create a pattern of variation very different from the pattern produced after all variation is removed. For example, in the case of a brief reduction in population size, there is a striking difference between residual patterns of variation found immediately after a bottleneck and the excess of new mutations that subsequently accumulate (TAJIMA 1989A ; FAY and WU 1999 ). To understand the effect of positive selection on linked variation, we describe the hitchhiking effect on the frequency spectrum as a function of the rate of recombination. We then adopt a statistical test that can detect hitchhiking even when there are few segregating sites. Application of the test to published surveys of variation in regions of low and high recombination in Drosophila indicates a departure from neutrality in the direction predicted by hitchhiking.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

Polymorphism data:
We restricted our analysis to published polymorphism surveys of Drosophila species. Recombination rates across the D. melanogaster genome have been estimated from the frequency of recombination per cytological band (COMERON et al. 1999 ). Regions of low recombination were designated as cytological divisions with rates of recombination < 5 x 10^-9/bp: 1A–2D, 20C–20F, 21A, 38A–44C, 60C–60F, 61A–61B, and 75F–84F (excluding the fourth chromosome where recombination is absent). This cutoff value is ~10 times lower than the division with the highest estimated rate of recombination. There are six published polymorphism surveys in the designated regions: achaete (ac, 1B1; MARTIN-CAMPOS et al. 1992 ), asense (ase, 1B4; HILTON et al. 1994 ), Phosphogluconate dehydrogenase (Pgd, 2D6; BEGUN and AQUADRO 1994 ), suppressor of forked [su(f), 20E; LANGLEY et al. 1993], fushi tarazu (ftz, 84B1; JENKINS et al. 1995 ), and Glucose dehydrogenase (Gld, 84D1; HAMBLIN and AQUADRO 1997 ). The Accessory gland protein 26Aa (Acp26Aa) is in a region of high recombination in D. melanogaster (26A, c estimated to be 5 x 10^-8/bp) and its intra- and interspecific variation has been extensively studied (AGUADE et al. 1992 ; TSAUR and WU 1997 ; AGUADE 1998 ; TSAUR et al. 1998 ). vermilion (v; STEPHAN and LANGLEY 1989 ; STEPHAN et al. 1998 ) is in a region of reduced recombination in D. ananassae.

The frequency spectrum of derived variants, including synonymous, nonsynonymous, and noncoding variants, was constructed for each survey using an outgroup. When multiple locations were sampled, the data were pooled for the determination of the frequency spectrum. For achaete, the European sample was scaled from 192 to 50, which is similar in size to the African and North American samples. All data are DNA polymorphism within D. melanogaster and divergence from D. simulans, except for achaete [which is a four-cutter restriction fragment length polymorphism (RFLP) study] and vermilion, which is a survey of polymorphism in D. ananassae with divergence from D. pallidosa.

Simulation methods:
A standard coalescence algorithm was used to generate neutral genealogies (HUDSON 1990 ). A modified coalescence algorithm, which tracks the genealogy of a neutral locus linked to an advantageous mutation, was used to generate hitchhiking genealogies. The algorithm is a simplification of that detailed in BRAVERMAN et al. 1995 , such that the recombination rate (c), selection coefficient (s), and time at which an advantageous mutation arises are set parameters instead of random variables. The selected allele confers a fitness of 1 + s in the heterozygous state and 1 + 2s in the homozygous state and decreases deterministically from a frequency of 1 - 1/2N to 1/2N in a population size (N) of 10⁶ individuals. Recombination events occur between the selected and neutral locus, and coalescent events occur between advantageous alleles or between nonadvantageous alleles, depending on the frequency of the advantageous mutation and the number of advantageous and nonadvantageous alleles. Once all alleles have coalesced, mutations are randomly distributed on the genealogy of the neutral locus. All results were generated from at least 1000 iterations of the algorithm.

Statistics:
Tajima's D statistic (TAJIMA 1989B ) is based on the difference between two commonly used estimators of = 4Nµ, where N is the effective population size and µ is the mutation rate: is calculated from average heterozygosity (TAJIMA 1983 ) and _W is calculated from the number of segregating sites (WATTERSON 1975 ),

(1)

(2)

where S_i is the number of derived variants found i times in a sample of n chromosomes. Let _H be an estimator of weighted by the homozygosity of the derived variants, as opposed to the ancestral variants,

(3)

where _H is an unbiased estimator of and is a special case, notated as L(-1), of a general class of recently derived estimators that vary in their weighting schemes (FU 1995 ). Let H be the difference between and _H. The null distributions of the D and H statistics were generated using 1000 iterations of a neutral coalescence algorithm conditioning on the observed number of segregating sites. P values were calculated as the probability of a neutral D or H value being less than the observed D or H value. The D and H statistics are probably conservative since critical values were generated in the absence of recombination (WALL 1999 ), unless otherwise noted. The power of the D and H statistics was calculated as the probability of rejecting a neutral genealogy given a hitchhiking genealogy. The critical D and H values ( = 0.05, one-sided) were generated for any given number of segregating sites. For calculating the power of the H statistic, the probability of misinference was 0.00375/site (see below).

An outgroup was used to infer the derived and ancestral states for all polymorphism data analyzed. However, a backmutation would result in the incorrect inference of the derived polymorphic state. As compensation, the probability of misinference was incorporated into the null distribution of the H statistic by exchanging the frequency of the derived and ancestral state, with probability equal to that of misinference, for each segregating site. For nucleotide sequences, the probability of a backmutation is d/3, where d is the net divergence or the average number of fixed differences between the two species. The observation that transitions occur at twice the rate of transversions (MORIYAMA and POWELL 1996 ) was incorporated and resulted in a probability of backmutation equal to 3d/8 (see Appendix at http://home.uchicago.edu/~jfay/appendix.html). The probability of misinference was calculated separately for synonymous, nonsynonymous, and noncoding sites and then weighted by the number of synonymous, nonsynonymous, and noncoding polymorphic sites. For restriction site data, where the outgroup sequence and one of the ingroup sequences are known, the probability of misinference differs, depending on the state of the outgroup (presence or absence of the restriction site). If the site is present in the outgroup, then the derived state is inferred to be a loss and the probability of misinference is d/(3 - 2d), whereas if the site is absent in the outgroup the probability of misinference is 4d/(1 + 3d) (see Appendix at http://home.uchicago.edu/~jfay/appendix.html). The probability of misinference was calculated assuming that only one backmutation could account for the observed data.

	THEORY

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

Deterministic:
The hitchhiking effect can be described by the change in allele frequency at a neutral polymorphic locus (B, b) due to the spread of a linked advantageous mutation, A, through a population. Clearly, B will increase or decrease in frequency depending on whether the B or b allele is originally linked to the advantageous mutation. In the absence of recombination, the variant originally linked to the advantageous mutation is fixed, but when recombination is present but rare, all variants that remain segregating are most likely present at either very low or high frequencies.

The frequency spectrum after hitchhiking can be found, knowing the frequency of variations before selection and their change in frequency due to selection. To distinguish between low- and high-frequency variants, ancestral (old) and derived (new) states must be distinguished (this can be done empirically using an outgroup). Let (x) be the frequency spectrum of derived variants in a population. Under neutrality, the expected number of sites where the derived variant is between a frequency of x and x + dx is given by

(4)

In Equation 4, = 4Nµ, where N is the effective population size and µ is the mutation rate of the region (WATTERSON 1975 ). Deterministically, if a neutral linked variant, B, hitchhikes with A, its frequency is transformed from x to 1 - c* + xc*; otherwise it is reduced to xc*, where c* is a scaled measure of recombination. c* is given in its exact form in MAYNARD SMITH and HAIGH 1974 but can be approximated by (c/s)ln(1/p₀), where s is the haploid selection coefficient, c is the rate of recombination between the two sites, and p₀ is the initial frequency of the advantageous mutation. The probability that B hitchhikes with A is equal to its initial frequency, x, so a uniform number of derived variants over the entire frequency spectrum, x(x)dx = dx, are transformed to high frequencies, while the rest, (/x - )dx, are transformed to low frequencies (inset of Fig 1). The frequency spectrum after hitchhiking can be approximated after a linear transformation of expectations (subscripted L and H are for low and high frequencies, respectively),

(5)

(6)

and zero otherwise (inset of Fig 1).

View larger version (21K): In this window In a new window Download PPT slide   Figure 1. The frequency spectrum of derived variants expected under neutrality (open areas) and subsequent to a hitchhiking event (shaded areas) with parameters = 10, n = 100, = 0.01, and s = 10-3 for 1000 iterations of a coalescence algorithm. The inset shows the proportion of variants in a population that are expected to increase or decrease in frequency during a hitchhiking event and the resulting frequency spectrum found under the deterministic model.

Stochastic:
Coalescent simulations of single hitchhiking events, where the advantageous mutation has just reached fixation, produce a skew in the frequency spectrum similar to that predicted by the deterministic theory (Fig 1). The hitchhiking algorithm incorporates stochastic mutation and recombination events but assumes that the advantageous allele follows a deterministic increase in frequency (MATERIALS AND METHODS). An exact treatment of the spread of an advantageous mutation through the population incorporates stochastic fluctuations in the frequency of the selected allele when it is at low or high frequencies but would produce similar results so long as selection is stronger than recombination (BARTON 1998 ).

The strength of hitchhiking is determined by the ratio of the recombination rate to the selection coefficient. As the strength of the hitchhiking effect increases, the ratio of low- to intermediate- to high-frequency variants becomes skewed toward low and high frequencies (Table 1). The frequency spectrum approaches that predicted under neutrality as the strength of hitchhiking decreases. Most importantly, there is a striking difference between a genealogy produced by complete (no recombination events) and incomplete hitchhiking. In the extreme case when all samples except one have coalesced during hitchhiking, nearly all mutations are present at a frequency of either (n - 1)/n or 1/n, with equal probability. In contrast, when there is no recombination a star genealogy is produced and any new mutation would be found once in a sample (Fig 2). Therefore, when a neutral variant is tightly, but not completely, linked to a site that has just been selectively driven to fixation, the salient feature of hitchhiking is a bipartite spectrum of allele frequencies. Subsequent to a hitchhiking event, variation will be regained, first as low-frequency variants, and will eventually return to equilibrium.

View larger version (10K): In this window In a new window Download PPT slide   Figure 2. Genealogies produced by positive selection (A) in the absence of recombination and (B) when there is a single recombinant in a sample. All variation is lost in the genealogy of A but in B mutations a and b produce low-frequency variants and mutations c and d produce high-frequency variants. Dotted lines indicate lineages lost due to hitchhiking.

Statistics:
A standard statistical test for hitchhiking is Tajima's D statistic, which compares the number of rare to intermediate-frequency variants (TAJIMA 1989B ). Tajima's D is the standardized difference between two commonly used measures of variability: _W, which is based on the total number of segregating sites and is influenced most by low-frequency variation (WATTERSON 1975 ), and , which is based on average heterozygosity and is influenced most by intermediate-frequency variants (TAJIMA 1989B ). While the D statistic should be useful in the detection of a partial hitchhiking event, many other processes, such as background selection or a change in population size (CHARLESWORTH et al. 1993 ), may generate D values that depart from neutrality due to an excess of low-frequency variants. On the other hand, strong deterministic forces are needed to drive variants to high frequencies. A statistic based on an excess of high-frequency variants would thus be uniquely sensitive to hitchhiking.

We define the test statistic, H, as the difference between , which is influenced most by variants at intermediate frequencies, and _H (MATERIALS AND METHODS), which is influenced most by high-frequency variants. This test is analogous to and was motivated by the F(0, -1) test described in FU 1995 , the only difference being F(0, -1) is the difference between and _H divided by the variance of the difference. Under neutrality the expected difference between two estimators of is zero, but following a hitchhiking event _H and _W should be >. This prediction is met over a range of c/s values (Fig 3).

View larger version (28K): In this window In a new window Download PPT slide   Figure 3. (A) The mean for each of the three different estimators of as a function of recombination rate where each point is calculated from 5000 coalescence simulations for s = 0.005 and = 5 (MATERIALS AND METHODS). In the absence of hitchhiking, the means of H, W, and from the simulations are 4.94, 4.99, and 4.97, respectively. The number of recombinants is the number of alleles in a sample of 50 that escaped coalescence during hitchhiking by recombining onto an advantageous chromosome. (B) The corresponding power of the D and H statistics. The solid squares indicate the proportion of cases where both statistics reject neutrality and the two statistics overlap.

The power of both the D and H tests in rejecting neutrality is similar and greatest for intermediate values of c/s. The power of these statistics compared with the percentage of simulated hitchhiking events where both the D and H tests reject neutrality shows that the two tests are complementary in their ability to detect a single hitchhiking event. For example, when = 0.01, 51.7 and 61.1% of the hitchhiking simulations were rejected by the D and H tests, respectively, while 38.6% were rejected by both. Both D and H reject neutrality when there is an excess of low- and high-frequency variants, but only D can reject neutrality when there is an excess of low- but not high-frequency variants and vice versa for H. While Fig 3 gives an indication of slight differences between the power of the two tests, this difference should be interpreted with caution since it depends on a number of parameters not explored here, such as the number of recombination events, level of variation before hitchhiking, and time since the hitchhiking event. For instance, when one recombinant is in a sample and there are few segregating sites, only very low- or high-frequency variants would be expected since there is nearly an equal probability that a derived variant is found at low (a and b) or high (c and d) frequency (Fig 2). Under this scenario the H test but not Tajima's D may be significant (Table 2). Regardless of the relative powers of the two tests, it is the excess of high- but not low-frequency variants that is the unique signature of positive selection.

	APPLICATION TO DATA

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

In regions of low recombination, the H test may distinguish hitchhiking from a neutral or background selection model, which is expected to have the same relative abundance of intermediate- to high-frequency variants (CHARLESWORTH et al. 1995 ; FU 1997 ). We identified seven published polymorphism surveys in four regions with greatly reduced rates of recombination (MATERIALS AND METHODS). Achaete, asense, and Pgd are located on the tip of the X chromosome and a significant excess of high-frequency variants is found at achaete (Table 2). An excess of high-frequency variants is not found at su(f) at the base of the X chromosome or at ftz or Gld, near the centromere of the third chromosome. In D. ananassae, there is a nonsignificant excess of high-frequency variants at vermilion near the centromere of the X chromosome.

The levels of variation and excess of high-frequency variants found at achaete, asense, and Pgd suggest a weakening hitchhiking effect with their distance from the tip of the X chromosome. Achaete has a significant excess of high-frequency variants and the largest reduction in levels of variation. A weaker hitchhiking pattern is found at asense and it is weaker still at Pgd. Out of eight comparisons one significant association was found between variation at achaete and Pgd (BEGUN and AQUADRO 1991 ). This is consistent with the hitchhiking effect, which produces linkage disequilibrium (THOMSON 1977 ). In the absence of hitchhiking, no genealogical association would be expected between loci separated by a recombinational distance, 4Nc, of >100 (HUDSON 1983 ). The three loci should be independent, given that the rate of recombination between achaete and asense is estimated to be 1.5 x 10^-4 and between achaete and Pgd, 6 x 10^-3, approximately (FLYBASE 1999 ). Treating the three surveys as independent, the three P values can be combined (SOKAL and ROHLF 1981 ) to show a significant departure from neutrality for the entire ac-ase-Pgd region (P = 0.022).

In regions of normal recombination, the hitchhiking effect may be limited to a small region so the location of the sequence under selection must be known. The accessory gland protein gene (Acp26Aa) has been shown to be under positive selection by interspecific comparisons (AGUADE et al. 1992 ; TSAUR and WU 1997 ; AGUADE 1998 ; TSAUR et al. 1998 ). On the basis of the expected ratio of synonymous and nonsynonymous differences within and between species (MCDONALD and KREITMAN 1991 ), only 29 nonsynonymous substitutions are expected between D. melanogaster and D. simulans, whereas 75 are observed (TSAUR et al. 1998 ). Thus, it is likely that a hitchhiking event occurred recently enough to influence patterns of intraspecific variation in D. melanogaster.

As previously noted (TSAUR et al. 1998 ), there is an excess of high-frequency variants in the Acp26Aa gene (Fig 4). Because the locus is in a region of high recombination, c is estimated to be 4.4 x 10^-8 (COMERON et al. 1999 ), and there is ample evidence of recombination between sites as close as 20 bp apart, a hitchhiking effect could be limited to a small region within the gene and each segregating site would experience a different strength of the hitchhiking effect. A sliding window plot of polymorphism vs. divergence shows an excess of high-frequency variants on both sides of a region where there is an absence of variation when compared to divergence (Fig 5). This is the expected pattern produced by positive selection over a continuous range of linked variation (Fig 3). Of the 13 variants found at a frequency of >50%, all lie within the region affected by hitchhiking and six lie on one side and seven lie on the other side of the region presumably affected by complete hitchhiking, which is indicated by the arrows. No such pattern is found in the Acp26Ab gene.

View larger version (23K): In this window In a new window Download PPT slide   Figure 4. The observed () and expected () numbers of derived variants that occur in a worldwide sample of 101 Acp26Aa D. melanogaster genes (AGUADE et al. 1992 ; AGUADE 1998 ; TSAUR et al. 1998 ).

View larger version (24K): In this window In a new window Download PPT slide   Figure 5. A sliding window plot of polymorphism and divergence (d, thick line) between D. melanogaster line NC1 and D. simulans (AGUADE et al. 1992 ). Window size is 100 bp and step size is 10 bp. On the abscissa, solid boxes are exons and shaded boxes are introns. The two double arrows indicate the size of the region affected by hitchhiking: the smaller one (~350 bp) has lost nearly all variation and the larger one (~700 bp) encompasses the region where there is an excess of high-frequency variants.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

In addition to reducing levels of variation (MAYNARD SMITH and HAIGH 1974 ), positive selection skews the frequency spectrum of linked variation to low and high frequencies. This hitchhiking effect is a function of recombinational distance from the site under selection and time since the hitchhiking event(s) (KAPLAN et al. 1989 ).

Recombination:
Hitchhiking is complete or incomplete in the removal of variation depending on whether or not recombination has occurred. The size of the region affected by partial hitchhiking is at least two orders of magnitude larger than that affected by complete hitchhiking. For a c/s value of 2 x 10^-4, hitchhiking is nearly complete, heterozygosity is reduced by 99%, and there is at least 1 recombinant out of a sample of 50 in 10% of hitchhiking events (Fig 3). In contrast, for a c/s value of 0.02, hitchhiking is partial and heterozygosity is reduced by 57%. For humans and flies, where the rate of recombination between adjacent nucleotides is on the order of 10^-8 (ASHBURNER 1989 ; BROMAN et al. 1998 ), complete hitchhiking ( = 2 x 10^-4) is limited to a region of 20 (s = 0.001) or 200 bp (s = 0.01).

For incomplete hitchhiking, the degree to which variation is reduced corresponds to the degree to which the frequency spectrum is skewed. This is the case for achaete, asense, and Pgd, which span a 2-Mb region on the tip of the X chromosome (FLYBASE 1999 ), and for the second exon of the Acp26Aa gene, where nearly all variation is eliminated from an ~350-bp region and the entire hitchhiking pattern is limited to an ~700-bp region (Fig 5). Taking c to be 5 x 10^-8/bp (COMERON et al. 1999 ) at the Acp26Aa gene, the selection coefficient can be estimated from the size of the region affected by complete hitchhiking. Assuming hitchhiking is complete (c/s < 2 x 10^-3), ~175 bp away from the site under selection the selection coefficient is at least 0.004. The actual selection coefficient may be even larger since 4Nc estimated from polymorphism is 0.0725/bp (HUDSON et al. 1987 ; ROZAS and ROZAS 1999 ) and hitchhiking produces linkage disequilibrium. In addition, population subdivision may weaken the hitchhiking effect (see below).

Time of hitchhiking:
Previous studies have shown that in the absence of recombination hitchhiking can be detected for a brief period, ~0.5N generations, some time after a hitchhiking event (SIMONSEN et al. 1995 ; FU 1997 ). We have shown that in the presence of recombination, hitchhiking can be detected immediately after the fixation of an advantageous mutation. Regardless of whether recombination is present or absent, there is only a transient period after hitchhiking when the deviation from the neutral frequency spectrum is detectable. As a population recovers from hitchhiking, D can be expected to be more powerful than H. This is because D but not H is sensitive to the influx of new mutations during the recovery period. As the excess of high-frequency variants is less likely to be influenced by other forces such as population history or background selection, monitoring this portion of the frequency spectrum represents a trade-off between power and distinguishing hitchhiking from several different hypotheses.

The significant excess of high-frequency variants in regions of low recombination as reported in Table 2 therefore suggests hitchhiking may be quite frequent in these regions, an inference supported by uniformly low levels of variation in regions of low recombination. As a result, recovery to equilibrium levels of variation may never occur and the effects of multiple or overlapping hitchhiking events must be considered. Models of strong recurrent hitchhiking events (BRAVERMAN et al. 1995 ) and numerous sites concurrently under weak selection (SANTIAGO and CABALLERO 1998 ) predict an excess of rare variants and negative D values. Although recurrent hitchhiking is not modeled in this study, a single strong hitchhiking event will push all variation to either very low or high frequencies, regardless of the frequency spectrum before selection. However, it is possible that a second round of hitchhiking may result in high-frequency variants being the only detectable polymorphism, as is observed at achaete. If the hitchhiking effect in the first round is strong and only low-frequency variants remain, during the second round of hitchhiking, low-frequency variants will either increase to high frequencies or decrease to nondetectable frequencies.

Population subdivision:
While an excess of low-frequency variants can be explained by positive selection, background selection, a change in population size, and population subdivision (CHARLESWORTH et al. 1993 ), an excess of high-frequency variants can only be explained by positive selection and specific demographic models. For example, if there are many fixed differences between populations and if only a few migrants are found in a sample, most variants would be at either low or high frequency. A more extreme scenario is that of introgression of alleles from one species into another, in which case all substitutions between the two species would appear to be either low- or high-frequency variants. Neither scenario is applicable to the data used in this article since ancient population subdivision is not supported in either D. melanogaster or D. ananassae (STEPHAN 1989 ; STEPHAN and LANGLEY 1989 ; MORIYAMA and POWELL 1996 ). However, population subdivision does affect the variance of all tests that assume an equilibrium nondifferentiated population (but see CHARLESWORTH 1998 ; STEPHAN et al. 1998 ).

The pooling of subpopulation data makes the H test conservative and may obscure patterns of hitchhiking in a single subpopulation. For instance, at su(f) there is one low (2/47) and one high (43/47) frequency-derived variant found in a North American sample. Neither of the two variants are found in a sample of 49 African flies (BEGUN and AQUADRO 1995 ). Thus, the pooled samples (2/96, 43/96) have no derived variants at a frequency of >50%. Population subdivision may also have obscured a clear hitchhiking pattern at vermilion. There is a nonsignificant excess of high-frequency variants when four subpopulations are pooled. Using a method independent of the frequency spectrum, which compares levels of variation within and between populations to divergence, hitchhiking was inferred to have homogenized two of the four populations, an observation incompatible with background selection (STEPHAN et al. 1998 ).

	CONCLUSIONS

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

An excess of very-high-frequency variants is a distinct consequence of hitchhiking. Background selection cannot easily account for this pattern, which is observed in some regions of low recombination, as revealed by the H statistic. Although it has been previously suggested that the lack of significant D values in these regions is inconsistent with a hitchhiking explanation (BRAVERMAN et al. 1995 ; CHARLESWORTH et al. 1995 ), cases of hitchhiking that yield significant H but not D and vice versa are not uncommon. The absence of significant D values can also be explained by the unexplored range of parameter space of the recurrent hitchhiking model, population subdivision, or by alternative models of selection, such as frequency-dependent selection (GILLESPIE 1991 ; BRAVERMAN et al. 1995 ). Further improvement of our ability to detect positive selection and distinguish its effects from other evolutionary processes will be made by incorporating population structure (CHARLESWORTH 1998 ; SLATKIN and WIEHE 1998 ; STEPHAN et al. 1998 ) and linkage disequilibrium (THOMSON 1977 ; GROTE et al. 1998 ) into models of selection.

	ACKNOWLEDGMENTS

We thank B. Charlesworth, R. Hudson, M. Jensen, C. Langley, T. Nagylaki, J. Spofford, C. Ting, J. Wall, and two reviewers for comments and suggestions. We thank S. Sun for help in deriving the probability of a back mutation. This work was supported by National Institutes of Health and National Science Foundation grants to C.-I Wu and a Genetics Training Grant to J. C. Fay.

Manuscript received August 7, 1999; Accepted for publication March 20, 2000.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS THEORY APPLICATION TO DATA DISCUSSION CONCLUSIONS LITERATURE CITED

AGUADÉ, M., 1998  Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex. Genetics 150:1079-1089[Abstract/Free Full Text].

AGUADÉ, M., N. MIYASHITA, and C. H. LANGLEY, 1989  Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster.. Genetics 122:607-615[Abstract/Free Full Text].

AGUADÉ, M., N. MIYASHITA, and C. H. LANGLEY, 1992  Polymorphism and divergence in the Mst26a male accessory gland gene region. Genetics 132:755-770[Abstract].

ASHBURNER, M., 1989 Drosophila: A Laboratory Handbook. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.

BARTON, N. H., 1998  The effect of hitch-hiking on neutral genealogies. Genet. Res. 72:123-133.

BEGUN, D. J. and C. F. AQUADRO, 1991  Molecular population genetics of the distal portion of the X chromosome in Drosophila: evidence for genetic hitchhiking of the yellow-achaete region. Genetics 129:1147-1158[Abstract].

BEGUN, D. J. and C. F. AQUADRO, 1992  Levels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogaster.. Nature 356:519-520[Medline].

BEGUN, D. J. and C. F. AQUADRO, 1994  Evolutionary inferences from DNA variation at the 6-phosphogluconate dehydrogenase locus in natural populations of Drosophila: selection and geographic differentiation. Genetics 136:155-171[Abstract].

BEGUN, D. J. and C. F. AQUADRO, 1995  Evolution at the tip and base of the X chromosome in an African population of Drosophila melanogaster.. Mol. Biol. Evol. 12:382-390[Abstract].

BRAVERMAN, J. M., R. R. HUDSON, N. L. KAPLAN, C. H. LANGLEY, and W. STEPHAN, 1995  The hitchhiking effect on the site frequency spectrum of DNA polymorphisms. Genetics 140:783-796[Abstract].

BROMAN, K. W., J. C. MURRAY, V. C. SHEFFIELD, R. L. WHITE, and J. L. WEBER, 1998  Comprehensive human genetic maps: individual and sex-specific variation in recombination. Am. J. Hum. Genet. 63:861-869[Medline].

CHARLESWORTH, B., 1998  Measures of divergence between populations and the effect of forces that reduce variability. Mol. Biol. Evol. 15:538-543[Abstract].

CHARLESWORTH, B., M. T. MORGAN, and D. CHARLESWORTH, 1993  The effect of deleterious mutations on neutral molecular variation. Genetics 134:1289-1303[Abstract].

CHARLESWORTH, D., B. CHARLESWORTH, and M. T. MORGAN, 1995  The pattern of neutral molecular variation under the background selection model. Genetics 141:1619-1632[Abstract].

COMERÓN, J. M., M. KREITMAN, and M. AGUADÉ, 1999  Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.. Genetics 151:239-249[Abstract/Free Full Text].

FAY, J. C. and C.-I WU, 1999  A human population bottleneck is compatible with the discordance between patterns of mitochondrial vs. nuclear DNA variation. Mol. Biol. Evol. 16:1003-1005[Medline].

FLYBASE,, 1999  The FlyBase Database of the Drosophila Genome Projects and community literature. Nucleic Acids Res. 27:85-88[Abstract/Free Full Text].

FU, Y. X., 1995  Statistical properties of segregating sites. Theor. Popul. Biol. 48:172-197[Medline].

FU, Y. X., 1997  Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147:915-925[Abstract].

GILLESPIE, J. H., 1991 The Causes of Molecular Evolution. Oxford University Press, Oxford.

GROTE, M. N., W. KLITZ, and G. THOMSON, 1998  Constrained disequilibrium values and hitchhiking in a three-locus system. Genetics 150:1295-1307[Abstract/Free Full Text].

HAMBLIN, M. T. and C. F. AQUADRO, 1997  Contrasting patterns of nucleotide sequence variation at the Glucose Dehydrogenase (Gld) locus in different populations of Drosophila melanogaster. Genetics 145:1053-1062[Abstract].

HILTON, H., R. M. KLIMAN, and J. HEY, 1994  Using hitchhiking genes to study adaptation and divergence during speciation within the Drosophila melanogaster species complex. Evolution 48:1900-1913.

HUDSON, R. R., 1983  Properties of a neutral allele model with intragenic recombination. Theor. Popul. Biol. 23:183-201[Medline].

HUDSON, R. R., 1990 Gene genealogies and the coalescent process, pp. 1–44 in Oxford Surveys in Evolutionary Biology, edited by D. FUTUYAMA and J. ANTONOVICS. Oxford University Press, New York.

HUDSON, R. R., M. KREITMAN, and M. AGUADÉ, 1987  A test of neutral molecular evolution based on nucleotide data. Genetics 116:153-159[Abstract/Free Full Text].

JENKINS, D. L., C. A. ORTORI, and J. F. Y. BROOKFIELD, 1995  A test for adaptive change in DNA sequences controlling transcription. Proc. R. Soc. Lond. Ser. B 261:203-207[Medline].

KAPLAN, N. L., R. R. HUDSON, and C. H. LANGLEY, 1989  The "hitchhiking effect" revisited. Genetics 123:887-899[Abstract/Free Full Text].

KIMURA, M., 1983 The Neutral Theory of Molecular Evolution. Cambridge University Press, Cambridge, UK.

LANGLEY, C. H., J. MCDONALD, N. MIYASHITA, and M. AGUADÉ, 1993  Lack of correlation between interspecific divergence and intraspecific polymorphism at the suppressor of forked region in Drosophila melanogaster and D. simulans.. Proc. Natl. Acad. Sci. USA 90:1800-1803[Abstract/Free Full Text].

MARTÍN-CAMPOS, J., J. M. COMERÓN, N. MIYASHITA, and M. AGUADÉ, 1992  Intraspecific and interspecific variation at the y-ac-sc region of Drosophila simulans and Drosophila melanogaster.. Genetics 130:805-816[Abstract].

MAYNARD SMITH, J. and J. HAIGH, 1974  The hitch-hiking effect of a favorable gene. Genet. Res. 23:23-35[Medline].

MCDONALD, J. H., 1998  Improved tests for heterogeneity across a region of DNA sequence in the ratio of polymorphism to divergence. Mol. Biol. Evol. 15:377-384[Abstract].

MCDONALD, J. H. and M. KREITMAN, 1991  Adaptive protein evolution at the Adh locus in Drosophila.. Nature 351:652-654[Medline].

MORIYAMA, E. N. and J. R. POWELL, 1996  Intraspecific nuclear DNA variation in Drosophila. Mol. Biol. Evol. 13:261-277[Abstract].

NACHMAN, M. W., 1997  Patterns of DNA variability at X-linked loci in Mus domesticus.. Genetics 147:1303-1316[Abstract].

NACHMAN, M. W., V. L. BAUER, S. L. CROWELL, and C. F. AQUADRO, 1998  DNA variability and recombination rates at X-linked loci in humans. Genetics 150:1133-1141[Abstract/Free Full Text].

NEI, M., 1987 Molecular Evolutionary Genetics. Columbia University Press, New York.

ROZAS, J. and R. ROZAS, 1999  DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 15:174-175[Abstract/Free Full Text].

SANTIAGO, E. and A. CABALLERO, 1998  Effective size and polymorphism of linked neutral loci in populations under directional selection. Genetics 149:2105-2117[Abstract/Free Full Text].

SIMONSEN, K. L., G. A. CHURCHILL, and C. F. AQUADRO, 1995  Properties of statistical tests of neutrality for DNA polymorphism data. Genetics 141:413-429[Abstract].

SLATKIN, M. and T. WIEHE, 1998  Genetic hitch-hiking in a subdivided population. Genet. Res. 71:155-160[Medline].

SOKAL, R. R., and F. J. ROHLF, 1981 Biometry. W. H. Freeman and Co., New York.

STEPHAN, W., 1989  Molecular genetic variation in the centromeric region of the X chromosome in three Drosophila ananassae populations. II. The Om(1D) locus. Mol. Biol. Evol. 6:624-635[Abstract].

STEPHAN, W. and C. H. LANGLEY, 1989  Molecular genetic variation in the centromeric region of the x chromosome in three Drosophila ananassae populations. I. Contrasts between vermilion and forked loci. Genetics 121:89-99[Abstract/Free Full Text].

STEPHAN, W. and C. H. LANGLEY, 1998  DNA polymorphism in Lycopersicon and crossing-over per physical length. Genetics 150:1585-1593[Abstract/Free Full Text].

STEPHAN, W. and S. J. MITCHELL, 1992  Reduced levels of DNA polymorphism and fixed between-population differences in the centromeric region of Drosophila ananassae.. Genetics 132:1039-1045[Abstract].

STEPHAN, W., T. H. E. WIEHE, and M. W. LENZ, 1992  The effect of strongly selected substitutions on neutral polymorphism-analytical results based on diffusion theory. Theor. Popul. Biol. 41:237-254.

STEPHAN, W., L. XING, D. A. KIRBY, and J. M. BRAVERMAN, 1998  A test of background selection hypothesis based on nucleotide data from Drosophila ananassae.. Proc. Natl. Acad. Sci. USA 95:5649-5654[Abstract/Free Full Text].

TAJIMA, F., 1983  Evolutionary relationship of DNA sequences in finite populations. Genetics 123:437-460.

TAJIMA, F., 1989a  The effect of change in population size on DNA polymorphism. Genetics 123:597-601[Abstract/Free Full Text].

TAJIMA, F., 1989b  Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585-595[Abstract/Free Full Text].

THOMSON, G., 1977  The effect of a selected locus on linked neutral loci. Genetics 85:753-788[Abstract/Free Full Text].

TSAUR, S. C. and C.-I WU, 1997  Positive selection and the molecular evolution of a gene of male reproduction, Acp26Aa of Drosophila.. Mol. Biol. Evol. 14:544-549[Abstract].

TSAUR, S. C., C. T. TING, and C.-I WU, 1998  Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila. II. Divergence versus polymorphism. Mol. Biol. Evol. 15:1040-1046[Abstract].

WALL, J. D., 1999  Recombination and the power of statistical tests of neutrality. Genet. Res. 74:65-79.

WATTERSON, G. A., 1975  On the number of segregating sites. Theor. Popul. Biol. 7:256-276[Medline].

This article has been cited by other articles:

				D. Prada Molecular population genetics and agronomic alleles in seed banks: searching for a needle in a haystack? J. Exp. Bot., July 1, 2009; 60(9): 2541 - 2552. [Abstract] [Full Text] [PDF]

				N. Svetec, P. Pavlidis, and W. Stephan Recent Strong Positive Selection on Drosophila melanogaster HDAC6, a Gene Encoding a Stress Surveillance Factor, as Revealed by Population Genomic Analysis Mol. Biol. Evol., July 1, 2009; 26(7): 1549 - 1556. [Abstract] [Full Text] [PDF]

				L. B. Scheinfeldt, S. Biswas, J. Madeoy, C. F. Connelly, E. E. Schadt, and J. M. Akey Population Genomic Analysis of ALMS1 in Humans Reveals a Surprisingly Complex Evolutionary History Mol. Biol. Evol., June 1, 2009; 26(6): 1357 - 1367. [Abstract] [Full Text] [PDF]

				R. Nielsen, M. J. Hubisz, I. Hellmann, D. Torgerson, A. M. Andres, A. Albrechtsen, R. Gutenkunst, M. D. Adams, M. Cargill, A. Boyko, et al. Darwinian and demographic forces affecting human protein coding genes Genome Res., May 1, 2009; 19(5): 838 - 849. [Abstract] [Full Text] [PDF]

				T. Stadler, B. Haubold, C. Merino, W. Stephan, and P. Pfaffelhuber The Impact of Sampling Schemes on the Site Frequency Spectrum in Nonequilibrium Subdivided Populations Genetics, May 1, 2009; 182(1): 205 - 216. [Abstract] [Full Text] [PDF]

				S. Boitard, C. Schlotterer, and A. Futschik Detecting Selective Sweeps: A New Approach Based on Hidden Markov Models Genetics, April 1, 2009; 181(4): 1567 - 1578. [Abstract] [Full Text] [PDF]

				T. Miyake, N. Takebayashi, and D. E. Wolf Possible Diversifying Selection in the Imprinted Gene, MEDEA, in Arabidopsis Mol. Biol. Evol., April 1, 2009; 26(4): 843 - 857. [Abstract] [Full Text] [PDF]

				J. L. Kelley, K. Turkheimer, M. Haney, and W. J. Swanson Targeted resequencing of two genes, RAGE and POLL, confirms findings from a genome-wide scan for adaptive evolution and provides evidence for positive selection in additional populations Hum. Mol. Genet., February 15, 2009; 18(4): 779 - 784. [Abstract] [Full Text] [PDF]

				S. Tang and D. C. Presgraves Evolution of the Drosophila Nuclear Pore Complex Results in Multiple Hybrid Incompatibilities Science, February 6, 2009; 323(5915): 779 - 782. [Abstract] [Full Text] [PDF]

				A. C. Carnaval, M. J. Hickerson, C. F. B. Haddad, M. T. Rodrigues, and C. Moritz Stability Predicts Genetic Diversity in the Brazilian Atlantic Forest Hotspot Science, February 6, 2009; 323(5915): 785 - 789. [Abstract] [Full Text] [PDF]

				M. Fumagalli, R. Cagliani, U. Pozzoli, S. Riva, G. P. Comi, G. Menozzi, N. Bresolin, and M. Sironi Widespread balancing selection and pathogen-driven selection at blood group antigen genes Genome Res., February 1, 2009; 19(2): 199 - 212. [Abstract] [Full Text] [PDF]

				A. Ramirez-Soriano and R. Nielsen Correcting Estimators of {theta} and Tajima's D for Ascertainment Biases Caused by the Single-Nucleotide Polymorphism Discovery Process Genetics, February 1, 2009; 181(2): 701 - 710. [Abstract] [Full Text] [PDF]

				W. Zhai, R. Nielsen, and M. Slatkin An Investigation of the Statistical Power of Neutrality Tests Based on Comparative and Population Genetic Data Mol. Biol. Evol., February 1, 2009; 26(2): 273 - 283. [Abstract] [Full Text] [PDF]

				Y. Kim and T. Wiehe Simulation of DNA sequence evolution under models of recent directional selection Brief Bioinform, January 1, 2009; 10(1): 84 - 96. [Abstract] [Full Text] [PDF]

				K. H. Takahashi, K. Tanaka, M. Itoh, and T. Takano-Shimizu Reduced X-Linked Rare Polymorphism in Males in Comparison to Females of Drosophila melanogaster J. Hered., January 1, 2009; 100(1): 97 - 105. [Abstract] [Full Text] [PDF]

				A. E. Palme, M. Wright, and O. Savolainen Patterns of Divergence among Conifer ESTs and Polymorphism in Pinus sylvestris Identify Putative Selective Sweeps Mol. Biol. Evol., December 1, 2008; 25(12): 2567 - 2577. [Abstract] [Full Text] [PDF]

				W. M. Muir, G. K.-S. Wong, Y. Zhang, J. Wang, M. A. M. Groenen, R. P. M. A. Crooijmans, H.-J. Megens, H. Zhang, R. Okimoto, A. Vereijken, et al. Genome-wide assessment of worldwide chicken SNP genetic diversity indicates significant absence of rare alleles in commercial breeds PNAS, November 11, 2008; 105(45): 17312 - 17317. [Abstract] [Full Text] [PDF]

				Y. Yu, T. Tang, Q. Qian, Y. Wang, M. Yan, D. Zeng, B. Han, C.-I Wu, S. Shi, and J. Li Independent Losses of Function in a Polyphenol Oxidase in Rice: Differentiation in Grain Discoloration between Subspecies and the Role of Positive Selection under Domestication PLANT CELL, November 1, 2008; 20(11): 2946 - 2959. [Abstract] [Full Text] [PDF]

				A. R. Omilian, D. G. Scofield, and M. Lynch Intron Presence-Absence Polymorphisms in Daphnia Mol. Biol. Evol., October 1, 2008; 25(10): 2129 - 2139. [Abstract] [Full Text] [PDF]

				C. Fan, Y. Zhang, Y. Yu, S. Rounsley, M. Long, and R. A. Wing The Subtelomere of Oryza sativa Chromosome 3 Short Arm as a Hot Bed of New Gene Origination in Rice Mol Plant, September 1, 2008; 1(5): 839 - 850. [Abstract] [Full Text] [PDF]

				L.-M. Chevin, S. Billiard, and F. Hospital Hitchhiking Both Ways: Effect of Two Interfering Selective Sweeps on Linked Neutral Variation Genetics, September 1, 2008; 180(1): 301 - 316. [Abstract] [Full Text] [PDF]

				P. K. Ingvarsson Multilocus Patterns of Nucleotide Polymorphism and the Demographic History of Populus tremula Genetics, September 1, 2008; 180(1): 329 - 340. [Abstract] [Full Text] [PDF]

				D. Zivkovic and T. Wiehe Second-Order Moments of Segregating Sites Under Variable Population Size Genetics, September 1, 2008; 180(1): 341 - 357. [Abstract] [Full Text] [PDF]

				G. Achaz Testing for Neutrality in Samples With Sequencing Errors Genetics, July 1, 2008; 179(3): 1409 - 1424. [Abstract] [Full Text] [PDF]

				C.-L. Huang, S.-Y. Hwang, Y.-C. Chiang, and T.-P. Lin Molecular Evolution of the Pi-ta Gene Resistant to Rice Blast in Wild Rice (Oryza rufipogon) Genetics, July 1, 2008; 179(3): 1527 - 1538. [Abstract] [Full Text] [PDF]

				C. Fan, Y. Chen, and M. Long Recurrent Tandem Gene Duplication Gave Rise to Functionally Divergent Genes in Drosophila Mol. Biol. Evol., July 1, 2008; 25(7): 1451 - 1458. [Abstract] [Full Text] [PDF]

				G. Coop, K. Bullaughey, F. Luca, and M. Przeworski The Timing of Selection at the Human FOXP2 Gene Mol. Biol. Evol., July 1, 2008; 25(7): 1257 - 1259. [Abstract] [Full Text] [PDF]

				A. Llopart and J. M. Comeron Recurrent Events of Positive Selection in Independent Drosophila Lineages at the Spermatogenesis Gene roughex Genetics, June 1, 2008; 179(2): 1009 - 1020. [Abstract] [Full Text] [PDF]

				J. M. Macpherson, J. Gonzalez, D. M. Witten, J. C. Davis, N. A. Rosenberg, A. E. Hirsh, and D. A. Petrov Nonadaptive Explanations for Signatures of Partial Selective Sweeps in Drosophila Mol. Biol. Evol., June 1, 2008; 25(6): 1025 - 1042. [Abstract] [Full Text] [PDF]

				P. Pfaffelhuber, A. Lehnert, and W. Stephan Linkage Disequilibrium Under Genetic Hitchhiking in Finite Populations Genetics, May 1, 2008; 179(1): 527 - 537. [Abstract] [Full Text] [PDF]

				A. Ramirez-Soriano, S. E. Ramos-Onsins, J. Rozas, F. Calafell, and A. Navarro Statistical Power Analysis of Neutrality Tests Under Demographic Expansions, Contractions and Bottlenecks With Recombination Genetics, May 1, 2008; 179(1): 555 - 567. [Abstract] [Full Text] [PDF]

				K. M. Olsen, S.-C. Hsu, and L. L. Small Evidence on the Molecular Basis of the Ac/ac Adaptive Cyanogenesis Polymorphism in White Clover (Trifolium repens L.) Genetics, May 1, 2008; 179(1): 517 - 526. [Abstract] [Full Text] [PDF]

				N. Warthmann, S. Das, C. Lanz, and D. Weigel Comparative Analysis of the MIR319a MicroRNA Locus in Arabidopsis and Related Brassicaceae Mol. Biol. Evol., May 1, 2008; 25(5): 892 - 902. [Abstract] [Full Text] [PDF]

				J. Lu, Y. Fu, S. Kumar, Y. Shen, K. Zeng, A. Xu, R. Carthew, and C.-I Wu Adaptive Evolution of Newly Emerged Micro-RNA Genes in Drosophila Mol. Biol. Evol., May 1, 2008; 25(5): 929 - 938. [Abstract] [Full Text] [PDF]

				I. M. Ehrenreich and M. D. Purugganan Sequence Variation of MicroRNAs and Their Binding Sites in Arabidopsis Plant Physiology, April 1, 2008; 146(4): 1974 - 1982. [Abstract] [Full Text] [PDF]

				J. L. Kelley and W. J. Swanson Dietary Change and Adaptive Evolution of enamelin in Humans and Among Primates Genetics, March 1, 2008; 178(3): 1595 - 1603. [Abstract] [Full Text] [PDF]

				A. Wong, M. C. Turchin, M. F. Wolfner, and C. F. Aquadro Evidence for Positive Selection on Drosophila melanogaster Seminal Fluid Protease Homologs Mol. Biol. Evol., March 1, 2008; 25(3): 497 - 506. [Abstract] [Full Text] [PDF]

				J. J. Cai PGEToolbox: A Matlab Toolbox for Population Genetics and Evolution J. Hered., February 29, 2008; (2008) esm127v1. [Abstract] [Full Text] [PDF]

				L. M. Matzkin The Molecular Basis of Host Adaptation in Cactophilic Drosophila: Molecular Evolution of a Glutathione S-Transferase Gene (GstD1) in Drosophila mojavensis Genetics, February 1, 2008; 178(2): 1073 - 1083. [Abstract] [Full Text] [PDF]

				H. Kuittinen, A. Niittyvuopio, P. Rinne, and O. Savolainen Natural Variation in Arabidopsis lyrata Vernalization Requirement Conferred by a FRIGIDA Indel Polymorphism Mol. Biol. Evol., February 1, 2008; 25(2): 319 - 329. [Abstract] [Full Text] [PDF]

				E. Eveno, C. Collada, M. A. Guevara, V. Leger, A. Soto, L. Diaz, P. Leger, S. C. Gonzalez-Martinez, M. T. Cervera, C. Plomion, et al. "Contrasting Patterns of Selection at Pinus pinaster Ait. Drought Stress Candidate Genes as Revealed by Genetic Differentiation Analyses" Mol. Biol. Evol., February 1, 2008; 25(2): 417 - 437. [Abstract] [Full Text] [PDF]

				N. Derome, E. Baudry, D. Ogereau, M. Veuille, and C. Montchamp-Moreau Selective Sweeps in a 2-Locus Model for Sex-Ratio Meiotic Drive in Drosophila simulans Mol. Biol. Evol., February 1, 2008; 25(2): 409 - 416. [Abstract] [Full Text] [PDF]

				J. D. Jensen, K. R. Thornton, and C. F. Aquadro Inferring Selection in Partially Sequenced Regions Mol. Biol. Evol., February 1, 2008; 25(2): 438 - 446. [Abstract] [Full Text] [PDF]

				V. Nolte and C. Schlotterer African Drosophila melanogaster and D. simulans Populations Have Similar Levels of Sequence Variability, Suggesting Comparable Effective Population Sizes Genetics, January 1, 2008; 178(1): 405 - 412. [Abstract] [Full Text] [PDF]

				S. De Mita, J. Ronfort, H. I. McKhann, C. Poncet, R. El Malki, and T. Bataillon Investigation of the Demographic and Selective Forces Shaping the Nucleotide Diversity of Genes Involved in Nod Factor Signaling in Medicago truncatula Genetics, December 1, 2007; 177(4): 2123 - 2133. [Abstract] [Full Text] [PDF]

				D. Garrigan, S. B. Kingan, M. M. Pilkington, J. A. Wilder, M. P. Cox, H. Soodyall, B. Strassmann, G. Destro-Bisol, P. de Knijff, A. Novelletto, et al. Inferring Human Population Sizes, Divergence Times and Rates of Gene Flow From Mitochondrial, X and Y Chromosome Resequencing Data Genetics, December 1, 2007; 177(4): 2195 - 2207. [Abstract] [Full Text] [PDF]

				T. Salcedo, A. Geraldes, and M. W. Nachman Nucleotide Variation in Wild and Inbred Mice Genetics, December 1, 2007; 177(4): 2277 - 2291. [Abstract] [Full Text] [PDF]

				T. Pyhajarvi, M. R. Garcia-Gil, T. Knurr, M. Mikkonen, W. Wachowiak, and O. Savolainen Demographic History Has Influenced Nucleotide Diversity in European Pinus sylvestris Populations Genetics, November 1, 2007; 177(3): 1713 - 1724. [Abstract] [Full Text] [PDF]

				W. Y. Low, H. L. Ng, C. J. Morton, M. W. Parker, P. Batterham, and C. Robin Molecular Evolution of Glutathione S-Transferases in the Genus Drosophila Genetics, November 1, 2007; 177(3): 1363 - 1375. [Abstract] [Full Text] [PDF]

				A. L. Evans, P. A. Mena, and B. F. McAllister Positive Selection Near an Inversion Breakpoint on the Neo-X Chromosome of Drosophila americana Genetics, November 1, 2007; 177(3): 1303 - 1319. [Abstract] [Full Text] [PDF]

				P. Wang, R. F. Lyman, S. A. Shabalina, T. F. C. Mackay, and R. R. H. Anholt Association of Polymorphisms in Odorant-Binding Protein Genes With Variation in Olfactory Response to Benzaldehyde in Drosophila Genetics, November 1, 2007; 177(3): 1655 - 1665. [Abstract] [Full Text] [PDF]

				K. R. Thornton The Neutral Coalescent Process for Recent Gene Duplications and Copy-Number Variants Genetics, October 1, 2007; 177(2): 987 - 1000. [Abstract] [Full Text] [PDF]

				J. D. Jensen, V. L. Bauer DuMont, A. B. Ashmore, A. Gutierrez, and C. F. Aquadro Patterns of Sequence Variability and Divergence at the diminutive Gene Region of Drosophila melanogaster: Complex Patterns Suggest an Ancestral Selective Sweep Genetics, October 1, 2007; 177(2): 1071 - 1085. [Abstract] [Full Text] [PDF]

				U. Arunyawat, W. Stephan, and T. Stadler Using Multilocus Sequence Data to Assess Population Structure, Natural Selection, and Linkage Disequilibrium in Wild Tomatoes Mol. Biol. Evol., October 1, 2007; 24(10): 2310 - 2322. [Abstract] [Full Text] [PDF]

				X. Wang and J. Zhang Rapid evolution of primate ESX1, an X-linked placenta- and testis-expressed homeobox gene Hum. Mol. Genet., September 1, 2007; 16(17): 2053 - 2060. [Abstract] [Full Text] [PDF]

				J. D. Jensen, K. R. Thornton, C. D. Bustamante, and C. F. Aquadro On the Utility of Linkage Disequilibrium as a Statistic for Identifying Targets of Positive Selection in Nonequilibrium Populations Genetics, August 1, 2007; 176(4): 2371 - 2379. [Abstract] [Full Text] [PDF]

				R. D. Hernandez, S. H. Williamson, and C. D. Bustamante Context Dependence, Ancestral Misidentification, and Spurious Signatures of Natural Selection Mol. Biol. Evol., August 1, 2007; 24(8): 1792 - 1800. [Abstract] [Full Text] [PDF]

				K. Zeng, S. Shi, and C.-I Wu Compound Tests for the Detection of Hitchhiking Under Positive Selection Mol. Biol. Evol., August 1, 2007; 24(8): 1898 - 1908. [Abstract] [Full Text] [PDF]

				T. Nora, C. Charpentier, O. Tenaillon, C. Hoede, F. Clavel, and A. J. Hance Contribution of Recombination to the Evolution of Human Immunodeficiency Viruses Expressing Resistance to Antiretroviral Treatment J. Virol., July 15, 2007; 81(14): 7620 - 7628. [Abstract] [Full Text] [PDF]

				G. Wen, J. Wessel, W. Zhou, G. B. Ehret, F. Rao, M. Stridsberg, S. K. Mahata, P. M. Gent, M. Das, R. S. Cooper, et al. An ancestral variant of Secretogranin II confers regulation by PHOX2 transcription factors and association with hypertension Hum. Mol. Genet., July 15, 2007; 16(14): 1752 - 1764. [Abstract] [Full Text] [PDF]

				K. Zeng, S. Mano, S. Shi, and C.-I Wu Comparisons of Site- and Haplotype-Frequency Methods for Detecting Positive Selection Mol. Biol. Evol., July 1, 2007; 24(7): 1562 - 1574. [Abstract] [Full Text] [PDF]

				A. De and R. Durrett Stepping-Stone Spatial Structure Causes Slow Decay of Linkage Disequilibrium and Shifts the Site Frequency Spectrum Genetics, June 1, 2007; 176(2): 969 - 981. [Abstract] [Full Text] [PDF]

				D. J. Orengo and M. Aguade Genome Scans of Variation and Adaptive Change: Extended Analysis of a Candidate Locus Close to the phantom Gene Region in Drosophila melanogaster Mol. Biol. Evol., May 1, 2007; 24(5): 1122 - 1129. [Abstract] [Full Text] [PDF]

				F. Yu, R. S. Hill, S. F. Schaffner, P. C. Sabeti, E. T. Wang, A. A. Mignault, R. J. Ferland, R. K. Moyzis, C. A. Walsh, and D. Reich Comment on "Ongoing Adaptive Evolution of ASPM, a Brain Size Determinant in Homo sapiens" Science, April 20, 2007; 316(5823): 370b - 370b. [Abstract] [Full Text] [PDF]

				A. Sanchez-Gracia and J. Rozas Unusual Pattern of Nucleotide Sequence Variation at the OS-E and OS-F Genomic Regions of Drosophila simulans Genetics, April 1, 2007; 175(4): 1923 - 1935. [Abstract] [Full Text] [PDF]

				S. Jagadeeshan and R. S. Singh Rapid Evolution of Outer Egg Membrane Proteins in the Drosophila melanogaster Subgroup: A Case of Ecologically Driven Evolution of Female Reproductive Traits Mol. Biol. Evol., April 1, 2007; 24(4): 929 - 938. [Abstract] [Full Text] [PDF]

				J. Gojobori, H. Tang, J. M. Akey, and C.-I Wu Adaptive evolution in humans revealed by the negative correlation between the polymorphism and fixation phases of evolution PNAS, March 6, 2007; 104(10): 3907 - 3912. [Abstract] [Full Text] [PDF]

				C. A. Machado, T. S. Haselkorn, and M. A. F. Noor Evaluation of the Genomic Extent of Effects of Fixed Inversion Differences on Intraspecific Variation and Interspecific Gene Flow in Drosophila pseudoobscura and D. persimilis Genetics, March 1, 2007; 175(3): 1289 - 1306. [Abstract] [Full Text] [PDF]

				G. McVean The Structure of Linkage Disequilibrium Around a Selective Sweep Genetics, March 1, 2007; 175(3): 1395 - 1406. [Abstract] [Full Text] [PDF]

				J. A. Shapiro, W. Huang, C. Zhang, M. J. Hubisz, J. Lu, D. A. Turissini, S. Fang, H.-Y. Wang, R. R. Hudson, R. Nielsen, et al. Adaptive genic evolution in the Drosophila genomes PNAS, February 13, 2007; 104(7): 2271 - 2276. [Abstract] [Full Text] [PDF]

				V. Mustonen and M. Lassig Adaptations to fluctuating selection in Drosophila PNAS, February 13, 2007; 104(7): 2277 - 2282. [Abstract] [Full Text] [PDF]

				D. C. Presgraves and W. Stephan Pervasive Adaptive Evolution among Interactors of the Drosophila Hybrid Inviability Gene, Nup96 Mol. Biol. Evol., January 1, 2007; 24(1): 306 - 314. [Abstract] [Full Text] [PDF]

				V. L. Bauer DuMont, H. A. Flores, M. H. Wright, and C. F. Aquadro Recurrent Positive Selection at Bgcn, a Key Determinant of Germ Line Differentiation, Does Not Appear to be Driven by Simple Coevolution with Its Partner Protein Bam Mol. Biol. Evol., January 1, 2007; 24(1): 182 - 191. [Abstract] [Full Text] [PDF]

				P. Pfaffelhuber, B. Haubold, and A. Wakolbinger Approximate Genealogies Under Genetic Hitchhiking Genetics, December 1, 2006; 174(4): 1995 - 2008. [Abstract] [Full Text] [PDF]

				D. Bachtrog and P. Andolfatto Selection, Recombination and Demographic History in Drosophila miranda Genetics, December 1, 2006; 174(4): 2045 - 2059. [Abstract] [Full Text] [PDF]

				A. Kawabe, S. Nasuda, and D. Charlesworth Duplication of Centromeric Histone H3 (HTR12) Gene in Arabidopsis halleri and A. lyrata, Plant Species With Multiple Centromeric Satellite Sequences Genetics, December 1, 2006; 174(4): 2021 - 2032. [Abstract] [Full Text] [PDF]

				J. de Meaux, A. Pop, and T. Mitchell-Olds Cis-regulatory Evolution of Chalcone-Synthase Expression in the Genus Arabidopsis Genetics, December 1, 2006; 174(4): 2181 - 2202. [Abstract] [Full Text] [PDF]

				M. Heuertz, E. De Paoli, T. Kallman, H. Larsson, I. Jurman, M. Morgante, M. Lascoux, and N. Gyllenstrand Multilocus Patterns of Nucleotide Diversity, Linkage Disequilibrium and Demographic History of Norway Spruce [Picea abies (L.) Karst] Genetics, December 1, 2006; 174(4): 2095 - 2105. [Abstract] [Full Text] [PDF]

				K. Zeng, Y.-X. Fu, S. Shi, and C.-I Wu Statistical Tests for Detecting Positive Selection by Utilizing High-Frequency Variants Genetics, November 1, 2006; 174(3): 1431 - 1439. [Abstract] [Full Text] [PDF]

				E. Betran, Y. Bai, and M. Motiwale Fast Protein Evolution and Germ Line Expression of a Drosophila Parental Gene and Its Young Retroposed Paralog Mol. Biol. Evol., November 1, 2006; 23(11): 2191 - 2202. [Abstract] [Full Text] [PDF]

				J. E. Pool and C. F. Aquadro History and Structure of Sub-Saharan Populations of Drosophila melanogaster Genetics, October 1, 2006; 174(2): 915 - 929. [Abstract] [Full Text] [PDF]

				S. Glinka, D. De Lorenzo, and W. Stephan Evidence of Gene Conversion Associated with a Selective Sweep in Drosophila melanogaster Mol. Biol. Evol., October 1, 2006; 23(10): 1869 - 1878. [Abstract] [Full Text] [PDF]

				T. M. Panhuis and W. J. Swanson Molecular Evolution and Population Genetic Analysis of Candidate Female Reproductive Genes in Drosophila Genetics, August 1, 2006; 173(4): 2039 - 2047. [Abstract] [Full Text] [PDF]

				I. M. Ehrenreich and M. D. Purugganan The molecular genetic basis of plant adaptation Am. J. Botany, July 1, 2006; 93(7): 953 - 962. [Abstract] [Full Text] [PDF]

				J. M. Akey, W. J. Swanson, J. Madeoy, M. Eberle, and M. D. Shriver TRPV6 exhibits unusual patterns of polymorphism and divergence in worldwide populations Hum. Mol. Genet., July 1, 2006; 15(13): 2106 - 2113. [Abstract] [Full Text] [PDF]

				H. Innan Modified Hudson-Kreitman-Aguade Test and Two-Dimensional Evaluation of Neutrality Tests Genetics, July 1, 2006; 173(3): 1725 - 1733. [Abstract] [Full Text] [PDF]

				P. C. Sabeti, S. F. Schaffner, B. Fry, J. Lohmueller, P. Varilly, O. Shamovsky, A. Palma, T. S. Mikkelsen, D. Altshuler, and E. S. Lander Positive natural selection in the human lineage. Science, June 16, 2006; 312(5780): 1614 - 1620. [Abstract] [Full Text] [PDF]

				K. M. Teshima, G. Coop, and M. Przeworski How reliable are empirical genomic scans for selective sweeps? Genome Res., June 1, 2006; 16(6): 702 - 712. [Abstract] [Full Text] [PDF]

				M. T. Hamblin, A. M. Casa, H. Sun, S. C. Murray, A. H. Paterson, C. F. Aquadro, and S. Kresovich Challenges of Detecting Directional Selection After a Bottleneck: Lessons From Sorghum bicolor Genetics, June 1, 2006; 173(2): 953 - 964. [Abstract] [Full Text] [PDF]

				J. M. Comeron Weak selection and recent mutational changes influence polymorphic synonymous mutations in humans PNAS, May 2, 2006; 103(18): 6940 - 6945. [Abstract] [Full Text] [PDF]

				A. J. Greenberg and C.-I Wu Molecular Genetics of Natural Populations Mol. Biol. Evol., May 1, 2006; 23(5): 883 - 886. [Abstract] [Full Text] [PDF]

				S. Ihle, I. Ravaoarimanana, M. Thomas, and D. Tautz An Analysis of Signatures of Selective Sweeps in Natural Populations of the House Mouse Mol. Biol. Evol., April 1, 2006; 23(4): 790 - 797. [Abstract] [Full Text] [PDF]

				K. Thornton and P. Andolfatto Approximate Bayesian Inference Reveals Evidence for a Recent, Severe Bottleneck in a Netherlands Population of Drosophila melanogaster Genetics, March 1, 2006; 172(3): 1607 - 1619. [Abstract] [Full Text] [PDF]