 Methodology article
 Open Access
 Published:
Ancestry inference using principal component analysis and spatial analysis: a distancebased analysis to account for population substructure
BMC Genomics volume 18, Article number: 789 (2017)
Abstract
Background
Accurate inference of genetic ancestry is of fundamental interest to many biomedical, forensic, and anthropological research areas. Genetic ancestry memberships may relate to genetic disease risks. In a genome association study, failing to account for differences in genetic ancestry between cases and controls may also lead to falsepositive results. Although a number of strategies for inferring and taking into account the confounding effects of genetic ancestry are available, applying them to large studies (tens thousands samples) is challenging. The goal of this study is to develop an approach for inferring genetic ancestry of samples with unknown ancestry among closely related populations and to provide accurate estimates of ancestry for application to largescale studies.
Methods
In this study we developed a novel distancebased approach, Ancestry Inference using Principal component analysis and Spatial analysis (AIPS) that incorporates an Inverse Distance Weighted (IDW) interpolation method from spatial analysis to assign individuals to population memberships.
Results
We demonstrate the benefits of AIPS in analyzing population substructure, specifically related to the four most commonly used tools EIGENSTRAT, STRUCTURE, fastSTRUCTURE, and ADMIXTURE using genotype data from various intraEuropean panels and EuropeanAmericans. While the aforementioned commonly used tools performed poorly in inferring ancestry from a large number of subpopulations, AIPS accurately distinguished variations between and within subpopulations.
Conclusions
Our results show that AIPS can be applied to largescale data sets to discriminate the modest variability among intracontinental populations as well as for characterizing intercontinental variation. The method we developed will protect against spurious associations when mapping the genetic basis of a disease. Our approach is more accurate and computationally efficient method for inferring genetic ancestry in the largescale genetic studies.
Background
During the last decade, genomewide association studies (GWAS) have helped identify a large number of allelic variants for common complex traits and diseases. Because many of the associations from these studies show small to modest effects in nature with a very strict alphalevel of statistical significance, robust conclusions from them require careful analysis to exclude falsepositive results. Population stratification, the presence of systematic allele frequency differences between populations or subpopulations, can cause spurious associations and distortions in effect estimates between genetic variants and disease [1,2,3,4,5]. Closely related individuals may have a more similar disease risk than distantly related individuals. This risk homogeneity among individuals of similar ancestries may result from lifestyle similarities or the presence of one or more riskconferring alleles [5]. However, several alleles may differ between ancestry groups that do not confer risk. Thus, some level of correlation with shared ancestry in GWAS can introduce bias leading to excess falsepositives unless a proper correction of population stratification is performed [2]. To detect whether there is confounding due to population stratification, genomic control and structured association applications are used. Several publications have described the selection of ancestry informative markers (AIMs), used to infer genetic ancestry [4, 6,7,8,9,10,11,12,13]. Basing analysis on AIMs rather than all markers that might have been analyzed in a GWAS allows a more parsimonious use of the data and the markers are typically selected to avoid strong linkage disequilibrium among the markers.
There are two commonly used types of analytical approaches to describe genetic similarities: distancebased and modelbased approaches. The distancebased approach adopts a pairwise distance matrix computed among each pair of individuals and the modelbased approach uses parametric models such as maximumlikelihood or Bayesian methods.
Menozzi et al. constructed synthetic maps of human gene frequencies in Europeans using genetic distance among population pairs [14]. They used principal component analysis (PCA) to generate a single geographic map from individual allele frequencies. The most commonly used software packages for accurately analyzing admixture population structures are EIGENSTRAT [15, 16], STRUCTURE [17] and fastStructure [18]. Price et al. developed EIGENSTRAT to detect and correct for population stratification using principal component analysis (PCA) of genotyped data to extract linear combinations of individuals that share the greatest similarities. EIGENSTRAT calculates the pattern of individual similarity in relation to markers. In the case of data with very large numbers of individuals in relation to markers, it is computationally demanding to compute the eigenvectors. Also, this does not provide any inference of ancestry membership. Pritchard introduced STRUCTURE, a Bayesian modelbased clustering method, to estimate population structure and assign individuals into population membership groups based on their genotypes under the assumption that the marker loci are unlinked and at linkage equilibrium with one another within populations [17]. With STRUCTURE, a variational Bayesian inference method was applied to compute approximate ancestry inference using the logmarginal likelihood of the data by proposing a family of tractable parametric posterior distributions over the hidden variables in the model. Inferring population structures in larger data sets with this method is computationally challenging because it requires intensive computation time and resources and may have convergence problems in fitting Markov Chain Monte Carlo based posterior samplings. In 2014, Raj proposed fastSTRUCTURE to reduce the computational time and complexity while attempting to achieve accuracy comparable to STRUCTURE [18]. ADMIXTURE is an additional popular program and uses a likelihoodbased approach [19, 20].
A distancebased approach such as multidimensional scaling could also be applied, but the groups identified from evaluating a pairwise distance (similarity) matrix may be heavily dependent on both the distance measure and the graphical representation. A challenge in largescale genetic studies is to understand the underlying data structure so as to identify whether individuals are from a homogeneous population or from heterogeneous subpopulations. When samples become larger and detected effects of genetic loci on disease phenotype become smaller, confounding with ancestry may introduce a greater number of falsepositive results. Guan et al. proposed a genetic similarity score matching method (GSM) to correct population stratification using individualbased matching [21]. GSM matches casecontrol subjects based on the average proportion of alleles using identitybystate (IBS) measures that indicate the degree of similarity over tens of thousands of SNPs. A different approach was taken by Lee et al. who developed a variation of genetic matching (GEM) called SpectralGEM that replaces the PCA used in GEM with significant ancestry components derived from the spectral graph theory [22].
More recently, Li et al. introduced an algorithmic approach, FastPop to infer the ancestry membership for the intercontinental study [23]. It is a distancebased method that reflects the clines of intermarriage among continental groups using a triangle connecting the known ancestry centroids. It could be easily applied to three or four intercontinental origins using triangle or tetrahedron shapes, respectively.
In this study we introduce a novel distancebased inference of ancestry membership with commonly used ancestry informative markers (AIMs). This novel approach can accurately infer ancestry memberships from a pairwise distance matrix calculated between individuals and centroids of the known populations using HapMap or Human Genome Diversity Project (HGDP) samples. The main aim of this method is to identify the unrevealed substructures and to infer the correct inference of ancestry memberships for samples with unknown ethnicity.
Methods
Principal component analysis
Principal component analysis (PCA) is one of the most useful statistical tools for analyzing multivariate data and has been widely applied to highdimensional genetics or genomics data. PCA uses spectral (eigenvalue) decomposition to transform a number of correlated variables into a smaller number of uncorrelated variables, which are called principal components (PCs) with a minimum loss of information. The reduced numbers of top ranked PCs are calculated by projecting samples onto spaces spanned by the eigenvectors of the sample covariance matrix and selecting the eigenvectors that comprise the largest contribution of sample variation [24]. To perform PCA, there are two approaches using eigenvalue decomposition (Pmode) and singular value decomposition (Qmode). The eigenvalue decomposition method uses the covariance relationships between markers and the singular value decomposition method uses covariance among individuals.
Initially genome wide association studies have a larger number of SNPs (p) compared to the size of samples (n), in which case principal components analysis is performed in the Qmode and can be obtained by calculating the eigenvectors and eigenvalues of a covariance matrix whose rank is at most n 1. The axes of the eigenvectors with the largest eigenvalues are usually important in describing withincontinent genetic variations and can correct for the confounding effects of population substructure. The eigenvectors so derived provide insights into variability among individuals but are specific to the specific population studied and cannot be applied to future populations. As N increases, the computational burden of computing the intraindividual correlation matrix increases exponentially. EIGENSTRAT was developed for analysis when the number of samples is far less than the number of markers, but more recent studies such as the Oncoarray [25] and the UK Biobank [26] present scenarios with very large sample sizes for which Qmode analysis is not practical.
Eigenvectors between markers derived by Pmode in a population can be used as the SNPweights (loadings) that enable researchers in a large consortium to compute the new variance components (scores) in new data with the nature of the similarity in the markers. The SNPweights so derived can help reduce the time to compute principal components by omitting the computational step of deriving the correlation matrix in markers with a specified AIMs. Because only a selected set of markers are informative about the population substructure, the number of markers that need to be included in a Pmode analysis can be limited to under ~25,000, which is computationally feasible for deriving correlation structures (Additional file 1: Supplementary Methods).
Spatial analysis; inverse distance weighted interpolation approach
Spatial analysis is used to manipulate spatial information to extract distance relationship information. Spatial interpolation is the application of spatial analysis to estimate values at unknown points with known values. As a common example, to predict precipitation in a certain area when not given entire weather information, spatial interpolation enables one to estimate precipitation in locations without recorded data using known weather information. In the Inverse Distance Weighted (IDW) interpolation method, the sample points are weighted during interpolation such that the influence of one point relative to another decreases with distance from the unknown point.
Ancestry inference using PCA scores and spatial interpolation
IDW interpolation assumes that points that are close to one another are more alike than those that are farther apart. To infer ancestry membership proportion for an unknown sample, IDW computes the distance metrics from each centroid of each known population. Those estimated values closest to the centroid of a known population will be assigned a higher proportion of ancestry that diminishes with distance and will be weighted greater than those populations that are farther away. To identify centroids of known populations, we obtained data from samples that had known European ancestry, as further described in the results.
For admixture membership, we sorted all distances among individuals by each population centroid, chose the number of admixtures denoted by s, and then found the first s closest population centroids to each individual. We then computed the distances from the centroid of the closest population to the centroids of other nearby populations. Next, we compared the distance from the second closest population centroid to each individual in relation to the distance between the two closest population centroids to each other. If the distance between the two closest population centroids was longer than the distance between the individual and the second closest population centroid, the second closest population was considered in admixture model, and so on as shown in Fig. 1.
Inverse distance weighted (IDW) interpolation in spatial analysis was then used to infer individual genetic ancestry. We applied two different spatial weights based on the centroid distances: powerdistance (PD) weights and exponentialdistance (ED) weights. Formally, let x _{ ik } be the kth score for the ith individual, x _{ jk } the kth centroid in the jth subpopulation. The powerdistance weights function, \( {\Delta}_{{\boldsymbol{p}}_{\boldsymbol{ijk}}}^{\boldsymbol{PD}} \) is a negative power function of distance given by,
The exponentialdistance weights function, \( {\Delta }_{{\boldsymbol{p}}_{\boldsymbol{ijk}}}^{\boldsymbol{ED}} \) with the negative exponential function is given by,
We developed a novel approach where eigenvalues contribute additional weights. The size of eigenvalues reflects the proportion of total variance explained by the eigenvector and larger eigenvalues should be upweighted to allow for greater variance compared to smaller eigenvalues. Formally, an ancestry inference, \( {\Delta}_{{\boldsymbol{p}}_{\boldsymbol{ijk}}}^{\boldsymbol{EVD}} \) is computed and normalized by the inverse distance weighted on each eigenvalue:
where EV _{ k } is the eigenvalue of kth score and EVD _{ ijk } is the weighted distance from each centroid of the known subpopulation to an individual. When we add eigenvalues as weights for inferring ancestry origin, the larger eigenvalue that has more weight can reveal which cluster may be closer and more appropriate to each individual (Additional file 2: Figure S1).
To utilize this methodology, we created the Rpackage AIPS that allows one to calculate SNP weights and scores from PCA, predict scores from SNP weights computed on the same predefined AIMs and infer genetic ancestry using predefined ancestry clustering information. AIPS can be performed on samples larger than markers and vice versa. After generating a matrix of SNP weights from large enough samples of AIMs, AIPS predicts a score matrix projected from the largest variance components. For missing genotype values, it computes the mean SNP value and replaces a missing genotype value with the mean SNP value. The eigenvectors and eigenvalues were calculated from correlation matrix based on standardizing each SNP column with zero mean and unit standard deviation.
Results
Application in European subpopulations and European AIMs
To demonstrate the application of AIPS, we performed an intraEuropean analysis involving 4376 individuals of European descent with a set of 25,732 preselected known IntraEuropean AIMs. For European genetic substructure studies presented in Fig. 2 (a) and Additional file 2: Table S1, we used data from the Human Genome Diversity Panel (HGDP), HapMap, Italian, Spanish, Swedish, and European Americans along with subpopulation unknown individuals from the New York Cancer Project and the Children’s Hospital of Philadelphia from the Illuminacontrol database (IControlDB). The approach to selecting subpopulations for characterizing European ancestry has been previously presented [13]. Of the 4376 individuals, 3424 participants from the New York Cancer Project and the Children’s Hospital of Philadelphia were selfidentified as Europeans and had unknown subpopulation ancestry among intraEuropean and closely related population clusters, while 952 individuals from 22 ancestryknown subpopulations (Adygei, Ashkenazi Jewish American, Basque, Bedouin, Druze, Palestine, CEPH European American, Eastern European American, German American, Greek American, Hungarian American, Irish, Italian American, Tuscan, Netherland, Orcadian, Russian, Sardinian, Scandinavian, Swedish, Spanish, United Kingdom American) were chosen to compute centroids from each European subpopulation as the known ancestry clusters [20]. In addition, we also collapsed 22 subpopulations into 7 major ethnic groups of Europeans, Arab, and Jews based on geographical and genetic differences (Fig. 2 b). To clarify the genetic heterogeneity among 7 major ethnic groups consisting of Northern European, Southern European, Great Britain, Russian, Basque, Arab, and Jew, we performed Hotelling’s T^{2} tests among 7 different centroids of subpopulations, as presented in Table 1.
We compared eigenvalues and principal components between AIPS and EIGENSTRAT. The correlation rates of eigenvalues and eigenvectors between AIPS and EIGENSTRAT are close to 1 even though the PC scales between the two approaches, scores projected from SNP weights and eigenvectors as the pattern of individual dissimilarity are different (Additional file 2: Figure S2). For the initial assessment, 952 individuals in either 22 or 7 collapsed Eurosubpopulations were selected to compare the proportions of ancestry population memberships among AIPS, STRUCTURE and fastSTRUCTURE. These samples and AIMs have been analyzed in many population studies [1, 4, 19, 20, 27]. Since the ancestry memberships for 952 samples are known, it is easy to identify the ethnic agreement between each individual and subpopulation cluster. STRUCTURE using Bayesian methods to differentiate population structures is feasible for limited sample sizes and small marker numbers [28]. Among 22 European and near Eastern subpopulations, AIPS performed better in inferring the proportions of ancestry memberships under the option in which each individual can be a descendent of between 2 and 4 populations (Additional file 2: Figure S3), while STRUCTURE with and without population labels and fastSTRUCTURE could not elucidate the heterogeneity and admixture among many of the populations in 22 clusters presented in Additional file 2: Table S1 (Additional file 2: Figure S4). In fastSTRUCTURE, we applied two types of priors; simple and logistic. fastSTRUCTURE could not recognize the differences between individuals of the 22 subpopulations. By default the number of eigenvalues in AIPS is five that are significant from the plot of the eigenvalues. AIPS allows one to have admixtures from up to number of populations. For 7 subpopulation study, we can assume at most 7 admixtures and AIPS computes 3 admixtures by default. AIPS using the top five ranked eigenvalues and the different number of admixtures displayed distinguishable population structures for inferring ethnic memberships whereas fastSTRUCTURE was unable to discriminate the ethnic heterogeneity among different population clusters.
We computed the pairwise difference of distances between centroids in two populations and ranked them based on the closeness among 22 European subpopulations (Additional file 2: Table S2 and S3). The ranks based on the geogenetic distances between them provide the clear interpretation between the geographical and population structures. We reduced the number of subpopulations using distancebased analysis and geographical relatedness. After grouping 22 subpopulations into geographically closer clusters based on PCA and distancebased analysis, we assigned 7 subpopulations: Northern European, Southern European, Great Britain, Russian, Arab, Basque, and Jews. To check whether 7 out of 22 subpopulations are substantially distinct from each other, we performed Hotelling’s T^{2} test, which compares the difference in two multivariate means. The reassigned clusters were clearly distinguishable in terms of genetic and geographical differences.
To assess the accuracy level of inferring ethnic membership in a large number of populations, we reanalyzed 952 individuals with identified ethnicities using AIPS, STRUCTURE, fastSTRUCTURE, and ADMIXTURE among 7 collapsed subpopulations. As shown in Fig. 3 AIPS in the different number of admixtures represented clear discrimination in ancestry memberships among 7 different clusters while STRUCTURE and fastSTRUCTURE performed very poorly in calculating these ancestry inferences and distinguishing all seven clusters. The graphical figure from STRUCTURE without prespecified Population IDs seemed to find largely homogeneous population structures; {Great Britain, Russian, Basque, and Northern European}; {Jews and Southern European}; and {Arab}, as further described in Table 2. To quantitate the average of accuracy in assigned groups, the average proportion for correct inference of each assigned ethnic group is computed using
where n _{ p } is a number of individuals in each population group, p a number of population group, and Avg % Correct _{ i } the correct population proportion for each individual. The average of Avg%Correct using AIPS is 0.81 among 7 population groups. As presented in Table 2, AIPS has correct classification rates between 0.68–0.90 for inferring the correct ancestry memberships whereas STRUCTURE without population information identifies correct classification between 0.21–0.64 and ADMIXTURE identifies correct classification between 0.37–0.92 among 7 European and closely related subpopulation clusters. It is not easy to identify and match true ancestry clusters in results from STRUCTURE when there are no distinct patterns between similarity and dissimilarity. Although STRUCTURE with prior population assignment identities has higher correct classification rates in the range of 0.46–1.00 than STRUCTURE without prior assignment in 0.21–0.64, STRUCTURE with prior assignment assigns only three major clusters; each cluster consists of 4, 2, and 1 subpopulations, respectively; the first cluster includes Northern European (NE, 0.93), Great Britain (GB, 1.00), Russian (1.00), and Basque (0.64); the second one includes Southern European (SE, 0.46) and Jew (1.00); the last one includes Arab (0.89). Europeans are commonly considered as a largely homogeneous population by STRUCTURE. AIPS can detect the distinction among NE, GB, Russia, and Basque while STRUCTURE is unable to distinguish among them. ADMIXTURE using supervised learning mode requires an additional file, specifying the ancestries of the reference (known) individuals. ADMIXTURE assigns 100% ancestry membership to all reference samples without further computation. We performed further comparison including 3424 Europeans with unknown subpopulation information between AIPS and ADMIXTURE. AIPS enables one to identify the ethnic heterogeneity whereas ADMIXTURE cannot recognize genetic dissimilarity between SE and Basque. According to Fig. 2 (b), there are a very small number of unknown samples of apparent Arab descent (in pink); AIPS assigned 2% out of 3424 samples into Arab subpopulation, while ADMIXTURE with and without reference information inferred 12 and 7% as Arab, respectively. In addition, we compared the average of Avg%Correct with AIPS and ADMIXTURE among 22 European subpopulations. AIPS assuming 3 admixtures has correct classification rates within 0.18–0.89 for inferring the correct ancestry memberships whereas ADMIXTURE without reference information identifies within 0.09–0.56 among 22 European subpopulations (Additional file 2: Table S4). AIPS identifies the genetic heterogeneity among 20 populations except CEU (0.34) and Italian (0.18) while ADMIXTURE clusters 22 subpopulations into 10 subpopulations, presenting no distinction in small genetic differences. Therefore, the assignment to subpopulations by AIPS outperforms the commonly used approaches, STRUCTURE and ADMIXTURE with or without prior (reference) subpopulation information.
Discussion
Population stratification in genomewide association studies can result in many falsepositive discoveries and mask the true associations [21]. Sometimes, genetic ancestry may not be available to the researchers and even though available, it may not be accurate for the underlying population genetic structure from selfreported questionnaire. It is important to confirm if selfreported ethnicity is correct and to infer the correct genetic ancestry of uncategorized individuals in many scientific studies.
The most common tool for accounting for the confounding effects of population stratification is principal component analysis (PCA). When the sample size is small, applying PCA is simple. However, because genomic highthroughput technologies are advancing, we now have larger data sets that are more difficult to analyze, especially related to inferring genetic ancestry. The widely used tool for detecting and adjusting population stratification is EIGENSOFT including two features; EIGENSTRAT and smartpca. The downside to EIGENSOFT is unable to provide correct ancestral origins while AIPS enables one to predict ancestry memberships with PCA scores as an input. The scores from PCA explain the similar patterns between samples and the eigenvectors called SNP weights (loadings) similarity between variables. Thus, the PCA scores can be used to adjust for population structures and identify ethnic origins in GWAS.
There are two types of ancestry inference approaches; distancebased and modelbased approaches. STRUCTURE and fastSTRUCTURE are the typical example of modelbased approach. Modelbased approach adapts parametric model; Bayesian or maximum likelihood method. For example, STRUCTURE uses the characteristic set of allele frequencies, HardyWeinberg equilibrium and complete linkage equilibrium between loci within populations to compute the ancestry inference in MCMC algorithm. Alternative approaches based on distance (similarity) matrix are GSM, SpectralGEM, and FastPop. GSM and SpectralGEM calculated the similarity matrix based on IBS measures and distance between two subjects that require computational intensity when the sample size is very large. FastPop results in complex computation and has not been established when inferring genetic ancestry among more than 4 population substructures. AIPS is a distancebased approach and very straightforward to infer ancestry origins. It combines two widely used statistical methods that are principal component analysis and spatial analysis. First, we compute scores of individuals and the centroid of each population in PCA and manipulate spatial information to extract distance relationship information in spatial analysis. The simplest spatial interpolation method, the inverse distance weighted interpolation is applied. This reveals the closeness between each centroid and score of individual. The calculation is very simple and straightforward and consequently the computational speed is faster. AIPS is a similar method compared to other existing population inference tools for estimating global ancestry membership like fastSTRUCTURE, because the eigenvectors from the covariance matrix are maximum likelihood estimators [24, 29]. Nevertheless, AIPS is comparably faster and achieves more accurate validation. For 952 samples using 25,732 ancestry informative markers, AIPS finished the ancestry inferences in less than 5 min to get principal components and less than 1 min to infer ancestry memberships compared with 19–23 h required by STRUCTURE, 3–4 h by fastSTRUCTURE for both 7 and 22 subpopulations, and about 20 min for 7 subpopulations and 6.5 h for 22 subpopulations by ADMIXTURE. For 4376 samples, AIPS took about 1 h to compute principal components and less than 2 min to make ancestry inferences for both 7 and 22 subpopulations while ADMIXTURE required about 28 min with reference information and about 5.5 h without reference information for 7 populations and 136.35 h without reference information for 22 populations. Furthermore, the heuristic ranks to closeness among each centroid of subpopulation provide a reasonably geogenetic relationship map to assign the given large subpopulations into the smaller clusters.
In this paper, we provide a distinct and reasonable population inference framework that achieves better accuracy comparable to STRUCTURE and fastSTRUCTURE with faster computational speed. While STRUCTURE and fastSTRUCTURE take quite long time to infer individual’s ancestry membership, AIPS takes about an hour to calculate the distance metrics of substructures for ancestry inference among 4376 individuals on 25,732 AIMs. In addition, AIPS allows one to choose the number of admixtures and top ranked eigenvalues that reflect the proportion of total variance explained by the eigenvectors. Plotting eigenvalues indicates how many top ranked eigenvalues should be included in the analysis.
If consortiums generate the large number of samples and would like to perform consistent approach, computing and sharing SNP weights (loading) consisting of the similarity in the markers (SNPs) on the specific set of AIMs are recommended. SNP weights on specific AIMs enable to predict the new variance components (scores) in new data that improves the computational efficiency and provide the consistent approach to perform multiple independent analyses in the large consortia. We recommend that the number of samples should be greater than one of markers due to shrinkage issue. In the case of analyzing genotyped data generated from same platform, AIPS can predict scores of new samples projected from SNP weights, which are eigenvectors, on the same predefined AIMs. This is an efficient computational framework to account for the confounding effects of population stratification and infer individual genetic ancestry in large consortiums.
For illustration, we selected population substructures in Europe. Europeans including EuropeanAmericans are considered as a single ethnic group such as “White” or “Caucasian” in many surveys [30]. In reality, Europeans have historically diverse ancestry and their genetic structure is strongly correlated with their geographical location [31]^{.} We demonstrated intraEuropean analysis involving 4376 individuals on 25,732 intraEuropean AIMs. Among them, 952 samples represented 22 ancestryknown subpopulations. We presented the comparisons among AIPS, STRUCTURE, and fastSTRUCTURE in graphical displays. In addition, we reduced the number of subpopulations to check the accuracy of ancestry classification. The reduced 7 clusters from 22 subpopulations within Europe are clearly distinct as suggested by Hotelling’s T^{2} test. We evaluated them with average of correctly inferred proportions. AIPS improves the level of accuracy for inferring ancestry memberships.
Better implementation of AIPS benefits from the choice of publicly available subpopulations. A pairwise distance matrix obtained between each subject and centroids of the known population substructures provides more accurate and clearer interpretation of the underlying substructures.
Conclusions
Genomewide association studies in the highdensity singlenucleotide polymorphism genotyping data have identified thousands of common variants associated to complex disease risks and traits. Because the frequency difference in genetic population structure between cases and controls due to systematic ancestry difference can lead to falsepositive results, an accurate inference of genetic ethnic membership is extremely important in many biomedical research areas. Although a few applications for detecting stratification and estimating genetic ancestry in population genetics have been developed, applying them to large genetic studies is challenging in computational time and cost. Analyzing large genotyped samples, which are becoming increasingly available, with selfreported or unknown ancestry labels, AIPS can improve accuracy in estimating ancestry memberships as well as computation efficiency. The Rpackage AIPS will be available for downloading at https://morgan.dartmouth.edu/~f000q4v/html/aips.html.
Abbreviations
 AIMs:

