Something old, something borrowed: admixture and adaptation in human evolution

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

Bayesian hybrid index and genomic cline estimation with the R package gghybrid

Admixture, the interbreeding of individuals from differentiated source populations, is now known to be a widespread phenomenon. Genomic studies of natural hybridisation can help to answer many questions on the impacts of admixture on adaptive evolution, reproductive isolation, and speciation. When a large variety of admixture proportions between two source populations exist, both geographic and genomic cline analysis are suitable methods for inferring biased, restricted or excessive gene flow at individual loci into the foreign genomic background, providing evidence for reproductive isolation, selection across an environmental transition, balancing selection, and adaptive introgression. Genomic cline analysis replaces geographic location with genome-wide hybrid index and is therefore useable in circumstances that violate geographic cline assumptions. Here, I introduce gghybrid, an R package for simple and flexible Bayesian estimation of Buerkle's hybrid index and Fitzpatrick's logit-logistic genomic clines using bi-allelic data, suitable for both small and large datasets. gghybrid allows any ploidy and uses Structure input file format. It has separate functions for hybrid index and cline estimation, treating each individual and locus respectively as an independent analysis, making it highly parallelisable. Admixture proportions from other software can alternatively be used in cline analysis, alongside parental allele frequencies. Parameters can be fixed and samples pooled for statistical model comparison with AIC or waic. Here, I describe the functions, pipeline, and statistical properties of gghybrid. Simulations reveal that model comparison with waic is preferred, and use of Bayesian posterior distributions and p values to select candidate non-null loci is problematic and should be avoided.

Ancestry and self-reported race in Brazilian breast cancer women

OBJECTIVE

This study aimed to evaluate the association between self-reported race/color and ancestry in Brazilian patients with breast cancer.

METHODS

This was an observational, transversal, epidemiological study, evaluating race and ancestry in 1,127 patients with breast cancer. For genetic ancestry, a 46-AIM-INDEL panel was used. The ancestral profile was evaluated with the Structure v.2.3.3 software. Descriptive statistics were performed. To assess differences between race and ancestry, an analysis of variance with Bonferoni adjustment was used.

RESULTS

