Allele frequency spectrum

In population genetics, the allele frequency spectrum, sometimes called the site frequency spectrum, is the distribution of the allele frequencies of a given set of loci (often SNPs) in a population or sample. ^{[1] [2] [3] [4]} Because an allele frequency spectrum is often a summary of or compared to sequenced samples of the whole population, it is a histogram with size depending on the number of sequenced individual chromosomes. Each entry in the frequency spectrum records the total number of loci with the corresponding derived allele frequency. Loci contributing to the frequency spectrum are assumed to be independently changing in frequency. Furthermore, loci are assumed to be biallelic (that is, with exactly two alleles present), although extensions for multiallelic frequency spectra exist. ^[5]

Many summary statistics of observed genetic variation are themselves summaries of the allele frequency spectrum, including estimates of ${\displaystyle \theta }$ such as Watterson's ${\displaystyle \theta _{W}}$ and Tajima's ${\displaystyle \theta _{\pi }}$, Tajima's D, Fay and Wu's H and the fixation index ${\displaystyle F_{ST}}$. ^[6]

Example

The allele frequency spectrum from a sample of ${\displaystyle n}$ chromosomes is calculated by counting the number of sites with derived allele frequencies ${\displaystyle 1\leq i\leq n-1}$. For example, consider a sample of ${\displaystyle n=6}$ individuals with eight observed variable sites. In this table, a 1 indicates that the derived allele is observed at that site, while a 0 indicates the ancestral allele was observed.

	SNP 1	SNP 2	SNP 3	SNP 4	SNP 5	SNP 6	SNP 7	SNP 8
Sample 1	0	1	0	0	0	0	1	0
Sample 2	1	0	1	0	0	0	1	0
Sample 3	0	1	1	0	0	1	0	0
Sample 4	0	0	0	0	1	0	1	1
Sample 5	0	0	1	0	0	0	1	0
Sample 6	0	0	0	1	0	1	1	0
Total	1	2	3	1	1	2	5	1

The allele frequency spectrum can be written as the vector ${\displaystyle \mathbf {x} =(x_{1},x_{2},x_{3},x_{4},x_{5})}$, where ${\displaystyle x_{i}}$ is the number of observed sites with derived allele frequency ${\displaystyle i}$. In this example, the observed allele frequency spectrum is ${\displaystyle (4,2,1,0,1)}$, due to four instances of a single observed derived allele at a particular SNP loci, two instances of two derived alleles, and so on.

Calculation

The expected allele frequency spectrum may be calculated using either a coalescent or diffusion approach. ^{[7] [8]} The demographic history of a population and natural selection affect allele frequency dynamics, and these effects are reflected in the shape of the allele frequency spectrum. For the simple case of selective neutral alleles segregating in a population that has reached demographic equilibrium (that is, without recent population size changes or gene flow), the expected allele frequency spectrum ${\displaystyle \mathbf {x} =(x_{1},\ldots ,x_{n-1})}$ for a sample of size ${\displaystyle n}$ is given by

${\displaystyle x_{i}=\theta {\frac {1}{i}},}$

where ${\displaystyle \theta =2N\mu }$ is the population scaled mutation rate. Deviations from demographic equilibrium or neutrality will change the shape of the expected frequency spectrum.

Calculating the frequency spectrum from observed sequence data requires one to be able to distinguish the ancestral and derived (mutant) alleles, often by comparing to an outgroup sequence. For example in human population genetic studies, the homologous chimpanzee reference sequence is typically used to estimate the ancestral allele. However, sometimes the ancestral allele cannot be determined, in which case the folded allele frequency spectrum may be calculated instead. The folded frequency spectrum stores the observed counts of the minor (most rare) allele frequencies. The folded spectrum can be calculated by binning together the ${\displaystyle i}$th and ${\displaystyle (n-i)}$th entries from the unfolded spectrum, where ${\displaystyle n}$ is the number of sampled individuals.

Multi-population allele frequency spectrum

The joint allele frequency spectrum (JAFS) is the joint distribution of allele frequencies across two or more related populations. The JAFS for ${\displaystyle d}$ populations, with ${\displaystyle n_{j}}$ sampled chromosomes in the ${\displaystyle j}$th population, is a ${\displaystyle d}$-dimensional histogram, in which each entry stores the total number of segregating sites in which the derived allele is observed with the corresponding frequency in each population. Each axis of the histogram corresponds to a population, and indices run from ${\displaystyle 0\leq i\leq n_{j}}$ for the ${\displaystyle j}$th population. ^{[9] [10]}

Example

Suppose we sequence diploid individuals from two populations, 4 individuals from population 1 and 2 individuals from population 2. The JAFS would be a ${\displaystyle 9\times 5}$ matrix, indexed from zero. The ${\displaystyle [3,2]}$ entry would record the number of observed polymorphic loci with derived allele frequency 3 in population 1 and frequency 2 in population 2. The ${\displaystyle [1,0]}$ entry would record those loci with observed frequency 1 in population 1, and frequency 0 in population 2. The ${\displaystyle [8,3]}$ entry would record those loci with the derived allele fixed in population 1 (seen in all chromosomes), and with frequency 3 in population 2.

Applications

