Evaluating the accuracy of AIM panels at quantifying genome ancestry
BMC Genomics volume 15, Article number: 543 (2014)
There is a growing interest among geneticists in developing panels of Ancestry Informative Markers (AIMs) aimed at measuring the biogeographical ancestry of individual genomes. The efficiency of these panels is commonly tested empirically by contrasting self-reported ancestry with the ancestry estimated from these panels.
Using SNP data from HapMap we carried out a simulation-based study aimed at measuring the effect of SNP coverage on the estimation of genome ancestry. For three of the main continental groups (Africans, East Asians, Europeans) ancestry was first estimated using the whole HapMap SNP database as a proxy for global genome ancestry; these estimates were subsequently compared to those obtained from pre-designed AIM panels. Panels that consider >400 AIMs capture genome ancestry reasonably well, while those containing a few dozen AIMs show a large variability in ancestry estimates. Curiously, 500-1,000 SNPs selected at random from the genome provide an unbiased estimate of genome ancestry and perform as well as any AIM panel of similar size. In simulated scenarios of population admixture, panels containing few AIMs also show important deficiencies to measure genome ancestry.
The results indicate that the ability to estimate genome ancestry is strongly dependent on the number of AIMs used, and not primarily on their individual informativeness. Caution should be taken when making individual (medical, forensic, or anthropological) inferences based on AIMs.
With the publication of the Human Genome Project (http://www.ornl.gov/sci/techresources/Human_Genome/home.shtml) in 2000 and the pioneering high-throughput single nucleotide polymorphism (SNP) genotyping projects (such as HapMap; http://www.hapmap.org) our perception of human genome has changed, as well as our understanding of human evolution and genome ancestry. The term ancestry refers to “the origin or background of something” (http://oxforddictionaries.com). Accordingly, in human genetics, ancestry is generally understood as the origin or background of our genomes. However, the question is far from trivial. Considering the way in which the DNA material is inherited through generations, most of it from both parents (the exception being the uniparental markers), entire blocks of our genome can have different ancestral origins. In the words of Svante Pääbo, “to understand what make us unique, both as individuals and as a species, we need to consider the genome as a mosaic of discrete segments, each with its own unique history and relatedness to different contemporary and ancestral individuals” .
Although genetic variation in humans shows gradients of allele frequencies extending over the entire world (within and among continents or among groups of individuals ), there is empirical evidence indicating that the most contrasting genomic patterns of diversity in humans occur at an inter-continental level; e.g. Africa, Europe, and Asia (often erroneously interpreted as genetic support for “races” ). The best way to characterize these continental patterns (as discrete clusters of variation) is by examining them to a genomic scale, given that single locus could not necessarily capture the global genomic scenario. With the arrival of new genotyping technologies and large-scale genomic projects, it is now possible to measure genomic ancestry using large genome-wide SNP panels or, more recently, next generation sequencing data (NGS; e.g. http://www.1000genomes.org) [3, 4]. However, this genomic approach is not always cost-effective and it can also represent a handicap in particular scenarios (low amount and/or degraded DNA; e.g. population and forensic routine casework). Alternatively, ancestry can be estimated using a selected number of SNPs ranging rom a few dozens to several hundreds; this option has been favored in different areas of biomedical research, including case-control association studies of complex disease (e.g. admixture mapping) [5–8], human population studies [9–11], and forensic genetics and police investigation [12–14]. The selected SNPs are commonly known as Ancestry Informative Markers (AIMs) and received this name because they exhibit large differences in allele frequencies between populations from different geographical or ethnic groups. By genotyping a number of AIMs, it seems possible to estimate the most likely geographical or ethnic origin of a given genomic profile, or to ascertain what proportion of ancestry in this profile is derived from different geographical regions or source populations.
Measuring ancestry is important in biomedical studies for a number of reasons. For instance, it has been demonstrated that population stratification represents an important confounding effect in case-control association studies of complex and multi-factorial diseases [15–19]. Estimating ancestry using AIMs panels can be used in these studies to control for population sub-structure in medical studies. Some companies have developed commercial kits (http://www.illumina.com/products/dna_test_panel.ilmn) aimed at measuring the ancestry of samples as a screening method before proceeding with their high-throughput genotyping or massive parallel sequencing.
The search for autosomal ancestry has also been a focus of attention in the forensic community [13, 20, 21]. Forensic geneticists have to deal with evidentiary samples containing little amounts of, and/or poorly preserved, DNA. In these cases, the limited amount of DNA available often allows a single PCR reaction only or, in cases where more DNA is available, it is generally preferred to preserve it in order to allow a second and independent test in a different laboratory. Forensic geneticists have also designed their own panels of AIMs allowing estimation of ancestry based on single-plex assays [13, 21].
At the same time, many private companies offer direct-to-consumer-tests (DTCT) specifically designed to measure ancestry [22, 23]. Although most of these tests do not aim to provide specific information about disease conditions, in reality they could reveal information relevant for the customer’s health. This is due to the fact that there are health disorders that can be more highly correlated with certain ancestries than others. However, the accuracy of DTCT has been questioned on several grounds. For instance, these companies often offer only to genotype the uniparental markers (the mitochondrial DNA [mtDNA] and/or the Y-chromosome ); however, these markers behave as single locus and therefore can only reflect a very tiny portion of the genomic individual ancestry [24, 25]. Ancestry inferences made using autosomal markers have been conflicting too .
Most of the AIM panels available in the literature have been designed by way of selecting SNPs from large genomic databases (e.g. HapMap) showing skewed population frequencies between the ancestral populations targeted. Usually, researchers do not evaluate the amount of genetic informativeness provided individually by the selected SNPs. An exception is the study by Galenter et al. , who used a multi-step algorithm that weighs the amount of information provided by their AIMs regarding the ancestral populations being considered.
The ability of an AIM panel to measure ancestry is generally evaluated empirically, that is, by examining its performance on a given set of DNA samples for which a given ancestry is already assumed. Several statistical techniques (such as principal component analysis [PCA], and admixture analysis) are then used to evaluate their efficiency. An AIM panel is generally considered to be efficient if e.g. it can differentiate the targeted populations in the Euclidian space represented by two or three principal components (PC) or if the inferred ancestry is consistent with some expectation (e.g. self-reported ancestry). The majority of the panels are designed with the aim of distinguishing main continental groups (e.g. Africans, Asians, Europeans, Americans) owing to the known difficulties of using small SNP panels to classify individuals when they belong to closely related populations.
The number of SNPs incorporated into AIM panels varies from a few dozens to a few hundreds (Table 1); and this number is generally constrained by the genotyping technique employed. For instance, most of the techniques allow genotyping only a few dozen SNPs in a single PCR reaction (e.g. SNaPshot , mass array spectrometry ). Although it seems reasonable to consider that the number of AIMs in a panel could be relevant when estimating ancestry, a comprehensive evaluation of this factor has not been carried out to date. Questions arise too about how many autosomal markers would be needed in scenarios of population admixture, where the use of a dense panel of AIMs could be even more important in order to better represent the admixed profiles. Difficulties in correctly estimating genome ancestry could also derive from the application of an AIM panel to a set of samples different from that used during the SNP selection process and the training tests (‘lack of portability’ ). At the same time, using mtDNA and Y-chromosome markers to measure genome ancestry could be justified in some genetic contexts , but not when trying to infer global individual genome ancestry .
The HapMap SNP database was retrieved from its repository (http://hapmap.ncbi.nlm.nih.gov). This database contains 1,440,616 SNPs genotyped in a total of 1,218 individual samples belonging to the following main continental groups: 472 Africans, 58 Americans, 101 Central-South Asians, 364 East Asians and 223 Europeans. Unless specified, for most of the simulation experiments, only three populations representing the main continental groups were taken from the full HapMap data, namely CEU (European ancestry), CHB (East Asian ancestry) and YRI (African ancestry), with 50 individuals in each group. This decision was based on the fact that most of the AIM panels available were designed to identify ancestry from main population groups.
For some simulation experiments, we created artificial scenarios of admixture by mixing at random the same proportion of SNPs from the following three HapMap datasets: CEU, CHB, and YRI. Therefore, the expected genome admixture in these artificially created hybrid genomic profiles (henceforth referred to as the “AA-genomes”) is 1/3 of ancestry from each of the main continental groups (Asia, Europe, and Africa).
Sample size and ancestry estimates
The dependence of ancestry inference on sample size was estimated through simulation experiments using a similar procedure to that in Heinz et al. . In brief, for each of the three main continental populations, we randomly selected 1,000 sub-samples of variable sizes (from five to 40 profiles; in stepwise increments of five and taken without replacement). Thus, for example, we obtained 1,000 sub-samples of size five, 1,000 sub-samples of size ten, and so on until a maximum sample size of 40. For each of the sub-samples we computed ancestry proportions as indicated below. Continental ancestry was estimated as the mean value obtained for the 1,000 sub-samples in each sample window, and bootstrapping intervals were built accordingly.
The software Admixture v. 1.22  was used to estimate individual and population ancestries. This software was run using default parameters. Cross validation errors were obtained from Admixture in order to determine the most likely K value (K indicating the number of inferred clusters showing the lowest cross validation error).
PLINK v.1.07  was used to obtain Identity-By-State (IBS) values between individuals, and IBS values were used to carry out two-dimensional PCA. PLINK was used with default settings. Only when calculating the effect of population sample size on the estimation of ancestry, individual profiles with missing data >10% were filtered out (call rates <90% could be critical when dealing with AIM panels containing low number of SNPs).
Locus specific branch length (LSBL) statistics was estimated using pairwise FST distances as carried out in Shriver et al. . LSBL aims to assist in the selection of AIMs in panels taking into account their level of individual informativeness with regards to the classification population sets. FST values were taken from SPSmart and ENGINES [34, 35].
Pre-designed AIM panels
Ancestry of the selected HapMap datasets was estimated using different AIM panels (Table 1): Corach et al.  (COR), Galenter et al.  (GAL), Halder et al.  (HAL), Kosoy et al.  (KOS), Lao et al.  (LAO), Nassir et al.  (NAS), Phillips et al.  (PHI), and the commercial DNA Test Panel from Illumina (ILU; http://www.illumina.com/products/dna_test_panel.ilmn). All of these panels were originally designed to differentiate the three main continental groups (Europe, East Asia and Sub-Saharan Africa), some of them also including other ancestral groups (see Table 1).
Estimating genome ancestry
The genome-wide set of SNPs in HapMap can be used to estimate (global) genome ancestry of continental groups. Admixture analysis shows an optimum value of K = 3 when considering CEU, CHB, and YRI. Each individual profile received virtually 100% of the expected genomic ancestry in one cluster (in agreement with their geographic origin/self-declared ancestry); that is, for example, a Yoruban profile receives a ~100% membership in a cluster that groups all samples of African ancestry (Additional file 1).A simulation experiment was carried out by selecting subsets of random SNPs (henceforth rSNPs) from the whole HapMap database and obtaining estimates of ancestry from these subsets. This procedure allows us to investigate the extent to which the estimation of genome ancestry degenerates when using decreasing amounts of SNPs. Panels of 10,000, 1,000, and 500 rSNPs where randomly selected; each panel size was sampled 500 times each in order to account for sampling variability. The simulations indicate that the inferred genome ancestry degenerates slightly as fewer rSNPs enter in a set (Figure 1). However, ancestry estimates using panels of 500 rSNPs approach well the genome ancestry (although the estimates show a moderate dispersion).
The estimates above were obtained considering the three main continental groups: Europe, Asia, and Africa. However, the number of SNPs needed to infer ancestry strongly depends on the evolutionary relatedness of the populations being considered: the closer the population under study, the larger the number of SNPs needed. The PCA plot in Additional file 2 indicates that the whole set of SNPs in HapMap clearly separates East Asian populations CHB + CHD (Chinese) from JPT (Japanese). However, using panels of 100,000, 50,000,10,000, 1,000 and 500 rSNPs, differentiating these two populations groups becomes increasingly difficult; see e.g. the overlapping patterns of profiles in the PCAs of Additional file 2 when using 500 rSNPs. In population scenarios considering very closely related groups, the whole power of a genome-wide dataset would be needed in order to differentiate populations; e.g. see the case for European populations in Novembre et al. .
Pre-designed AIM panels
Ancestry estimates were obtained using different pre-designed panels and compared to the genome ancestry inferred using sets of 10,000, 1,000 and 500 rSNPs. As shown in the bar-plots in Figure 2, inference of ancestry degenerates as fewer AIMs are considered in the panels. The PCA plots of Figure 3 mirror the same scenario. It is noticeable that the panel of 500 rSNPs yields broadly similar results to the GAL and ILL panels, that is, the panels containing a similar number of AIMs. PCA clearly shows that the panels containing fewer AIMs show more disperse patterns. The LAO panel, which considers only 10 AIMs, represents the most extreme case, showing an extremely large variability (Table 1, Figure 3).
Measuring informativeness of SNPs in AIM panels
LSBL can be used to measure the informativeness of specific AIMs in the panels and their potential to measure different ancestries when applied to the HapMap populations. To the best of our knowledge, only the GAL panel was designed using LSBL as a criterion to balance the informativeness of the AIMs incorporated in the panels (Table 2). As shown in Additional file 3, all the panels show unbalanced accumulated LSBL values, including GAL. However, both of the panels that contain the largest number of AIMs (GAL and ILL) perform reasonably well when estimating genome ancestry (see above). This evidence suggests that, even though measuring the informativeness of the AIMs would seem a logical way to proceed when designing a panel, ensuring that a large number of AIMs is considered constitutes at least an equally important parameter.
We also followed the LSBL criteria for selecting the best AIMs from the HapMap database (as done in Galanter et al. ); in each case, the number of simulated AIMs selected was the same as the number used in the different panels tested. Additional file 3 shows that these test panels work better than their pre-designed counterpart panel using analogous continental populations (TSI as representative of Europe; CHD from East Asia and LWK from Africa). However, those containing a higher number of HapMap-AIMs perform much better than those considering lower numbers of SNPs.
Effect of population sample size on the estimation of ancestry
Simulations were carried out in order to estimate the effect of population sample size when inferring genome ancestry using different AIM panels. These analyses are significant because geneticists are often interested in measuring the average genetic ancestry of a given population (see e.g. Heinz et al. ), but sample sizes are very low in several studies.As expected, simulations indicate that as the sample size increases the variability of the ancestry estimations decreases (Figure 4). The most remarkable finding is that, for equal sample sizes, the largest AIM panels, GAL and ILL, show lower variability than smaller panels.
Inferring ancestry in admixed genomes
The hybrid AA-genomes were created in order to allow some simulations to investigate the number of SNPs needed to capture genome ancestry in admixture scenarios. Figure 5 shows the patterns of ancestry estimated using the different panels and compared to the expected ancestry in these genomes (equal membership in Africa, East Asia and Europe; ~33%). The estimates of individual ancestry are stable and close to expectation when using rSNPs (>1000) and also using GAL (and slightly worse with ILL); while other panels such as PHI, COR and particularly LAO show more arbitrary patterns.
Quantifying errors in ancestry estimates
The standard deviation can be used to measure the error of the estimated ancestry obtained from panels, compared to global genome ancestry (Figure 6). The errors in ancestry estimates are more pronounced when the number of rSNPs falls to ~500 SNPs (Figure 6). The ability of panels to capture different proportions of ancestry varies significantly from panel to panel, but the panels with more AIMs (specially GAL and ILL) perform much better than those containing fewer SNPs. For instance, panels PHI and COR show very large variability in ancestry estimates, and this variability is extremely large in the case of LAO (Figure 6). Additional file 4 shows that there is a negative correlation between the number of SNPs in a panel and the error associated to the estimates of ancestry. Furthermore, the error differs regarding the kind of ancestry that is measured: for instance, the panels with the larger numbers of AIMs have more difficulties to measure the European ancestry than the African one. Curiously, the error in ancestry estimates seems to have a more balanced behavior regarding the inferences of the different ancestries in scenarios of admixture (using AA-genomes).
Measuring genome ancestry is an issue of interest in different fields of biomedical research, including case-control association studies, forensic casework and police investigation, and anthropological studies. It is also of interest for private companies, given the growing social interest in knowing more about ancestry coupled with the progressive reduction of the cost of DNA tests.
The present study aims to estimate the number of SNPs needed to reliably infer genome ancestry using unbiased sets of SNPs (rSNPs) and sets of pre-fabricated AIM panels.
The results indicate that 10,000 SNPs selected at random from an individual can be used to infer genome ancestry with negligible error when considering the three HapMap populations CEU, CHB, and YRI. Even so, panels of 500 rSNPs perform reasonably well in this population scenario. Below this number, errors in the inference of ancestry increase noticeably as the number of rSNPs is reduced. As expected, the number of rSNPs needed to infer ancestry strongly depends on the evolutionary proximity of the populations under study. For instance, we made simulations to test the number of rSNPs needed to differentiate ancestry in two different East Asian populations, Chinese and Japanese. Here the number of rSNPs needed to differentiate these populations increases significantly more than one order of magnitude; therefore, the need for searching panels of highly discriminating AIMs is more justified. The distinction between individual ancestries within Asian populations (or other closely related groups) would require genome-wide screenings  or very large panels of AIMs (probably containing thousands of SNPs).
During the last few years, several panels of SNPs have been designed in order to estimate ancestry using only a few markers (AIM panels). Analyses were carried out in the present study in order to assess the performance of these panels when applied to three main HapMap continental populations, CEU, CHB, and YRI. The results indicate that inference of ancestry can be seriously compromised when using panels containing small numbers of AIMs. For instance, out of the panels tested in the present study, those showing the best performance are ILU and GAL, that is, those that have more AIMs, while the ones including only a few dozen AIMs show higher errors and variability (Additional file 4).
It is interesting to note that neither GAL nor ILU were specifically conceived to discriminate exactly between the three tested populations from HapMap (Table 1). In fact, the cumulative LSBL value for the three HapMap populations indicates that these AIMs are not balanced for these population groups (Figure 5). Therefore, the good performance of these panels is based to a great extent on the large number of AIMs contained in these panels, and not exclusively on the individual discriminatory power of the selected SNPs.
Of the different AIM panels tested in the present study, only GAL  was initially designed using a criterion of SNP informativeness; thus, markers were selected on the basis of balanced cumulative LSBL values in the targeted populations. The present study reveals two limitations in this procedure. First, this method does not specify how many SNPs should enter the ancestry SNP panel; thus, different amounts of SNPs could fit the criteria of similar cumulative LSBL values . Second, it is hard to predict the extent to which the good LSBL characteristics of the AIMs in a panel (in training set populations) can be extrapolated to other population sample sets (which may or may not belong to a closely related geographic/ethnic group). The results of the present study indicate that the best way to ensure the good performance of a panel is to incorporate the largest possible number of AIMs (at least >400 when considering main continental groups).
Our results allow further relevant conclusions. First, inferences related to population demography (e.g. molecular anthropological studies) could be biased if using panels containing a small number of AIMs. Second, DTCT should consider employing panels containing large amounts of markers in order to provide the most accurate service to the public . Third, one of the most important limitations in forensic casework and police DNA investigation is the amount and quality of DNA available from evidentiary samples; here the use of AIM panels could play an important role given that only a limited number of SNPs can enter a single PCR reaction. However, forensic specialists and police investigators should be aware of the limitations of the approach; where possible, a large number of AIMs should be analyzed in order to provide the most precise inferences on the ancestry of evidentiary samples. Inferences of ancestry could be particularly compromised in scenarios of admixture. In such scenarios, SNP coverage can be more crucial given the need to represent the genome more densely than in scenarios of non-admixture (where only one main component has to be measured). The arrival of NGS technologies may help overcome these limitations; see however some caveats in Bandelt and Salas .
Caution should be exercised when inferring ancestry using AIM panels. The concept of ancestry is a complex one and although it can be operational for particular purposes, it can lead to erroneous perceptions of human variability. As stated by Sankar and Cho : “the appearance of clustering is a function of how populations are sampled, of how criteria for boundaries between clusters are set, and of the level of resolution used. In the same way that the earth can be described by many different kinds of maps —from topological to economic—so, too, can the naturally occurring genetic variation among populations be divided in numerous ways and be made to highlight any chosen similarity or difference”. This conclusion is particularly important for the general public, who is often not aware of the limitations of ancestry DNA tests; and also in police investigation, where over-interpretation of an ancestry test could have important consequences on the investigation of forensic DNA evidence.
Pääbo S: The mosaic that is our genome. Nature. 2003, 421 (6921): 409-412. 10.1038/nature01400.
Serre D, Pääbo S: Evidence for gradients of human genetic diversity within and among continents. Genome Res. 2004, 14 (9): 1679-1685. 10.1101/gr.2529604.
Reich D, Patterson N, Campbell D, Tandon A, Mazieres S, Ray N, Parra MV, Rojas W, Duque C, Mesa N, Mesa N, Garcia LF, Triana O, Blair S, Maestre A, Dib JC, Bravi CM, Bailliet G, Corach D, Hunemeier T, Bortolini MC, Salzano FM, Petzl-Erler ML, Acuna-Alonzo V, Aguilar-Salinas C, Canizales-Quinteros S, Tusie-Luna T, Riba L, Rodriguez-Cruz M, Lopez-Alarcon M, et al: Reconstructing Native American population history. Nature. 2012, 488 (7411): 370-374. 10.1038/nature11258.
Novembre J, Johnson T, Bryc K, Kutalik Z, Boyko AR, Auton A, Indap A, King KS, Bergmann S, Nelson MR, Stephens M, Bustamante CD: Genes mirror geography within Europe. Nature. 2008, 456 (7218): 98-101. 10.1038/nature07331.
Galanter JM, Fernández-López JC, Gignoux CR, Barnholtz-Sloan J, Fernandez-Rozadilla C, Via M, Hidalgo-Miranda A, Contreras AV, Figueroa LU, Raska P, Jimenez-Sanchez G, Silva Zolezzi I, Torres M, Ponte CR, Ruiz Y, Salas A, Nguyen E, Eng C, Borjas L, Zabala W, Barreto G, Rondon González F, Ibarra A, Taboada P, Porras L, Moreno F, Bigham A, Gutierrez G, Brutsaert T, Leon-Velarde F, et al: Development of a panel of genome-wide ancestry informative markers to study admixture throughout the Americas. PLoS Genet. 2012, 8 (3): e1002554-10.1371/journal.pgen.1002554.
Montana G, Pritchard JK: Statistical tests for admixture mapping with case-control and cases-only data. Am J Hum Genet. 2004, 75 (5): 771-789. 10.1086/425281.
Rosenberg NA, Li LM, Ward R, Pritchard JK: Informativeness of genetic markers for inference of ancestry. Am J Hum Genet. 2003, 73 (6): 1402-1422. 10.1086/380416.
Rodrigues AC, Perin PM, Purim SG, Silbiger VN, Genvigir FD, Willrich MA, Arazi SS, Luchessi AD, Hirata MH, Bernik MM, Dorea EL, Santos C, Faludi AA, Bertolami MC, Salas A, Freire A, Lareu MV, Phillips C, Porras-Hurtado L, Fondevila M, Carracedo A, Hirata RD: Pharmacogenetics of OATP transporters reveals that SLCO1B1 c.388A> G variant is determinant of increased atorvastatin response. Int J Mol Sci. 2011, 12 (9): 5815-5827.
Taboada-Echalar P, Álvarez-Iglesias V, Heinz T, Vidal-Bralo L, Gómez-Carballa A, Catelli L, Pardo-Seco J, Pastoriza A, Carracedo Á, Torres-Balanza A, Rocabado O, Vullo C, Salas A: The genetic legacy of the pre-colonial period in contemporary Bolivians. PLoS One. 2013, 8 (3): e58980-10.1371/journal.pone.0058980.
Cerezo M, Achilli A, Olivieri A, Perego UA, Gómez-Carballa A, Brisighelli F, Lancioni H, Woodward SR, López-Soto M, Carracedo Á, Capelli C, Torroni A, Salas A: Reconstructing ancient mitochondrial DNA links between Africa and Europe. Genome Res. 2012, 22 (5): 821-826. 10.1101/gr.134452.111.
Lao O, van Duijn K, Kersbergen P, de Knijff P, Kayser M: Proportioning whole-genome single-nucleotide-polymorphism diversity for the identification of geographic population structure and genetic ancestry. Am J Hum Genet. 2006, 78 (4): 680-690. 10.1086/501531.
Phillips C, Prieto L, Fondevila M, Salas A, Gómez-Tato A, Álvarez-Dios J, Alonso A, Blanco-Verea A, Brión M, Montesino M, Carracedo A, Lareu MV: Ancestry analysis in the 11-M Madrid bomb attack investigation. PLoS One. 2009, 4 (8): e6583-10.1371/journal.pone.0006583.
Phillips C, Salas A, Sánchez JJ, Fondevila M, Gómez-Tato A, Álvarez-Dios J, Calaza M, de Cal MC, Ballard D, Lareu MV, Carracedo A: Inferring ancestral origin using a single multiplex assay of ancestry-informative marker SNPs. Forensic Sci Int Genet. 2007, 1 (3–4): 273-280.
Salas A, Phillips C, Carracedo A: Ancestry vs physical traits: the search for ancestry informative markers (AIMs). Int J Legal Med. 2006, 120 (3): 188-189. 10.1007/s00414-005-0032-3.
Pritchard JK, Rosenberg NA: Use of unlinked genetic markers to detect population stratification in association studies. Am J Hum Genet. 1999, 65 (1): 220-228. 10.1086/302449.
Reich DE, Goldstein DB: Detecting association in a case-control study while correcting for population stratification. Genet Epidemiol. 2001, 20 (1): 4-16. 10.1002/1098-2272(200101)20:1<4::AID-GEPI2>3.0.CO;2-T.
Bacanu SA, Devlin B, Roeder K: Association studies for quantitative traits in structured populations. Genet Epidemiol. 2002, 22 (1): 78-93. 10.1002/gepi.1045.
Devlin B, Roeder K, Wasserman L: Genomic control, a new approach to genetic-based association studies. Theor Popul Biol. 2001, 60 (3): 155-166. 10.1006/tpbi.2001.1542.
Marchini J, Cardon LR, Phillips MS, Donnelly P: The effects of human population structure on large genetic association studies. Nat Genet. 2004, 36 (5): 512-517. 10.1038/ng1337.
Sánchez JJ, Borsting C, Balogh K, Berger B, Bogus M, Butler JM, Carracedo Á, Court DS, Dixon LA, Filipovic B, Fondevila M, Gill P, Harrison CD, Hohoff C, Huel R, Ludes B, Parson W, Parsons TJ, Petkovski E, Phillips C, Schmitter H, Schneider PM, Vallone PM, Morling N: Forensic typing of autosomal SNPs with a 29 SNP-multiplex–results of a collaborative EDNAP exercise. Forensic Sci Int Genet. 2008, 2 (3): 176-183. 10.1016/j.fsigen.2007.12.002.
Sánchez JJ, Phillips C, Børsting C, Balogh K, Bogus M, Fondevila M, Harrison CD, Musgrave-Brown E, Salas A, Syndercombe-Court D, Schneider PM, Carracedo A: A multiplex assay with 52 single nucleotide polymorphisms for human identification. Electrophoresis. 2006, 27: 13-24.
Goldsmith L, Jackson L, O’Connor A, Skirton H: Direct-to-consumer genomic testing: systematic review of the literature on user perspectives. Eur J Hum Genet. 2012, 20 (8): 811-816. 10.1038/ejhg.2012.18.
Bandelt H-J, Yao Y-G, Richards MB, Salas A: The brave new era of human genetic testing. Bioessays. 2008, 30 (11–12): 1246-1251.
Egeland T, Bøvelstad HM, Storvik GO, Salas A: Inferring the most likely geographical origin of mtDNA sequence profiles. Ann Hum Genet. 2004, 68 (Pt 5): 461-471.
Egeland T, Salas A: Estimating haplotype frequency and coverage of databases. PLoS One. 2008, 3 (12): e3988-10.1371/journal.pone.0003988.
Ng PC, Murray SS, Levy S, Venter JC: An agenda for personalized medicine. Nature. 2009, 461 (7265): 724-726. 10.1038/461724a.
Cerezo M, Černý V, Carracedo Á, Salas A: Applications of MALDI-TOF MS to large-scale human mtDNA population-based studies. Electrophoresis. 2009, 30 (21): 3665-3673. 10.1002/elps.200900294.
Ribas G, González-Neira A, Salas A, Milne RL, Vega A, Carracedo B, González E, Barroso E, Fernández LP, Yankilevich P, Robledo M, Carracedo A, Benitez J: Evaluating HapMap SNP data transferability in a large-scale genotyping project involving 175 cancer-associated genes. Hum Genet. 2006, 118 (6): 669-679. 10.1007/s00439-005-0094-9.
Salas A, Jaime JC, Álvarez-Iglesias V, Carracedo Á: Gender bias in the multi-ethnic genetic composition of Central Argentina. J Hum Genet. 2008, 53: 662-674. 10.1007/s10038-008-0297-8.
Heinz T, Alvarez-Iglesias V, Taboada-Echalar P, Gomez-Carballa A, Torres-Balanza A, Rocabado O, Carracedo A, Vullo C, Salas A: Ancestry analysis reveals a predominant Native American component with moderate European admixture in Bolivians. Forensic Sci Int Genet. 2013, 7 (2013): 537-542. in press
Alexander DH, Novembre J, Lange K: Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19 (9): 1655-1664. 10.1101/gr.094052.109.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.
Shriver MD, Kennedy GC, Parra EJ, Lawson HA, Sonpar V, Huang J, Akey JM, Jones KW: The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs. Hum Genomics. 2004, 1 (4): 274-286. 10.1186/1479-7364-1-4-274.
Amigo J, Salas A, Phillips C, Carracedo A: SPSmart: adapting population based SNP genotype databases for fast and comprehensive web access. BMC Bioinformatics. 2008, 9: 6-10.1186/1471-2105-9-6.
Amigo J, Salas A, Phillips C: ENGINES: exploring single nucleotide variation in entire human genomes. BMC Bioinformatics. 2011, 12: 6-10.1186/1471-2105-12-6.
Corach D, Lao O, Bobillo C, van Der Gaag K, Zuniga S, Vermeulen M, van Duijn K, Goedbloed M, Vallone PM, Parson W, de Knijff P, Kayser M: Inferring continental ancestry of argentineans from Autosomal, Y-chromosomal and mitochondrial DNA. Ann Hum Genet. 2010, 74 (1): 65-76. 10.1111/j.1469-1809.2009.00556.x.
Halder I, Shriver M, Thomas M, Fernandez JR, Frudakis T: A panel of ancestry informative markers for estimating individual biogeographical ancestry and admixture from four continents: utility and applications. Hum Mutat. 2008, 29 (5): 648-658. 10.1002/humu.20695.
Kosoy R, Nassir R, Tian C, White PA, Butler LM, Silva G, Kittles R, Alarcon-Riquelme ME, Gregersen PK, Belmont JW, De La Vega FM, Seldin MF: Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America. Hum Mutat. 2009, 30 (1): 69-78. 10.1002/humu.20822.
Nassir R, Kosoy R, Tian C, White PA, Butler LM, Silva G, Kittles R, Alarcon-Riquelme ME, Gregersen PK, Belmont JW, De La Vega FM, Seldin MF: An ancestry informative marker set for determining continental origin: validation and extension using human genome diversity panels. BMC Genet. 2009, 10: 39-
Bandelt H-J, Salas A: Current next generation sequencing technology may not meet forensic standards. Forensic Sci Int Genet. 2012, 6 (1): 143-145. 10.1016/j.fsigen.2011.04.004.
Sankar P, Cho MK: Genetics. Toward a new vocabulary of human genetic variation. Science. 2002, 298 (5597): 1337-1338. 10.1126/science.1074447.
The research leading to these results has received funding from the “Ministerio de Ciencia e Innovación” (SAF2008-02971) and from the Plan Galego IDT, Xunta de Galicia (EM 2012/045) (A.S.) and Consellería de Sanidade/Xunta de Galicia (RHI07/2-intensificación actividad investigadora and 10PXIB918184PR), Instituto Carlos III (Intensificación de la actividad investigadora) and Fondo de Investigación Sanitaria (FIS; PI070069 and PI1000540) del Plan Nacional de I + D + I and ‘fondos FEDER’ (F.M.T.), and the grant from the Sistema Universitario Gallego- Modalidad REDES (2012-PG226) of the Consellería de Cultura, Educación e Ordenación Universitaria of the Xunta de Galicia (A.S., F.M.T.). There are no conflicts of interest in this study.
The authors declare that they have no competing interests.
JPS, FMT and AS conceived the study. JPS and AS and analysed the data. AS drafted the article and all the authors have critically revised the manuscript and given final approval of the version to be published. All authors read and approved the final manuscript.
Jacobo Pardo-Seco and Antonio Salas contributed equally to this work.
Electronic supplementary material
Additional file 2: PCA plots of East Asian populations using the whole set of SNPs in HapMap and 500 re-samples of rSNPs taking at random 100,00, 50,000, 10,000, 1,000 and 500 SNPs. As in Figure 2, the profiles that correspond to one of the re-samples are shown in color, while the remaining re-samples are shown in grey. (TIFF 7 MB)
Additional file 3: The distributions on the left show the cumulative LSBL values for the different AIM panels compared to the distributions generated by the best AIMs in the HapMap datasets (indicated in the figure as “present study”). The bar-plots on the right mirror the ancestry inferred using these panels on the HapMap populations TSI (representing Europe), CHD (representing East Asia) and LWK (representing sub-Saharan Africa), in order to reflect the ‘portability’ of the different SNP panels in other population groups. (TIFF 8 MB)
Additional file 4: Correlation values between the standard errors computed when comparing genomic ancestry with the estimated ancestries using different AIM panels, and considering admixed and non-admixed genomes. These values correspond to the distributions in Figure 6. (DOC 33 KB)
Authors’ original submitted files for images
About this article
Cite this article
Pardo-Seco, J., Martinón-Torres, F. & Salas, A. Evaluating the accuracy of AIM panels at quantifying genome ancestry. BMC Genomics 15, 543 (2014). https://doi.org/10.1186/1471-2164-15-543