The race distribution was 77.7% white, 17.6% brown, 4.1% black, 0.4% yellow, and 0.3% cafuse. The African ancestry proportion was significantly (p<0.001) more evident in black [0.63±0.21 (0.17-0.96)], followed by brown [0.25±0.16 (0.02-0.70)], and less frequent in white skin color. The European ancestry proportion was significantly (p<0.001) higher in white [0.72±0.17 (0.02-0.97)], followed by brown [0.57±0.19 (0.12-0.92)], yellow [0.27±0.31 (0.12-0.620], and black [0.24±0.19 (0.02-0.72)]. The Asiatic ancestry proportion is significantly (p<0.001) higher in yellow [0.48±0.51 (0.04-0.93)]. The Amerindian ancestry proportion frequency was the least frequent in all groups, and cafuse patients did not express differences between all race groups. The brown race group presented differences in African and European ancestry.

CONCLUSION

Although we found many similarities between white European ancestry, black African ancestry, and yellow Asian ancestry, there is great miscegenation between patients. Although they can be labeled as having one race, they do present many ancestral genes that would allow their inclusion in another race group.

Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans

Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.

Human evolution: When admixture met selection

Admixture has been a major force during human evolution. Two new studies using ancient DNA now show how two key admixture events in the evolutionary history of Europeans altered their adaptive trajectories and facilitated rapid evolution.

Ecology and age, but not genetic ancestry, predict fetal loss in a wild baboon hybrid zone

OBJECTIVES

Pregnancy failure represents a major fitness cost for any mammal, particularly those with slow life histories such as primates. Here, we quantified the risk of fetal loss in wild hybrid baboons, including genetic, ecological, and demographic sources of variance. We were particularly interested in testing the hypothesis that hybridization increases fetal loss rates. Such an effect would help explain how baboons may maintain genetic and phenotypic integrity despite interspecific gene flow.

MATERIALS AND METHODS

We analyzed outcomes for 1020 pregnancies observed over 46 years in a natural yellow baboon-anubis baboon hybrid zone. Fetal losses and live births were scored based on records of female reproductive state and the appearance of live neonates. We modeled the probability of fetal loss as a function of a female's genetic ancestry (the proportion of her genome estimated to be descended from anubis [vs. yellow] ancestors), age, number of previous fetal losses, dominance rank, group size, climate, and habitat quality using binomial mixed effects models.

RESULTS

Female genetic ancestry did not predict fetal loss. Instead, the risk of fetal loss is elevated for very young and very old females. Fetal loss is most robustly predicted by ecological factors, including poor habitat quality prior to a home range shift and extreme heat during pregnancy.

DISCUSSION

Our results suggest that gene flow between yellow and anubis baboons is not impeded by an increased risk of fetal loss for hybrid females. Instead, ecological conditions and female age are key determinants of this component of female reproductive success.

Sequential hybridization may have facilitated ecological transitions in the Southwestern pinyon pine syngameon

Multispecies interbreeding networks, or syngameons, have been increasingly reported in natural systems. However, the formation, structure, and maintenance of syngameons have received little attention. Through gene flow, syngameons can increase genetic diversity, facilitate the colonization of new environments, and contribute to hybrid speciation. In this study, we evaluated the history, patterns, and consequences of hybridization in a pinyon pine syngameon using morphological and genomic data to assess genetic structure, demographic history, and geographic and climatic data to determine niche differentiation. We demonstrated that Pinus edulis, a dominant species in the Southwestern US and a barometer of climate change, is a core participant in the syngameon, involved in the formation of two drought-adapted hybrid lineages including the parapatric and taxonomically controversial fallax-type. We found that species remain morphologically and genetically distinct at range cores, maintaining species boundaries while undergoing extensive gene flow in areas of sympatry at range peripheries. Our study shows that sequential hybridization may have caused relatively rapid speciation and facilitated the colonization of different niches, resulting in the rapid formation of two new lineages. Participation in the syngameon may allow adaptive traits to be introgressed across species barriers and provide the changes needed to survive future climate scenarios.

A Comprehensive Overview of Baboon Phylogenetic History

Baboons (genus Papio) are an intriguing study system to investigate complex evolutionary processes and the evolution of social systems. An increasing number of studies over the last 20 years has shown that considerable incongruences exist between phylogenies based on morphology, mitochondrial, and nuclear sequence data of modern baboons, and hybridization and introgression have been suggested as the main drivers of these patterns. Baboons, therefore, present an excellent opportunity to study these phenomena and their impact on speciation. Advances both in geographic and genomic coverage provide increasing details on the complexity of the phylogeography of baboons. Here, we compile the georeferenced genetic data of baboons and review the current knowledge on baboon phylogeny, discuss the evolutionary processes that may have shaped the patterns that we observe today, and propose future avenues for research.

Moving beyond the adaptationist paradigm for human evolution, and why it matters

The Journal of Human Evolution (JHE) was founded 50 years ago when much of the foundation for how we think about human evolution was in place or being put in place, providing the main framework for how we consider our origins today. Here, we will explore historical developments, including early JHE outputs, as they relate to our understanding of the relationship between phenotypic variation and evolutionary process, and use that as a springboard for considering our current understanding of these links as applied to human evolution. We will focus specifically on how the study of variation itself has shifted us away from taxonomic and adaptationist perspectives toward a richer understanding of the processes shaping human evolutionary history, using literature searches and specific test cases to highlight this. We argue that natural selection, gene exchange, genetic drift, and mutation should not be considered individually when considering the production of hominin diversity. In this context, we offer suggestions for future research directions and reflect on this more complex understanding of human evolution and its broader relevance to society. Finally, we end by considering authorship demographics and practices in the last 50 years within JHE and how a shift in these demographics has the potential to reshape the science of human evolution going forward.

Neandertal introgression partitions the genetic landscape of neuropsychiatric disorders and associated behavioral phenotypes

Despite advances in identifying the genetic basis of psychiatric and neurological disorders, fundamental questions about their evolutionary origins remain elusive. Here, introgressed variants from archaic humans such as Neandertals can serve as an intriguing research paradigm. We compared the number of associations for Neandertal variants to the number of associations of frequency-matched non-archaic variants with regard to human CNS disorders (neurological and psychiatric), nervous system drug prescriptions (as a proxy for disease), and related, non-disease phenotypes in the UK biobank (UKBB). While no enrichment for Neandertal genetic variants were observed in the UKBB for psychiatric or neurological disease categories, we found significant associations with certain behavioral phenotypes including pain, chronotype/sleep, smoking and alcohol consumption. In some instances, the enrichment signal was driven by Neandertal variants that represented the strongest association genome-wide. SNPs within a Neandertal haplotype that was associated with smoking in the UKBB could be replicated in four independent genomics datasets.Our data suggest that evolutionary processes in recent human evolution like admixture with Neandertals significantly contribute to behavioral phenotypes but not psychiatric and neurological diseases. These findings help to link genetic variants in a population to putative past beneficial effects, which likely only indirectly contribute to pathology in modern day humans.

Hybrid fitness effects modify fixation probabilities of introgressed alleles

Hybridization is a common occurrence in natural populations, and introgression is a major source of genetic variation. Despite the evolutionary importance of adaptive introgression, classical population genetics theory does not take into account hybrid fitness effects. Specifically, heterosis (i.e. hybrid vigor) and Dobzhansky–Muller incompatibilities influence the fates of introgressed alleles. Here, we explicitly account for polygenic, unlinked hybrid fitness effects when tracking a rare introgressed marker allele. These hybrid fitness effects quickly decay over time due to repeated backcrossing, enabling a separation-of-timescales approach. Using diffusion and branching process theory in combination with computer simulations, we formalize the intuition behind how hybrid fitness effects affect introgressed alleles. We find that hybrid fitness effects can significantly hinder or boost the fixation probability of introgressed alleles, depending on the relative strength of heterosis and Dobzhansky–Muller incompatibilities effects. We show that the inclusion of a correction factor (α, representing the compounded effects of hybrid fitness effects over time) into classic population genetics theory yields accurate fixation probabilities. Despite having a strong impact on the probability of fixation, hybrid fitness effects only subtly change the distribution of fitness effects of introgressed alleles that reach fixation. Although strong Dobzhansky–Muller incompatibility effects may expedite the loss of introgressed alleles, fixation times are largely unchanged by hybrid fitness effects.

The genetic and evolutionary determinants of COVID-19 susceptibility

Devastating pandemics, such as that due to COVID-19, can provide strong testimony to our knowledge of the genetic and evolutionary determinants of infectious disease susceptibility and severity. One of the most remarkable aspects of such outbreaks is the stunning interindividual variability observed in the course of infection. In recent decades, enormous progress has been made in the field of the human genetics of infectious diseases, and an increasing number of human genetic factors have been reported to explain, to a great extent, the observed variability for a large number of infectious agents. However, our understanding of the cellular, molecular, and immunological mechanisms underlying such disparities between individuals and ethnic groups, remains very limited. Here, we discuss recent findings relating to human genetic predisposition to infectious disease, from an immunological or population genetic perspective, and show how these and other innovative approaches have been applied to deciphering the genetic basis of human susceptibility to COVID-19 and the severity of this disease. From an evolutionary perspective, we show how past demographic and selection events characterizing the history of our species, including admixture with archaic humans, such as Neanderthals, facilitated modern human adaptation to the threats imposed by ancient pathogens. In the context of emerging infectious diseases, these past episodes of genetic adaptation may contribute to some of the observed population differences in the outcome of SARS-CoV-2 infection and the severity of COVID-19 illness.

Paleogenomics: The demographic past of prehistoric Europeans

Ancient DNA provides answers to long-standing debates about past human history. New work using demographic modeling on ancient genomes documents the nature and timing of the demographic processes - population size changes, divergences and admixture - that took place in prehistoric Europe.

Mitochondrial Pseudogenes Suggest Repeated Inter-Species Hybridization among Direct Human Ancestors

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

OUP accepted manuscript

Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve “ghosts,” that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA–BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.]

