 BRIEF REPORTS
 Open Access
 Published:
AdmixKJump: identifying population structure in recently diverged groups
Source Code for Biology and Medicine volume 10, Article number: 2 (2015)
Abstract
Motivation
Correctly modeling population structure is important for understanding recent evolution and for association studies in humans. While preexisting knowledge of population history can be used to specify expected levels of subdivision, objective metrics to detect population structure are important and may even be preferable for identifying groups in some situations. One such metric for genomic scale data is implemented in the crossvalidation procedure of the program ADMIXTURE, but it has not been evaluated on recently diverged and potentially cryptic levels of population structure. Here, I develop a new method, AdmixKJump, and test both metrics under this scenario.
Findings
I show that AdmixKJump is more sensitive to recent population divisions compared to the crossvalidation metric using both realistic simulations, as well as 1000 Genomes Project European genomic data. With two populations of 50 individuals each, AdmixKJump is able to detect two populations with 100% accuracy that split at least 10KYA, whereas crossvalidation obtains this 100% level at 14KYA. I also show that AdmixKJump is more accurate with fewer samples per population. Furthermore, in contrast to the crossvalidation approach, AdmixKJump is able to detect the population split between the Finnish and Tuscan populations of the 1000 Genomes Project.
Conclusion
AdmixKJump has more power to detect the number of populations in a cohort of samples with smaller sample sizes and shorter divergence times.
Availability
A java implementation can be found at https://sites.google.com/site/igsevolgenomicslab/home/downloads
Introduction
Correctly identifying population structure is important both to understand population history and to mitigate potential confounding signals in association analyzes in molecular epidemiology [1]. Recent population divisions can be statistically difficult to recognize as there has not been substantial time for the groups to differentiate. Objective methods to identify recent population divisions are needed. STRUCTURE [2] was the first to do this, but its Bayesian framework is not computationally efficient with whole genome sequence data. ADMIXTURE [3] has implemented a crossvalidation approach to select the correct number of K (i.e. clusters or putative populations) [4], but how this statistic performs on recent population divisions with realistic simulations has yet to be evaluated. Also, alternative approaches may be better suited to recent demographic events. In this paper, I present an implementation and adaptation of the “jump method" of Sugar and James [5] for the problem of identifying populations in genomic sequence data and termed this method the AdmixKJump approach and compare its performance with the crossvalidation approach.
Method
To better facilitate ease of use and comparison with ADMIXTURE, I have implemented AdmixKJump in java with input files that can be taken directly from ADMIXTURE’s output. Parameters can be estimated from different sources (e.g. STRUCTURE or ADMIXTURE) and then used here to identify the number of clusters.
The method makes use of an informationtheoretical approach where the distortions for any given K (d _{ K }) can be calculated using the mean squared error between the genotypes and those predicted by the admixture model:
where for individual i we sum across all M markers (typically single nucleotide variants) indexed by l, and K clusters indexed by k. g represents the genotypes in the form 0, 1, or 2; $\hat {p}$ represents the estimated allele frequency for a specific cluster; and $\hat {q}$ represents the modeled proportion of each individual to each cluster and is usually interpreted as the ancestry percentage [2,3]. $\hat {d}_{K}$ is calculated as the average of $\hat {d}_{\textit {Ki}}$ for all N individuals. Note, Sugar and James [5] originally formulate the d _{ K } value with the Mahalanobis distance, but simplified it to the mean squared error because of the complexity of calculating the covariance matrix. I have made the same simplification in this implementation.
The jump statistic (J _{ K }) is a measure to identify the “elbow” in the monotonically decreasing values of $\hat {d}_{K}$ as the value of K increases. As per Sugar and James [5], it is calculated using a transformation value Y as:
and the estimated number of clusters (K ^{∗}) is selected by:
The transformation value shifts the focus to a particular part of the $\hat {d}_{K}$ distortion curve, e.g. smaller values bias towards a lower K ^{∗}. To mitigate the subjectivity of selecting Y, I modified the selection of K ^{∗} by estimating the lowest possible value of Y that would select for a given K (Y ^{∗}[k]). K ^{∗} is then selected by taking the largest value of K where Y ^{∗}[K+1]−Y ^{∗}[K]>0, or in other words the largest K where some value of Y supports its selection. This is the metric I evaluate in simulation and with real data.
Testing
An important parameter space for these methods is recent population splits such as within continent population divisions. Accordingly, I test how each method fares in identifying the correct number of clusters as a function of time. I use a coalescent simulation framework based on realistic parameters from the Exome Sequencing Project [6,7] to generate whole genome sequences (i.e. 3,000 megabases) with two populations. I used these parameters with the coalescent simulator MSMS [8] to generate the data. The perl script that runs these parameters and imputes them to MSMS can be found in the program distribution. One extra parameter is added which allows me to vary the split time between two populations. This value ranged from 0 to 50K with 50 replicates for every 2K years. Sample size per population was also set to 10, 30, and 50 (20 and 40 not shown, but are consistent). The generated data was subsequently filtered with linkagedisequilibrium pruning and the removal of all singletons as is typical in admixture analysis [2,3]. I evaluated the accuracy of AdmixKJump and crossvalidation by the number of times it correctly identified K ^{∗}=2.
I then apply both metrics to pairwise European populations of the 1000 Genomes Project [9]. I filter individuals, sequenced using SOLiD technology which add additional structure to the data from technical artifacts (see Figure S4 of [9] for evidence of this effect). I also excluded the Spanish population, as after filtering, only six samples remained. I then LD pruned the remaining data. This resulted in over 600K single nucleotide variants for 347 samples from 4 populations.
Results and discussion
In simulation, I find that population structure signals evaporate after the exponential expansion in human population size at about 5 thousand years ago (KYA). One potential explanation for this lack of signal is the reduced effect of genetic drift due to increased population size. With a sample size (N) of 50 for each of two populations, AdmixKJump reaches 100% accuracy at 10KYA, whereas the crossvalidation metric obtains 100% power at about 14KYA. The data generated with the test demographic model produce data with an average F _{ ST } of 0.009 for 10KYA and 0.015 for 14KYA.
I also find that the new measure has more power with smaller sample sizes, for instance N=30 is 100% at 12KYA for AdmixKJump (see Figure 1).
With the 1000 Genomes populations, the crossvalidation approach identifies K ^{∗}=1 for all pairs. The AdmixKJump algorithm finds similar results for all comparisons except the Finnish/Tuscan pair, where it finds K ^{∗}=2, consistent with the known population bottleneck of the Finnish population [10], the greatest geographic separation within these populations, and larger F _{ ST } values (see Table 1). Thus, for closely related populations AdmixKJump is more sensitive.
Conclusions
I have developed a new and powerful approach to classify population structure. I evaluated the current standard metric for large scale data sets, crossvalidation, and found both in simulations and 1000 Genomes data that AdmixKJump is more powerful in recently diverged populations and with smaller sample sizes.
Abbreviations
 KYA:

Thousand years ago
References
 1
O’Connor TD, Kiezun A, Bamshad M, Rich SS, Smith JD, Turner E, et al.Finescale patterns of population stratification confound rare variant association tests. PloS one. 2013; 8(7):65834.
 2
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155(2):945–59.
 3
Alexander DH, Novembre J, Lange K. Fast modelbased estimation of ancestry in unrelated individuals. Genome Res. 2009; 19(9):1655–64.
 4
Alexander DH, Lange K. Enhancements to the admixture algorithm for individual ancestry estimation. BMC Bioinformatics. 2011; 12(1):246.
 5
Sugar CA, James GM. Finding the number of clusters in a dataset. J Am Stat Assoc. 2003; 98(463):750–63.
 6
Tennessen JA, Bigham AW, O’Connor TD, Fu W, Kenny EE, Gravel S, et al.Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2012; 337(6090):64–9.
 7
Fu W, O’Connor TD, Jun G, Kang HM, Abecasis G, Leal SM, et al.Analysis of 6,515 exomes reveals the recent origin of most human proteincoding variants. Nature. 2012; 493(7431):216–20.
 8
Ewing G, Hermisson J. MSMS: a coalescent simulation program including recombination, demographic structure and selection at a single locus. Bioinformatics. 2064; 26(16).
 9
1000 Genomes Project Consortium. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012; 491:56–65.
 10
Workman P, Mielke J, Nevanlinna H. The genetic structure of finland. Am J Phys Anthropol. 1976; 44(2):341–67.
Acknowledgments
I want to thank Joshua Akey and Jessica Chong for helpful comments, and Wenqing Fu for providing the initial MSMS commands for the simulations.
Author information
Additional information
Competing interests
The author declares that he has no competing interests.
Authors’ contributions
Study design, implementation, analysis, and writing of the manuscript were completed by TDO.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Admixture
 Population genetics
 1000 Genomes project
 Fine scale population structure