Inferring positive selection in humans from genomic data
© Wollstein and Stephan; licensee BioMed Central. 2015
Received: 21 November 2014
Accepted: 23 February 2015
Published: 1 April 2015
Adaptation can be described as an evolutionary process that leads to an adjustment of the phenotypes of a population to their environment. In the classical view, new mutations can introduce novel phenotypic features into a population that leave footprints in the genome after fixation, such as selective sweeps. Alternatively, existing genetic variants may become beneficial after an environmental change and increase in frequency. Although they may not reach fixation, they may cause a shift of the optimum of a phenotypic trait controlled by multiple loci. With the availability of polymorphism data from various organisms, including humans and chimpanzees, it has become possible to detect molecular evidence of adaptation and to estimate the strength and target of positive selection. In this review, we discuss the two competing models of adaptation and suitable approaches for detecting the footprints of positive selection on the molecular level.
KeywordsSelective sweep Evolution Polygenic adaptation
Understanding the genetic architecture and evolution of phenotypes that are present in populations adapting to heterogeneous environments has been a long-standing interest in evolutionary biology [1-3]. This question has been studied by means of quantitative genetics and population genetics. Quantitative genetics provides the methods to describe differences in the distribution of phenotypes, determine their heritability and map relevant regions controlling the phenotype in the genome . In contrast, population genetics provides a framework to describe changes of allele frequencies that are known to be mostly determined by genetic drift  and selection . The latter field produced a profound theory about the process of adaptation on the sequence level , which gave rise to an abundance of population genetic tools that can be applied to genetic data independent of phenotypes [8-14].
In the early years of the study of molecular adaptation, candidate genes with regard to certain phenotypes were conducted . The progress in sequencing and genotyping methods, however, allowed researchers to produce genome-wide single nucleotide polymorphism data in humans and various other organisms [15-20]. This has motivated many genome-wide scans to search for signatures of positive selection [21-30] and quantify fitness effects of various classes of segregating sites [31-34].
Results from genome-wide scans, however, were often inconclusive [8,35,36]. The lack of reproducibility has been attributed to an insufficient power of the tests , the presence of masking signals of positive selection by purifying selection (for example, ) or complex demographic histories (for example, ). Furthermore, the classical model of adaptation in which single new mutations with large effects are favoured by recent positive selection has been questioned [40,41]. And the problem was raised whether evidence for more general models of adaptation (in particular those involving quantitative genetic variation) could be detected on the genomic level [39,40]. This latter issue became particularly interesting in the face of an influx of huge amounts of data from genome-wide association studies [42,43].
In this review, we summarize the population genetic and quantitative genetic models of adaptation and describe the methods to detect the footprints of adaptation in the genome. Furthermore, we provide examples of adaptation in humans that illustrate these theoretical accomplishments.
Population genetic models of adaptation
Genetic adaptation is the result of fitness differences of alleles. Consider the alleles a and A at a bi-allelic locus in diploid organisms as mutant and wild type, respectively. A fitness value may be assigned to each possible genotype aa, aA and AA. Mutations are neutral if the fitness effects are equal (that is, w aa = w aA = w AA ), which is the case for most of the genetic variation observed in humans . In the classical model, positive selection occurs when the derived allele has a higher fitness than the ancestral allele, and negative (or purifying) selection, when the derived allele is detrimental to the organism. Balancing selection occurs in the case of heterozygote advantage and in situations of spatial and temporal heterogeneity of selection. Nucleotide changes in the DNA sequence may have some direct or indirect effect on the phenotype of the individual that generates a fitness advantage or disadvantage and hence are assumed to occur in coding regions of genes or regulatory sequences .
Expected patterns of positive selection in the genome: A beneficial mutation may rise quickly under positive selection. If the beneficial allele is going to fixation, genetic hitchhiking [46,47] results in depletion of variation around the selected site, also termed as selective sweep [46,47]. If the beneficial allele has not yet reached fixation, a sweep is called incomplete, partial, or ongoing. Sweep patterns that arise from a newly introduced mutation or migrant are considered as a ‘hard sweep’. If, however, the beneficial allele arises from standing variation, that is, after an environment change, the pattern of nucleotide polymorphism generated after fixation of the beneficial allele is called ‘soft sweep’ [48-50]. In this latter model, adaptation is not limited to the occurrence of new mutations and can therefore occur more rapidly after an environmental change . The resulting pattern of variation of a soft sweep becomes very similar to that of a hard sweep in case the initial frequency of the beneficial allele is low. This situation may occur if the allele is initially in a mutation-selection balance and becomes positively selected after an environmental shift [46,47].
The genomic signatures of recent adaptation can be measured by means of the site frequency spectrum (SFS), which summarizes the counts of derived variants in a region. Under the action of positive directional selection, the SFS exhibits an excess of both rare and high-frequency derived variants around the selected site that are present in the population at the time of fixation of the beneficial allele [51,52]. The size of the region with depleted variation is expected to be larger when recombination is low and/or selection is strong [47,53], if hitchhiking has started from a selected allele with low frequency. The transient phase, until the beneficial mutation reaches fixation in the population, is inversely proportional to the population size .
Furthermore, in a subdivided population a frequency shift of a beneficial allele may lead to increased genetic differentiation between subpopulations in comparison to a population that has not been subjected to selection . In its extreme, fixed differences between subpopulations may be observed.
The signature of linkage disequilibrium (LD) around the selected site is another characteristic of the hitchhiking process. LD emerges between pairs of sites due to non-random association of alleles. When selection is strong and a sweep is in progress, LD among hitchhiking alleles will strongly increase [56,57], due to limited time for recombination events to occur. However, after the beneficial allele driving hitchhiking has reached an intermediate frequency around 50%, LD between variants across the selected site decreases rapidly and eventually disappears when fixation has occurred. In contrast, the LD between polymorphisms on either side remains high and decreases only slowly. The establishment of the well-known long-range haplotypes in a population  is a consequence of the strong LD around the selected site in the first half of the selective phase (until the beneficial allele reaches intermediate frequency). Therefore, these extended haplotypes can be used to detect incomplete (ongoing) sweeps that are typical for humans [9,36]. The use of LD has the advantage that it is relatively robust against purifying selection .
Once a beneficial mutation has been fixed in a population the signature of linkage disequilibrium decreases and the pattern of polymorphism in the neighbourhood can be restored. The time range to detect these LD signatures of recent adaptation in a single population is rather limited (for example, in the scale of 10,000 years in the case of humans ) and measurable only when adaptation is still ongoing or has only recently ceased. The fixed differences between populations or species remain evident much longer (millions of years, humans compared to chimp, ).
These latter genomic signatures of positive selection, however, may not be unique. It is known that population expansion as well as sudden decreases in population size (bottlenecks) can result in similar genomic patterns, such as an excess of rare and intermediate-frequency derived variants, respectively [58,59]. For example, while human populations were migrating out of Africa, consecutive population bottlenecks followed by population expansion occurred [60,61]. Such a cascade of demographic events is expected to leave patterns in the genome that are very similar to the one of selective sweeps .
Furthermore, population structure can mimic the signature of balancing selection . The inference of the demographic history of a population can in addition be confounded by the genotyping technology leading to single nucleotide polymorphism (SNP) ascertainment bias [64,65]. Choosing SNPs from a too small discovery sample for genotyping can skew the resulting site frequency spectrum toward intermediate frequencies.
Alternative modes of selection may also result in similar signatures as those produced by positive directional selection. In particular, background selection may also lead to a depletion of variation [66,67], yet without causing shifts of low- and high-frequency-derived variants in the SFS. This signature may resemble that of multiple selective sweeps (recurrent sweeps; ) and may result in a lack of high-frequency-derived variants . Selective sweeps may also be difficult to distinguish from recombination hotspots . If recombination is strong, the region of depleted variation may become too small to be recognized. In contrast, a recombination cold spot can generate a pattern of increased LD that is similar as the pattern of a sweep in progress . Furthermore, varying recombination rate on a fine scale may also confound the long-range haplotype signature of sweeps.
A common statistical approach dealing with these difficulties is to derive a likelihood by comparing a statistical null model that includes all the aforementioned non-selective effects to an alternative model that in addition contains positive selection. Many of the confounding factors, however, are difficult to model jointly in a likelihood framework. In an alternative approach, summary statistics are constructed that quantify specific patterns of selective forces and are applied genome-wide. Regions with the strongest signals are considered as outliers. Statistical significance is then assessed by simulating a null model using the coalescent . In the following, we review statistical approaches and their applications taking these confounding effects into account.
Statistical tests to detect deviations from neutrality: Several tests have been developed that make use of the aforementioned signatures of hitchhiking, that is, the reduction of genetic variation, the skew in the frequency spectrum and the pattern of linkage disequilibrium. These tests may be broadly categorized into three classes: (i) tests that use only data from one population, (ii) tests that compare genetic signatures among multiple populations and (iii) comparative tests that use a closely related species as an out-group. The tests can be further classified into model-free and model-based methods. The latter use the neutral theory  to build the null hypothesis and can be applied to compare single candidate regions to a neutral expectation, when full genome data is not available. In contrast, model-free methods try to quantify the characteristic signatures of hitchhiking and are usually applied in an outlier approach to genome-wide data. Regions that show the strongest signals are assumed to be candidates for sweeps [8,10,73].
The most widely used method in the first class of tests is Tajima’s D statistic  that compares the number of segregating sites to an expected value when the population size is assumed to be constant over time (standard neutral model). Large positive values indicate an excess of variation in the tested region that could be due to balancing selection, whereas negative values indicate a depletion of variation due to positive directional selection. The interpretation of the Tajima’s D statistic, however, may be ambiguous as the demographic history of a population needs to be taken into account. Therefore, several more recent developments corrected Tajima’s D statistic, for instance, by including population size changes  or SNP ascertainment bias  that can arise from genotyping technology .
The Fay and Wu’s H test  uses, in addition, data from an out-group species to get information of the ancestral state of a polymorphism and detect selective sweeps by an excess of high-frequency-derived polymorphisms. In contrast, the Fu and Li’s D statistic  takes advantage of low-frequency variation that is enriched in regions that recently underwent genetic hitchhiking. The maximum frequency of derived mutations (MFDM) test  utilizes the MFDM to estimate the presence of an unbalanced tree topology in a given sample that is thought to arise in the adjacency of a locus that is under positive selection due to hitchhiking [46,52]. In line with coalescent theory, the tree topology is independent of changes in population size, which makes the MFDM statistics evidently robust against demographic events, such as bottlenecks or expansions . To obtain good estimates for the MFDM statistics, large sample sizes of at least 42 chromosomes (21 diploids) are necessary  that have to be unaffected by migration, admixture or any hidden population substructure.
A statistic that uses the full site frequency spectrum has been introduced by Kim and Stephan . Here, a composite likelihood ratio (CLR) is calculated by multiplying the probabilities of all polymorphic sites of a genomic region, which makes it possible to estimate the strength and location of a selective sweep. The method returns a likelihood of a complete sweep compared to a population that evolves under standard neutrality, and an estimate about the selection parameter and the target of selection. This test has been further developed by Nielsen et al.  to detect deviations from a background spectrum that includes deviations from neutrality due to demographic history and SNP ascertainment bias under the assumption that the selective sweep has been completed. A demographic model consisting of two epochs of population sizes has been incorporated into the CLR approach by Williamson et al. . Finally, LD has been combined with this composite likelihood framework by Pavlidis et al. , which is reducing the number of false positives. Currently, the most advanced CLR-based test is SweeD  that includes a demographic model with an arbitrary number of instantaneous changes in population size . The power of this test increases with up to a sample size of about 500.
A large fraction of model-free tests are also based on the patterns of LD. Many tests take advantage of the haplotype homozygosity as introduced by Sabeti et al. , which is a measure of genetic diversity regarding multiple polymorphic sites . The decay of the extended haplotype homozygosity (EHH) as calculated step by step from a defined core haplotype was designed as a test for positive selection. This test, however, cannot easily distinguish between complete and incomplete sweeps. Several modifications of the EHH test statistic have been introduced that account for the confounding effect of varying recombination rates. The relative extended haplotype homozygosity (REHH) is defined to be the extended homozygosity of a core haplotype divided by the homozygosity of the remaining core haplotypes combined . The integrated haplotype score (iHS) as proposed by Voight et al.  compares the decay of the ancestral allele against the derived allele. If the derived allele is beneficial, its underlying haplotype will take longer to decay than the ancestral one. While this test cannot be applied to sites that are already fixed, it is useful to detect recent sweeps that are still in progress (that is, incomplete sweeps). As the latter mentioned tests do not compare the observation with a theoretical expectation, they are mostly used in a statistical outlier approach.
The second class of tests compares recently diverging populations under the assumption that adaptation was acting differently on the populations. A test for detecting differentiation in allele frequencies between populations by means of Wright’s fixation index F ST  has been first formulated by Lewontin and Krakauer . This idea has been incorporated into various frequency- and LD-based test statistics. The CLR approach has been extended by Chen et al. . It models population structure by multi-locus allele frequency differentiation between two populations (XP-CLR). However, population size changes and associations between polymorphic sites were not considered in the model. The model-free Rsb measure proposed by Tang et al.  compares the haplotype homozygosity decay at homologous sites between two populations that diverged recently. Similarly, the XP-EHH method  compares the homozygosity decay among different populations. The latter tests take advantage of the assumption that local adaptation increases population differentiation compared to neutrally evolving subpopulations. Another extension of measuring population differentiation between populations on a haplotype level is a method proposed by Fariello et al.  and Ferrer-Admetlla et al.  that has been shown to have more power to detect soft sweeps over SFS-based methods . A combination of class one and class two tests has been proposed in . The composite of multiple signals (CMS) test combines the different priors of detecting extended haplotypes (XP-EHH, iHS), high-frequency-derived alleles (iHS), and polymorphic sites that exhibit population differentiation and results in a score that represents a posterior probability that a certain variant is under selection .
The third class of tests uses the information of an out-group species to detect selection. Most widely used is the dN/dS ratio, also known as Ka/Ks statistic . The basic idea is that the ratio of non-synonymous and synonymous substitutions is close to one under neutrality. The Hudson-Kreitmann-Aquadé test (HKA, ) compares polymorphisms within species by means of Watterson’s estimator  and divergence between species across two or more loci. Under neutrality, they are expected to be identical, which is tested by means of a goodness of fit test. In contrast, the McDonald-Kreitman test compares polymorphism within populations and divergence between species at single loci for two classes of sites (for example, synonymous and non-synonymous sites) .
Quantitative genetic models of adaptation
Quantitative genetic models of adaptation date back to the time before the genetic mechanisms of inheritance were fully discovered [1,94]. Quantitative phenotypes in a population are characterized by a distribution of gradual differences among individuals that are controlled by a multitude of genes. In varying environments, different phenotypes may be favoured. This leads to a change in the population mean phenotype that is known to depend on the additive genetic variation present in the population. When a population deviates from its optimum, mutations are favoured according to their effect size and distance to the optimum. The mean step size of such an adaptive walk has been shown to be approximately exponentially distributed ; that is, alleles with larger effects are favoured when the population resides far from the optimum, whereas alleles with smaller effects are favoured during the adaptive fine-tuning close to the population optimum.
The impact of beneficial mutations in the process of adaptation depends on the mutation rate and population size . In humans, the most non-synonymous mutations have been shown to be neutral (27% to 29% ) or mildly deleterious (30% to 42% [31,33]). In comparison with chimpanzees, 10% to 20% of the fixations appear to be adaptive . However, beneficial mutations that lead to fixation in recent time have been shown to be rare (1% ), so that adaptation from standing variation may be the most important mode of recent adaptation.
In this scenario, classical selective sweeps play only a role if the beneficial alleles are driven to fixation from low frequency by strong selection [40,97]. Instead, small frequency shifts of selected alleles at the quantitative trait loci driving a trait value towards its optimum may occur predominantly.
In case the trait optima of populations are ordered along clines [98,99], effective alleles are expected to change in frequency accordingly . This may be detected by means of the Lewontin and Krakauer test  and other F ST-based statistics (for example, ). To be able to distinguish these adaptive frequency changes from drift, Coop et al.  proposed a model that analyses whether allele frequencies correlate with environmental variables along a population gradient. A test for polygenic adaptation that also incorporates estimates of phenotypic values from genome-wide association data and compares those with environmental variables has been recently introduced by Berg and Coop . However, phenotypic and genotypic data for many populations are required for this test.
Evidence for adaptation in humans
As the migration out of Africa  and the settlement around the world exposed humans to different environmental conditions with regard to temperature, amount of light, humidity, oxygen levels, and agriculture , many adaptations in non-African populations must have occurred in the recent past . In line with this, positive selection has been shown to be a less important determinant in various African populations . The most accepted examples from different genome scans show human adaptations to (i) agriculture , (ii) environmental variables, such as amount of light, temperature, or oxygen levels, and (iii) pathogen resistance [107-109].
The most prominent example of adaptation in humans to agriculture is the ability to digest lactose from milk products in adulthood . Indeed, an extended haplotype homozygosity as a signature for a selective sweep around the LCT gene was observed [22,83]. The activity of the LCT gene is usually reduced in adult mammals . However, the presence of the beneficial mutation provides a selective advantage of about 1.4% to 19% . The most likely explanation for the evolutionary advantage of the mutation is the additional caloric and calcium source it produces because it reduces the risk for diseases related to bone mineralization caused by a lack of vitamin D [110,112]. The frequency of the allele associated with lactase persistence has been shown to decline from Northwest Europe to the southern populations  and the mutation is absent in African populations. In African rural, populations show strong evidence for parallel adaptation to digest lactose from milk products. Other alleles have been associated with lactase persistence  that show similar LD patterns and high selection coefficients of 4% to 9% [113,114].
Skin pigmentation is another example of adaptation to environmental conditions in humans. It is known to be controlled by the amount of eumelanin and pheomelanin that are produced in the melanosomes [115,116]. The dark pigmented skin is assumed to be ancestral, whereas lighter pigmented skin has emerged after the migration out of Africa . Skin colour has long been speculated to evolve under positive selection and is another example for convergent evolution [115,117,118]. Many genes have been shown to be associated with variation in skin colour in different human populations [115,116]. The MC1R gene is a main switch in the production of the lighter pheomelanin and darker eumelanin pigments in the melanosomes . Strong selection for the persistence of the dark pigment has been found in African  and southern European populations . The gene SLC24A5 regulates calcium levels in melanosomes and has been associated with lighter pigmentation in Europeans . In genome-wide scans, it has been shown that SLC24A5 is surrounded by a region of decreased variability and increased LD levels [22,23,83,117] and is substantially differentiated among different populations [23,105,122]. In East Asian populations, another candidate gene, OCA2, has been shown to be subject to positive selection . Furthermore, there are several other candidate genes, such as UGT1A and BNC2 that are associated with skin pigmentation . However, an adaptive signature has not been observed for these genes yet, most likely due to lower effect sizes so that the establishment of a sweep signature and/or frequency changes become too small to be identified.
Human height is a classical quantitative trait that has been studied since the beginning of the last century [124-126] and shows evidence for phenotypic adaptation to different environmental factors, such as temperature (for example, Bergmans rule ), with extreme differences among populations of up to 30 cm . More than 180 loci have been associated with it , with no evidence of selective sweeps so far. Turchin et al.  demonstrated that alleles that contribute to a tall stature are enriched in northern European populations, which is better explained by small selection coefficients of 0.001% to 0.1% than drift. Since human height can be expected to be under stabilizing selection [115,116], the probability of observing selective sweeps is rather low .
Another example of parallel adaptation to low oxygen levels in high altitude has been described in Tibetan, Andean, and Ethiopian populations. Tibetans and Ethiopians adapted differently to the low oxygen levels compared to Andeans . Andeans show an increased haemoglobin blood concentration that elevates the oxygen transport in blood, whereas Tibetans and Ethiopians exhibit an increased lung capacity and breathing rate . The EPAS1 and EGLN1 genes show strong signatures of selective sweeps in Tibetans; that is, an increased differentiation in allele frequency compared to East Asian populations and an increased LD [131-133]. Variants of the EPAS1 and EGLN1 genes have been associated with haemoglobin concentration levels in blood [134,135]. It has been shown that the EPAS1 gene has likely been introgressed from an archaic human, the Denisovans into the Tibetans . In the Andean population, different genes (NOS2A and PRKAA1) have been identified as targets of adaptation [131,137]. Ethiopian high-altitude populations that have a similar phenotype as Tibetan populations, also show a different set of genes (CBARA1, VAV3, ARNT2 and THRB), with evidence for positive selection . Variants associated with haemoglobin variation in Tibetans do not overlap with variants associated in Ethiopians .
Classical sweeps have been shown to be rare in humans [13,96,105] and, if they exist, they occur around loci with large-effects alleles. As selective sweeps are rare in humans (in contrast to species with large effective sizes such as Drosophila), the emphasis of human population genetics in the near future must be to identify adaptive signatures for polygenic phenotypes. There is an urgent need for more theoretical modelling and better statistical methods to analyse the evolution of polygenic traits for populations of varying environments and demographies.
composite likelihood ratio
composite of multiple signals
extended haplotype homozygosity
integrated haplotype score
maximum frequency of derived mutations
relative extended haplotype homozygosity
site frequency spectrum
single nucleotide polymorphism
cross population composite likelihood ratio
cross population extended haplotype homozygosity
We would like to thank Oscar Lao for helpful comments on the manuscript. AW was funded by the Volkswagen Stiftung (grant 86042) and WS by DFG grant 325/14-1 from the priority program 1590.
- Orr HA. The genetic theory of adaptation: a brief history. Nat Rev Genet. 2005;6:119–27.PubMedGoogle Scholar
- Mackay TF. The genetic architecture of quantitative traits. Annu Rev Genet. 2001;35:303–39.PubMedGoogle Scholar
- Fu W, O’Connor TD, Akey JM. Genetic architecture of quantitative traits and complex diseases. Curr Opin Genet Dev. 2013;23:678–83.PubMedGoogle Scholar
- Hill WG. Understanding and using quantitative genetic variation. Philos Trans R Soc Lond B Biol Sci. 2010;365:73–85.PubMed CentralPubMedGoogle Scholar
- Kimura M. The neutral theory of molecular evolution. Cambridge: Cambridge University Press; 1983.Google Scholar
- Orr HA. Fitness and its role in evolutionary genetics. Nat Rev Genet. 2009;10:531–9.PubMed CentralPubMedGoogle Scholar
- Crow JF, Kimura M. An introduction to population genetics theory. 1970.Google Scholar
- Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39:197–218.PubMedGoogle Scholar
- Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG. Recent and ongoing selection in the human genome. Nat Rev Genet. 2007;8:857–68.PubMed CentralPubMedGoogle Scholar
- Pavlidis P, Hutter S, Stephan W. A population genomic approach to map recent positive selection in model species. Mol Ecol. 2008;17:3585–98.PubMedGoogle Scholar
- Stephan W. Detecting strong positive selection in the genome. Mol Ecol Resour. 2010;10:863–72.PubMedGoogle Scholar
- Vitti JJ, Grossman SR, Sabeti PC. Detecting natural selection in genomic data. Annu Rev Genet. 2013;47:97–120.PubMedGoogle Scholar
- Fu W, Akey JM. Selection and adaptation in the human genome. Annu Rev Genomics Hum Genet. 2013;14:467–89.PubMedGoogle Scholar
- Jeong C, Di Rienzo A. Adaptations to local environments in modern human populations. Curr Opin Genet Dev. 2014;29C:1–8.Google Scholar
- International HapMap Consortium. The International HapMap Project. Nature. 2003;426:789–96.Google Scholar
- Chimpanzee Sequencing and Analysis Consortium. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005;437:69–87.Google Scholar
- International HapMap Consortium. A haplotype map of the human genome. Nature. 2005;437:1299–320.Google Scholar
- 1000 Genomes Project Consortium, Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, et al. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–73.Google Scholar
- Cann HM, de Toma C, Cazes L, Legrand M-F, Morel V, Piouffre L, et al. A human genome diversity cell line panel. Science. 2002;296:261–2.PubMedGoogle Scholar
- Hinds DA, Stuve LL, Nilsen GB, Halperin E, Eskin E, Ballinger DG, et al. Whole-genome patterns of common DNA variation in three human populations. Science. 2005;307:1072–9.PubMedGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.PubMedGoogle Scholar
- Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:e72.PubMed CentralPubMedGoogle Scholar
- Lao O, de Gruijter JM, van Duijn K, Navarro A, Kayser M. Signatures of positive selection in genes associated with human skin pigmentation as revealed from analyses of single nucleotide polymorphisms. Ann Hum Genet. 2007;71(Pt 3):354–69.PubMedGoogle Scholar
- Carlson CS, Thomas DJ, Eberle MA, Swanson JE, Livingston RJ, Rieder MJ, et al. Genomic regions exhibiting positive selection identified from dense genotype data. Genome Res. 2005;15:1553–65.PubMed CentralPubMedGoogle Scholar
- Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, et al. Natural selection on protein-coding genes in the human genome. Nature. 2005;437:1153–7.PubMedGoogle Scholar
- Akey JM, Zhang G, Zhang K, Jin L, Shriver MD. Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002;12:1805–14.PubMed CentralPubMedGoogle Scholar
- Shriver MD, Kennedy GC, Parra EJ, Lawson HA, Sonpar V, Huang J, et al. The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs. Hum Genomics. 2004;1:274–86.PubMed CentralPubMedGoogle 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:980–9.PubMed CentralPubMedGoogle Scholar
- Tang K, Thornton KR, Stoneking M. A new approach for using genome scans to detect recent positive selection in the human genome. PLoS Biol. 2007;5:e171.PubMed CentralPubMedGoogle Scholar
- Kimura R, Ohashi J, Matsumura Y, Nakazawa M, Inaoka T, Ohtsuka R, et al. Gene flow and natural selection in oceanic human populations inferred from genome-wide SNP typing. Mol Biol Evol. 2008;25:1750–61.PubMedGoogle Scholar
- Williamson SH, Hernandez R, Fledel-Alon A, Zhu L, Nielsen R, Bustamante CD. Simultaneous inference of selection and population growth from patterns of variation in the human genome. Proc Natl Acad Sci U S A. 2005;102:7882–7.PubMed CentralPubMedGoogle Scholar
- Eyre-Walker A, Keightley PD. High genomic deleterious mutation rates in hominids. Nature. 1999;397:344–7.PubMedGoogle Scholar
- Boyko AR, Williamson SH, Indap AR, Degenhardt JD, Hernandez RD, Lohmueller KE, et al. Assessing the evolutionary impact of amino acid mutations in the human genome. PLoS Genet. 2008;4:e1000083.PubMed CentralPubMedGoogle Scholar
- Eyre-Walker A, Woolfit M, Phelps T. The distribution of fitness effects of new deleterious amino acid mutations in humans. Genetics. 2006;173:891–900.PubMed CentralPubMedGoogle Scholar
- Akey JM. Constructing genomic maps of positive selection in humans: where do we go from here? Genome Res. 2009;19:711–22.PubMed CentralPubMedGoogle Scholar
- Sabeti PC, Schaffner SF, Fry B, Lohmueller J, Varilly P, Shamovsky O, et al. Positive natural selection in the human lineage. Science. 2006;312:1614–20.PubMedGoogle Scholar
- Zhai W, Nielsen R, Slatkin M. An investigation of the statistical power of neutrality tests based on comparative and population genetic data. Mol Biol Evol. 2009;26:273–83.PubMed CentralPubMedGoogle Scholar
- Enard D, Messer PW, Petrov DA. Genome-wide signals of positive selection in human evolution. Genome Res. 2014;24:885–95.PubMed CentralPubMedGoogle Scholar
- Teshima KM, Coop G, Przeworski M. How reliable are empirical genomic scans for selective sweeps? Genome Res. 2006;16:702–12.PubMed CentralPubMedGoogle Scholar
- Pritchard JK, Pickrell JK, Coop G. The genetics of human adaptation: hard sweeps, soft sweeps, and polygenic adaptation. Curr Biol. 2010;20:R208–15.PubMed CentralPubMedGoogle Scholar
- Pritchard JK, Di Rienzo A. Adaptation - not by sweeps alone. Nat Rev Genet. 2010;11:665–7.PubMedGoogle Scholar
- McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JPA, et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008;9:356–69.PubMedGoogle Scholar
- Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, et al. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009;106:9362–7.PubMed CentralPubMedGoogle Scholar
- Hellmann I, Ebersberger I, Ptak SE, Pääbo S, Przeworski M. A neutral explanation for the correlation of diversity with recombination rates in humans. Am J Hum Genet. 2003;72:1527–35.PubMed CentralPubMedGoogle Scholar
- Zhen Y, Andolfatto P. Methods to detect selection on noncoding DNA. Methods Mol Biol. 2012;856:141–59.PubMed CentralPubMedGoogle Scholar
- Kaplan NL, Hudson RR, Langley CH. The “hitchhiking effect” revisited. Genetics. 1989;123:887–99.PubMed CentralPubMedGoogle Scholar
- Stephan W, Wiehe THE, Lenz MW. The effect of strongly selected substitutions on neutral polymorphism: analytical results based on diffusion theory. Theor Popul Biol. 1992;41:237–54.Google Scholar
- Innan H, Kim Y. Pattern of polymorphism after strong artificial selection in a domestication event. Proc Natl Acad Sci U S A. 2004;101:10667–72.PubMed CentralPubMedGoogle Scholar
- Hermisson J, Pennings PS. Soft sweeps: molecular population genetics of adaptation from standing genetic variation. Genetics. 2005;169:2335–52.PubMed CentralPubMedGoogle Scholar
- Przeworski M, Coop G, Wall JD. The signature of positive selection on standing genetic variation. Evolution. 2005;59:2312–23.PubMedGoogle Scholar
- Braverman JM, Hudson RR, Kaplan NL, Langley CH, Stephan W. The hitchhiking effect on the site frequency spectrum of DNA polymorphisms. Genetics. 1995;140:783–96.PubMed CentralPubMedGoogle Scholar
- Fay JC, Wu CI. Hitchhiking under positive Darwinian selection. Genetics. 2000;155:1405–13.PubMed CentralPubMedGoogle Scholar
- Maynard Smith J, Haigh J. The hitch-hiking effect of a favourable gene. Genet Res. 1974;23:23–35.Google Scholar
- Kim Y, Stephan W. Detecting a local signature of genetic hitchhiking along a recombining chromosome. Genetics. 2002;160:765–77.PubMed CentralPubMedGoogle Scholar
- Kim Y. Stochastic patterns of polymorphism after a selective sweep over a subdivided population. Genet Res. 2013;95:57–67.Google Scholar
- Kim Y, Nielsen R. Linkage disequilibrium as a signature of selective sweeps. Genetics. 2004;167:1513–24.PubMed CentralPubMedGoogle Scholar
- Stephan W, Song YS, Langley CH. The hitchhiking effect on linkage disequilibrium between linked neutral loci. Genetics. 2006;172:2647–63.PubMed CentralPubMedGoogle Scholar
- Tajima F. The effect of change in population size on DNA polymorphism. Genetics. 1989;123:597–601.PubMed CentralPubMedGoogle Scholar
- Evans SN, Shvets Y, Slatkin M. Non-equilibrium theory of the allele frequency spectrum. Theor Popul Biol. 2007;71:109–19.PubMedGoogle Scholar
- DeGiorgio M, Jakobsson M, Rosenberg NA. Out of Africa: modern human origins special feature: explaining worldwide patterns of human genetic variation using a coalescent-based serial founder model of migration outward from Africa. Proc Natl Acad Sci U S A. 2009;106:16057–62.PubMed CentralPubMedGoogle Scholar
- Veeramah KR, Hammer MF. The impact of whole-genome sequencing on the reconstruction of human population history. Nat Rev Genet. 2014;15:149–62.PubMedGoogle Scholar
- Sousa V, Peischl S, Excoffier L. Impact of range expansions on current human genomic diversity. Curr Opin Genet Dev. 2014;29:22–30.PubMedGoogle Scholar
- Muirhead CA. Consequences of population structure on genes under balancing selection. Evolution. 2001;55:1532–41.PubMedGoogle Scholar
- Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R. Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005;15:1496–502.PubMed CentralPubMedGoogle Scholar
- Nielsen R. Population genetic analysis of ascertained SNP data. Hum Genomics. 2004;1:218–24.PubMed CentralPubMedGoogle Scholar
- Stephan W. Genetic hitchhiking versus background selection: the controversy and its implications. Philos Trans R Soc Lond B Biol Sci. 2010;365:1245–53.PubMed CentralPubMedGoogle Scholar
- Charlesworth B, Morgan MT, Charlesworth D. The effect of deleterious mutations on neutral molecular variation. Genetics. 1993;134:1289–303.PubMed CentralPubMedGoogle Scholar
- Stephan W. An improved method for estimating the rate of fixation of favorable mutations based on DNA polymorphism data. Mol Biol Evol. 1995;12:959–62.PubMedGoogle Scholar
- Kim Y. Allele frequency distribution under recurrent selective sweeps. Genetics. 2006;172:1967–78.PubMed CentralPubMedGoogle Scholar
- Reed FA, Tishkoff SA. Positive selection can create false hotspots of recombination. Genetics. 2006;172:2011–4.PubMed CentralPubMedGoogle Scholar
- McVean G. The structure of linkage disequilibrium around a selective sweep. Genetics. 2007;175:1395–406.PubMed CentralPubMedGoogle Scholar
- Hudson RR. Others: gene genealogies and the coalescent process. Oxford Surveys Evol Biol. 1990;7:44.Google Scholar
- Scheinfeldt LB, Tishkoff SA. Recent human adaptation: genomic approaches, interpretation and insights. Nat Rev Genet. 2013;14:692–702.PubMedGoogle Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Zivković D, Wiehe T. Second-order moments of segregating sites under variable population size. Genetics. 2008;180:341–57.PubMed CentralPubMedGoogle Scholar
- Ramírez-Soriano A, Nielsen R. Correcting estimators of theta and Tajima’s D for ascertainment biases caused by the single-nucleotide polymorphism discovery process. Genetics. 2009;181:701–10.PubMed CentralPubMedGoogle Scholar
- Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.PubMed CentralPubMedGoogle 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:365–75.PubMedGoogle Scholar
- Pavlidis P, Jensen JD, Stephan W. Searching for footprints of positive selection in whole-genome SNP data from nonequilibrium populations. Genetics. 2010;185:907–22.PubMed CentralPubMedGoogle Scholar
- Pavlidis P, Živkovic D, Stamatakis A, Alachiotis N. SweeD: likelihood-based detection of selective sweeps in thousands of genomes. Mol Biol Evol. 2013;30:2224–34.PubMed CentralPubMedGoogle Scholar
- Živković D, Stephan W. Analytical results on the neutral non-equilibrium allele frequency spectrum based on diffusion theory. Theor Popul Biol. 2011;79:184–91.PubMedGoogle Scholar
- Sabatti C, Risch N. Homozygosity and linkage disequilibrium. Genetics. 2002;160:1707–19.PubMed CentralPubMedGoogle Scholar
- Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–8.PubMed CentralPubMedGoogle Scholar
- Wright S. Genetical structure of populations. Nature. 1950;166:247–9.PubMedGoogle Scholar
- Lewontin RC, Krakauer J. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics. 1973;74:175–95.PubMed CentralPubMedGoogle Scholar
- Chen H, Patterson N, Reich D. Population differentiation as a test for selective sweeps. Genome Res. 2010;20:393–402.PubMed CentralPubMedGoogle Scholar
- Fariello MI, Boitard S, Naya H, SanCristobal M, Servin B. Detecting signatures of selection through haplotype differentiation among hierarchically structured populations. Genetics. 2013;193:929–41.PubMed CentralPubMedGoogle Scholar
- Ferrer-Admetlla A, Liang M, Korneliussen T, Nielsen R. On detecting incomplete soft or hard selective sweeps using haplotype structure. Mol Biol Evol. 2014;31:1275–91.PubMed CentralPubMedGoogle Scholar
- Grossman SR, Shlyakhter I, Shylakhter I, Karlsson EK, Byrne EH, Morales S, et al. A composite of multiple signals distinguishes causal variants in regions of positive selection. Science. 2010;327:883–6.PubMedGoogle Scholar
- Miyata T, Yasunaga T. Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application. J Mol Evol. 1980;16:23–36.PubMedGoogle Scholar
- Hudson RR, Kreitman M, Aguadé M. A test of neutral molecular evolution based on nucleotide data. Genetics. 1987;116:153–9.PubMed CentralPubMedGoogle Scholar
- Watterson GA. On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975;7:256–76.PubMedGoogle Scholar
- McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991;351:652–4.PubMedGoogle Scholar
- Barton NH. The geometry of adaptation. Nature. 1998;395:751–2.PubMedGoogle Scholar
- Elena SF, Wilke CO, Ofria C, Lenski RE. Effects of population size and mutation rate on the evolution of mutational robustness. Evolution. 2007;61:666–74.PubMedGoogle Scholar
- Hernandez RD, Kelley JL, Elyashiv E, Melton SC, Auton A, McVean G. Classic selective sweeps were rare in recent human evolution. Science. 2011;331:920–4.PubMed CentralPubMedGoogle Scholar
- Wollstein A, Stephan W. Adaptive fixation in two-locus models of stabilizing selection and genetic drift. Genetics. 2014;198(2):685–97.PubMedGoogle Scholar
- Haldane JBS. The theory of a cline. J Genet. 1948;48:277–84.PubMedGoogle Scholar
- Barton NH. Gene flow past a cline. Heredity. 1979;43:333–9.Google Scholar
- Rosenberg NA, Li LM, Ward R, Pritchard JK. Informativeness of genetic markers for inference of ancestry. Am J Hum Genet. 2003;73:1402–22.PubMed CentralPubMedGoogle Scholar
- Coop G, Witonsky D, Di Rienzo A, Pritchard JK. Using environmental correlations to identify loci underlying local adaptation. Genetics. 2010;185:1411–23.PubMed CentralPubMedGoogle Scholar
- Berg JJ, Coop G. A population genetic signal of polygenic adaptation. PLoS Genet. 2014;10:e1004412.PubMed CentralPubMedGoogle Scholar
- Stoneking M, Krause J. Learning about human population history from ancient and modern genomes. Nat Rev Genet. 2011;12:603–14.PubMedGoogle Scholar
- Beja-Pereira A, Caramelli D, Lalueza-Fox C, Vernesi C, Ferrand N, Casoli A, et al. The origin of European cattle: evidence from modern and ancient DNA. Proc Natl Acad Sci U S A. 2006;103:8113–8.PubMed CentralPubMedGoogle Scholar
- Coop G, Pickrell JK, Novembre J, Kudaravalli S, Li J, Absher D, et al. The role of geography in human adaptation. PLoS Genet. 2009;5:e1000500.PubMed CentralPubMedGoogle Scholar
- Granka JM, Henn BM, Gignoux CR, Kidd JM, Bustamante CD, Feldman MW. Limited evidence for classic selective sweeps in African populations. Genetics. 2012;192:1049–64.PubMed CentralPubMedGoogle Scholar
- Daub JT, Hofer T, Cutivet E, Dupanloup I, Quintana-Murci L, Robinson-Rechavi M, et al. Evidence for polygenic adaptation to pathogens in the human genome. Mol Biol Evol. 2013;30:1544–58.PubMedGoogle Scholar
- Siddle KJ, Quintana-Murci L. The Red Queen’s long race: human adaptation to pathogen pressure. Curr Opin Genet Dev. 2014;29:31–8.PubMedGoogle Scholar
- Barreiro LB, Quintana-Murci L. From evolutionary genetics to human immunology: how selection shapes host defence genes. Nat Rev Genet. 2010;11:17–30.PubMedGoogle Scholar
- Swallow DM. Genetics of lactase persistence and lactose intolerance. Annu Rev Genet. 2003;37:197–219.PubMedGoogle Scholar
- Bersaglieri T, Sabeti PC, Patterson N, Vanderploeg T, Schaffner SF, Drake JA, et al. Genetic signatures of strong recent positive selection at the lactase gene. Am J Hum Genet. 2004;74:1111–20.PubMed CentralPubMedGoogle Scholar
- Gerbault P, Liebert A, Itan Y, Powell A, Currat M, Burger J, et al. Evolution of lactase persistence: an example of human niche construction. Philos Trans R Soc Lond B Biol Sci. 2011;366:863–77.PubMed CentralPubMedGoogle Scholar
- Tishkoff SA, Reed FA, Ranciaro A, Voight BF, Babbitt CC, Silverman JS, et al. Convergent adaptation of human lactase persistence in Africa and Europe. Nat Genet. 2007;39:31–40.PubMed CentralPubMedGoogle Scholar
- Breton G, Schlebusch CM, Lombard M, Sjödin P, Soodyall H, Jakobsson M. Lactase persistence alleles reveal partial East African ancestry of Southern African Khoe pastoralists. Curr Biol. 2014;24:852–8.PubMedGoogle Scholar
- Sturm RA. Molecular genetics of human pigmentation diversity. Hum Mol Genet. 2009;18:R9–17.PubMedGoogle Scholar
- Liu F, Wen B, Kayser M. Colorful DNA polymorphisms in humans. Semin Cell Dev Biol. 2013;24:562–75.PubMedGoogle Scholar
- Norton HL, Kittles RA, Parra E, McKeigue P, Mao X, Cheng K, et al. Genetic evidence for the convergent evolution of light skin in Europeans and East Asians. Mol Biol Evol. 2007;24:710–22.PubMedGoogle Scholar
- Sturm RA, Duffy DL. Human pigmentation genes under environmental selection. Genome Biol. 2012;13:248.PubMed CentralPubMedGoogle Scholar
- Harding RM, Healy E, Ray AJ, Ellis NS, Flanagan N, Todd C, et al. Evidence for variable selective pressures at MC1R. Am J Hum Genet. 2000;66:1351–61.PubMed CentralPubMedGoogle Scholar
- Martínez-Cadenas C, López S, Ribas G, Flores C, García O, Sevilla A, et al. Simultaneous purifying selection on the ancestral MC1R allele and positive selection on the melanoma-risk allele V60L in south Europeans. Mol Biol Evol. 2013;30:2654–65.PubMedGoogle Scholar
- Lamason RL, Mohideen M-APK, Mest JR, Wong AC, Norton HL, Aros MC, et al. SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science. 2005;310:1782–6.PubMedGoogle Scholar
- Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009;19:826–37.PubMed CentralPubMedGoogle Scholar
- Jacobs LC, Wollstein A, Lao O, Hofman A, Klaver CC, Uitterlinden AG, et al. Comprehensive candidate gene study highlights UGT1A and BNC2 as new genes determining continuous skin color variation in Europeans. Hum Genet. 2013;132:147–58.PubMedGoogle Scholar
- Galton F. Hereditary stature. Nature. 1886;33:317.Google Scholar
- Blakeslee AF. Corn and men. The interacting influence of heredity and environment - movements for betterment of men, or corn, or any other living thing, one-sided unless they take both fact. J Hered. 1914;5(11):511–8.Google Scholar
- Fisher RA. The correlation between relatives on the supposition of Mendelian inheritance. Trans R Soc Edinb. 1919;52:399–433.Google Scholar
- Mendizabal I, Marigorta UM, Lao O, Comas D. Adaptive evolution of loci covarying with the human African Pygmy phenotype. Hum Genet. 2012;131:1305–17.PubMed CentralPubMedGoogle Scholar
- Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010;467:832–8.PubMed CentralPubMedGoogle Scholar
- Turchin MC, Chiang CWK, Palmer CD, Sankararaman S, Reich D, Genetic Investigation of Anthropometric Traits (GIANT) Consortium, et al. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nat Genet. 2012;44:1015–9.PubMed CentralPubMedGoogle Scholar
- Beall CM. Two routes to functional adaptation: Tibetan and Andean high-altitude natives. Proc Natl Acad Sci U S A. 2007;104 Suppl 1:8655–60.PubMed CentralPubMedGoogle Scholar
- Bigham A, Bauchet M, Pinto D, Mao X, Akey JM, Mei R, et al. Identifying signatures of natural selection in Tibetan and Andean populations using dense genome scan data. PLoS Genet. 2010;6:e1001116.PubMed CentralPubMedGoogle Scholar
- Bigham AW, Wilson MJ, Julian CG, Kiyamu M, Vargas E, Leon-Velarde F, et al. Andean and Tibetan patterns of adaptation to high altitude. Am J Hum Biol. 2013;25:190–7.PubMedGoogle Scholar
- Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, et al. Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010;329:75–8.PubMed CentralPubMedGoogle Scholar
- Beall CM, Cavalleri GL, Deng L, Elston RC, Gao Y, Knight J, et al. Natural selection on EPAS1 (HIF2α) associated with low hemoglobin concentration in Tibetan highlanders. Proc Natl Acad Sci U S A. 2010;107:11459–64.PubMed CentralPubMedGoogle Scholar
- Simonson TS, Yang Y, Huff CD, Yun H, Qin G, Witherspoon DJ, et al. Genetic evidence for high-altitude adaptation in Tibet. Science. 2010;329:72–5.PubMedGoogle Scholar
- Huerta-Sánchez E, Jin X, Asan, Bianba Z, Peter BM, Vinckenbosch N, et al. Altitude adaptation in Tibetans caused by introgression of Denisovan-like DNA. Nature. 2014;512:194–7.PubMed CentralPubMedGoogle Scholar
- Bigham AW, Mao X, Mei R, Brutsaert T, Wilson MJ, Julian CG, et al. Identifying positive selection candidate loci for high-altitude adaptation in Andean populations. Hum Genomics. 2009;4:79–90.PubMed CentralPubMedGoogle Scholar
- Scheinfeldt LB, Soi S, Thompson S, Ranciaro A, Woldemeskel D, Beggs W, et al. Genetic adaptation to high altitude in the Ethiopian highlands. Genome Biol. 2012;13:R1.PubMed CentralPubMedGoogle Scholar
- Alkorta-Aranburu G, Beall CM, Witonsky DB, Gebremedhin A, Pritchard JK, Di Rienzo A. The genetic architecture of adaptations to high altitude in Ethiopia. PLoS Genet. 2012;8:e1003110.PubMed CentralPubMedGoogle Scholar
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.