Analysis of whole exome sequencing in severe mental illness hints at selection of brain development and immune related genes

Evolutionary trends may underlie some aspects of the risk for common, non-communicable disorders, including psychiatric disease. We analyzed whole exome sequencing data from 80 unique individuals from India coming from families with two or more individuals with severe mental illness. We used Population Branch Statistics (PBS) to identify variants and genes under positive selection and identified 74 genes as candidates for positive selection. Of these, 20 were previously associated with Schizophrenia, Alzheimer’s disease and cognitive abilities in genome wide association studies. We then checked whether any of these 74 genes were involved in common biological pathways or related to specific cellular or molecular functions. We found that immune related pathways and functions related to innate immunity such as antigen binding were over-represented. We also evaluated for the presence of Neanderthal introgressed segments in these genes and found Neanderthal introgression in a single gene out of the 74 candidate genes. However, the introgression pattern indicates the region is unlikely to be the source for selection. Our findings hint at how selection pressures in individuals from families with a history of severe mental illness may diverge from the general population. Further, it also provides insights into the genetic architecture of severe mental illness, such as schizophrenia and its link to immune factors.

Contemporary biomedical engineering perspective on volitional evolution for human radiotolerance enhancement beyond low-earth orbit

A primary objective of the National Aeronautics and Space Administration (NASA) is expansion of humankind's presence outside low-Earth orbit, culminating in permanent interplanetary travel and habitation. Having no inherent means of physiological detection or protection against ionizing radiation, humans incur capricious risk when journeying beyond low-Earth orbit for long periods. NASA has made large investments to analyze pathologies from space radiation exposure, emphasizing the importance of characterizing radiation's physiological effects. Because natural evolution would require many generations to confer resistance against space radiation, immediately pragmatic approaches should be considered. Volitional evolution, defined as humans steering their own heredity, may inevitably retrofit the genome to mitigate resultant pathologies from space radiation exposure. Recently, uniquely radioprotective genes have been identified, conferring local or systemic radiotolerance when overexpressed in vitro and in vivo. Aiding in this process, the CRISPR/Cas9 technique is an inexpensive and reproducible instrument capable of making limited additions and deletions to the genome. Although cohorts can be identified and engineered to protect against radiation, alternative and supplemental strategies should be seriously considered. Advanced propulsion and mild synthetic torpor are perhaps the most likely to be integrated. Interfacing artificial intelligence with genetic engineering using predefined boundary conditions may enable the computational modeling of otherwise overly complex biological networks. The ethical context and boundaries of introducing genetically pioneered humans are considered.

