# HaploPOP: a software that improves population assignment by combining markers into haplotypes

- Nicolas Duforet-Frebourg
^{1, 2, 3}Email author, - Lucie M. Gattepaille
^{4}, - Michael G.B Blum
^{1, 2}and - Mattias Jakobsson
^{4, 5}Email author

**Received: **12 November 2014

**Accepted: **3 July 2015

**Published: **31 July 2015

## Abstract

### Background

In ecology and forensics, some population assignment techniques use molecular markers to assign individuals to known groups. However, assigning individuals to known populations can be difficult if the level of genetic differentiation among populations is small. Most assignment studies handle independent markers, often by pruning markers in Linkage Disequilibrium (LD), ignoring the information contained in the correlation among markers due to LD.

### Results

To improve the accuracy of population assignment, we present an algorithm, implemented in the *HaploPOP* software, that combines markers into haplotypes, without requiring independence. The algorithm is based on the Gain of Informativeness for Assignment that provides a measure to decide if a pair of markers should be combined into haplotypes, or not, in order to improve assignment. Because complete exploration of all possible solutions for constructing haplotypes is computationally prohibitive, our approach uses a greedy algorithm based on windows of fixed sizes. We evaluate the performance of *HaploPOP* to assign individuals to populations using a split-validation approach. We investigate both simulated SNPs data and dense genotype data from individuals from Spain and Portugal.

### Conclusions

Our results show that constructing haplotypes with *HaploPOP* can substantially reduce assignment error. The *HaploPOP* software is freely available as a command-line software at www.ieg.uu.se/Jakobsson/software/HaploPOP/.

## Background

Molecular markers provide powerful approaches in forensic science and ecology to assign individuals into predefined populations [1, 2]. With the advent of new sequencing technologies, the number of available molecular markers in different species is rapidly increasing. At the same time, dense datasets tend to contain increasingly correlated markers because Single Nucleotide Polymorphisms (SNPs) that are physically close on a chromosome, often are in Linkage Disequilibrium (LD). Such correlations are usually perceived as a nuisance factor in statistical analyses since it violates a common assumption of independence among markers. This statistical nuisance can be overcome by pruning SNPs using for example the software *PLINK* [3]. However much information may be lost because of the pruning process. Another approach is to explicitly model the correlation between markers to control for LD [4–6], or to include the pruning process in the statistical analysis [7]. In addition, it has been shown that it can be useful to combine correlated markers into haplotypes to augment the information about population structure at a finer scale [8]. Such an approach is valuable for assignment methods when the level of genetic differentiation among groups is low [9].

Gattepaille and Jakobsson [10] introduced the Gain of Informativeness for Assignment (GIA), which is a statistic measuring the gain in information for population assignment by combining two markers into haplotypes. GIA is based upon an ancestry information criterion that measures to what extent a molecular marker is informative about population assignment [11]. GIA is defined as the difference between the ancestry information carried by two markers and the ancestry information carried by the haplotypes resulting from the combination of the two markers. Building haplotypes with GIA increases correct assignment to predefined populations [10]. However, a major combinatorial challenge arises when using GIA because of the prohibitively large number of pairs of markers that can be combined into haplotypes.

In this article, we present a new algorithm that efficiently uses GIA to build informative haplotypes for population assignment. The algorithm needs *reference* individuals whose population of origin is known. Based on these *reference* individuals, the algorithm uses GIA to construct informative haplotypes. To handle large numbers of markers, we provide a heuristic approach where only markers located within the same genomic region can be combined to form haplotypes. Combining markers into haplotypes is a recursive process so that haplotypes can result from the combination of two or more markers. The raw genotype data are recoded into multi-allelic haplotype data and the new data file containing both genotypic and haplotypic information can be used to assign individuals to populations based on for instance Principal Component Analysis, or model-based assignment approaches [12, 13].

