Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Software

AWclust: point-and-click software for non-parametric population structure analysis

Xiaoyi Gao1* and Joshua D Starmer23

Author Affiliations

1 Miami Institute for Human Genomics, University of Miami Miller School of Medicine, Miami, FL 33136, USA

2 Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA

3 Curriculum in Toxicology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA

For all author emails, please log on.

BMC Bioinformatics 2008, 9:77  doi:10.1186/1471-2105-9-77


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/9/77


Received:14 August 2007
Accepted:31 January 2008
Published:31 January 2008

© 2008 Gao and Starmer; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Population structure analysis is important to genetic association studies and evolutionary investigations. Parametric approaches, e.g. STRUCTURE and L-POP, usually assume Hardy-Weinberg equilibrium (HWE) and linkage equilibrium among loci in sample population individuals. However, the assumptions may not hold and allele frequency estimation may not be accurate in some data sets. The improved version of STRUCTURE (version 2.1) can incorporate linkage information among loci but is still sensitive to high background linkage disequilibrium. Nowadays, large-scale single nucleotide polymorphisms (SNPs) are becoming popular in genetic studies. Therefore, it is imperative to have software that makes full use of these genetic data to generate inference even when model assumptions do not hold or allele frequency estimation suffers from high variation.

Results

We have developed point-and-click software for non-parametric population structure analysis distributed as an R package. The software takes advantage of the large number of SNPs available to categorize individuals into ethnically similar clusters and it does not require assumptions about population models. Nor does it estimate allele frequencies. Moreover, this software can also infer the optimal number of populations.

Conclusion

Our software tool employs non-parametric approaches to assign individuals to clusters using SNPs. It provides efficient computation and an intuitive way for researchers to explore ethnic relationships among individuals. It can be complementary to parametric approaches in population structure analysis.

Background

Population structure analysis is important to genetic association studies [1-4] and evolutionary investigations [5-9]. Many statistical methods have been proposed to infer population structure and to assign individuals to ethnically similar clusters using multilocus genotype data, among which there are two major categories: parametric and non-parametric approaches.

Parametric approaches usually need to estimate population parameters such as allele frequencies and genotype frequencies and calculate likelihood, assuming Hardy-Weinberg equilibrium (HWE) and linkage equilibrium (LE) among loci for each population [10,11]. Two representative programs for parametric approaches are: STRUCTURE, a Bayesian method which uses a Markov chain Monte Carlo (MCMC) algorithm based on the Gibbs sampler algorithm [10], and L-POP, a frequentist method which uses the Expectation-Maximization (EM) algorithm [11]. In the extended version of STRUCTURE (version 2.1), the program can account for loose linkage between loci, but not high background linkage disequilibrium (LD) [12,13]. High background LD increases the chance of spurious clusters [13]. There are many other parametric Bayesian methods [14-20] and frequentist methods [21,22], which require similar or more complicated model assumptions. Two major challenges for the parametric approaches are the accuracy of allele frequencies estimates with small sample sizes, and the model assumptions that may not hold for some data sets. Moreover, assumptions of LE or loosely linked loci put a restriction on the number of genome-wide SNP loci that can be used.

In contrast to parametric approaches, non-parametric approaches do not rely on model assumptions about the properties of the sub-populations, nor do they require allele frequency estimates. In situations where parametric model assumptions can not be verified, or there is only a limited number of individuals from a single sub-population, non-parametric methods are more powerful for inference. However, when the model assumptions do hold and allele frequencies can be accurately estimated, then parametric methods provide more information. Thus, the two approaches are complementary in that one method is stronger where the other is weaker.

As stated by Liu and Zhao [23], non-parametric methods use a two-stage design. They start by calculating pair-wise distances [6,7,9], or some other form of dimension reduction, e.g. singular value decomposition (SVD) [23], and then rely on statistical clustering methods, e.g. neighbor joining (NJ) [6,7], K-means method [23], principal coordinates analysis (PCoA) [9,24] or multidimensional scaling (MDS) [25,26], to separate individuals. Recently, Gao and Starmer proposed a non-parametric method for population structure analysis and showed its advantages when genome-wide SNPs are available [27]. Liu and Zhao also proposed a non-parametric approach [23], but it requires missing genotypes be imputed explicitly and the software is not widely available. In recent publications, researchers tend to use both parametric and non-parametric approaches in their reports [24,25,28].

