Genetic variation in taste receptor pseudogenes provides evidence for a dynamic role in human evolution
© Risso et al.; licensee BioMed Central Ltd. 2014
Received: 20 May 2014
Accepted: 1 September 2014
Published: 13 September 2014
Human bitter taste receptors are encoded by a gene family consisting of 25 functional TAS2R loci. In addition, humans carry 11 TAS2R pseudogenes, some of which display evidence for substantial diversification among species, showing lineage-specific loss of function. Since bitter taste is thought to help prevent the intake of toxic substances, diversity at TAS2R genes could reflect the action of natural selection on the ability to recognize some bitter compounds rather than others. Whether species-specific variation in TAS2R pseudogenes is solely the result of genetic drift or whether it may have been influenced by selection due to different feeding behaviors has been an open question.
In this study, we analyzed patterns of variation at human TAS2R pseudogenes in both African and non-African populations, and compared them to those observable in nonhuman primates and archaic human species. Our results showed a similar worldwide distribution of allelic variation for most of the pseudogenes, with the exception of the TAS2R6P and TAS2R18P loci, both of which presented an unexpected higher frequency of derived alleles outside Africa. At the TAS2R6P locus, two SNPs were found in strong linkage disequilibrium (r2 > 0.9) with variants in the functional TAS2R5 gene, which showed signatures of selection. The human TAS2R18P carried a species-specific stop-codon upstream of four polymorphic insertions in the reading frame. SNPs at this locus showed significant positive values in a number of neutrality statistics, and age estimates indicated that they arose after the homo-chimp divergence.
The similar distribution of variation of many human bitter receptor pseudogenes among human populations suggests that they arose from the ancestral forms by a unidirectional loss of function. However we explain the higher frequency of TAS2R6P derived alleles outside Africa as the effect of the balancing selection acting on the closely linked TAS2R5 gene. In contrast, TAS2R18P displayed a more complex history, suggesting an acquired function followed by a recent pseudogenization that predated the divergence of human modern and archaic species, which we hypothesize was associated with adaptions to dietary changes.
The perception of bitter taste is thought to have evolved as a protection from toxic and harmful foods ,. However, not all bitter compounds evoke the same aversive reaction in different species and, in particular, among primates . Therefore, it has been hypothesized that consumption or rejection response depends upon the relative occurrence of bitter and potentially toxic compounds in an animal’s diet . The number of different bitter taste receptors varies greatly across species , but all of them belong to a family of seven-transmembrane G protein coupled receptors (GPCRs) known as TAS2Rs or T2Rs. In humans, this gene family is encoded by 25 functional TAS2R loci, which reside on chromosomes 5, 7 and 12 . In addition to these genes, humans also carry 11 TAS2R pseudogenes . Polymorphisms in these genes have been shown to modulate the taste response to different compounds, both natural and synthetic -, and showed evidence of evolutionary pressures -, highlighting the importance of studying both TAS2R genes and pseudogenes.
Pseudogenes are considered genomic fossils, classically defined as genomic loci with sequence similarity to functional genes, but lacking coding potential ,, often due to disruptive mutations such as frameshifts, premature stop codons and deletions. However, recent studies have demonstrated that some pseudogenes may likely have a function, providing evidence for their plasticity and a dynamic role in evolution -. A number of authors have described these events with a “less-is-more” hypothesis, suggesting that gene loss, or pseudogenization, may serve as an engine of evolutionary change, especially in human evolution ,. In fact, previous phylogenetic analyses demonstrated that remarkable variation exists in both TAS2Rs genes and pseudogenes in different lineages, suggesting that the ability to perceive bitterness may be mostly due to the repertoire of TAS2R genes ,. For instance, the mouse genome contains 33 functional TAS2Rs loci but only three pseudogenes, and the evolutionary relationships between human and mouse genes were shown to fall into three categories, depending on their orthology: 1) one-to-one orthology; 2) one-to-multiple orthology; 3) multiple-to-one orthology ,,. In contrast, zebrafish and chicken have only four and three TAS2R loci, respectively ,. Comparative analyses showed that primate TAS2Rs had a higher ratio of nonsynonymous/synonymous substitutions and a lower selective pressure on this gene family compared to rodents ,.
To better understand the evolutionary mechanisms underlying taste receptor pseudogenization events in the human lineage, we performed a survey of the human genetic variation at nine human TAS2RP loci. We then compared these to the homologous pseudogenes present in archaic human forms (i.e. Neandertal and Denisovan), as well as in other five primate species, Chimpanzee, Gorilla, Orangutan, Gibbon and Mouse Lemur. We examined the intra- and inter-species variation patterns to search for possible footprints of natural selection at these loci, with the goal of understanding the potential adaptive role of taste receptor pseudogenes in the evolutionary history of modern humans.
Genes and data sets
TAS2R pseudogenes chromosomal positions (GRCh37)
Haploview  was used to identify Linkage Disequilibrium (LD) patterns and haplotype blocks at the selected loci using the solid spine approach. Genealogical relationships among inferred haplotypes were constructed using the median-joining algorithm implemented in the Network 4.5 program . Haplotypes distribution across human populations were investigated using PLINK v.1.07 .
Population genetics analyses
Arlequin v.3.5  was used to compute summary statistics, such as nucleotide diversity (π), estimated heterozygosity (EH) and number of polymorphic sites (PS). Chi square-tests were performed to compare allele frequencies among different populations and the adopted significance threshold was adjusted using the Bonferroni correction (i.e. adjusted p = p value X number of individual tests). Metric Multidimensional Scaling (MDS) analyses based on the obtained FST values, were performed with STATISTICA v. 6.0 (Stat-Soft Inc, Tulsa, OK).
MUSCLE v.3.3  (Edgar 2004) was used to perform multi-alignments of the examined sequences. A maximum likelihood tree was constructed with MEGA v.6.0  using the Tamura-Nei substitution model. To assess the relative support for each clade, bootstrap values were calculated from 10.000 analysis replicates, and the cut-off point for bootstrap replication was 50%.
Neutrality test and age estimates
To test whether patterns of allele/haplotype frequencies and tree topology were consistent with neutral expectations, we performed three neutrality tests. For this purpose, the DNASP package  was used to calculate Tajima’s D and Fu’s FS values at each locus. Since deviations found using these tests could be caused by selection and/or demographic processes (e.g. population expansion and/or bottlenecks), we also performed the Li’s MFDM test , which is more robust in the presence of population size changes. Finally, GENETREE , was used to infer the estimated age of selected variants and the time of the most recent common ancestor (TMRCA).
A total of 47 single nucleotide polymorphisms (SNPs) (Additional file 1: Table S1), annotated according to dbSNP Build 137, were observed in the nine selected pseudogenes. 32 of them had a minor allele frequency (MAF) above 5%. Unexpected patterns of genetic variation were observed for two pseudogenes, TAS2R18P and TAS2R6P, in populations from different continents. In particular, the derived alleles at four TAS2R18P SNPs (rs2290318, rs2290319, rs61928604 and rs61928603) showed significantly increased frequency (0.5) in non-African populations compared to African populations (0.053) (Fisher’s exact test, adjusted p < 0.01). Similarly, the derived alleles of TAS2R6P polymorphisms rs1859645 and rs11761380 showed a significantly different distribution (0.24 vs. 0.57) between Africa and other continents (Fisher’s exact test, adjusted p < 0.01). Both Neandertal and Denisovan genomes carried the derived alleles of the TAS2R18P SNPs rs2290318, rs2290319, rs61928604 and rs61928603. For the TAS2R6P SNPs rs1859645 and rs11761380, the Denisovan genome showed ancestral alleles, while the Neandertal genome was heterozygous at both SNPs. In addition, derived alleles of two SNPs located in TAS2R67P (rs319269 and rs34648613) and one in TAS2R64P (rs68071847) were more frequent in African populations with respect to non-African ones (0.68 vs. 0.11, Fisher’s exact test, adjusted p < 0.01), and were also present in both Neandertal and Denisovan genomes. Finally, the derived allele of one SNP located in TAS2R63P (rs2597986) was present only in a few African individuals (1.58%), as well as in the two archaic species.
Comparisons of TAS2R pseudogene haplotype blocks among the studied populations
To test the significance of the observed structure, the distribution of the H5 and H6 haplotypes in human populations was also investigated by means of Analysis of Molecular Variance (AMOVA). For both TAS2R6P and TAS2R18P, most of variation was accounted for by differences within populations (86.67% and 72.07%, respectively), with a smaller percentage attributed to differences among (13.2% and 27.12%) and within groups (0.13% and 0.81%). The global FST values were 0.13 for TAS2R6P and 0.28 for TAS2R18P, while those related to H5 and H6 haplotypes were 0.14 and 0.31, respectively.
We constructed median-joining networks for the TAS2R6 and TAS2R16 pseudogenes in order to better understand the relationships between the inferred haplotypes. The resultant topologies (Additional file 3: Figure S2 A-B) identified of two major clusters defined by the presence of either the derived (GC for TAS2R6P, CA for TAS2R18P) or ancestral (AA for TAS2R6P, GC for TAS2R18P) alleles at the H5 and H6 haplotypes.
Summary statistics and population structure
The pattern of diversity shown in most of these pseudogenes was in accordance with the usual distribution of human genetic variation , where the diversity is higher in African populations than in non-African ones. However, a different situation was observed in TAS2R6P and TAS2R18P. African groups (i.e. ASW, LWK and YRI) showed lower values of both nucleotide diversity (π) and estimated heterozigosity (EH) at theTAS2R18P locus. With the exception of the ASW group, the same pattern was observed for the TAS2R6P gene. The number of polymorphic sites (PS) was similar in all the studied populations (Additional file 4: Table S2 A-B).
However, DNA sequence similarities confirmed these results, showing a high percentage of TAS2R6P identity between Hominidae and Hylobatidae (>95%) and Prosimians, Hominidae and Hylobatidae (>75%). TAS2R18P sequences showed a lower degree of similarities (<70%) between primates and Prosimians (Additional file 6: Table S3 A-B).
Relationships with functional genes
A 200 kb (100 kb in both 5’ and 3’ directions) region surrounding both TAS2R6P and TAS2R18P was explored for linkage disequilibrium patterns. The TAS2R6P rs11761380 and rs1859645 polymorphisms were found to be in strong LD (r2 > 0.9) with SNPs rs62477710, rs10952507 and rs6962558 that reside in the functional TAS2R5 gene (Additional file 7: Table S4). In contrast, TAS2R18P showed no r2 values above 0.3 between variants. In addition, we analyzed the sequences of all the genes found in this enlarged genomic interval and compared them to the pseudogene sequences, in order to document the sequence similarities between pseudogenes and functional genes. TAS2R6P shares a common ancestor with the functional TAS2R5 gene (62.59% identity), while TAS2R18P has a high level of identity with the functional TAS2R42 gene and with the TAS2R67P pseudogene (66.59 and 67.02%, respectively) (Additional file 8: Figure S4). Pseudogenes derived from a very recent pseudogenization event (i.e. TAS2R64P) showed a much higher sequence similarity (91.32%) with their functional forms (i.e. TAS2R48), indicating ancient pseudogenization events at the TAS2R6P and TAS2R18P loci.
Test of neutrality and mutation age estimates
Finally, the ages of these variants were estimated using a coalescent-based method implemented in the GENETREE package. Estimated ages turned out to be considerably ancient. For SNPs in TAS2R18P, rs2290318 was estimated to be 1,553,750 +/− 632,500 years old, while rs2290319 was estimated to be 1,113,750 +/− 495,000 years old, with a TMRCA of 3,520,000 +/−1,198,750 years. Similar values were also obtained for the TAS2R6P SNPs rs11761380 and rs1859645, with age estimates of 1,382,462 +/− 326,856 and 1,600,637 +/− 421,558 years, and a TMRCA of 3,457,712 +/− 1,326,856 years.
A small number of studies have examined genes encoding functional bitter taste receptors and their corresponding pseudogenes from an evolutionary perspective ,,,,. In the present study, we focused on the recent evolutionary history of bitter taste receptor pseudogenes in modern humans. We found a modest level of variation and few haplotype blocks, as expected for pseudogenes in general. However, African populations showed a significantly reduced frequency of derived forms at TAS2R6P and TAS2R18P polymorphisms, at both the allelic and haplotypic level. Genomes of the archaic Neandertal and Denisova species also differed in the distribution of these derived alleles compared to modern human populations. In contrast to the usual pattern of human variation, the highest genetic diversity was found in non-African populations. Similarly, heterozygosity was higher in non-African populations. In addition, results from both population structure analyses and AMOVA demonstrated that most of variation was due to differences within populations. In fact, Europeans and Latin Americans fell into the same cluster and did not differ in allelic distribution, suggesting a similar pattern of SNP distribution in these two populations and highlighting the admixed nature of the 1000 Genomes American populations.
Previous work ,, indicated that the pseudogenization event of TAS2R6P is very old, predating the divergence of Antrhopoidea from Prosimians. The TAS2R18P pseudogenization event occurred later in the primate genealogy, but still prior to the appearance of the Hominidae family.
Our results were consistent with these findings: TAS2R6P variation patterns enabled the identification of a unique cluster grouping Prosimians, Hylobatidae and Hominidae families and high level of sequence similarity among primates (>95%), as well as between primates and Prosimians (>75%). TAS2R18P sequence similarity was lower between the mouse lemur and other primates (<70%) and the reconstructed phylogenetic trees showed longer branch lengths after the divergence of the two lineages. In addition, sequence similarity confirmed this view, showing a lower percentage (<70%) of identity between these two pseudogenes and their functional forms, when compared to the similarities between very recent pseudogenes and their functional forms (>90%).
The analysis of the genomic regions surrounding these pseudogenes showed that TAS2R6P polymorphisms rs11761380 and rs1859645 were in strong linkage disequilibrium (r2 > 0.9) with rs62477710, rs6962558 and rs6962558 SNPs of the functional TAS2R5 gene, suggesting that evolutionary forces acting on this gene could have driven the allele differentiation at the TAS2R6P locus. To test this assumption, Tajima’s D and Fu’s FS tests were performed: these analyses showed that TAS2R6P did not show evidence of departure from neutral expectations while the neighboring TAS2R5 gene, which shared more than 60% sequence similarity with TAS2R6P, showed significant positive values for both these statistics. Since positive values of these tests may indicate balancing selection or may be the result of the confounding effect due to demographic history ,, the maximum frequency of derived mutations (MFDM) test was applied to further investigate the evolutionary history of these regions. The results of this test indicate that TAS2R6P has not undergone recent selective pressures, whereas TAS2R5 showed significant signatures of selection. This is consistent with balancing selection maintaining multiple alleles for long evolutionary times, and extending to maintain allele frequencies at closely linked neutral sites ,. The application of this test to TAS2R18 also showed significant signatures of selection and both Tajima’s D and Fu’s FS values were significantly positive at this locus in European and Latin American populations. This suggests that balancing selection potentially acted on the genes in these groups. In addition, these populations showed increased values of nucleotide diversity and heterozygosity in the TAS2R18P pseudogene compared to neighboring loci. Moreover, global FST values were unusually high compared to the typical range of 0.10 to 0.16 for estimated FST values in global populations ,, suggesting a high level of genetic differentiation among worldwide populations. These data are consistent with a scenario of balancing selection maintaining TAS2R18P alleles and enhancing genetic diversity at this locus in European and Latin American populations.
Variants at TAS2R18P showed no evidence of association with adjacent functional genes. This pseudogene carries a human-specific stop-codon that is shared among H. sapiens Neandertal and Denisova. In addition, the exonic region of this pseudogene carried four polymorphic insertions in the reading frame, upstream of the shared stop-codons and downstream of the human-specific one. These data suggest that these insertions may have shifted the human TAS2R18P reading frame, with a consequent acquired function of this gene, shortly after the homo-chimp divergence. Such an event, followed by balancing selection operating outside Africa, would produce the observed TAS2R18P rs2290318 and rs2290318 different allele distributions in human populations. The estimated ages of these mutations indicate that they arose after the divergence of humans from chimpanzee that occurred 7–8 million years ago . We hypothesize that the human-specific stop codon located upstream of all these insertions represents a second inactivation, which happened before the split between H. sapiens, Neandertal and Denisova (from 400,000 to 800,000 years ago).
Our results provide evidence for a dynamic role for TAS2R18P in primate evolution, suggesting that this locus may have acquired its function during the evolution of the human lineage, shortly after the homo-chimp divergence. This was followed by a much more recent deactivation due to the stop-gain mutation which was shared among modern humans, Neandertals and Denisovans. We speculate that this event may have been due to the disappearance of some bitter compound only found in Africa, which was specifically recognized by the product of this pseudogene, and that therefore the functionality of TAS2R18P was useful only in a given stage of human evolution.
Availability of supporting data
The article does not report new empirical data since the analyzed sequences were already deposited at public databases, including the 1000 Genomes Project, the UCSC Genome Browser and the Ensemble Genome Browser (see the Methods section for further details).
This research was supported by the National Institute on Deafness and Other Communication Disorders Intramural grant number Z1A-000046-14. We would like to thank Dr. Marco Sazzini for help in preparing the manuscript and Drs. Igor Rogozin and Rob Morell for their helpful suggestions.
- Garcia J, Hankins WG: The evolution of bitter and the acquisition of toxiphobia. Olfaction and Taste. 1975, 5: 39-45.Google Scholar
- Chandrashekar J, Mueller KL, Hoon MA, Adler E, Feng L, Guo W, Zuker CS, Ryba NJ: T2Rs function as bitter taste receptors. Cell. 2000, 100 (6): 703-711. 10.1016/S0092-8674(00)80706-0.PubMedView ArticleGoogle Scholar
- Imai H, Suzuki N, Ishimaru Y, Sakurai T, Yin L, Pan W, Abe K, Misaka T, Hirai H: Functional diversity of bitter taste receptor TAS2R16 in primates. Biol Lett. 2012, 8 (4): 652-656. 10.1098/rsbl.2011.1251.PubMedPubMed CentralView ArticleGoogle Scholar
- Glendinning JI: Is the bitter rejection response always adaptive?. Physiol Behav. 1994, 56 (6): 1217-1227. 10.1016/0031-9384(94)90369-7.PubMedView ArticleGoogle Scholar
- Shi P, Zhang J: Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol Biol Evol. 2006, 23 (2): 292-300. 10.1093/molbev/msj028.PubMedView ArticleGoogle Scholar
- Adler E, Hoon MA, Mueller KL, Chandrashekar J, Ryba NJ, Zuker CS: A novel family of mammalian taste receptors. Cell. 2000, 100 (6): 693-702. 10.1016/S0092-8674(00)80705-9.PubMedView ArticleGoogle Scholar
- Go Y, Satta Y, Takenaka O, Takahata N: Lineage-specific loss of function of bitter taste receptor genes in humans and nonhuman primates. Genetics. 2005, 170 (1): 313-326. 10.1534/genetics.104.037523.PubMedPubMed CentralView ArticleGoogle Scholar
- Bufe B, Breslin PA, Kuhn C, Reed DR, Tharp CD, Slack JP, Kim UK, Drayna D, Meyerhof W: The molecular basis of individual differences in phenylthiocarbamide and propylthiouracil bitterness perception. Curr Biol. 2005, 15 (4): 322-327. 10.1016/j.cub.2005.01.047.PubMedPubMed CentralView ArticleGoogle Scholar
- Wooding S, Gunn H, Ramos P, Thalmann S, Xing C, Meyerhof W: Genetics and bitter taste responses to goitrin, a plant toxin found in vegetables. Chem Senses. 2010, 35 (8): 685-692. 10.1093/chemse/bjq061.PubMedView ArticleGoogle Scholar
- Risso D, Morini G, Pagani L, Quagliariello A, Giuliani C, De Fanti S, Sazzini M, Luiselli D, Tofanelli S: Genetic signature of differential sensitivity to stevioside in the Italian population. Genes Nutr. 2014, 9 (3): 401-10.1007/s12263-014-0401-y.PubMedPubMed CentralView ArticleGoogle Scholar
- Wooding S, Kim UK, Bamshad MJ, Larsen J, Jorde LB, Drayna D: Natural selection and molecular evolution in PTC, a bitter-taste receptor gene. Am J Hum Genet. 2004, 74 (4): 637-646. 10.1086/383092.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell MC, Ranciaro A, Froment A, Hirbo J, Omar S, Bodo JM, Nyambo T, Lema G, Zinshteyn D, Drayna D, Breslin PA, Tishkoff SA: Evolution of functionally diverse alleles associated with PTC bitter taste sensitivity in Africa. Mol Biol Evol. 2012, 29 (4): 1141-1153. 10.1093/molbev/msr293.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell MC, Ranciaro A, Zinshteyn D, Rawlings-Goss R, Hirbo J, Thompson S, Woldemeskel D, Froment A, Rucker JB, Omar SA, Bodo JM, Nyambo T, Belay G, Drayna D, Breslin PA, Tishkoff SA: Origin and differential selection of allelic variation at TAS2R16 associated with salicin bitter taste sensitivity in Africa. Mol Biol Evol. 2014, 31 (2): 288-302. 10.1093/molbev/mst211.PubMedPubMed CentralView ArticleGoogle Scholar
- Balakirev ES, Ayala FJ: Pseudogenes: are they “junk” or functional DNA?. Annu Rev Genet. 2003, 37: 123-151. 10.1146/annurev.genet.37.040103.103949.PubMedView ArticleGoogle Scholar
- Zhang Z, Gerstein M: Large-scale analysis of pseudogenes in the human genome. Curr Opin Genet Dev. 2004, 14 (4): 328-335. 10.1016/j.gde.2004.06.003.PubMedView ArticleGoogle Scholar
- Bekpen C, Marques-Bonet T, Alkan C, Antonacci F, Leogrande MB, Ventura M, Kidd JM, Siswara P, Howard JC, Eichler EE: Death and resurrection of the human IRGM gene. PLoS Genet. 2009, 5 (3): e1000403-10.1371/journal.pgen.1000403.PubMedPubMed CentralView ArticleGoogle Scholar
- Korrodi-Gregório L, Abrantes J, Muller T, Melo-Ferreira J, Marcus K, da Cruz e Silva QA, Fardilha M, Esteves PJ: Not so pseudo: the evolutionary history of protein phosphatase 1 regulatory subunit 2 and related pseudogenes. BMC Evol Biol. 2013, 13: 242-10.1186/1471-2148-13-242.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Q: Using pseudogene database to identify lineage-specific genes and pseudogenes in humans and chimpanzees. J Hered. 2014, 105 (3): 436-443. 10.1093/jhered/est097.PubMedView ArticleGoogle Scholar
- Olson MV: When less is more: gene loss as an engine of evolutionary change. Am J Hum Genet. 1999, 64 (1): 18-23. 10.1086/302219.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang X, Grus WE, Zhang J: Gene losses during human origins. PLoS Biol. 2006, 4 (3): e52-10.1371/journal.pbio.0040052.PubMedPubMed CentralView ArticleGoogle Scholar
- Shi P, Zhang J, Yang H, Zhang YP: Adaptive diversification of bitter taste receptor genes in Mammalian evolution. Mol Biol Evol. 2003, 20 (5): 805-814. 10.1093/molbev/msg083.PubMedView ArticleGoogle Scholar
- Dong D, Jones G, Zhang S: Dynamic evolution of bitter taste receptor genes in vertebrates. BMC Evol Biol. 2009, 9: 12-10.1186/1471-2148-9-12.PubMedPubMed CentralView ArticleGoogle Scholar
- Conte C, Ebeling M, Marcuz A, Nef P, Andres-Barquin PJ: Evolutionary relationships of the Tas2r receptor gene families in mouse and human. Physiol Genomics. 2003, 14 (1): 73-82.PubMedView ArticleGoogle Scholar
- Go Y: Proceedings of the SMBE Tri-National Young InvestigatorsInvestigators’ Workshop. Lineage-specific expansions and contractions of the bitter taste receptor gene repertoire in vertebrates. Mol Biol Evol. 2005, 23 (5): 964-972. 10.1093/molbev/msj106.View ArticleGoogle Scholar
- Wang X, Thomas SD, Zhang J: Relaxation of selective constraint and loss of function in the evolution of human bitter taste receptor genes. Hum Mol Genet. 2004, 13 (21): 2671-2678. 10.1093/hmg/ddh289.PubMedView ArticleGoogle Scholar
- Fischer A, Gilad Y, Man O, Pääbo S: Evolution of bitter taste receptors in humans and apes. Mol Biol Evol. 2005, 22 (3): 432-436. 10.1093/molbev/msi027.PubMedView ArticleGoogle Scholar
- Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, Kang HM, Marth GT, McVean GA: An integrated map of genetic variation from 1,092 human genomes. Nature. 2012, 491 (7422): 56-65. 10.1038/nature11632.PubMedView ArticleGoogle Scholar
- 1000 Genomes Browser.. , [http://browser.1000genomes.org/index.html]
- Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, Kent WJ: The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004, 32: D493-D496. 10.1093/nar/gkh103.PubMedPubMed CentralView ArticleGoogle Scholar
- UCSC Table Browser.. , [http://genome.ucsc.edu/cgi-bin/hgTables]
- Flicek P, Ahmed I, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, García-Girón C, Gordon L, Hourlier T, Hunt S, Juettemann T, Kähäri AK, Keenan S, Komorowska M, Kulesha E, Longden I, Maurel T, McLaren WM, Muffato M, Nag R, Overduin B, Pignatelli M, Pritchard B, Pritchard E, et al: Ensembl 2013.Nucleic Acids Res41(D1):D48–D55.,Google Scholar
- Ensembl Genome Browser.. , [http://www.ensembl.org/index.html]
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.PubMedView ArticleGoogle Scholar
- Bandelt HJ, Forster P, RRöhlhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48. 10.1093/oxfordjournals.molbev.a026036.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Excoffier L, Lischer HE: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10 (3): 564-567. 10.1111/j.1755-0998.2010.02847.x.PubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013, 30 (12): 2725-2729. 10.1093/molbev/mst197.PubMedPubMed CentralView ArticleGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.PubMedView ArticleGoogle Scholar
- Li H: A new test for detecting recent positive selection that is free from the confounding impacts of demography. Mol Biol Evol. 2011, 28 (1): 365-375. 10.1093/molbev/msq211.PubMedView ArticleGoogle Scholar
- Bahlo M, Griffiths RC: Inference from gene trees in a subdivided population. Theor Popul Biol. 2000, 57 (2): 79-95. 10.1006/tpbi.1999.1447.PubMedView ArticleGoogle Scholar
- GENETREE software.. , [http://www.stats.ox.ac.uk/~griff/software.html]
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMedPubMed CentralGoogle Scholar
- Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133 (3): 693-709.PubMedPubMed CentralGoogle Scholar
- Schierup MH, Charlesworth D, Vekemans X: The effect of hitch-hiking on genes linked to a balanced polymorphism in a subdivided population. Genet Res. 2000, 76 (1): 63-73. 10.1017/S0016672300004547.PubMedView ArticleGoogle Scholar
- Navarro A, Barton NH: The effects of multilocus balancing selection on neutral variability. Genetics. 2002, 161 (2): 849-863.PubMedPubMed CentralGoogle Scholar
- Weir BS, Cardon LR, Anderson AD, Nielsen DM, Hill WG: Measures of human population structure show heterogeneity among genomic regions. Genome Res. 2005, 15 (11): 1468-1476. 10.1101/gr.4398405.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheng F, Chen W, Richards E, Deng L, Zeng C: SNP@Evolution: a hierarchical database of positive selection on the human genome. BMC Evol Biol. 2009, 9: 221-10.1186/1471-2148-9-221.PubMedPubMed CentralView ArticleGoogle Scholar
- Langergraber KE, Prüfer K, Rowney C, Boesch C, Crockford C, Fawcett K, Inoue E, Inoue-Muruyama M, Mitani JC, Muller MN, Robbins MM, Schubert G, Stoinski TS, Viola B, Watts D, Wittig RM, Wrangham RW, Zuberbühler K, Pääbo S, Vigilant L: Generation times in wild chimpanzees and gorillas suggest earlier divergence times in great ape and human evolution. Proc Natl Acad Sci U S A. 2012, 109 (39): 15716-15721. 10.1073/pnas.1211740109.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.