Biology:Tajima's D

Short description: Population genetic test statistic

Tajima's D is a population genetic test statistic created by and named after the Japanese researcher Fumio Tajima. Tajima's D is computed as the difference between two measures of genetic diversity: the mean number of pairwise differences and the number of segregating sites, each scaled so that they are expected to be the same in a neutrally evolving population of constant size.

The purpose of Tajima's D test is to distinguish between a DNA sequence evolving randomly ("neutrally") and one evolving under a non-random process, including directional selection or balancing selection, demographic expansion or contraction, genetic hitchhiking, or introgression. A randomly evolving DNA sequence contains mutations with no effect on the fitness and survival of an organism. The randomly evolving mutations are called "neutral", while mutations under selection are "non-neutral". For example, a mutation that causes prenatal death or severe disease would be expected to be under selection. In the population as a whole, the frequency of a neutral mutation fluctuates randomly (i.e. the percentage of individuals in the population with the mutation changes from one generation to the next, and this percentage is equally likely to go up or down) through genetic drift.

The strength of genetic drift depends on population size. If a population is at a constant size with constant mutation rate, the population will reach an equilibrium of gene frequencies. This equilibrium has important properties, including the number of segregating sites [math]\displaystyle{ S }[/math], and the number of nucleotide differences between pairs sampled (these are called pairwise differences). To standardize the pairwise differences, the mean or 'average' number of pairwise differences is used. This is simply the sum of the pairwise differences divided by the number of pairs, and is often symbolized by [math]\displaystyle{ \pi }[/math].

The purpose of Tajima's test is to identify sequences which do not fit the neutral theory model at equilibrium between mutation and genetic drift. In order to perform the test on a DNA sequence or gene, you need to sequence homologous DNA for at least 3 individuals. Tajima's statistic computes a standardized measure of the total number of segregating sites (these are DNA sites that are polymorphic) in the sampled DNA and the average number of mutations between pairs in the sample. The two quantities whose values are compared are both method of moments estimates of the population genetic parameter theta, and so are expected to equal the same value. If these two numbers only differ by as much as one could reasonably expect by chance, then the null hypothesis of neutrality cannot be rejected. Otherwise, the null hypothesis of neutrality is rejected.

Scientific explanation

Under the neutral theory model, for a population at constant size at equilibrium:

[math]\displaystyle{ E[\pi]=\theta=E\left[\frac{S}{\sum_{i=1}^{n-1} \frac{1}{i}}\right]=4N\mu }[/math]

for diploid DNA, and

[math]\displaystyle{ E[\pi]=\theta=E\left[\frac{S}{\sum_{i=1}^{n-1} \frac{1}{i}}\right]=2N\mu }[/math]

for haploid.

In the above formulas, S is the number of segregating sites, n is the number of samples, N is the effective population size, [math]\displaystyle{ \mu }[/math] is the mutation rate at the examined genomic locus, and i is the index of summation. But selection, demographic fluctuations and other violations of the neutral model (including rate heterogeneity and introgression) will change the expected values of [math]\displaystyle{ S }[/math] and [math]\displaystyle{ \pi }[/math], so that they are no longer expected to be equal. The difference in the expectations for these two variables (which can be positive or negative) is the crux of Tajima's D test statistic.

[math]\displaystyle{ D\, }[/math] is calculated by taking the difference between the two estimates of the population genetics parameter [math]\displaystyle{ \theta\, }[/math]. This difference is called [math]\displaystyle{ d\, }[/math], and D is calculated by dividing [math]\displaystyle{ d\, }[/math] by the square root of its variance [math]\displaystyle{ \sqrt{\hat{V}(d)} }[/math] (its standard deviation, by definition).

[math]\displaystyle{ D=\frac {d} {\sqrt {\hat{V}(d)} } }[/math]

Fumio Tajima demonstrated by computer simulation that the [math]\displaystyle{ D\, }[/math] statistic described above could be modeled using a beta distribution. If the [math]\displaystyle{ D\, }[/math] value for a sample of sequences is outside the confidence interval then one can reject the null hypothesis of neutral mutation for the sequence in question. However, in real world uses, one must be careful as past population changes (for instance, a population bottleneck) can bias the value of the [math]\displaystyle{ D\, }[/math] statistic.