Since its publication in 2000, the freely available program STRUCTURE has become quite popular and dominated population structure analysis, while the non-parametric methods have not received much attention. However, with the vast amount of genotype data available, non-parametric approaches may be preferred because of their robustness to model assumptions and fast calculation. Recently, it was shown in an empirical study that non-parametric methods can give accurate results in fine-scale population structure detection and even separated Chinese and Japanese individuals using genome-wide random SNPs [27]. The separation of Chinese and Japanese individuals was also observed by Purcell et al. using MDS [26].

R is a convenient fast growing statistical computing environment with considerable popularity in the research community. It is freely available on a wide range of platforms, comes with implementations of many standard statistical methods, and can be easily extended through packages. We borrowed the strength of R and developed an add-on package that specifically focused on non-parametric population structure analysis. The motivation behind the package is to make recent developments in non-parametric population structure analysis available to researchers with an easy to use and intuitive graphical interface.

Implementation

We have developed point-and-click software for non-parametric population structure exploration. The program was written in R and Tk and can be installed as an R package. The package is named AWclust (

    A
llele sharing distance and
    W
ard's minimum variance hierarchical clustering). The major modules enclosed are: allele sharing distance (ASD) calculation, MDS 2D plot, MDS 3D plot, hierarchical plot by Ward's minimum variance algorithm, gap statistic calculation for inferring the optimal number of clusters and saving the clustered results.

Results

We automated the non-parametric population structure analysis procedures and packed all the routine steps in an intuitive graphical interface. The software can save researchers' time in data exploration. The outputs from AWclust are ready for publication and further analysis. It is distributed as R installation packages available for all popular operating systems: MS Windows, Mac OS X and Linux/Unix. The software comes with two example data sets, is fully documented, and includes a tutorial to help users become familiar with how to use it.

After installing and loading the AWclust package, users are presented with a GUI interface (see Figure 1). To experiment with the program, users can load one of the two sample data sets: hapmap500 or perlegen500. The hapmap500 dataset contains 500 genome-wide random SNPs from 209 unrelated individuals from the HapMap project, specifically 60 Yoruba from Ibadan, Nigeria (YRI), 60 CEPH Utah residents with ancestry from northern and western Europe (CEU), 45 Han Chinese from Beijing, China (CHB), and 44 Japanese from Tokyo, Japan (JPT). The perlegen500 dataset contains 500 genome-wide random SNPs from 71 unrelated individuals from the Perlegen project, including 23 African Americans (AA), 24 European Americans (EA) and 24 Han Chinese (HC). In each dataset, the SNP information is encoded as numeric values (i.e. 0, 1, or 2) to represent the number of variant SNP alleles in genotypes, and -1 is used to represent missing values.

thumbnailFigure 1. The AWclust interface.

After loading a SNP data set, the next step is to calculate the allele sharing distance (ASD) matrix. Once this is done the user can perform non-parametric exploration with the SNP data set. The user can generate multidimensional scaling (MDS) 2D/3D plots to get a general idea of how the data clusters and to detect any outliers in the dataset. MDS is a statistical technique for allowing differences and similarities to be visualized. The differences are represented by distances between points on a graph. Elements in the ASD matrix that are close together will tend to cluster together.

If the MDS plot does not reveal any outliers in the dataset, then it makes sense to create a hierarchical plot of the data. For large datasets, the screen size and resolution may cause the IDs for the individuals to become difficult to read, however, this problem can be solved by saving and viewing the PDF output separately. The hierarchical plot can help users identify clusters and general relationships among individuals. The cluster tree can be cut at any level of similarity according to researchers' need and be saved to a text file.

One of the important features of AWclust is that it will calculate the gap statistic, a method for estimating the number of clusters in a data set [29]. It compares the pooled within-cluster sum of squares with its expectation from a null reference distribution. The precision of this method requires multiple simulations from the null reference distribution and thus, can be computationally intensive. For example, a dataset with 209 individuals, each with 1000 SNPs and running 60 simulations (K = 1 to 6) takes slightly longer than two minutes to run on an Intel Core2, 2.4 GHz CPU with 2 GB of RAM. The data points are then plotted for a range of cluster sizes and the optimal size maximizes the distance between the observed and expected pooled within-cluster sum of squares. Information from MDS and hierarchical plots may also help interpret of the gap statistic plots.

When using AWclust on datasets other than the provided samples, it is important to note that the more closely related the subpopulations are, the larger the number of genome-wide random SNP loci needed for good separation. Fortunately, AWclust can quickly process large datasets. For the number of SNP loci required for major human populations, users can refer the empirical studies by Gao and Starmer [27].

Discussion

Genome-wide SNPs can be easily obtained using Affymetrix or Illumina chips and thus, making full use of the vast amount of genetic markers available is an open issue for parametric approaches since the LE assumption does not hold when SNPs are densely genotyped. Without complicated LD modeling, there is a bottleneck for the amount of data that can be used by parametric approaches. Moreover, the allele frequency estimation in parametric approaches suffers from high variation if the number of individuals for a sub-population is extremely small, e.g. less than ten, and this leads to questionable results.

Besides STRUCTURE, there are several other parametric software packages that are compared by Wu et al. [22], such as PARTITION [15], BAPS 2 [19], GENELAND [20] and PSMIX [22]. However, all parametric approaches suffer similar limitations because they all require verification of their model assumptions. When they are not justified, and for densely genotyped SNPs or when there is only a limited number of individuals from a single sub-population this is likely to be the case, non-parametric approaches provide better results.

Recently researchers have began to use both parametric and non-parametric approaches to aid population structure analysis [24,25,28]. However, non-parametric approaches have been limited to the utilization of NJ, K-means method, PCoA or MDS, all of which can not give the optimal number of clusters objectively. Furthermore, to our knowledge, we are not aware of an easy to use software package for non-parametric population structure analysis. Therefore, it is imperative to promote a non-parametric software package that automates routine steps and can take advantage of the vast amount of SNP markers available while inferring the optimal K objectively.

A limitation for the AWclust software is that it does not estimate the proportion of genome that belongs to each sub-population, which would necessarily require allele frequency estimation and many other assumptions about HWE and LD, while AWclust assigns each individual to one and only one cluster. Furthermore, AWclust only handles SNPs, currently the most popular genetic marker. Also, we did not aim to cover as many features in population genetics analysis as possible, like Arlequin does [30] since we did not want to replicate functions available in other software. There is also algorithm/software targeting on admixture mapping, i.e. ADMIXMAP [14,16,17] and ANCESTRYMAP [31], which requires more complicated model assumptions and is out of the scope of non-parametric methods. AWclust performs classification rather than test statistic adjustment as genomic control [32-34] and EigenStrat [35] do. However, investigators can conduct further analysis conditional on the cluster information provided by AWclust.

A general challenge in population structure analysis is to infer the optimal number of populations, K, and this is no different for the AWclust software. Before using the gap statistic to infer the optimal K, we suggest to plot the data using MDS first in addition to the hierarchical plot. In our experience, hierarchical plots give general information about the clusters embedded in the data, while MDS is suitable for outlier detection in exploration of individual relationships. To give an extreme situation, if there is just one individual for a particular sub-population, it is unlikely this individual can form a stand-alone cluster among all of the hierarchical clusters. Therefore, it is better to take outlier individuals out before applying the hierarchical plot and the gap statistic. The optimal K can also be explained in combination with other prior information about populations in the data sets and experience of the field.

Conclusion

In summary, we have developed new software, AWclust, for non-parametric population structure exploration. The software does not require HWE and LE for sample individuals. Nor does it need to estimate allele frequency. Most importantly, it can identify the optimal number of populations objectively. AWclust provides a user friendly GUI interface and does not require any prior programming skill from users. This non-parametric software is complementary to the parametric population structure programs because it is useful when HWE and LE can not be assumed to hold and when the accuracy of allele frequency estimation is questionable.

Availability and requirements

Project name: AWclust

Project home page: http://awclust.sourceforge.net/ webcite

Online users' manual: http://awclust.sourceforge.net/docs/index.html webcite

Operating system(s): Platform independent

Programming language: R, Tk

Other requirements: R 2.5 or higher (with Tk Widgets)

License: GPL

Any restrictions to use by non-academics: none

Authors' contributions

XG designed and programmed the software. JS participated in programming and wrote the manual. XG and JS co-drafted the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

This work was supported in part by NIH grants NS39764 and NIEHS T32 ES007126. We are grateful for Drs. Eden Martin, Robert Tibshirani, Amy Kapp and Alan Dabney's help, which made the software possible.

References

  1. Lander ES, Schork NJ: Genetic dissection of complex traits.

    Science 1994, 265:2037-2048. PubMed Abstract | Publisher Full Text OpenURL

  2. Risch NJ: Searching for genetic determinants in the new millennium.

    Nature 2000, 405:847-856. PubMed Abstract | Publisher Full Text OpenURL

  3. Marchini J, Cardon L, Phillips M, Donnelly P: The effects of human population structure on large genetic association studies.

    Nat Genet 2004, 36(5):512-517. PubMed Abstract | Publisher Full Text OpenURL

  4. Freedman M, Reich D, Penney K, McDonald G, Mignault A, Patterson N, Gabriel S, Topol E, Smoller J, Pato C, Pato M, Petryshen T, Kolonel L, Lander E, Sklar P, Henderson B, Hirschhorn J, Altshuler D: Assessing the impact of population stratification on genetic association studies.

    Nat Genet 2004, 36:388-393. PubMed Abstract | Publisher Full Text OpenURL

  5. Cavalli-Sforza LL, Menozzi P, Piazza A: The history and geography of human genes. Princeton University Press. Princeton, NJ; 1994. OpenURL

  6. Bowcock A, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd J, Cavalli-Sforza L: High resolution of human evolutionary trees with polymorphic microsatellites.

    Nature 1994, 368:455-457. PubMed Abstract | Publisher Full Text OpenURL

  7. Mountain J, Cavalli-Sforza L: Multilocus genotypes, a tree of individuals, and human evolutionary history.

    Am J Hum Genet 1997, 61:705-718. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Rosenberg N, Pritchard J, Weber J, Cann H, Kidd K, Zhivotovsky L, Feldman M: Genetic structure of human populations.

    Science 2002, 298:2381-2385. PubMed Abstract | Publisher Full Text OpenURL

  9. Shriver M, Kennedy G, Parra E, Lawson H, Sonpar V, Huang J, Akey J, Jones K: The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs.

    Hum Genomics 2004, 1:274-286. PubMed Abstract OpenURL

  10. Pritchard JK, Stephens M, Donelly P: Inference of population structure using multilocus genotype data.

    Am J Hum Genet 2000, 67:945-959. Publisher Full Text OpenURL

  11. Purcell S, Sham P: Properties of structured association approaches to detecting population stratification.

    Hum Hered 2004, 58:93-107. PubMed Abstract | Publisher Full Text OpenURL

  12. Falush D, Stephens M, Pritchard J: Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies.

    Genetics 2003, 164:1567-1587. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Kaeuffer R, Réale D, Coltman DW, Pontier D: Detecting population structure using STRUCTURE software: effect of background linkage disequilibrium.

    Heredity 2007, 99:374-380. PubMed Abstract | Publisher Full Text OpenURL

  14. McKeigue PM, Carperter JR, Parra EJ, Shriver MD: Estimation of admixture and detection of linkage in admixed populations by a Bayesian approach: application to African-American populations.

    Ann Hum Genet 2000, 64:171-186. PubMed Abstract | Publisher Full Text OpenURL

  15. Dawson KJ, Belkhir K: A Bayesian approach to the identication of panmictic populations and the assignment of individuals.

    Genet Res 2001, 78:59-77. PubMed Abstract | Publisher Full Text OpenURL

  16. Hoggart CJ, Parra EJ, Shriver MD, Bonilla C, Kittles RA, Clayton DG, McKeigue PM: Control of confounding of genetic associations in stratified populatinos.

    Am J Hum Genet 2003, 72:1492-1504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Hoggart CJ, Shriver MD, Kittles RA, Clayton DG, McKeigue PM: Design and analysis of admixture mapping studies.

    Am J Hum Genet 2004, 74:965-78. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Corander J, Waldmann P, Sillanpaa MJ: Bayesian analysis of genetic differentiation between populations.

    Genetics 2003, 163:367-374. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Corander J, Waldmann P, Marttinen P, Sillanpaa MJ: BAPS 2: enhanced possibilities for the analysis of genetic population structure.

    Bioinformatics 2004, 20:2363-2369. PubMed Abstract | Publisher Full Text OpenURL

  20. Guillot G, Mortier F, Estoup A: Geneland: A program for landscape genetics.

    Molecular Ecology Notes 2005, 5:712-715. Publisher Full Text OpenURL

  21. Tang H, Peng J, Wang P, Rish N: Estimation of individual admixture: analytical and study design considerations.

    Genet Epi 2005, 28:289-301. Publisher Full Text OpenURL

  22. Wu B, Liu N, Zhao H: PSMIX: an R package for population structure inference via maximum likelihood method.

    BMC Bioinformatics 2006, 7:317. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Liu N, Zhao H: A non-parametric approach to population structure inference using multilocus genotypes.

    Human Genomics 2006, 2:353-364. PubMed Abstract OpenURL

  24. Bauchet M, McEvoy B, Pearson LN, Quillen EE, Sarkisian T, Hovhannesyan K, Deka R, Bradley DG, Shriver MD: Measuring European Population Stratification with Microarray Genotype Data.

    Am J Hum Genet 2007, 80:948-956. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. 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:680-690. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. 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:559-575. PubMed Abstract | Publisher Full Text OpenURL

  27. Gao X, Starmer J: Human population structure detection via multilocus genotype clustering.

    BMC Genet 2007, 8(1):34. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Guthery SL, Salisbury BA, Pungliya MS, Stephens JC, Bamshad M: The Structure of Common Genetic Variation in United States Populations.

    Am J Hum Genet 2007, 81:1221-1231. PubMed Abstract | Publisher Full Text OpenURL

  29. Tibshirani R, Walther G, Hastie T: Estimating the number of clusters in a data set via the gap statistic.

    J R Statist Soc B 2001, 63:411-423. Publisher Full Text OpenURL

  30. Excoffer L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis.

    Evolutionary Bioinformatics Online 2005, 1:47-50. OpenURL

  31. Patterson N, Hattangadi N, Lane B, Lohmueller KE, Hafler DA, Oksenberg JR, Hauser SL, Smith MW, O'Brien SJ, Altshuler D, Daly MJ, Reich D: Methods for high-density admixture mapping of disease genes.

    Am J Hum Genet 2004, 74:979-1000. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Devlin B, Roeder K: Genomic control for association studies.

    Biometrics 1999, 55:997-1004. PubMed Abstract | Publisher Full Text OpenURL

  33. Devlin B, Roeder K, Wasserman L: Genomic control, a new approach to genetic-based association studies.

    Theoretical Population Biology 2001, 60:155-166. PubMed Abstract | Publisher Full Text OpenURL

  34. Devlin B, Roeder K, Bacanu S: Unbiased methods for population-based association studies.

    Genet Epi 2001, 21:273-284. Publisher Full Text OpenURL

  35. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies.

    Nat Genet 2006, 38:904-9. PubMed Abstract | Publisher Full Text OpenURL