Palaeogenetics: Dirt, what is it good for? Everything

Researchers in the ancient DNA community have suspected for well over a decade that ancient whole genomes can be found in sediments. Three new studies now provide such evidence and, with it, endless possibilities for future studies of sedimentary ancient DNA.

New insights into human immunity from ancient genomics

Population genetic studies have clearly indicated that immunity and host defense are among the functions most frequently subject to natural selection, and increased our understanding of the biological relevance of the corresponding genes and their contribution to variable immune traits and diseases. Herein, we will focus on some recently studied forms of human adaptation to infectious agents, including hybridization with now-extinct hominins, such as Neanderthals and Denisovans, and admixture between modern human populations. These studies, which are partly enabled by the technological advances in the sequencing of DNA from ancient remains, provide new insight into the sources of immune response variation in contemporary humans, such as the recently reported link between Neanderthal heritage and susceptibility to severe COVID-19 disease. Furthermore, ancient DNA analyses, in both humans and pathogens, allow to measure the action of natural selection on immune genes across time and to reconstruct the impact of past epidemics on the evolution of human immunity.

Inferring archaic introgression from hominin genetic data

Questions surrounding the timing, extent, and evolutionary consequences of archaic admixture into human populations have a long history in evolutionary anthropology. More recently, advances in human genetics, particularly in the field of ancient DNA, have shed new light on the question of whether or not Homo sapiens interbred with other hominin groups. By the late 1990s, published genetic work had largely concluded that archaic groups made no lasting genetic contribution to modern humans; less than a decade later, this conclusion was reversed following the successful DNA sequencing of an ancient Neanderthal. This reversal of consensus is noteworthy, but the reasoning behind it is not widely understood across all academic communities. There remains a communication gap between population geneticists and paleoanthropologists. In this review, we endeavor to bridge this gap by outlining how technological advancements, new statistical methods, and notable controversies ultimately led to the current consensus.

The cis-regulatory effects of modern human-specific variants

The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.

Harnessing pluripotent stem cells as models to decipher human evolution

The study of human evolution, long constrained by a lack of experimental model systems, has been transformed by the emergence of the induced pluripotent stem cell (iPSC) field. iPSCs can be readily established from noninvasive tissue sources, from both humans and other primates; they can be maintained in the laboratory indefinitely, and they can be differentiated into other tissue types. These qualities mean that iPSCs are rapidly becoming established as viable and powerful model systems with which it is possible to address questions in human evolution that were until now logistically and ethically intractable, especially in the quest to understand humans' place among the great apes, and the genetic basis of human uniqueness. In this review, we discuss the key lessons and takeaways of this nascent field; from the types of research, iPSCs make possible to lingering challenges and likely future directions. We provide a comprehensive overview of how the seemingly unlikely combination of iPSCs and explicit evolutionary frameworks is transforming what is possible in our understanding of humanity's past and present.

The Population-Specific Impact of Neandertal Introgression on Human Disease

Since the discovery of admixture between modern humans and Neandertals, multiple studies investigated the effect of Neandertal-derived DNA on human disease and nondisease phenotypes. These studies have linked Neandertal ancestry to skin- and hair-related phenotypes, immunity, neurological, and behavioral traits. However, these inferences have so far been limited to cohorts with participants of European ancestry. Here, I analyze summary statistics from 40 disease GWAS (genome-wide association study) cohorts of ∼212,000 individuals provided by the Biobank Japan Project for phenotypic effects of Neandertal DNA. I show that Neandertal DNA is associated with autoimmune diseases, prostate cancer and type 2 diabetes. Many of these disease associations are linked to population-specific Neandertal DNA, highlighting the importance of studying a wider range of ancestries to characterize the phenotypic legacy of Neandertals in people today.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Gene flow as a simple cause for an excess of high-frequency-derived alleles