Mathematical details

[math]\displaystyle{ D=\frac {d} {\sqrt {\hat{V}(d)} } = \frac {\hat{k} - \frac{S}{a_1} } {\sqrt {[e_1S+e_2S(S-1)]} } }[/math]

where

[math]\displaystyle{ e_1 = \frac {c_1}{a_1} }[/math]	[math]\displaystyle{ e_2 = \frac{c_2}{a_1^2+a_2} }[/math]
[math]\displaystyle{ c_1 = b_1 - \frac {1}{a_1} }[/math]	[math]\displaystyle{ c_2 = b_2 - \frac{n+2}{a_{1}n} +\frac{a_2}{a_{1}^{2}} }[/math]
[math]\displaystyle{ b_1 = \frac {n+1}{3(n-1)} }[/math]	[math]\displaystyle{ b_2 = \frac{2(n^{2}+n+3)}{9n(n-1)} }[/math]
[math]\displaystyle{ a_1 = \sum_{i=1}^{n-1} \frac{1}{i} }[/math]	[math]\displaystyle{ a_2 = \sum_{i=1}^{n-1} \frac{1}{i^2} }[/math]

[math]\displaystyle{ \hat{k}\, }[/math] and [math]\displaystyle{ \frac{S}{a_1} }[/math] are two estimates of the expected number of single nucleotide polymorphisms (SNPs) between two DNA sequences under the neutral mutation model in a sample size [math]\displaystyle{ n\, }[/math] from an effective population size [math]\displaystyle{ N }[/math].

The first estimate is the average number of SNPs found in (n choose 2) pairwise comparisons of sequences [math]\displaystyle{ (i,j) }[/math] in the sample,

[math]\displaystyle{ \hat{k}= \frac { \sum\sum_{i\lt j} k_{ij} } { \binom{n}{2} }. }[/math]

The second estimate is derived from the expected value of [math]\displaystyle{ S }[/math], the total number of polymorphisms in the sample

[math]\displaystyle{ E(S)=a_1M. }[/math]

Tajima defines [math]\displaystyle{ M=4N\mu }[/math], whereas Hartl & Clark use a different symbol to define the same parameter [math]\displaystyle{ \theta=4N\mu }[/math].

Example

Suppose you are a geneticist studying an unknown gene. As part of your research you get DNA samples from four random people (plus yourself). For simplicity, you label your sequence as a string of zeroes, and for the other four people you put a zero when their DNA is the same as yours and a one when it is different. (For this example, the specific type of difference is not important.)

                    1           2
Position  12345 67890 12345 67890
Person Y  00000 00000 00000 00000
Person A  00100 00000 00100 00010
Person B  00000 00000 00100 00010
Person C  00000 01000 00000 00010
Person D  00000 01000 00100 00010

Notice the four polymorphic sites (positions where someone differs from you, at 3, 7, 13 and 19 above). Now compare each pair of sequences and get the average number of polymorphisms between two sequences. There are "five choose two" (ten) comparisons that need to be done.

Person Y is you!

You vs A: 3 polymorphisms

Person Y     00000 00000 00000 00000
Person A     00100 00000 00100 00010

You vs B: 2 polymorphisms

Person Y     00000 00000 00000 00000
Person B     00000 00000 00100 00010

You vs C: 2 polymorphisms

Person Y     00000 00000 00000 00000
Person C     00000 01000 00000 00010

You vs D: 3 polymorphisms

Person Y     00000 00000 00000 00000
Person D     00000 01000 00100 00010

A vs B: 1 polymorphism

Person A     00100 00000 00100 00010
Person B     00000 00000 00100 00010

A vs C: 3 polymorphisms

Person A     00100 00000 00100 00010
Person C     00000 01000 00000 00010

A vs D: 2 polymorphisms

Person A     00100 00000 00100 00010
Person D     00000 01000 00100 00010

B vs C: 2 polymorphisms

Person B     00000 00000 00100 00010
Person C     00000 01000 00000 00010

B vs D: 1 polymorphism

Person B     00000 00000 00100 00010
Person D     00000 01000 00100 00010

C vs D: 1 polymorphism

Person C     00000 01000 00000 00010
Person D     00000 01000 00100 00010

