The fixation index (FST) is a measure of population differentiation due to genetic structure. It is frequently estimated from genetic polymorphism data, such as single-nucleotide polymorphisms (SNP) or microsatellites. Developed as a special case of Wright's F-statistics, it is one of the most commonly used statistics in population genetics. Its values range from 0 to 1, with 0 being no differentiation and 1 being complete differentiation.
This comparison of genetic variability within and between populations is frequently used in applied population genetics. The values range from 0 to 1. A zero value implies complete panmixia; that is, that the two populations are interbreeding freely. A value of one implies that all genetic variation is explained by the population structure, and that the two populations do not share any genetic diversity.
For idealized models such as Wright's finite island model, FST can be used to estimate migration rates. Under that model, the migration rate is
where m is the migration rate per generation, and is the mutation rate per generation. [1]
The interpretation of FST can be difficult when the data analyzed are highly polymorphic. In this case, the probability of identity by descent is very low and FST can have an arbitrarily low upper bound, which might lead to misinterpretation of the data. Also, strictly speaking FST is not a distance in the mathematical sense, as it does not satisfy the triangle inequality.
Two of the most commonly used definitions for FST at a given locus are based on 1) the variance of allele frequencies among populations, and on 2) the probability of identity by descent.
If is the average frequency of an allele in the total population, is the variance in the frequency of the allele among different subpopulations, weighted by the sizes of the subpopulations, and is the variance of the allelic state in the total population, FST is defined as [2]
Wright's definition illustrates that FST measures the amount of genetic variance that can be explained by population structure. This can also be thought of as the fraction of total diversity that is not a consequence of the average diversity within subpopulations, where diversity is measured by the probability that two randomly selected alleles are different, namely . If the allele frequency in the th population is and the relative size of the th population is , then
Alternatively, [3]
where is the probability of identity by descent of two individuals given that the two individuals are in the same subpopulation, and is the probability that two individuals from the total population are identical by descent. Using this definition, FST can be interpreted as measuring how much closer two individuals from the same subpopulation are, compared to the total population. If the mutation rate is small, this interpretation can be made more explicit by linking the probability of identity by descent to coalescent times: Let T0 and T denote the average time to coalescence for individuals from the same subpopulation and the total population, respectively. Then,
This formulation has the advantage that the expected time to coalescence can easily be estimated from genetic data, which led to the development of various estimators for FST.
In practice, none of the quantities used for the definitions can be easily measured. As a consequence, various estimators have been proposed. A particularly simple estimator applicable to DNA sequence data is: [4]
where and represent the average number of pairwise differences between two individuals sampled from different sub-populations () or from the same sub-population (). The average pairwise difference within a population can be calculated as the sum of the pairwise differences divided by the number of pairs. However, this estimator is biased when sample sizes are small or if they vary between populations. Therefore, more elaborate methods are used to compute FST in practice. Two of the most widely used procedures are the estimator by Weir & Cockerham (1984), [5] or performing an Analysis of molecular variance. A list of implementations is available at the end of this article.
FST values depend strongly on the choice of populations. Closely related ethnic groups, such as the Danes vs. the Dutch, or the Portuguese vs. the Spaniards show values significantly below 1%, indistinguishable from panmixia. Within Europe, the most divergent ethnic groups have been found to have values of the order of 7% (Sámi vs. Sardinians).
Larger values are found if highly divergent homogenous groups are compared: the highest such value found was at close to 46%, between Mbuti and Papuans. [6]
A genetic distance of 0.125 implies that kinship between unrelated individuals of the same ancestry relative to the world population is equivalent to kinship between half siblings in a randomly mating population. This also implies that if a human from a given ancestral population has a mixed half-sibling, that human is closer genetically to an unrelated individual of their ancestral population than to their mixed half-sibling. [7]
In their study The History and Geography of Human Genes (1994), Cavalli-Sforza, Menozzi and Piazza provide some of the most detailed and comprehensive estimates of genetic distances between human populations, within and across continents. Their initial database contains 76,676 gene frequencies (using 120 blood polymorphisms), corresponding to 6,633 samples in different locations. By culling and pooling such samples, they restrict their analysis to 491 populations.
They focus on aboriginal populations that were at their present location at the end of the 15th century when the great European migrations began. [8] When studying genetic difference at the world level, the number is reduced to 42 representative populations, aggregating subpopulations characterized by a high level of genetic similarity. For these 42 populations, Cavalli-Sforza and coauthors report bilateral distances computed from 120 alleles. Among this set of 42 world populations, the greatest genetic distance observed is between Mbuti Pygmies and Papua New Guineans, where the Fst distance is 0.4573, while the smallest genetic distance (0.0021) is between the Danish and the English.
When considering more disaggregated data for 26 European populations, the smallest genetic distance (0.0009) is between the Dutch and the Danes, and the largest (0.0667) is between the Lapps and the Sardinians. The mean genetic distance among the 861 available pairings of the 42 selected populations was found to be 0.1338.[ page needed ].
The following table shows Fst calculated by Cavalli-Sforza (1994) for some populations:
Bantu | Nio-Saharan | W. African | Mbuti | Japanese | Korean | Thai | Filipino | S. Chinese | Danish | English | Melanesian | New Guinean | Australian | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Bantu | 0.00 | |||||||||||||
Nio-Saharan | 0.01 | 0.00 | ||||||||||||
W. African | 0.02 | 0.02 | 0.00 | |||||||||||
Mbuti | 0.07 | 0.08 | 0.08 | 0.00 | ||||||||||
Japanese | 0.24 | 0.25 | 0.23 | 0.31 | 0.00 | |||||||||
Korean | 0.27 | 0.24 | 0.18 | 0.30 | 0.01 | 0.00 | ||||||||
Thai | 0.34 | 0.30 | 0.25 | 0.39 | 0.07 | 0.06 | 0.00 | |||||||
Filipino | 0.29 | 0.30 | 0.23 | 0.38 | 0.10 | 0.12 | 0.06 | 0.00 | ||||||
S. Chinese | 0.30 | 0.29 | 0.20 | 0.34 | 0.05 | 0.05 | 0.01 | 0.03 | 0.00 | |||||
Danish | 0.17 | 0.17 | 0.15 | 0.15 | 0.12 | 0.09 | 0.13 | 0.13 | 0.13 | 0.00 | ||||
English | 0.23 | 0.18 | 0.15 | 0.24 | 0.12 | 0.10 | 0.11 | 0.11 | 0.12 | 0.00 | 0.00 | |||
Melanesian | 0.34 | 0.31 | 0.27 | 0.40 | 0.11 | 0.11 | 0.08 | 0.05 | 0.06 | 0.14 | 0.16 | 0.00 | ||
New Guinean | 0.34 | 0.33 | 0.28 | 0.46 | 0.12 | 0.14 | 0.18 | 0.13 | 0.15 | 0.16 | 0.16 | 0.07 | 0.00 | |
Australian | 0.33 | 0.36 | 0.27 | 0.43 | 0.06 | 0.07 | 0.13 | 0.13 | 0.11 | 0.14 | 0.15 | 0.11 | 0.10 | 0.00 |
A 2012 study based on International HapMap Project data estimated FST between the three major "continental" populations of Europeans (combined from Utah residents of Northern and Western European ancestry from the CEPH collection and Italians from Tuscany), East Asians (combining Han Chinese from Beijing, Chinese from metropolitan Denver and Japanese from Tokyo, Japan) and Sub-Saharan Africans (combining Luhya of Webuye, Kenya, Maasai of Kinyawa, Kenya and Yoruba of Ibadan, Nigeria). It reported a value close to 12% between continental populations, and values close to panmixia (smaller than 1%) within continental populations. [9]
Europe (CEU) | Sub-Saharan Africa (Yoruba) | East-Asia (Japanese) | |
---|---|---|---|
Sub-Saharan Africa (Yoruba) | 0.153 | ||
East-Asia (Japanese) | 0.111 | 0.190 | |
East-Asia (Chinese) | 0.110 | 0.192 | 0.007 |
Italians | Palestinians | Swedish | Finns | Spanish | Germans | Russians | French | Greeks | |
---|---|---|---|---|---|---|---|---|---|
Palestinians | 0.0064 | ||||||||
Swedish | 0.0064-0.0090 | 0.0191 | |||||||
Finns | 0.0130-0.0230 | 0.0050-0.0110 | |||||||
Spanish | 0.0010-0.0050 | 0.0101 | 0.0040-0055 | 0.0110-0.0170 | |||||
Germans | 0.0029-0.0080 | 0.0136 | 0.0007-0.0010 | 0.0060-0.0130 | 0.0015-0.0030 | ||||
Russians | 0.0088-0.0120 | 0.0202 | 0.0030-0.0036 | 0.0060-0.0120 | 0.0070-0.0079 | 0.0030-0.0037 | |||
French | 0.0030-0.0050 | 0.0020 | 0.0080-0.0150 | 0.0010 | 0.0010 | 0.0050 | |||
Greeks | 0.0000 | 0.0057 | 0.0084 | 0.0035 | 0.0039 | 0.0108 | |||
Pairwise Fst values among several populations based on whole exome sequencing (WES) in 2016: [12]
The Cauchy distribution, named after Augustin Cauchy, is a continuous probability distribution. It is also known, especially among physicists, as the Lorentz distribution, Cauchy–Lorentz distribution, Lorentz(ian) function, or Breit–Wigner distribution. The Cauchy distribution is the distribution of the x-intercept of a ray issuing from with a uniformly distributed angle. It is also the distribution of the ratio of two independent normally distributed random variables with mean zero.
The median of a set of numbers is the value separating the higher half from the lower half of a data sample, a population, or a probability distribution. For a data set, it may be thought of as the “middle" value. The basic feature of the median in describing data compared to the mean is that it is not skewed by a small proportion of extremely large or small values, and therefore provides a better representation of the center. Median income, for example, may be a better way to describe the center of the income distribution because increases in the largest incomes alone have no effect on the median. For this reason, the median is of central importance in robust statistics.
In probability theory and statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real-valued random variable. The general form of its probability density function is The parameter is the mean or expectation of the distribution, while the parameter is the variance. The standard deviation of the distribution is . A random variable with a Gaussian distribution is said to be normally distributed, and is called a normal deviate.
In statistics, the standard deviation is a measure of the amount of variation of the values of a variable about its mean. A low standard deviation indicates that the values tend to be close to the mean of the set, while a high standard deviation indicates that the values are spread out over a wider range. The standard deviation is commonly used in the determination of what constitutes an outlier and what does not.
The weighted arithmetic mean is similar to an ordinary arithmetic mean, except that instead of each of the data points contributing equally to the final average, some data points contribute more than others. The notion of weighted mean plays a role in descriptive statistics and also occurs in a more general form in several other areas of mathematics.
The Allan variance (AVAR), also known as two-sample variance, is a measure of frequency stability in clocks, oscillators and amplifiers. It is named after David W. Allan and expressed mathematically as . The Allan deviation (ADEV), also known as sigma-tau, is the square root of the Allan variance, .
In statistics, maximum likelihood estimation (MLE) is a method of estimating the parameters of an assumed probability distribution, given some observed data. This is achieved by maximizing a likelihood function so that, under the assumed statistical model, the observed data is most probable. The point in the parameter space that maximizes the likelihood function is called the maximum likelihood estimate. The logic of maximum likelihood is both intuitive and flexible, and as such the method has become a dominant means of statistical inference.
In population genetics, F-statistics describe the statistically expected level of heterozygosity in a population; more specifically the expected degree of (usually) a reduction in heterozygosity when compared to Hardy–Weinberg expectation.
Directional statistics is the subdiscipline of statistics that deals with directions, axes or rotations in Rn. More generally, directional statistics deals with observations on compact Riemannian manifolds including the Stiefel manifold.
Estimation theory is a branch of statistics that deals with estimating the values of parameters based on measured empirical data that has a random component. The parameters describe an underlying physical setting in such a way that their value affects the distribution of the measured data. An estimator attempts to approximate the unknown parameters using the measurements. In estimation theory, two approaches are generally considered:
In Bayesian statistics, a maximum a posteriori probability (MAP) estimate is an estimate of an unknown quantity, that equals the mode of the posterior distribution. The MAP can be used to obtain a point estimate of an unobserved quantity on the basis of empirical data. It is closely related to the method of maximum likelihood (ML) estimation, but employs an augmented optimization objective which incorporates a prior distribution over the quantity one wants to estimate. MAP estimation can therefore be seen as a regularization of maximum likelihood estimation.
In probability theory and directional statistics, the von Mises distribution is a continuous probability distribution on the circle. It is a close approximation to the wrapped normal distribution, which is the circular analogue of the normal distribution. A freely diffusing angle on a circle is a wrapped normally distributed random variable with an unwrapped variance that grows linearly in time. On the other hand, the von Mises distribution is the stationary distribution of a drift and diffusion process on the circle in a harmonic potential, i.e. with a preferred orientation. The von Mises distribution is the maximum entropy distribution for circular data when the real and imaginary parts of the first circular moment are specified. The von Mises distribution is a special case of the von Mises–Fisher distribution on the N-dimensional sphere.
Genetic distance is a measure of the genetic divergence between species or between populations within a species, whether the distance measures time from common ancestor or degree of differentiation. Populations with many similar alleles have small genetic distances. This indicates that they are closely related and have a recent common ancestor.
In estimation theory and decision theory, a Bayes estimator or a Bayes action is an estimator or decision rule that minimizes the posterior expected value of a loss function. Equivalently, it maximizes the posterior expectation of a utility function. An alternative way of formulating an estimator within Bayesian statistics is maximum a posteriori estimation.
In population genetics, fixation is the change in a gene pool from a situation where there exists at least two variants of a particular gene (allele) in a given population to a situation where only one of the alleles remains. That is, the allele becomes fixed. In the absence of mutation or heterozygote advantage, any allele must eventually either be lost completely from the population, or fixed, i.e. permanently established at 100% frequency in the population. Whether a gene will ultimately be lost or fixed is dependent on selection coefficients and chance fluctuations in allelic proportions. Fixation can refer to a gene in general or particular nucleotide position in the DNA chain (locus).
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.
Population structure is the presence of a systematic difference in allele frequencies between subpopulations. In a randomly mating population, allele frequencies are expected to be roughly similar between groups. However, mating tends to be non-random to some degree, causing structure to arise. For example, a barrier like a river can separate two groups of the same species and make it difficult for potential mates to cross; if a mutation occurs, over many generations it can spread and become common in one subpopulation while being completely absent in the other.
In probability theory and directional statistics, a wrapped normal distribution is a wrapped probability distribution that results from the "wrapping" of the normal distribution around the unit circle. It finds application in the theory of Brownian motion and is a solution to the heat equation for periodic boundary conditions. It is closely approximated by the von Mises distribution, which, due to its mathematical simplicity and tractability, is the most commonly used distribution in directional statistics.
Isolation by distance (IBD) is a term used to refer to the accrual of local genetic variation under geographically limited dispersal. The IBD model is useful for determining the distribution of gene frequencies over a geographic region. Both dispersal variance and migration probabilities are variables in this model and both contribute to local genetic differentiation. Isolation by distance is usually the simplest model for the cause of genetic isolation between populations. Evolutionary biologists and population geneticists have been exploring varying theories and models for explaining population structure. Yoichi Ishida compares two important theories of isolation by distance and clarifies the relationship between the two. According to Ishida, Sewall Wright's isolation by distance theory is termed ecological isolation by distance while Gustave Malécot's theory is called genetic isolation by distance. Isolation by distance is distantly related to speciation. Multiple types of isolating barriers, namely prezygotic isolating barriers, including isolation by distance, are considered the key factor in keeping populations apart, limiting gene flow.
In probability theory, a log-Cauchy distribution is a probability distribution of a random variable whose logarithm is distributed in accordance with a Cauchy distribution. If X is a random variable with a Cauchy distribution, then Y = exp(X) has a log-Cauchy distribution; likewise, if Y has a log-Cauchy distribution, then X = log(Y) has a Cauchy distribution.