Most human populations exhibit an excess of high-frequency variants, leading to a U-shaped site-frequency spectrum (uSFS). This pattern has been generally interpreted as a signature of ongoing episodes of positive selection, or as evidence for a mis-assignment of ancestral/derived allelic states, but uSFS has also been observed in populations receiving gene flow from a ghost population, in structured populations, or after range expansions. In order to better explain the prevalence of high-frequency variants in humans and other populations, we describe here which patterns of gene flow and population demography can lead to uSFS by using extensive coalescent simulations. We find that uSFS can often be observed in a population if gene flow brings a few ancestral alleles from a well-differentiated population. Gene flow can either consist in single pulses of admixture or continuous immigration, but different demographic conditions are necessary to observe uSFS in these two scenarios. Indeed, an extremely low and recent gene flow is required in the case of single admixture events, while with continuous immigration, uSFS occurs only if gene flow started recently at a high rate or if it lasted for a long time at a low rate. Overall, we find that a neutral uSFS occurs under more restrictive conditions in populations having received single pulses of gene flow than in populations exposed to continuous gene flow. We also show that the uSFS observed in human populations from the 1000 Genomes Project can easily be explained by gene flow from surrounding populations without requiring past episodes of positive selection. These results imply that uSFS should be common in non-isolated populations, such as most wild or domesticated plants and animals.

Triangulating Neanderthal cognition: A tale of not seeing the forest for the trees

The inference of Neanderthal cognition, including their cultural and linguistic capabilities, has persisted as a fiercely debated research topic for decades. This lack of consensus is substantially based on inherent uncertainties in reconstructing prehistory out of indirect evidence as well as other methodological limitations. Further factors include systemic difficulties within interdisciplinary discourse, data artifacts, historic research biases, and the sheer scope of the relevant research. Given the degrees of freedom in interpretation ensuing from these complications, any attempt to find approximate answers to the yet unsettled pertinent discourse may not rest on single studies, but instead a careful and comprehensive interdisciplinary synthesis of findings. Triangulating Neanderthals' cognition by considering the plethora of data, diverse perspectives and aforementioned complexities present within the literature constitutes the currently most reliable pathway to tentative conclusions. While some uncertainties remain, such an approach paints the picture of an extensive shared humanity between anatomically modern humans and Neanderthals. This article is categorized under: Cognitive Biology > Evolutionary Roots of Cognition Linguistics > Evolution of Language.

Evolution of Human Brain Size-Associated NOTCH2NL Genes Proceeds toward Reduced Protein Levels

Ever since the availability of genomes from Neanderthals, Denisovans, and ancient humans, the field of evolutionary genomics has been searching for protein-coding variants that may hold clues to how our species evolved over the last ∼600,000 years. In this study, we identify such variants in the human-specific NOTCH2NL gene family, which were recently identified as possible contributors to the evolutionary expansion of the human brain. We find evidence for the existence of unique protein-coding NOTCH2NL variants in Neanderthals and Denisovans which could affect their ability to activate Notch signaling. Furthermore, in the Neanderthal and Denisovan genomes, we find unusual NOTCH2NL configurations, not found in any of the modern human genomes analyzed. Finally, genetic analysis of archaic and modern humans reveals ongoing adaptive evolution of modern human NOTCH2NL genes, identifying three structural variants acting complementary to drive our genome to produce a lower dosage of NOTCH2NL protein. Because copy-number variations of the 1q21.1 locus, encompassing NOTCH2NL genes, are associated with severe neurological disorders, this seemingly contradicting drive toward low levels of NOTCH2NL protein indicates that the optimal dosage of NOTCH2NL may have not yet been settled in the human population.

Phylogenomic test of mitochondrial clues to archaic ancestors in a group of hybridizing swallowtail butterflies

Genomics has revolutionized our understanding of hybridization and introgression, but most of the early evidence for these processes came from studies of mitochondrial introgression. To expand these evolutionary insights from mitochondrial patterns, we evaluate phylogenetic discordance across the nuclear genomes of a hybridizing system, the Papilio machaon group of swallowtail butterflies. This species group contains three hybrid lineages (P. brevicauda, P. joanae, and P. m. kahli) that are geographically disjunct across North America and have complete fixation of a mitochondrial lineage that is otherwise primarily found in P. m. hudsonianus, a boreal subspecies of the Holarctic P. machaon. Genome-wide nuclear markers place the three hybrid lineages as a monophyletic group that is sister to P. polyxenes/P. zelicaon rather than P. machaon, although ancient hybridization between a subspecies of P. machaon and the ancestor of these three lineages is also shown by their greater nuclear affinity to P. m. hudsonianus than to other subspecies of P. machaon. Individuals from contemporary hybrid swarms in Alberta, where mitochondrial DNA fixation has not occurred, were more intermediate between their respective parent species, demonstrating diversity in mito-nuclear discordance following hybrid interactions. Our new phylogenetic findings for the P. machaon species group also include: subspecific paraphyly within P. machaon itself across its Holarctic distribution; paraphyly of P. zelicaon relative to P. polyxenes; and more divergent placement of a Mediterranean species, P. hospiton. These results provide the first comprehensive genomic evaluation of relationships within this species group and provide insight into the evolutionary dynamics of hybridization and mitochondrial introgression.

