Academic journal article Genetics

Modified Hudson-Kreitman-Aguadé Test and Two-Dimensional Evaluation of Neutrality Tests

Academic journal article Genetics

Modified Hudson-Kreitman-Aguadé Test and Two-Dimensional Evaluation of Neutrality Tests

Article excerpt


There are a number of polymorphism-based statistical tests of neutrality, but most of them focus on either the amount or the pattern of polymorphism. In this article, a new test called the two-dimensional (2D) test is developed. This test evaluates a pair of summary statistics in a two-dimentional field. One statistic should summarize the pattern of polymorphism, while the other could be a measure of the level of polymorphism. For the latter summary statistic, the polymorphism-divergence ratio is used following the idea of the Hudson-Kreitman-Aguadé (HKA) test. To incorporate the HKA test in the 2D test, a summary statistic-based version of the HKA test is developed such that the polymorphism-divergence ratio at a particular region of interest is examined if it is consistent with the average of those in other independent regions.

(ProQuest Information and Learning: ... denotes formulae omitted.)

SINCE the development of the coalescent theory (KINGMAN 1982; HUDSON 1983; TAJIMA 1983), a number of polymorphism-based statistical tests have been developed to examine a neutral null model (i.e., neutrality tests). With increasing intraspecific variation data in various species, these tests have been ubiquitous tools inmolecular population genetic analysis (KREITMAN 2000).

Neutrality tests include the following two major categories, although there are other types of tests available such as haplotype tests (HUDSON et al. 1994; FU 1996; SABETI et al. 2002) (see INNAN et al. 2005, for a recent review of haplotype tests). The first category focuses on the amount of polymorphism. Balancing selection increases the level of polymorphism because multiple alleles are likely maintained for a long time (HUDSON and KAPLAN 1988), while the level of polymorphism is reduced shortly after a fixation of adaptive mutation. This event is called a selective sweep because the fixation of a beneficial allele could sweep out the variation in the surrounding region of the selection target site by the hitchhiking effect (KAPLAN et al. 1989). The Hudson- Kreitman-Aguadé (HKA) test (HUDSON et al. 1987) focuses on this effect of selection by comparing the levels of polymorphism and divergence from an outgroup (see also WRIGHT and CHARLESWORTH 2004).

The second major category of neutrality tests examines whether the observed frequency spectrum of nucleotide polymorphism is consistent with the neutral expectation. Tajima (1989) has devised a simple method that compares ... and ..., two unbiased estimators of θ, the population mutation rate. ... is identical to the average number of pairwise nucleotide differences, which can be a direct estimator of θ (TAJIMA 1983). ... is an estimator based on the number of segregating sites, S (WATTERSON 1975). Tajima's D is defined as ..., and its expectation under the standard neutral model of a constant-size population is ~0. Balancing selection creates an excess of alleles in intermediate frequencies so that Tajima's D is likely positive, while Tajima's D tends to be negative in a region shortly after a selective sweep or under the pressure of purifying selection due to an excess of variation in low frequencies. A number of tests similar to Tajima's D have been developed (FU and LI 1993; SIMONSEN et al. 1995; FAY and WU 2000).

Thus, most neutrality tests use either the amount or the allele frequency spectrum of polymorphism. That is, those tests do not use part of the important information, which could result in a loss of power to detect selection. For example, consider a gene that experienced a recent selective sweep so that no polymorphism is observed. The HKA test could work, but the second category of tests cannot be performed when the number of segregating sites is zero. This article introduces a simple algorithm to examine both the amount and the freqeuncy spectrum of polymorphism simultaneously, which is referred to as the two-dimensional (2D) test because it examines a pair of test statistics in a two-dimensional field. …

Search by... Author
Show... All Results Primary Sources Peer-reviewed


An unknown error has occurred. Please click the button below to reload the page. If the problem persists, please try again in a little while.