Because the construction of haplotypes uses predefined populations, there is a risk of overfitting. For example if the evaluation of population assignment is performed with the same individuals that were used to construct the haplotypes, the assignment errors may be underestimated. Additionally, the constructed haplotypes can generate artificial population structure although there is no true stratification among the predefined populations. Both problems arise because the construction of haplotypes can exaggerate the differentiation among populations. To get a fair evaluation of population assignment, we implement a split-validation approach where we use different individuals to construct the haplotypes and to evaluate assignment [14]. Haplotypes are built using a subset of the individuals, consisting of a training set. The quality of population assignment can then be assessed using the remaining individuals (the validation set). If the individuals in the validation set cluster with individuals in the training set, there is evidence for some level of population structure, which may not have been detected based only on genotype markers.

Our new algorithm for combining markers into haplotypes is implemented in the software *HaploPOP*. The software is a command-line program written in C. We give examples of how to use *haploPOP* to perform population assignment with SNP data that were simulated from a population divergence model. We also show that *HaploPOP* improves assignment of individuals from Spain and Portugal using the POPRES dataset that contains 447,245 SNPs [15].

## Methods

### Gain of informativeness for assignment

*Gain of Informativeness for Assignment*(

*GIA*) is a one-dimensional statistic that provides a criterion to decide whether markers should be combined into haplotypes in order to improve population assignment [10]. It is based on the

*Informativeness for Assignment*(

*IA*) statistic, which measures how informative a marker is for assigning individuals to different populations [11]. The more different the allele frequencies are in a set of predefined populations, the more informative the marker is to assign individuals of unknown origin to their source population, and the larger is the

*IA*statistic. Denoting by

*K*the number of populations, by

*N*the number of alleles of the marker under consideration, by \(p_{j}^{(i)}\) the frequency of allele

*j*in population

*i*, and by \(\overline {p}_{j}\) the average frequency of allele

*j*across all populations, the

*IA*statistic is computed as follows [11]

*M*

_{1}and

*M*

_{2}, the question is whether combining

*M*

_{1}with

*M*

_{2}into a haplotype marker

*H*improves the assignment of individuals to predefined populations.

*GIA*computes the difference between the informativeness for assignment of

*H*and the sum of the informativeness of

*M*

_{1}and

*M*

_{2}

If *GIA* is positive, it suggests that population assignment is improved by considering haplotype *H* instead of using the two markers separately. However, if *GIA* is negative, there is no advantage of combining the two markers into a haplotype. In particular, it can be shown that if the two markers are in linkage equilibrium, *GIA* is expected to be non-positive [10].

### Maximizing the informativeness for assignment

*n*individuals at

*l*molecular markers (

*M*

_{1},...

*M*

_{ l }). We also assume that the dataset has been phased, where all individuals have been phased together in one go to avoid introducing any haplotype difference due to phasing (note that there may still be switch errors from the phasing, but these should affect all individuals similarly). The approach implemented in the software

*HaploPOP*builds a set of haplotypes that increases the total informativeness for assignment contained in the genotype data. To find the optimal haplotype set

*Γ*

_{0}, we address the maximization problem

where *P*
*a*
*r*
*t*(*M*
_{1}...*M*
_{
l
}) is the set of all possible partitions of the *l* markers. The number of partitions in a dataset of *l* markers is given by Bell’s number [16]. Because this number is large, we cannot evaluate the objective function for all possible partitions. A commonly used heuristic is to apply a greedy strategy, although it can perform arbitrarily good. In the case of increasing Informativeness for Assignment, the resulting haplotypes always provide genetic data with augmented informativeness. Because the cost of this algorithm increases rapidly with the number of genetic markers, we limit potential combinations of markers within windows of fixed size.

### Algorithm

In a first step, the algorithm constructs haplotypes from the phased genotype file of individuals with known origin and returns a haplotype coding file that provides the correspondence between haplotypes and initial markers. This is the LEARN option of *HaploPOP*. The construction of haplotypes is constrained by a predefined window-size. The set (*M*
_{1},…*M*
_{
l
}) of markers is divided into subsets of contiguous markers corresponding to the genomic windows. Haplotypes are constrained to be combinations of markers of the same window. The window size is chosen by the user and can be defined based on number of markers, on base pairs, or genetic distance. By choosing genetic distances, one can account for non-uniform recombination rates. In every window, the *GIA* statistic is computed for all pairs of markers, and the pair with the greatest *GIA* value is merged to form a haplotype. Combinations proceed recursively until there is no pair of markers for which *GIA* >0 (or a certain positive user-defined threshold). A particular haplotype-loci formed by a combination of markers is thereafter treated as a (potentially multi-allelic) marker of the particular window and can be combined with other markers in a recursive manner.