VolcanoFinder: Genomic scans for adaptive introgression

Recent research shows that introgression between closely-related species is an important source of adaptive alleles for a wide range of taxa. Typically, detection of adaptive introgression from genomic data relies on comparative analyses that require sequence data from both the recipient and the donor species. However, in many cases, the donor is unknown or the data is not currently available. Here, we introduce a genome-scan method—VolcanoFinder—to detect recent events of adaptive introgression using polymorphism data from the recipient species only. VolcanoFinder detects adaptive introgression sweeps from the pattern of excess intermediate-frequency polymorphism they produce in the flanking region of the genome, a pattern which appears as a volcano-shape in pairwise genetic diversity. Using coalescent theory, we derive analytical predictions for these patterns. Based on these results, we develop a composite-likelihood test to detect signatures of adaptive introgression relative to the genomic background. Simulation results show that VolcanoFinder has high statistical power to detect these signatures, even for older sweeps and for soft sweeps initiated by multiple migrant haplotypes. Finally, we implement VolcanoFinder to detect archaic introgression in European and sub-Saharan African human populations, and uncovered interesting candidates in both populations, such as TSHR in Europeans and TCHH-RPTN in Africans. We discuss their biological implications and provide guidelines for identifying and circumventing artifactual signals during empirical applications of VolcanoFinder.

The process by which beneficial alleles are introduced into a species from a closely-related species is termed adaptive introgression. We present an analytically-tractable model for the effects of adaptive introgression on non-adaptive genetic variation in the genomic region surrounding the beneficial allele. The result we describe is a characteristic volcano-shaped pattern of increased variability that arises around the positively-selected site, and we introduce an open-source method VolcanoFinder to detect this signal in genomic data. Importantly, VolcanoFinder is a population-genetic likelihood-based approach, rather than a comparative-genomic approach, and can therefore probe genomic variation data from a single population for footprints of adaptive introgression, even from a priori unknown and possibly extinct donor species.

Excavating ghost footprints and tangled trees from modern genomes

Due to pervasive gene flow and admixture, simple bifurcating trees often do not provide an accurate representation of relationships among diverging lineages, but limited resolution in the available genomic data and the spatial distribution of samples has hindered detailed insights regarding the evolutionary and demographic history of many species and populations. In this issue of Molecular Ecology, Foote et al. (2019) combine a powerful sampling design with novel analytical methods adopted from human genetics to describe previously unrecognized patterns of recurrent vicariance and admixture among lineages in the globally distributed killer whale (Orcinus orca). Based on sequence data from modern samples alone, they discover clear signatures of ancient admixture with a now extinct "ghost" lineage, providing one of the first accounts of archaic introgression in a nonhominid species. Coupling a cost-effective sequencing strategy with novel analytical approaches, their paper provides a roadmap for advancing inference of evolutionary history in other nonmodel species, promising exciting times ahead for our field.

Adaptive introgression as a driver of local adaptation to climate in European white oaks

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

Evolutionary and population (epi)genetics of immunity to infection

Immune response is one of the functions that have been more strongly targeted by natural selection during human evolution. The evolutionary genetic dissection of the immune system has greatly helped to distinguish genes and functions that are essential, redundant or advantageous for human survival. It is also becoming increasingly clear that admixture between early Eurasians with now-extinct hominins such as Neanderthals or Denisovans, or admixture between modern human populations, can be beneficial for human adaptation to pathogen pressures. In this review, we discuss how the integration of population genetics with functional genomics in diverse human populations can inform about the changes in immune functions related to major lifestyle transitions (e.g., from hunting and gathering to farming), the action of natural selection to the evolution of the immune system, and the history of past epidemics. We also highlight the need of expanding the characterization of the immune system to a larger array of human populations-particularly neglected human groups historically exposed to different pathogen pressures-to fully capture the relative contribution of genetic, epigenetic, and environmental factors to immune response variation in humans.

Patterns of Genomic Differentiation in the Drosophila nasuta Species Complex

The Drosophila nasuta species complex contains over a dozen recently diverged species that are distributed widely across South-East Asia, and which shows varying degrees of pre- and postzygotic isolation. Here, we assemble a high-quality genome for D. albomicans using single-molecule sequencing and chromatin conformation capture, and draft genomes for 11 additional species and 67 individuals across the clade, to infer the species phylogeny and patterns of genetic diversity in this group. Our assembly recovers entire chromosomes, and we date the origin of this radiation ∼2 Ma. Despite low levels of overall differentiation, most species or subspecies show clear clustering into their designated taxonomic groups using population genetics and phylogenetic methods. Local evolutionary history is heterogeneous across the genome, and differs between the autosomes and the X chromosome for species in the sulfurigaster subgroup, likely due to autosomal introgression. Our study establishes the nasuta species complex as a promising model system to further characterize the evolution of pre- and postzygotic isolation in this clade.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6

