SNP@Evolution: a hierarchical database of positive selection on the human genome
© Cheng et al; licensee BioMed Central Ltd. 2009
Received: 25 November 2008
Accepted: 5 September 2009
Published: 5 September 2009
Positive selection is a driving force that has shaped the modern human. Recent developments in high throughput technologies and corresponding statistics tools have made it possible to conduct whole genome surveys at a population scale, and a variety of measurements, such as heterozygosity (HET), F ST , and Tajima's D, have been applied to multiple datasets to identify signals of positive selection. However, great effort has been required to combine various types of data from individual sources, and incompatibility among datasets has been a common problem. SNP@Evolution, a new database which integrates multiple datasets, will greatly assist future work in this area.
As part of our research scanning for evolutionary signals in HapMap Phase II and Phase III datasets, we built SNP@Evolution as a multi-aspect database focused on positive selection. Among its many features, SNP@Evolution provides computed F ST and HET of all HapMap SNPs, 5+ HapMap SNPs per qualified gene, and all autosome regions detected from whole genome window scanning. In an attempt to capture multiple selection signals across the genome, selection-signal enrichment strength (ES) values of HET, F ST , and P-values of iHS of most annotated genes have been calculated and integrated within one frame for users to search for outliers. Genes with significant ES or P-values (with thresholds of 0.95 and 0.05, respectively) have been highlighted in color. Low diversity chromosome regions have been detected by sliding a 100 kb window in a 10 kb step. To allow this information to be easily disseminated, a graphical user interface (GBrowser) was constructed with the Generic Model Organism Database toolkit.
Available at http://bighapmap.big.ac.cn, SNP@Evolution is a hierarchical database focused on positive selection of the human genome. Based on HapMap Phase II and III data, SNP@Evolution includes 3,619,226/1,389,498 SNPs with their computed HET and F ST , as well as qualified genes of 21,859/21,099 with ES values of HET and F ST . In at least one HapMap population group, window scanning for selection signals has resulted in 1,606/10,138 large low HET regions. Among Phase II and III geographical groups, 660 and 464 regions show strong differentiation.
Natural selection has played an essential role in the formation of the human genome and in the diversity of phenotypes. The identification of the functional targets of positive selection, however, has been a major challenge in understanding the evolution of human beings. Traditional investigations to localize loci that have undergone selection have focused on the coding sequences of individual genes. For example, since the finding of hemoglobin B gene under the selective pressure for malaria resistance, only a limited number of genes, including G6PD, LCT, and ASPM, have been functionally determined as targets of positive selection [1–3]. In recent years, the rapid development in high throughput DNA technologies, as well as in statistical analysis and bioinformatics tools, has promoted whole genome surveys in multiple aspects of genetic variation [4–6]. Since the construction of the human HapMap, many massive genome-wide projects aiming to search and scan for SNPs, indels, copy number variations (CNVs), functional DNA elements, DNA methylation sites, and expression quantitative trait loci (eQTL) have been accomplished or undertaken [7–17].
These enormous maps and datasets have made great contributions to trace our evolutionary history. Although only a handful of genome-wide measurements with limited marker density have been developed to detect selection signals, these initial studies suggest that large scale and highly detailed analyses will greatly illuminate our understanding of human evolution [18–21]. Using whole genome SNP data, selection signals may be demonstrated through the computation of classical measurements including heterozygosity (a measurement which is used to estimate the frequency of heterozygote in a population, also referred to as HET), Tajima's D (a statistical test which is used to determine whether a genetic locus is under neutral selection), and the fixation index (a measurement which is used to compare the genetic variability within and between populations, also referred to as F ST ) [18, 22, 23]. Moreover, haplotype based measurements--such as extended haplotype homozygosity (EHH) and relative EHH (REHH), or more complex values including integrated haplotype homozygosity (iHH) and integrated haplotype score (iHS)--have been successfully applied to find signals of recent positive selection across the human genome [24, 25]. A few databases, such as Haplotter and SNP@Ethnos, serve as public tools in population genetics . Built to find ethnically related SNPs, SNP@Ethnos provides 100,736 individual variants with large ethnic differences. With the determination of ancestral or derived allele state, Haplotter is efficient in finding recent positive selections and their affected genes by the frequency of haplotypes that extend from a core SNP.
To complement these and other publicly available databases, we have built an integrated data library on human evolution ("SNP@Evolution") which offers several novel features: (1) the inclusion of ES (enrichment strength of selection-signals) to estimate the selection strength on a specific gene by computation of outlier ratios; (2) a sliding window scanning method that uses measures of HET and F ST to locate genomic segments under selection (with the resulting regional signals reflecting geographical adaptation, founder effects, and fixed or unfixed selections); (3) HapMap Phase III data from samples of 11 populations; (4) the integration with a haplotype-based dataset of iHS that simplifies the process of comparing multiple datasets by providing a simple, easy-to-read table.
One of SNP@Evolution's primary functions is to find selection signals from both chromosome regions and individual SNPs. To achieve this aim, our major strategy has involved the computation and comparison of HET, F ST , and iHS of all SNPs and their regional values using a sliding window method, then followed by the demonstration of outliers of each measurement as the selection signals. ES values of various measurements in most annotated genes are also listed as an independent dataset. Selection signals can therefore be detected in the data query or visualization interface by (1) SNP outliers, (2) low HET regions merged from adjacent window outliers, and (3) ES values of individual genes. With the comparison and integration of multiple datasets, SNP@Evolution provides candidate genes and regions that will assist researchers in locating positive selection signals genome-wide.
Construction and content
We first computed HET and Akey's F ST to measure the polymorphism within each population and the differentiation among geographical groups . For gene analysis, we considered a genic region to be a gene with 2 kb flanking regions. Genic regions containing no less than 5 genotyped SNPs in the dataset were chosen for subsequent analysis. From 24,011 annotated genes in NCBI build36, we identified 21,859 and 21,099 qualified genic regions in HapMap Phase II and Phase III, respectively.
To search for outliers of HET and F ST , the SNPs with HET values larger than average were considered as outliers directly. For F ST , we defined an "outlier" as a SNP whose value extends beyond the nearest quartile by a length of at least 1.5 times the inter-quartile range. With simulated data, we first sought to determine whether this method (henceforth referred to as the "empirical distribution method") is efficient and reliable to obtain selection signals without having to consider the apriori distributions. Genotype data under different selection strength and neutral selection generated by SelSim were mixed together , then the empirical distribution method was applied to find outliers. The results indicated that the empirical distribution method separated data under positive selection from the mixed data pool with a high sensitivity, stability, and low error rate (manuscript in preparation). (To test the empirical distribution, the positive predictive value, or PPV, can be used to describe the proportion of outliers that have been the targets of selection . For instance, when the simulated proportion of selection target is 0.05 with a certain recombination environment--Rec rate = 10-8, the PPV of the empirical distribution ranges from 0.50 to 0.90 under the selection strength 20-200.) Therefore, we used the empirical distribution method to identify SNP outliers of F ST as well as the window outliers of HET (see below).
The ratios of outlier numbers of both HET and F ST in each genic region to the corresponding total SNPs were then calculated. Relevant enrichment value ES of both HET and F ST were determined as the percentile values on the distribution of outlier ratios for all qualified genes. Thus, the ES value of HET or F ST in a particular gene represents the enrichment strength of the outliers of each measure. (Distributing from 0 to 1, a larger ES value indicates a higher outlier ratio of the corresponding measure in a gene.)
To estimate regional genetic diversity, we set up a 100 kb window and slid it with a 10 kb step throughout the human autosomes to obtain the averages of HET and F ST . This window scanning resulted in 267,069 and 266,650 regions (with 90 kb overlapping each other) in HapMap Phase II and Phase III data, respectively. Regional HET was first normalized within each chromosome and then the empirical distribution method described above was used to obtain window outliers of F ST and normalized HET. Adjacent window outliers were merged to regions and the qualification of each region was determined by a bootstrap test. Briefly, we picked 1,000 regions randomly with the similar length to the testing region, and then calculated the HET values of these picked regions plus the testing region. The percentile values of these 1,001 HET values were defined as P boot , and only those with P boot < 0.01 were accepted as low HET regions.
Consequently, 1,606 low HET regions were obtained from Phase II Build36 data, including 434 in ASN, 576 in CEU, and 596 in YRI. In addition, merged window outliers of F ST identified 660 highly differentiated regions among three groups. Moreover, 10,138 low HET regions were found from Phase III data, including 807 in CHB, 802 in CHD, 803 in JPT, 995 in CEU, 1,000 in TSI, 881 YRI, 965 in ASW, 916 in LWK, 962 in MKK, 1,023 in GIH, 984 in MEX. Additionally, 464 regions show strong differentiation among four geographical groups.
Datasets containing evolutionary signals
The distribution of F ST outliers across autosomes in data of HapMap Phase II Build36
As for all individual 24,011 genes listed with annotations, 21,859 genes with 5 or more genotyped SNPs in their genic regions are provided with ES values of HET and F ST , resulting in 1,094 genes of low HET which exceed the ES threshold of 0.95. In these 1,094 genes, ASN and CEU share 370, ASN and YRI share 147, CEU and YRI share 169, and all three groups share 81. Such a pattern--non-Africans sharing twice the number of genes containing enriched selection signals in comparison with that between non-African and African populations--suggests either geographical selection or founder effects by ancestors of ASN and CEU after their migration out of Africa. (Accordingly, we performed simulations based on both selection and neutral models. Very little founder effect was revealed (data not shown), suggesting that the selection events most likely resulted in these signals in genes. Our analysis results are consistent with recent reports which showed also geographic selections rather than founder effects as the major force of evolution in large genomic regions [29, 30].) For the 535 genes of high F ST (ES threshold of 0.95), 369 were seen in at least one population which showing low HET (247 genes in ASN, 163 genes in CEU, and 42 in YRI respectively). In addition, P-values of iHS are included in 10,375 genes. For those with significant iHS (P < 0.05), 441, 527, and 433 genes are in ASN, CEU, and YRI, respectively, among which 32/441, 43/527, and 41/433 are also shown as low HET.
There are 1,389,498 genotyped SNPs in HapMap Phase III Build36, and among them, 351 outliers of F ST were found. ES values of HET and F ST in 21,099 qualified genic regions of all 11 populations are provided in our database. Finally, in our study we have found that for a fixed sample size and in regions with the same number of polymorphic loci, the regional HET is linear to Tajima's D [see Additional file 1]; therefore, we present HET data only in SNP@Evolution. Considering that iHS was derived from EHH and iHH, we have chosen iHS to be the measurement of haplotype diversity in our database.
Utility and Discussion
There are two user interfaces in SNP@Evolution, one for data query and another for data visualization. In the data query interface, users submit one or more SNP rs# labels, gene symbols, gene sequence accession IDs, or specific chromosome regions. Results are displayed in three tables, as shown in the EDAR example in Fig. 2. Measurements--including the highly differentiated region among three geographical groups, ES values of HET, F ST , and P-value of iHS of the gene, as well as the HET, F ST , and iHS of the SNPs--indicate that the EDAR gene is under positive selection in Asians as previously reported [19, 31]. Like all queries, these EDAR results can be downloaded in a table formatted file.
Furthermore, results in the data query interface can be easily linked to the visualization interface to see the features under the background of the chromosome. To help users obtain additional information, we also provide the SNP links to HapMap database and NCBI dbSNP database and the gene links to Entrez Nucleotide and OMIM. To facilitate the operation, a page of 'Complete Guide' is included for detailed introduction and efficient use of SNP@Evolution. The terms used in measurements and population groups in this database are hyperlinked to the 'Abbreviations' page for brief reference.
Generally speaking, SNP@Evolution allows users to access all data through the 'Download' link. In these datasets, most of computations were conducted based on HapMap and related projects. As the data came from no more than a hundred of individuals in each HapMap population group, one shall also consider the sample size while making conclusions. The aim of SNP@Evolution is to provide genome wide signals of positive selection on human being, to generate fine scale traces of natural selection in enlarged samples is a long term goal in our research and computation.
SNP@Evolution is a valuable and useful resource for finding and verifying signals of natural selection, and we will continue to update SNP@Evolution as research in positive selection progresses. At regions showing strong selection signals, we plan to add additional SNP information obtained from resequencing data from our work and from public datasets. Genotype data of new individuals from various sources will also be added to the database. By evaluating the effects of amino acid substitution with the method of Sorting Intolerant From Tolerant (SIFT) , we also plan to include protein functions at genic regions showing strong selection signals.
Availability and requirements
SNP@Evolution is freely available at http://bighapmap.big.ac.cn/. A minimum screen resolution of 1,152 × 864 is recommended. Please send all questions, comments, and suggestions to email@example.com.
We thank Guanghua Wu for editing this manuscript. We also appreciate the help from Dake Zhang and Yubin Sun on website construction. This work was supported by a grant of National Natural Science Foundation of China (30890031) to CZ.
- Verrelli BC, McDonald JH, Argyropoulos G, Destro-Bisol G, Froment A, Drousiotou A, Lefranc G, Helal AN, Loiselet J, Tishkoff SA: Evidence for balancing selection from nucleotide sequence analyses of human G6PD. American journal of human genetics. 2002, 71 (5): 1112-1128. 10.1086/344345.PubMed CentralView ArticlePubMedGoogle Scholar
- Bersaglieri T, Sabeti PC, Patterson N, Vanderploeg T, Schaffner SF, Drake JA, Rhodes M, Reich DE, Hirschhorn JN: Genetic signatures of strong recent positive selection at the lactase gene. American journal of human genetics. 2004, 74 (6): 1111-1120. 10.1086/421051.PubMed CentralView ArticlePubMedGoogle Scholar
- Mekel-Bobrov N, Gilbert SL, Evans PD, Vallender EJ, Anderson JR, Hudson RR, Tishkoff SA, Lahn BT: Ongoing adaptive evolution of ASPM, a brain size determinant in Homo sapiens. Science. 2005, 309 (5741): 1720-1722. 10.1126/science.1116815.View ArticlePubMedGoogle Scholar
- Biswas S, Akey JM: Genomic insights into positive selection. Trends Genet. 2006, 22 (8): 437-446. 10.1016/j.tig.2006.06.005.View ArticlePubMedGoogle Scholar
- Ronald J, Akey JM: Genome-wide scans for loci under selection in humans. Human genomics. 2005, 2 (2): 113-125.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang ET, Kodama G, Baldi P, Moyzis RK: Global landscape of recent inferred Darwinian selection for Homo sapiens. Proceedings of the National Academy of Sciences of the United States of America. 2006, 103 (1): 135-140. 10.1073/pnas.0509691102.PubMed CentralView ArticlePubMedGoogle Scholar
- The International HapMap Consortium: A haplotype map of the human genome. Nature. 2005, 437 (7063): 1299-1320. 10.1038/nature04226.PubMed CentralView ArticleGoogle Scholar
- Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449 (7164): 851-861. 10.1038/nature06258.View ArticlePubMedGoogle Scholar
- Conrad DF, Andrews TD, Carter NP, Hurles ME, Pritchard JK: A high-resolution survey of deletion polymorphism in the human genome. Nature genetics. 2006, 38 (1): 75-81. 10.1038/ng1697.View ArticlePubMedGoogle Scholar
- Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nature genetics. 2004, 36 (9): 949-951. 10.1038/ng1416.View ArticlePubMedGoogle Scholar
- Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, Thurman RE, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447 (7146): 799-816. 10.1038/nature05874.View ArticlePubMedGoogle Scholar
- Rakyan VK, Hildmann T, Novik KL, Lewin J, Tost J, Cox AV, Andrews TD, Howe KL, Otto T, Olek A, et al: DNA methylation profiling of the human major histocompatibility complex: a pilot study for the human epigenome project. PLoS Biol. 2004, 2 (12): e405-10.1371/journal.pbio.0020405.PubMed CentralView ArticlePubMedGoogle Scholar
- Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, et al: Global variation in copy number in the human genome. Nature. 2006, 444 (7118): 444-454. 10.1038/nature05329.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, Lee C, et al: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315 (5813): 848-853. 10.1126/science.1136678.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, Lyle R, Hunt S, Kahl B, Antonarakis SE, Tavare S, et al: Genome-wide associations of gene expression variation in humans. PLoS genetics. 2005, 1 (6): e78-10.1371/journal.pgen.0010078.PubMed CentralView ArticlePubMedGoogle Scholar
- Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al: A genome-wide association study of global gene expression. Nature genetics. 2007, 39 (10): 1202-1207. 10.1038/ng2109.View ArticlePubMedGoogle Scholar
- McCarroll SA, Hadnott TN, Perry GH, Sabeti PC, Zody MC, Barrett JC, Dallaire S, Gabriel SB, Lee C, Daly MJ, et al: Common deletion polymorphisms in the human genome. Nature genetics. 2006, 38 (1): 86-92. 10.1038/ng1696.View ArticlePubMedGoogle Scholar
- Akey JM, Zhang G, Zhang K, Jin L, Shriver MD: Interrogating a high-density SNP map for signatures of natural selection. Genome research. 2002, 12 (12): 1805-1814. 10.1101/gr.631202.PubMed CentralView ArticlePubMedGoogle Scholar
- Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne EH, McCarroll SA, Gaudet R, et al: Genome-wide detection and characterization of positive selection in human populations. Nature. 2007, 449 (7164): 913-918. 10.1038/nature06250.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. Human genomics. 2004, 1 (4): 274-286.PubMed CentralView ArticlePubMedGoogle Scholar
- Shriver MD, Mei R, Parra EJ, Sonpar V, Halder I, Tishkoff SA, Schurr TG, Zhadanov SI, Osipova LP, Brutsaert TD, et al: Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation. Human genomics. 2005, 2 (2): 81-89.PubMed CentralPubMedGoogle Scholar
- Carlson CS, Thomas DJ, Eberle MA, Swanson JE, Livingston RJ, Rieder MJ, Nickerson DA: Genomic regions exhibiting positive selection identified from dense genotype data. Genome research. 2005, 15 (11): 1553-1565. 10.1101/gr.4326505.PubMed CentralView ArticlePubMedGoogle Scholar
- Weir BS, Cardon LR, Anderson AD, Nielsen DM, Hill WG: Measures of human population structure show heterogeneity among genomic regions. Genome research. 2005, 15 (11): 1468-1476. 10.1101/gr.4398405.PubMed CentralView ArticlePubMedGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, et al: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837. 10.1038/nature01140.View ArticlePubMedGoogle Scholar
- Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol. 2006, 4 (3): e72-10.1371/journal.pbio.0040072.PubMed CentralView ArticlePubMedGoogle Scholar
- Park J, Hwang S, Lee YS, Kim SC, Lee D: SNP@Ethnos: a database of ethnically variant single-nucleotide polymorphisms. Nucleic acids research. 2007, D711-715. 10.1093/nar/gkl962. 35 Database
- Spencer CC, Coop G: SelSim: a program to simulate population genetic data with natural selection and recombination. Bioinformatics. 2004, 20 (18): 3673-3675. 10.1093/bioinformatics/bth417.View ArticlePubMedGoogle Scholar
- Kelley JL, Madeoy J, Calhoun JC, Swanson W, Akey JM: Genomic signatures of positive selection in humans and the limits of outlier approaches. Genome Res. 2006, 16 (8): 980-989. 10.1101/gr.5157306.PubMed CentralView ArticlePubMedGoogle Scholar
- Akey JM: Constructing genomic maps of positive selection in humans: where do we go from here?. Genome Res. 2009, 19 (5): 711-722. 10.1101/gr.086652.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, Srinivasan BS, Barsh GS, Myers RM, Feldman MW, et al: Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009, 19 (5): 826-837. 10.1101/gr.087577.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujimoto A, Kimura R, Ohashi J, Omi K, Yuliwulandari R, Batubara L, Mustofa MS, Samakkarn U, Settheetham-Ishida W, Ishida T, et al: A scan for genetic determinants of human hair morphology: EDAR is associated with Asian hair thickness. Hum Mol Genet. 2008, 17 (6): 835-843. 10.1093/hmg/ddm355.View ArticlePubMedGoogle Scholar
- Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, et al: The generic genome browser: a building block for a model organism system database. Genome research. 2002, 12 (10): 1599-1610. 10.1101/gr.403602.PubMed CentralView ArticlePubMedGoogle Scholar
- Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic acids research. 2003, 31 (13): 3812-3814. 10.1093/nar/gkg509.PubMed CentralView ArticlePubMedGoogle Scholar