*n*the number of

*reference*individuals whose population of origin is known, and by

*l*the total number of initial markers. The greedy algorithm proceeds as follows:

- 1.
divide the 2

*n*×*l*data matrix in contiguous windows. - 2.for every window do
- (a)
Calculate

*GIA*for all pair of markers. - (b)while for all markers
*M*and*M*’, \(\displaystyle {max}_{M, M^{\prime }} (\textit {GIA}(M, M^{\prime }) > 0)\), do- i.
\(\displaystyle (M_{0},M^{\prime }_{0}) = {argmax}_{M, M^{\prime }} \textit {GIA}(M,M^{\prime }) > 0\).

- ii.
Combine the markers

*M*_{0}and \(M^{\prime }_{0}\) to form a haplotype marker*H*_{0}. - iii.
Remove the

*GIA*statistics involving*M*_{0}and \(M^{\prime }_{0}\) and compute the new*GIA*statistics with pairs of markers that include*H*_{0}.

- i.

- (a)

At every end of the inner loop, the algorithm partitions the markers into a set of haplotypes that increases the score of the objective function (3). It stops when no additional pairwise combination improves the total score of the partition. A warning is raised when the number of haplotype-alleles reaches the number of chromosomes 2*n* making haplotypes useless because they become private to every individual and do not provide any useful information for assignment.

In a second step, *HaploPOP* combines the SNPs in the initial genotype file into haplotypes according to the combinations of markers constructed at the first step, and generates the haplotype data file. The genotype file can contain individuals of unknown origin that the user is trying to determine, as well as the individuals of known or suggested origin used to construct the haplotypes. This corresponds to the APPLY function of *HaploPOP*. When two markers are combined, the resulting haplotype-alleles are coded in a range from 0 to the number of haplotype-alleles minus one, in order of appearance in the list of individuals.

### Window size

*S*markers, the number of windows is

*n*

_{window}=

*l*/

*S*, and the cost

*C*of the algorithm in number of operations is

where *K* is the number of populations in the data. The algorithm scales very well for genome wide datasets, since for a given window size *S*, the cost of the algorithm is proportional to the number of markers *l* in the data. The term proportional to *S*
^{3} corresponds to the iterative maximum search in all possible pairs of the window. The term *S*
^{3} is an upper bound for a search that is done in the worst case *S* times in a matrix of size *S*×*S* or less.

In the event of choosing a large window size, there may be a large number of haplotype-alleles, which could fit closely to the distribution of haplotype-alleles of individuals in the training set. Such a set of haplotypes would likely perform poorly for other sets of individuals from the same reference population, and reduces accuracy of population assignment. We refer to this phenomenon as overfitting. Limiting the size of the window is one way to avoid overfitting. We demonstrate in the Results section that the window size has a strong impact on the performance of the created combinations of markers, and an optimal value generally exists. The optimal window size depends on multiple factors, including the effective population sizes of the investigated groups and the extent of Linkage Disequilibrium in the groups.

### Split-validation

To validate the gain in assignment accuracy provided by the constructed haplotypes, we implement a split-validation technique [14]. For each population, we randomly split the set of individuals into two subsets consisting of the training subset used to learn the haplotypes and the validation subset used to compute assignment accuracy. It is important that the division between validation and training set is done after phasing. Phasing performed on the two datasets separately could introduce haplotypic differences and weaken the informativeness for assignment of the haplotypes built by the algorithm. To assign individuals to populations, we use Principal Component Analysis (PCA) as implemented in the software EIGENSOFT [12]. For each of the constructed haplotype-loci, we enumerate all haplotype-alleles present in the dataset. We use a presence/absence coding for each haplotype-allele. In particular, we add one column per haplotype-allele and note 1 for a chromosome carrying the allele, and 0 otherwise. The number of principal components we consider equals the number of populations used for constructing the haplotypes minus one [12]. We determine the PC axes using individuals from the training and the validation set. For each individual of the validation set, we compute Euclidean distances on the PC space between this individual and the barycentric coordinates of each population computed from the training set of individuals. We assign individuals to the population that has the closest barycenter. Because the origins of individuals in the validation sets were known for all examples (see below), we can measure the number of incorrectly assigned individuals in these examples. Note that the assessment of individuals to populations depends on the assignment procedure itself (here we use PCA) and that different assignment procedures may lead to different assignment errors (see [10] for a comparison of different assignment strategies). However, since we are primarily interested in the comparison between assignment using the raw genotype data and assignment using combined markers found with *HaploPOP*, we focus on a single assignment approach based on PCA. With this assignment approach each haplotype-allele is treated as a unique allele with the same relationship to all other alleles.