BACKGROUND

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.

RESULTS

Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.

CONCLUSIONS

Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

The repertoire of family A-peptide GPCRs in archaic hominins

Given the importance of G-protein coupled receptors in the regulation of many physiological functions, deciphering the relationships between genotype and phenotype in past and present hominin GPCRs is of main interest to understand the evolutionary process that contributed to the present-day variability in human traits and health. Here, we carefully examined the publicly available genomic and protein sequence databases of the archaic hominins (Neanderthal and Denisova) to draw up the catalog of coding variations in GPCRs for peptide ligands, in comparison with living humans. We then searched in the literature the functional changes, phenotypes and risk of disease possibly associated with the detected variants. Our survey suggests that Neanderthal and Denisovan hominins were likely prone to lower risk of obesity, to enhanced platelet aggregation in response to thrombin, to better response to infection, to less anxiety and aggressiveness and to favorable sociability. While some archaic variants were likely advantageous in the past, they might be responsible for maladaptive disorders today in the context of modern life and/or specific regional distribution. For example, an archaic haplotype in the neuromedin receptor 2 is susceptible to confer risk of diabetic nephropathy in type 1 diabetes in present-day Europeans. Paying attention to the pharmacological properties of some of the archaic variants described in this study may be helpful to understand the variability of therapeutic efficacy between individuals or ethnic groups.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

The role of recombination on genome-wide patterns of local ancestry exemplified by supplemented brook charr populations

Assessing the immediate and long-term evolutionary consequences of human-mediated hybridization is of major concern for conservation biology. Several studies have documented how selection in interaction with recombination modulates introgression at a genome-wide scale, but few have considered the dynamics of this process within and among chromosomes. Here, we used an exploited freshwater fish, the brook charr (Salvelinus fontinalis), for which decades of stocking practices have resulted in admixture between wild populations and an introduced domestic strain, to assess both the temporal dynamics and local chromosomal variation in domestic ancestry. We provide a detailed picture of the domestic ancestry patterns across the genome using about 33,000 mapped single nucleotide polymorphisms genotyped in 611 individuals from 24 supplemented populations. For each lake, we distinguished early- and late-generation hybrids using information regarding admixture tracts. To assess the selective outcomes following admixture we then evaluated the relationship between recombination and admixture proportions at three different scales: the whole genome, chromosomes and within 2-Mb windows. This allowed us to detect a wide range of evolutionary mechanisms varying along the genome, as reflected by the finding of favoured or disfavoured introgression of domestic haplotypes. Among these, the main factor modulating local ancestry was probably the presence of deleterious recessive mutations in the wild populations, which can be efficiently hidden to selection in the presence of long admixture tracts. Overall, our results emphasize the relevance of taking into consideration local ancestry information to assess both the temporal and the chromosomal variation in local admixture ancestry toward better understanding post-hybridization evolutionary outcomes.

Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe

European populations display low genetic differentiation as the result of long-term blending of their ancient founding ancestries. However, it is unclear how the combination of ancient ancestries related to early foragers, Neolithic farmers, and Bronze Age nomadic pastoralists can explain the distribution of genetic variation across Europe. Populations in natural crossroads like the Italian peninsula are expected to recapitulate the continental diversity, but have been systematically understudied. Here, we characterize the ancestry profiles of Italian populations using a genome-wide dataset representative of modern and ancient samples from across Italy, Europe, and the rest of the world. Italian genomes capture several ancient signatures, including a non-steppe contribution derived ultimately from the Caucasus. Differences in ancestry composition, as the result of migration and admixture, have generated in Italy the largest degree of population structure detected so far in the continent, as well as shaping the amount of Neanderthal DNA in modern-day populations.

The presence and impact of reference bias on population genomic studies of prehistoric human populations