The average number of polymorphisms is [math]\displaystyle{ {3 + 2 + 2 + 3 + 1 + 3 + 2 + 2 + 1 + 1\over 10} = 2 }[/math].

The second estimate of the equilibrium is M=S/a1

Since there were n=5 individuals and S=4 segregating sites

a1=1/1+1/2+1/3+1/4=2.08

M=4/2.08=1.92

The lower-case d described above is the difference between these two numbers—the average number of polymorphisms found in pairwise comparison (2) and M. Thus [math]\displaystyle{ d = 2 - 1.92= .08 }[/math].

Since this is a statistical test, you need to assess the significance of this value. A discussion of how to do this is provided below.

Interpreting Tajima's D

A negative Tajima's D signifies an excess of low frequency polymorphisms relative to expectation, indicating population size expansion (e.g., after a bottleneck or a selective sweep). A positive Tajima's D signifies low levels of both low and high frequency polymorphisms, indicating a decrease in population size and/or balancing selection. However, calculating a conventional "p-value" associated with any Tajima's D value that is obtained from a sample is impossible. Briefly, this is because there is no way to describe the distribution of the statistic that is independent of the true, and unknown, theta parameter (no pivot quantity exists). To circumvent this issue, several options have been proposed.

Value of Tajima's D	Mathematical reason	Biological interpretation 1	Biological interpretation 2
Tajima's D=0	Theta-Pi equivalent to Theta-k (Observed=Expected). Average Heterozygosity= # of Segregating sites.	Observed variation similar to expected variation	Population evolving as per mutation-drift equilibrium. No evidence of selection
Tajima's D<0	Theta-Pi less than Theta-k (Observed<Expected). Fewer haplotypes (lower average heterozygosity) than # of segregating sites.	Rare alleles abundant (excess of rare alleles)	Recent selective sweep, population expansion after a recent bottleneck, linkage to a swept gene
Tajima's D>0	Theta-Pi greater than Theta-k (Observed>Expected). More haplotypes (more average heterozygosity)than # of segregating sites.	Rare alleles scarce (lack of rare alleles)	Balancing selection, sudden population contraction

However, this interpretation should be made only if the D-value is deemed statistically significant.

Determining significance

When performing a statistical test such as Tajima's D, the critical question is whether the value calculated for the statistic is unexpected under a null process. For Tajima's D, the magnitude of the statistic is expected to increase the more the data deviates from a pattern expected under a population evolving according to the standard coalescent model.

Tajima (1989) found an empirical similarity between the distribution of the test statistic and a beta distribution with mean zero and variance one. He estimated theta by taking Watterson's estimator and dividing it by the number of samples. Simulations have shown this distribution to be conservative, and now that the computing power is more readily available this approximation is not frequently used.

A more nuanced approach was presented in a paper by Simonsen et al. These authors advocated constructing a confidence interval for the true theta value, and then performing a grid search over this interval to obtain the critical values at which the statistic is significant below a particular alpha value. An alternative approach is for the investigator to perform the grid search over the values of theta which they believe to be plausible based on their knowledge of the organism under study. Bayesian approaches are a natural extension of this method.

A very rough rule of thumb to significance is that values greater than +2 or less than -2 are likely to be significant. This rule is based on an appeal to asymptotic properties of some statistics, and thus +/- 2 does not actually represent a critical value for a significance test.

Finally, genome wide scans of Tajima's D in sliding windows along a chromosomal segment are often performed. With this approach, those regions that have a value of D that greatly deviates from the bulk of the empirical distribution of all such windows are reported as significant. This method does not assess significance in the traditional statistical sense, but is quite powerful given a large genomic region, and is unlikely to falsely identify interesting regions of a chromosome if only the greatest outliers are reported.

See also

• Fay and Wu's H

Notes

• Hartl, Daniel L.; Clark, Andrew G. (2007). Principles of Population Genetics (4th ed.). Sinauer Associates. ISBN 978-0878933082. 

External links

Computational tools:

• DNAsp (Windows)
• Variscan (Mac OS X, Linux, Windows)
• Arlequin (Windows)
• Online view of Tajima's D values in human genome
• Python3 package for computation of Tajima's D
• MEGA4 or MEGA5
• Bio::PopGen::Statistics in BioPerl

• A video explanation of Tajima's D, and its application to DNA sequences, is available online.