## Results and discussion

We evaluate the performance of the approach and the *HaploPOP* software on both simulated and empirical data.

### Application to simulated data

We evaluated the assignment approach and the *HaploPOP* software for simulated data generated by the software *ms* [17]. We simulated 200 kb sequences from a 3-population divergence model. We set the effective population sizes of all populations to *N*
_{
e
}=1,000, the mutation rate to *μ*=0.012, so that *θ*=48, and we considered a sample of 100 individuals in each population. The population divergence between population 2 and population 3 was set to occur at *T*
_{1}=0.025 coalescent time units (or 100 generations) in the past and the population divergence between population 2 and population 1, was set to occur at *T*
_{2}=0.05 coalescent time units (or 200 generations) in the past. We generated four datasets for a hypothetical 200 kb region with effective recombination rates of *ρ*=30,60,120, and 240, and replicated this procedure 10 times for each value of the recombination rate.

*HaploPOP*, we used a split-validation technique. The training set and the validation set contained each 50 randomly chosen individuals in each population. Figure 1 a shows that assignment accuracy improves by constructing haplotypes. The assignment error decreases as the window size increases (up to a certain level). However, most of the improvement occurs when moving from genotypes to haplotypes spanning up to 50 kb (Fig. 1 a). Compared to the error of assignment obtained with genotype data, constructing haplotypes decreases the error by 20−70

*%*depending on the recombination rate. Constructing haplotypes for 200 kb windows compared to 50 kb windows only reduces the error by at most an additional 8

*%*.

Furthermore, we find that the mean incorrect assignment is lower with greater recombination rates. This emphasizes the fact that strongly correlated polymorphisms tend to carry less information for assignment than the same number of independant polymorphisms. Since simulated sequences have on average the same number of SNPs, sequences with a greater recombination rate carry more informativeness for assignment.

For a fixed window size of 50 kb, we construct the haplotypes for different sizes of training sets ranging from 2 to 60 individuals. When comparing the mean incorrect assignment of individuals from the validation set, we find a decay of MIA with increasing numbers of individuals in the training set for all recombination rates. However, when using a fraction of individuals greater than 10 *%* of the overall population, the change in MIA is minimal (Fig. 1
b). Hence, even a fairly small fraction of individuals in a sample can be used to accurately train the algorithm.

### Application on human data

*HaploPOP*improves population assignment of respectively 133 and 125 self-reported Spanish and Portuguese individuals from the POPRES dataset, which contains 447,245 SNPs [15]. We first phased the data using

*fastPhase*[10, 18]. No pruning of SNPs were performed because our aim is to capture the information for assignment for all markers, including markers in LD. All markers were therefore retained and used to build haplotypes. Considering the first two PCs based on all SNP-genotype data, we found that the two populations cannot be distinguished (Fig. 2). A thorough PCA exploration of all European individuals of the POPRES collection was further unable to distinguish between Spanish and Portuguese individuals [19]. Using

*HaploPOP*, we constructed the haplotypes that are informative to discriminate between Spanish and Portuguese individuals. We randomly selected 67 Spanish individuals and 63 Portuguese individuals for constructing the training set. We then performed PCA based on the haplotype markers generated by

*HaploPOP*and compute the PC scores for the individuals from the training set and the validation set (Fig. 3 a). On PC1, the Spanish and Portuguese samples from the training set are clearly separated. By contrast, the 66 Spanish individuals and 62 Portuguese individuals from the validation set overlap but a large majority of these individuals are pulled in the direction of their population of origin.