Haploid high quality reference genomes are an important resource in genomic research projects. A consequence is that DNA fragments carrying the reference allele will be more likely to map successfully, or receive higher quality scores. This reference bias can have effects on downstream population genomic analysis when heterozygous sites are falsely considered homozygous for the reference allele. In palaeogenomic studies of human populations, mapping against the human reference genome is used to identify endogenous human sequences. Ancient DNA studies usually operate with low sequencing coverages and fragmentation of DNA molecules causes a large proportion of the sequenced fragments to be shorter than 50 bp-reducing the amount of accepted mismatches, and increasing the probability of multiple matching sites in the genome. These ancient DNA specific properties are potentially exacerbating the impact of reference bias on downstream analyses, especially since most studies of ancient human populations use pseudo-haploid data, i.e. they randomly sample only one sequencing read per site. We show that reference bias is pervasive in published ancient DNA sequence data of prehistoric humans with some differences between individual genomic regions. We illustrate that the strength of reference bias is negatively correlated with fragment length. Most genomic regions we investigated show little to no mapping bias but even a small proportion of sites with bias can impact analyses of those particular loci or slightly skew genome-wide estimates. Therefore, reference bias has the potential to cause minor but significant differences in the results of downstream analyses such as population allele sharing, heterozygosity estimates and estimates of archaic ancestry. These spurious results highlight how important it is to be aware of these technical artifacts and that we need strategies to mitigate the effect. Therefore, we suggest some post-mapping filtering strategies to resolve reference bias which help to reduce its impact substantially.

Impact and Evolutionary Determinants of Neanderthal Introgression on Transcriptional and Post-Transcriptional Regulation

Archaic admixture is increasingly recognized as an important source of diversity in modern humans, and Neanderthal haplotypes cover 1%-3% of the genome of present-day Eurasians. Recent work has shown that archaic introgression has contributed to human phenotypic diversity, mostly through the regulation of gene expression. Yet the mechanisms through which archaic variants alter gene expression and the forces driving the introgression landscape at regulatory regions remain elusive. Here, we explored the impact of archaic introgression on transcriptional and post-transcriptional regulation. We focused on promoters and enhancers across 127 different tissues as well as on microRNA (miRNA)-mediated regulation. Although miRNAs themselves harbor few archaic variants, we found that some of these variants may have a strong impact on miRNA-mediated gene regulation. Enhancers were by far the regulatory elements most affected by archaic introgression: up to one-third of the tissues we tested presented significant enrichments. Specifically, we found strong enrichments of archaic variants in adipose-related tissues and primary T cells, even after accounting for various genomic and evolutionary confounders such as recombination rate and background selection. Interestingly, we identified signatures of adaptive introgression at enhancers of some key regulators of adipogenesis, raising the interesting hypothesis of a possible adaptation of early Eurasians to colder climates. Collectively, this study sheds new light on the mechanisms through which archaic admixture has impacted gene regulation in Eurasians and, more generally, increases our understanding of the contribution of Neanderthals to the regulation of acquired immunity and adipose homeostasis in modern humans.

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

Dating admixture events is unsolved problem in multi-way admixed populations

Advances in human sequencing technologies, coupled with statistical and computational tools, have fostered the development of methods for dating admixture events. These methods have merits and drawbacks in estimating admixture events in multi-way admixed populations. Here, we first provide a comprehensive review and comparison of current methods pertinent to dating admixture events. Second, we assess various admixture dating tools. We do so by performing various simulations. Third, we apply the top two assessed methods to real data of a uniquely admixed population from South Africa. Results reveal that current dating admixture models are not sufficiently equipped to estimate ancient admixtures events and to identify multi-faceted admixture events in complex multi-way admixed populations. We conclude with a discussion of research areas where further work on dating admixture-based methods is needed.

Time estimations by network of beta globin gene cluster haplotypes linked with Hb D-Los Angeles [β121 (GH4) Glu → Gln GAA → CAA] mutation in the world populations

Background

β‐Globin gene cluster haplotypes associated with the Hb D‐Los Angeles mutation have been reported in many different locations in different populations including Italy, Iran, Thailand, Belgium, Mexico, Holland, and Turkey. In this study, we have identified genetic relationships and formation periods between the haplotypes reported in the world regarding the Hb D‐Los Angeles.

Methods

We comparatively analyzed the RFLP (restriction fragment length polymorphism) data in Denizli region and world populations using Arlequin 3.5 statistical software program. The data obtained from the Arlequin software were then entered into the Phylogenetic Network software to calculate the age estimates and to discover possible links between the haplotypes.

Results

We observed the frequencies of the β‐globin gene haplotypes for the seven polymorphic restriction sites around the world and calculated the estimated time of haplotypes using Network software on the basis of ancestral haplotypes. We performed the phylogenetic network analysis of the haplotypes linked with Hb D‐Los Angeles mutation by processing the data of frequency and age estimations with Network software.

Conclusion

Our period of time results suggests that HAP1 was formed before modern human migration to Asia and/or independent origin of the Hb D‐Los Angeles mutation from other populations. Considering that the population in Denizli region started the Hb D‐Los Angeles mutation past about 40,000 years ago, it can be said that HAP1, HAP15, and HAP21 belong to the gene pool with an external effect. Our period of time results of HAP6 is compatible with published dating results.