Ancestry Informative Markers
 AIPS:

Ancestry Inference using Principal component analysis and Spatial analysis
 ED:

ExponentialDistance
 GB:

Great Britain
 GEM:

Genetic Matching
 GSM:

Genetic Similarity Score Matching
 GWAS:

GenomeWide Association Studies
 HGDP:

Human Genome Diversity Project
 IBS:

IdentityByState
 IDW:

Inverse Distance Weighted
 NE:

Northern European
 PCA:

Principal Component Analysis
 PD:

PowerDistance
 SE:

Southern European
References
 1.
Amos CI, Wang LE, Lee J, et al. Genomewide association study identifies novel loci predisposing to cutaneous melanoma. Hum Mol Gen. 2011;20:5012–23.
 2.
Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999;55:997–1004.
 3.
Price AL, Zaitlen NA, Reich D, Patterson N. New approaches to population stratification in genomewide association studies. Nature Rev Genet. 2010;11:459–63.
 4.
Tian C, Gregersen PK, Seldin MF. Accounting for ancestry: population substructure and genomewide association studies. Hum Mol Gen. 2008;17:143–50.
 5.
Wacholder S, Rothman N, Caporaso N. Counterpoint: bias from population stratification is not a major threat to the validity of conclusions from epidemiological studies of common polymorphisms and cancer. Cancer Epidemiol Biomarkers Prevent. 2002;11:513–20.
 6.