To show that the population labels, Portugal and Spain, correspond to true population differentiation, we generated a control training set. We arbitrarily assign each individual of the training set to a label A or B and by construction, individuals labeled by A (or B) contain both individuals from Spain and Portugal. Using this training set of half the individuals in the dataset, we learn the haplotypes that are informative for discriminating between A and B. Using the validation set, we find that the haplotypes learned with populations A and B cannot distinguish between Spanish and Portuguese ancestry (Fig. 3 b). This analysis shows that the above demonstrated separation between Spanish and Portuguese individuals corresponds to true population differentiation and that the separation is not a consequence of overfitting.

These two examples emphasize two important features of *HaploPOP*. First, the haplotypes constructed based on the training set are very efficient in separating individuals of the training set, regardless of any true stratification between candidate populations. Comparing only the individuals from the training set on a PC plot can either lead to the wrong conclusion that two populations can be distinguished (Fig. 3
b) or at least exaggerate the ability to distinguish between the two populations (Fig. 3
a) because of overfitting. Second, evidence for population structure comes from the ability of the constructed haplotypes to distinguish between individuals that were not used in the training process. If a validation set of individuals can be assigned to the candidate populations, it is a good indication of fine-scale level of stratification between the candidate populations that might be difficult to detect using SNPs only.

*HaploPOP*can be used to study samples where prior belief suggests that there is population structure but SNP-genotype data fail to detect it. When computing the error for assignment of the validation set for the POPRES data, we find that there is an optimal window size (Fig. 4) at which the assignment error can be reduced by 45 %. Intuitively, combining SNPs into haplotypes can only improve the power to assign individuals to groups up to some level: for too large window sizes, we run into overfitting problems where trained haplotypes are well suited to separate the particular individuals of the training set but not of the individuals in the validation set. The optimal value of the window size depends on many factors, such as the extent of linkage disequilibrium in the groups, or the degree of genetic differentiation between groups. The strategy we advocate for choosing the window size is to try different window sizes and to find the minimal assignment error as estimated with a split-validation approach. Such a strategy is computationally costly and requires for each of the chosen window size a run of

*HaploPOP*, where the cost will be dependent of the window size

*S*as described in equation (4).

Recently, many model-based methods have been proposed to assign unlabeled individuals to populations [13, 20–22]. These methods can be used together with *HaploPOP* to reduce the proportion of incorrect assignment, as it is shown with the software *Structure* [23] in a previous article [10]. In this article we focus on using Principal Component Analysis and, from a statistical point of view, model-based approaches and PCA are related [20]. In the case of assigning individuals to labeled populations, we expect that most of these methods will result in similar assignment accuracy.

## Conclusions

In this article, we present a new algorithm that uses the *GIA* statistic to construct haplotypes with a window-based approach. The algorithm is implemented in the command-line software *HaploPOP*. The software allows users to apply a 2-step procedure. First, *HaploPOP* constructs haplotypes that are informative about population assignment from a training set of individuals. Second, *HaploPOP* recodes the genotype data to haplotypes. These new haplotype data can then be used to assign unknown individuals to candidate populations or investigate fine-scale population structure using e.g. PCA. We have shown how constructing haplotypes with *HaploPOP* can substantially reduce mis-assignment of individuals to candidate populations. For SNP data simulated in a 3-population divergence model, the assignment-error was reduced by 20 *%* to 70 *%*. Using the 447,245 SNPs of the POPRES data, the assignment-error was reduced by 45 *%* when trying to distinguish Portuguese from Spanish individuals.

Constructing Haplotypes with *HaploPOP* is a promising approach to assign individuals into populations in forensic science and ecology. It can also confirm prior belief about fine-scale population structure which is a main confounding factor for association studies with rare variants [24].

## Availability and requirements

Linux and windows versions of the software are available at: www.ieg.uu.se/Jakobsson/software/HaploPOP/.

## Declarations

### Acknowledgements

This work was supported by a grant from the Swedish Foundation for International Cooperation in Research and Higher Education (STINT) awarded to Mattias Jakobsson and Michael Blum. A grant from the French national research agency provided support to Michael Blum and Nicolas Duforet-Frebourg (DATGEN project, ANR-2010-JCJC-1607-01) and a grant from the Swedish research council provided support to Mattias Jakobsson. The POPRES data were obtained from dbGaP (accession number phs000145.v1.p1).