The shape of the allele frequency spectrum is sensitive to demography, such as population size changes, migration, and substructure, as well as natural selection. By comparing observed data summarized in a frequency spectrum to the expected frequency spectrum calculated under a given demographic and selection model, one can assess the goodness of fit of that the model to the data, and use likelihood theory to estimate the best fit parameters of the model.

For example, suppose a population experienced a recent period of exponential growth and ${\displaystyle n}$ sample sequences were obtained from the population at the end of the growth and the observed (data) allele frequency spectrum was calculated using putatively neutral variation. The demographic model would have parameters for the exponential growth rate ${\displaystyle \rho }$, the time ${\displaystyle T}$ for which the growth occurred, and a reference population size ${\displaystyle N_{ref}}$, assuming that the population was at equilibrium when the growth began. The expected frequency spectrum for a given parameter set ${\displaystyle (\rho ,T,N_{ref})}$ can be obtained using either diffusion or coalescent theory, and compared to the data frequency spectrum. The best fit parameters can be found using maximum likelihood.

This approach has been used to infer demographic and selection models for many species, including humans. For example, Marth et al. (2004) used the single population allele frequency spectra for a group of Africans, Europeans, and Asians to show that population bottlenecks have occurred in the Asian and European demographic histories, but not in the Africans. ^[11] More recently, Gutenkunst et al. (2009) used the joint allele frequency spectrum for these same three populations to infer the time at which the populations diverged and the amount of subsequent ongoing migration between them (see out of Africa hypothesis). ^[10] Additionally, these methods may be used to estimate patterns of selection from allele frequency data. For example, Boyko et al. (2008) inferred the distribution of fitness effects for newly arising mutations using human polymorphism data that controlled for the effects of non-equilibrium demography. ^[12]

References

1. Fisher, Ronald A. (1930). "The distribution of gene ratios for rare mutations". Proceedings of the Royal Society of Edinburgh. 50: 205–220.
2. Wright, Sewall (1938). "The distribution of gene frequencies under irreversible mutation". Proc. Natl. Acad. Sci. USA. 24 (7): 253–259. Bibcode:1938PNAS...24..253W. doi: 10.1073/pnas.24.7.253 . PMC   1077089 . PMID   16577841.
3. Kimura, Motoo (1964). "Diffusion models in population genetics". J. Appl. Probab. 1 (2): 177–232. doi:10.2307/3211856. JSTOR   3211856. S2CID   86705023.
4. Evans, Steven N.; Shvets, Yelena; Slatkin, Montgomery (2007). "Non-equilibrium theory of the allele frequency spectrum". Theoretical Population Biology. 71 (1): 109–119. arXiv: q-bio/0604010 . doi:10.1016/j.tpb.2006.06.005. PMID   16887160. S2CID   924340.
5. Jenkins, Paul A.; Mueller, Jonas W.; Song, Yun S. (2014). "General triallelic frequency spectrum under demographic models with variable population size". Genetics. 196 (1): 295–311. arXiv: 1310.3444 . doi:10.1534/genetics.113.158584. PMC   3872192 . PMID   24214345.
6. Durrett, Rick (2008). Probability Models for DNA Sequence Evolution (PDF) (2 ed.).
7. Wakeley, John (22 April 2016). Coalescent Theory: An Introduction. Roberts & Company Publishers. ISBN   978-0974707754.
8. Crow, James F.; Kimura, Motoo (1970). An introduction to population genetics theory ([Reprint] ed.). New Jersey: Blackburn Press. ISBN   9781932846126.
9. Chen, H.; Green, R. E.; Paabo, S.; Slatkin, M. (29 July 2007). "The Joint Allele-Frequency Spectrum in Closely Related Species". Genetics. 177 (1): 387–398. doi:10.1534/genetics.107.070730. PMC   2013700 . PMID   17603120.
10. Gutenkunst, Ryan N.; Hernandez, Ryan D.; Williamson, Scott H.; Bustamante, Carlos D. (23 October 2009). "Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data". PLOS Genetics. 5 (10): e1000695. arXiv: 0909.0925 . doi: 10.1371/journal.pgen.1000695 . PMC   2760211 . PMID   19851460.
11. Marth, Gabor T.; Czabarka, Eva; Murvai, Janos; Sherry, Stephen T. (1 January 2004). "The Allele Frequency Spectrum in Genome-Wide Human Variation Data Reveals Signals of Differential Demographic History in Three Large World Populations". Genetics. 166 (1): 351–372. doi:10.1534/genetics.166.1.351. PMC   1470693 . PMID   15020430.
12. Boyko, Adam R.; Williamson, Scott H.; Indap, Amit R.; Degenhardt, Jeremiah D.; Hernandez, Ryan D.; Lohmueller, Kirk E.; Adams, Mark D.; Schmidt, Steffen; Sninsky, John J.; Sunyaev, Shamil R.; White, Thomas J.; Nielsen, Rasmus; Clark, Andrew G.; Bustamante, Carlos D. (30 May 2008). "Assessing the Evolutionary Impact of Amino Acid Mutations in the Human Genome". PLOS Genetics. 4 (5): e1000083. doi: 10.1371/journal.pgen.1000083 . PMC   2377339 . PMID   18516229.