Amirisetty S, Hershey GKK, Baye TM. AncestrySNPminer: a bioinformatics tool to retrieve and develop ancestry informative SNP panels. Genomics. 2012;100:57–63.
 7.
Baye TM, Tiwari HK, Allison DB, Go RC. Database mining for selection of SNP markers useful in admixture mapping. BMC BioData Mining. 2009; doi:10.1186/1756–0381–21.
 8.
Kodaman N, Aldrich MC, Smith JR, et al. A small number of candidate gene SNPs reveal continental ancestry in African Americans. Ann Hum Genet. 2013;77:56–66.
 9.
Kosoy R, Nassir R, Tian C, et al. Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America. Hum Mutat. 2009;30:69–78.
 10.
Ma J, Amos CI. Principal components analysis of population admixture. PLoS One. 2012;7(7):e40115. doi:10.1371/journal.pone.0040115.
 11.
PardoSeco J, MartinonTorres F, Salas A. Evaluating the accuracy of AIM panels at quantifying genome ancestry. BMC Genomics. 2014;15:543.
 12.
Tian C, Plenge RM, Ransom M, et al. Analysis and application of European genetic substructure using 300 K SNP information. PLoS Genet. 2008; doi: 10.1371/journal.pgen.0040004.
 13.
Tian C, Kosoy R, Nassir R, et al. European population genetic substructure: further definition of ancestry informative markers for distinguishing among diverse European ethnic groups. Mol Med. 2009;15:371–83.
 14.