## Authors’ Affiliations

## References

- Manel S, Gaggiotti OE, Waples RS. Assignment methods: matching biological questions with appropriate techniques. Trends Ecol Evol. 2005; 20(3):136–42.View ArticlePubMedGoogle Scholar
- Kayser M, de Knijff P. Improving human forensics through advances in genetics, genomics and molecular biology. Nat Rev Genet. 2011; 12(3):179–92.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81(3):559–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003; 164(4):1567–87.PubMedPubMed CentralGoogle Scholar
- Baran Y, Quintela I, Carracedo Á, Pasaniuc B, Halperin E. Enhanced localization of genetic samples through linkage-disequilibrium correction. Am J Hum Genet. 2013; 92(6):882–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Duforet-Frebourg N, Bazin E, Blum MGB. Genome scans for detecting footprints of local adaptation using a Bayesian factor model. Mol Biol Evol. 2014; 31(9):2483–2495.View ArticlePubMedPubMed CentralGoogle Scholar
- Intarapanich A, Shaw PJ, Assawamakin A, Wangkumhang P, Ngamphiw C, Chaichoompu K, et al. Iterative pruning PCA improves resolution of highly structured populations. BMC bioinformatics. 2009; 10:382.View ArticlePubMedPubMed CentralGoogle Scholar
- Lawson DJ, Hellenthal G, Myers S, Falush D. Inference of population structure using dense haplotype data. PLoS Genet. 2012; 8:e1002453.View ArticlePubMedPubMed CentralGoogle Scholar
- Latch EK, Dharmarajan G, Glaubitz JC, Rhodes Jr OE. Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation. Conserv Genet. 2006; 7(2):295–302.View ArticleGoogle Scholar
- Gattepaille LM, Jakobsson M. Combining markers into haplotypes can improve population structure inference. Genetics. 2012; 190:159–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosenberg NA, Li LM, Ward R, Pritchard JK. Informativeness of genetic markers for inference of ancestry. Am J Hum Genet. 2003; 73(6):1402–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS Genet. 2006; 2(12):e190.View ArticlePubMedPubMed CentralGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19(9):1655–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Hastie T, Friedman J, Tibshirani R, Vol. 2. The elements of statistical learning. Berlin: Springer; 2009.View ArticleGoogle Scholar
- Nelson M, Bryc K, King K, Indap A. The population reference sample, POPRES: a resource for population, disease, and pharmacological genetics research. Am J Hum Genet. 2008; 83:347–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Rota GC. The number of partitions of a set. Am Math Mon. 1964; 71(5):498–504.View ArticleGoogle Scholar
- Hudson RR. Generating samples under a Wright–Fisher neutral model of genetic variation. Bioinformatics. 2002; 18(2):337–8.View ArticlePubMedGoogle Scholar
- Scheet P, Stephens M. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006; 78(4):629–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Novembre J, Johnson T, Bryc K, Kutalik Z, Boyko AR, Auton A, et al. Genes mirror geography within Europe. Nature. 2008; 456(7218):98–101.View ArticlePubMedPubMed CentralGoogle Scholar
- Engelhardt BE, Stephens M. Analysis of population structure: a unifying framework and novel methods based on sparse factor analysis. PLoS Genet. 2010; 6(9):e1001117.View ArticlePubMedPubMed CentralGoogle Scholar
- Lao O, Liu F, Wollstein A, Kayser M. GAGA: A new algorithm for genomic inference of geographic ancestry reveals fine level population substructure in Europeans. PLoS Comput Biol. 2014; 10(2):e1003480.View ArticlePubMedPubMed CentralGoogle Scholar
- Frichot E, Mathieu F, Trouillon T, Bouchard G, François O. Fast and efficient estimation of individual ancestry coefficients. Genetics. 2014; 196(4):973–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155(2):945–59.PubMedPubMed CentralGoogle Scholar
- Mathieson I, McVean G. Differential confounding of rare and common variants in spatially structured populations. Nat Genet. 2012; 44(3):243–6.View ArticlePubMedPubMed CentralGoogle Scholar

## Copyright

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.