Menozzi P, Piazza A, CavalliSforza L. Synthetic maps of human gene frequencies in Europeans. Science. 1978;201:786–92.
 15.
Patterson N, Price AL, Reich D. Population Structure and Eigenanalysis. PLoS Genetics. 2006; doi: 10.1371/journal.pgen.0020190.
 16.
Price AL, Patterson NJ, Plenge RM, et al. Principal components analysis corrects for stratification in genomewide association studies. Nature Genet. 2006;38:904–9.
 17.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
 18.
Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: Variational inference of Populatioin structure in large SNP data sets. Genetics. 2014;197:573–89.
 19.
Alexander DH, Novembre J, Lange K. Fast modelbased estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.
 20.
Alexander DH, Lange K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinformatics. 2011;12:246.
 21.
Guan W, Liang L, Boehnke M, Abecasis GR. Genotypebased matching to correct for population stratification in largescale casecontrol genetic association studies. Genet Epidemiol. 2009;33:508–17.
 22.
Lee AB, Luca D, Klei L, et al. Discovering genetic ancestry using spectral graph theory. Genet Epidemiol. 2010;34:51–9.
 23.
Li Y, Byun J, Cai G, et al. FastPop: a rapid principal component derived method to infer intercontinental ancestry using genetic data. BMC Bioinformatics. 2016;17:122. doi: 10.1186/s1285901609651.
 24.
Lee S, Zou F, Wright FA. Convergence and prediction of principal component scores in highdimensional settings. Ann Stat. 2010;38:3605–29.
 25.
Amos CI, Dennis J, Wang Z, Byun J, et al. The OncoArray consortium: a network for understanding the genetic architecture of common cancers. Cancer Epidemiol Biomark Prev. 2017;26(1):126–35.
 26.
Sudlow C, Gallacher J, Allen N, Beral V, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.
 27.
Novembre J, Johnson T, Bryc K, et al. Genes mirror geography within Europe. Nature. 2008;456:98–101.
 28.
PorrasHurtado L, Ruiz Y, Santos C, Phillips C, et al. An overview of STRUCTURE: applications, parameter settings, and supporting software. Front Genet. 2013;4:98.
 29.
Girshick M. Principal components. JASA. 1936;31:519–28.
 30.
Price AL, Butler J, Patterson N, et al. Discerning the ancestry of European Americans in genetic association studies. PLoS Genet. 2008;4:e236. doi:10.1371/journal.pgen.0030236.
 31.
Nelis M, Esko T, Magi R, et al. Genetic structure of Europeans: a view from the northeast. PLoS One. 2009;4:e5472. doi:10.1371/journal.pone.0005472.
Funding
This work was supported by the National Institutes of Health grants from P20GM103534, P30CA0123108, and U19CA148127.
Availability of data and materials
HDGP: http://www.hagsc.org/hgdp/files.html
HapMap: ftp://ftp.ncbi.nlm.nih.gov/hapmap/samples_individuals/
dbGaP: https://www.ncbi.nlm.nih.gov/projects/gap/cgibin/study.cgi?study_id=phs000122.v1.p1
AIPS: https://morgan.dartmouth.edu/~f000q4v/html/aips.html
The curated genetic ancestry data can be downloaded from https://morgan.dartmouth.edu/~f000q4v/Public/IntraEuro/AIPS/Data/
Author information
Affiliations
Contributions
JY, YH, and CIA conceived and designed study. JY and YH developed methods. JY performed the computational analyses. JY wrote the original draft and JY, YH, IPG, JAB, MFS, and CIA reviewed and edited the manuscript. MFS provided the curated genetic ancestry data. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Supplementary Methods. Mathematical definition of principal component analysis. (DOCX 20 kb)
Additional file 2:
RSupplementary Materials. The attached file includes 4 supplementary figures and 4 supplementary tables. Figure S1. C_{N.Euro} and C_{Rus} present centroid of known ancestry samples from Northern Europeans and Russians, respectively. (a) With the first three scores from PCA, individual A seems to be closer to Russian group on the proportion of total variance explained by eigenvalues. (b) In the two dimensional plot with the top two principal component scores, individual A seems to be closer to N. European. Figure S2. Comparison of eigenvalues and top 3 principal components from AIPS and EIGENSTRAT. The options that were set in EIGENSTRAT were numoutlieriter = 0;outliermode = 2(no outlier removal) and in AIPS the option was method = eigen. Figure S3. Graphical Comparison of Population Structure using AIPS among 22 European subpopulations. Only 952 known ancestry individuals were used in 22 subpopulations within Europe. The scores from PCA were first calculated then InverseDistance Weighted Interpolation without and with eigenvalue weight were applied to infer the ancestry membership. The number of admixture indicated the definition of admixture in AIPS. Figure S4. Graphical Comparison of Population Structure using STRUCTURE and fastSTRUCTURE among 22 European subpopulations. The inferences of ancestry membership for 952 individuals were calculated by STRUCTURE and fastSTRUCTURE. (a) The prior population information was not given to compute inference of population membership using STRUCTURE. (b) To infer the population membership within 22 Europe countries, the prior population information was assigned in STRUCTURE. (c) fastSTRUCTURE was applied to infer 22 European subpopulations with simple model. (d) fastSTRUCTURE was used with logistic prior model. Table S1. Distancebased clustering among 952 known and 3426 unknown ancestry Europeans on 25,732 AIMs. Table S2. Rankbased on Closeness among 22 European subpopulations. Table S3. Distance between two centroids among 22 European subpopulations. Table S4. The Average Percent of Correctly Inferred Proportions from AIPS and ADMIXTURE without Population Information. (DOCX 922 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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
Cite this article
Byun, J., Han, Y., Gorlov, I.P. et al. Ancestry inference using principal component analysis and spatial analysis: a distancebased analysis to account for population substructure. BMC Genomics 18, 789 (2017). https://doi.org/10.1186/s1286401741668
Received:
Accepted:
Published:
Keywords
 Ancestry inference
 Principal component analysis
 Spatial analysis
 Inverse distance weighted interpolation