The date of interbreeding between Neandertals and modern humans

Aurignacian groups at Isturitz (France) adapted to a shifting environment upon their arrival in Western Europe ∼42,000 years ago

The Marine Isotope Stage 3 is a context of considerable climatic instability. Establishing the link between global climate changes and their impact on the local ecological contexts and prey exploited by human populations is challenging. Still, it is necessary to understand better the local conditions where humans lived to unravel how they adapted to fluctuating environmental conditions. Here, we address this question by studying 250 osteodental elements from animals hunted and consumed by human groups at Isturitz, a rich and well-documented French archaeological site and one of the earliest in Western Europe where the Aurignacian technoculture has been attested. To do so, we set up a multiproxy approach (archaeozoology, three-dimensional dental microwear texture analyses, and stable isotopic analyses of δ18O and δ13C in enamel bioapatite and δ13C, δ15N, and δ34S in bone collagen) that informs us on a timeline from the first years to the last few days of an animal's life. We reconstructed their ecologies and paleoenvironments during the different Aurignacian phases at Isturitz. Our findings indicate that the first human occupations at Isturitz occurred under cold and arid conditions, rapidly becoming even cooler and drier. Limited changes are observed in the human-environment-prey relationship despite this unstable climatic context where significant changes in rainfall, temperature, and a gradual opening of environments and some changes in the faunal assemblage occurred. Our findings suggest that human groups hunted in similar territories and utilized comparable strategies throughout the temporal sequence. Our multiproxy approach, combining complementary analyses, provides a better understanding of the adaptation strategies when the first phases of the Upper Paleolithic were emerging in Western Europe.

Tracing the Evolution of Human Immunity Through Ancient DNA

Infections have imposed strong selection pressures throughout human evolution, making the study of natural selection's effects on immunity genes highly complementary to disease-focused research. This review discusses how ancient DNA studies, which have revolutionized evolutionary genetics, increase our understanding of the evolution of human immunity. These studies have shown that interbreeding between modern humans and Neanderthals or Denisovans has influenced present-day immune responses, particularly to viruses. Additionally, ancient genomics enables the tracking of how human immunity has evolved across cultural transitions, highlighting strong selection since the Bronze Age in Europe (<4,500 years) and potential genetic adaptations to epidemics raging during the Middle Ages and the European colonization of the Americas. Furthermore, ancient genomic studies suggest that the genetic risk for noninfectious immune disorders has gradually increased over millennia because alleles associated with increased risk for autoimmunity and inflammation once conferred resistance to infections. The challenge now is to extend these findings to diverse, non-European populations and to provide a more global understanding of the evolution of human immunity.

Climate change in Europe between 90 and 50 kyr BP and Neanderthal territorial habitability

After having lived as the dominant human species in Europe for over 200 kyr, Homo neanderthalensis (the Neanderthals) disappeared around 40 kyr BP (Before Present) Higham T (2014). Competition with Homo sapiens, who arrived in Europe around the same time, is often invoked to explain this extinction. Others have argued that climate change may have reduced the living space of this population making its disappearance more rapid. In order to test the climate change hypothesis we modelled the Neanderthals' ecological niches in Europe between 90 and 50 kyr BP through paleoenvironmental reconstructions and Eco-Cultural Niche Modelling. We selected five environmental variables (orographic height, mean annual precipitation, mean temperature of the coldest month, carrying capacity and friction, see below) from climate model simulations of 5 periods between 90 and 50 kyr BP in Europe. We used Structural Similarity (SSIM) index to compare the probability maps of suitable niches to Neanderthals performed by Maxent. After a strong initial environmental change between the first (P1 = 90 to 83 kyr BP) and second (P2 = 83 to 69 kyr BP) periods, our results show that large areas highly suitable for Neanderthal occupation persisted across Europe. As our results show an increase/stability of the areas suitable to Neanderthals, the question of the cause of the decrease or displacement of the Neanderthal population towards southern Europe after this climatic change remains open.

Gene expression patterns of the developing human face at single cell resolution reveal cell type contributions to normal facial variation and disease risk

Craniofacial development gives rise to the complex structures of the face and involves the interplay of diverse cell types. Despite its importance, our understanding of human-specific craniofacial developmental mechanisms and their genetic underpinnings remains limited. Here, we present a comprehensive single-nucleus RNA sequencing (snRNA-seq) atlas of human craniofacial development from craniofacial tissues of 24 embryos that span six key time points during the embryonic period (4-8 post-conception weeks). This resource resolves the transcriptional dynamics of seven major cell types and uncovers distinct major cell types, including muscle progenitors and cranial neural crest cells (CNCCs), as well as dozens of subtypes of ectoderm and mesenchyme. Comparative analyses reveal substantial conservation of major cell types, alongside human biased differences in gene expression programs. CNCCs, which play a crucial role in craniofacial morphogenesis, exhibit the lowest marker gene conservation, underscoring their evolutionary plasticity. Spatial transcriptomics further localizes cell populations, providing a detailed view of their developmental roles and anatomical context. We also link these developmental processes to genetic variation, identifying cell type-specific enrichments for common variants associated with facial morphology and rare variants linked to orofacial clefts. Intriguingly, Neanderthal-introgressed sequences are enriched near genes with biased expression in cartilage and specialized ectodermal subtypes, suggesting their contribution to modern human craniofacial features. This atlas offers unprecedented insights into the cellular and genetic mechanisms shaping the human face, highlighting conserved and distinctly human aspects of craniofacial biology. Our findings illuminate the developmental origins of craniofacial disorders, the genetic basis of facial variation, and the evolutionary legacy of ancient hominins. This work provides a foundational resource for exploring craniofacial biology, with implications for developmental genetics, evolutionary biology, and clinical research into congenital anomalies.

Earliest modern human genomes constrain timing of Neanderthal admixture

Modern humans arrived in Europe more than 45,000 years ago, overlapping at least 5,000 years with Neanderthals1-4. Limited genomic data from these early modern humans have shown that at least two genetically distinct groups inhabited Europe, represented by Zlatý kůň, Czechia3 and Bacho Kiro, Bulgaria2. Here we deepen our understanding of early modern humans by analysing one high-coverage genome and five low-coverage genomes from approximately 45,000-year-old remains from Ilsenhöhle in Ranis, Germany4, and a further high-coverage genome from Zlatý kůň. We show that distant familial relationships link the Ranis and Zlatý kůň individuals and that they were part of the same small, isolated population that represents the deepest known split from the Out-of-Africa lineage. Ranis genomes harbour Neanderthal segments that originate from a single admixture event shared with all non-Africans that we date to approximately 45,000-49,000 years ago. This implies that ancestors of all non-Africans sequenced so far resided in a common population at this time, and further suggests that modern human remains older than 50,000 years from outside Africa represent different non-African populations.

Ignoring population structure in hominin evolutionary models can lead to the inference of spurious admixture events

Genomic and ancient DNA data have revolutionized palaeoanthropology and our vision of human evolution, with indisputable landmarks like the sequencing of Neanderthal and Denisovan genomes. Yet, using genetic data to identify, date and quantify evolutionary events-such as ancient bottlenecks or admixture-is not straightforward, as inferences may depend on model assumptions. In the last two decades, the idea that Neanderthals and members of the Homo sapiens lineage interbred has gained momentum. From the status of unlikely theory, it has reached consensus among human evolutionary biologists. This theory is mainly supported by statistical approaches that depend on demographic models minimizing or ignoring population structure, despite its widespread occurrence and the fact that, when ignored, population structure can lead to the inference of spurious demographic events. We simulated genomic data under a structured and admixture-free model of human evolution, and found that all the tested admixture approaches identified long Neanderthal fragments in our simulated genomes and an admixture event that never took place. We also observed that several published admixture models failed to predict important empirical diversity or admixture statistics, and that we could identify several scenarios from our structured model that better predicted these statistics jointly. Using a simulated time series of ancient DNA, the structured scenarios could also predict the trajectory of the empirical D statistics. Our results suggest that models accounting for population structure are fundamental to improve our understanding of human evolution, and that admixture between Neanderthals and H. sapiens needs to be re-evaluated in the light of structured models. Beyond the Neanderthal case, we argue that ancient hybridization events, which are increasingly documented in many species, including with other hominins, may also benefit from such re-evaluation.

Patterns of population structure and genetic variation within the Saudi Arabian population

The Arabian Peninsula is considered the initial site of historic human migration out of Africa. The modern-day indigenous Arabians are believed to be the descendants who remained from the ancient split of the migrants into Eurasia. Here, we investigated how the population history and cultural practices such as endogamy have shaped the genetic variation of the Saudi Arabians. We genotyped 3,352 individuals and identified twelve genetic sub-clusters that corresponded to the geographical distribution of different tribal regions, differentiated by distinct components of ancestry based on comparisons to modern and ancient DNA references. These sub-clusters also showed variation across ranges of the genome covered in runs of homozygosity, as well as differences in population size changes over time. Using 25,488,981 variants found in whole genome sequencing data (WGS) from 302 individuals, we found that the Saudi tend to show proportionally more deleterious alleles than neutral alleles when compared to Africans/African Americans from gnomAD (e.g. a 13% increase of deleterious alleles annotated by AlphaMissense between 0.5 - 5% frequency in Saudi, compared to 7% decrease of the benign alleles; P < 0.001). Saudi sub-clusters with greater inbreeding and lower effective population sizes showed greater enrichment of deleterious alleles as well. Additionally, we found that approximately 10% of the variants discovered in our WGS data are not observed in gnomAD; these variants are also enriched with deleterious annotations. To accelerate studying the population-enriched deleterious alleles and their health consequences in this population, we made available the allele frequency estimates of 25,488,981 variants discovered in our samples. Taken together, our results suggest that Saudi's population history impacts its pattern of genetic variation with potential consequences to the population health. It further highlights the need to sequence diverse and unique populations so to provide a foundation on which to interpret medical- and pharmaco- genomic findings from these populations.

High-resolution genomic history of early medieval Europe

Many known and unknown historical events have remained below detection thresholds of genetic studies because subtle ancestry changes are challenging to reconstruct. Methods based on shared haplotypes1,2 and rare variants3,4 improve power but are not explicitly temporal and have not been possible to adopt in unbiased ancestry models. Here we develop Twigstats, an approach of time-stratified ancestry analysis that can improve statistical power by an order of magnitude by focusing on coalescences in recent times, while remaining unbiased by population-specific drift. We apply this framework to 1,556 available ancient whole genomes from Europe in the historical period. We are able to model individual-level ancestry using preceding genomes to provide high resolution. During the first half of the first millennium CE, we observe at least two different streams of Scandinavian-related ancestry expanding across western, central and eastern Europe. By contrast, during the second half of the first millennium CE, ancestry patterns suggest the regional disappearance or substantial admixture of these ancestries. In Scandinavia, we document a major ancestry influx by approximately 800 CE, when a large proportion of Viking Age individuals carried ancestry from groups related to central Europe not seen in individuals from the early Iron Age. Our findings suggest that time-stratified ancestry analysis can provide a higher-resolution lens for genetic history.

Neanderthal ancestry through time: Insights from genomes of ancient and present-day humans

Gene flow from Neanderthals has shaped genetic and phenotypic variation in modern humans. We generated a catalog of Neanderthal ancestry segments in more than 300 genomes spanning the past 50,000 years. We examined how Neanderthal ancestry is shared among individuals over time. Our analysis revealed that the vast majority of Neanderthal gene flow is attributable to a single, shared extended period of gene flow that occurred between 50,500 to 43,500 years ago, as evidenced by ancestry correlation, colocalization of Neanderthal segments across individuals, and divergence from the sequenced Neanderthals. Most natural selection-positive and negative-on Neanderthal variants occurred rapidly after the gene flow. Our findings provide new insights into how contact with Neanderthals shaped modern human origins and adaptation.

Haplotype-Based Approach Represents Locus Specificity in the Genomic Diversification Process in Humans (Homo sapiens)

BACKGROUND/OBJECTIVES

Recent progress in evolutionary genomics on human (Homo sapiens) populations has revealed complex demographic events and genomic changes. These include population expansion with complicated migration, substantial population structure, and ancient introgression from other hominins, as well as human characteristics selections. Nevertheless, the genomic regions in which such evolutionary events took place have remained unclear.

METHODS

Here, we focused on eight loci containing the haplotypes that were previously presented as atypical for the mutation pattern in sequence and/or geographic distribution pattern with the model of recent African origin, which constitute two major clusters: African only, and global. This was the consensus model before information regarding introgression from Neanderthal (Homo neanderthalensis) was available. We compared diversity in identical datasets of the modern human population genome, with the 1000 Genomes project among them.

RESULTS/CONCLUSIONS

This study identified representative genomic regions that show traces of various demographic events and genomic changes that modern humans have undergone by categorizing the relationships in sequence similarity and in worldwide geographic distribution among haplotypes.

Towards Reliable Detection of Introgression in the Presence of Among-Species Rate Variation

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression-the exchange of genetic material through hybridization and backcrossing-are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression"- referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to old systems, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites along the genome and implemented this test in the program Dsuite.

Life history and ancestry of the late Upper Palaeolithic infant from Grotta delle Mura, Italy

The biological aspects of infancy within late Upper Palaeolithic populations and the role of southern refugia at the end of the Last Glacial Maximum are not yet fully understood. This study presents a multidisciplinary, high temporal resolution investigation of an Upper Palaeolithic infant from Grotta delle Mura (Apulia, southern Italy) combining palaeogenomics, dental palaeohistology, spatially-resolved geochemical analyses, direct radiocarbon dating, and traditional anthropological studies. The skeletal remains of the infant - Le Mura 1 - were directly dated to 17,320-16,910 cal BP. The results portray a biological history of the infant's development, early life, health and death (estimated at ~72 weeks). They identify, several phenotypic traits and a potential congenital disease in the infant, the mother's low mobility during gestation, and a high level of endogamy. Furthermore, the genomic data indicates an early spread of the Villabruna-like components along the Italian peninsula, confirming a population turnover around the time of the Last Glacial Maximum, and highlighting a general reduction in genetic variability from northern to southern Italy. Overall, Le Mura 1 contributes to our better understanding of the early stages of life and the genetic puzzle in the Italian peninsula at the end of the Last Glacial Maximum.

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

The ear of the Sima de los Huesos hominins (Atapuerca, Spain)

Previous studies on the morphology of the inner ear (semicircular canals and cochlea) in the Sima de los Huesos hominin sample have provided important results on the evolution of these structures in the Neandertal lineage. Similarly, studies of the anatomy of the external and middle ear cavities of the Sima de los Huesos hominins have also provided important data on the auditory capacities of this European Middle Pleistocene population. The present contribution provides unpublished data on three new middle ear variables from the Sima de los Huesos fossils and compares these data with values from samples of Pan troglodytes, Homo neanderthalensis and Homo sapiens. The results of this analysis are combined with those obtained in previous studies to characterize the anatomy of the outer, middle and inner ear in the Sima de los Huesos fossils, as well as to establish the order of appearance of the features that characterize Neandertal ears. As in other cranial structures, the ear region in the Sima de los Huesos show a mosaic evolutionary pattern that includes primitive traits, others shared exclusively with Neandertals, and others that are specific to the Sima de los Huesos hominins. Neandertals and Sima de los Huesos hominins share two exclusive features of the middle ear that are among the first characteristics of the Neandertal lineage: a long tympanic cavity and a large entrance and exit of the mastoid antrum. Along with these traits, the Sima de los Huesos hominins present two specialized features: large volumes of the tympanic cavity and the mastoid antrum. Finally, the middle ear of the Neandertals is characterized by the presence of small angles between the tympanic axis and the plane of the oval window.

Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations

Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.

Neandertal ancestry through time: Insights from genomes of ancient and present-day humans

Gene flow from Neandertals has shaped the landscape of genetic and phenotypic variation in modern humans. We identify the location and size of introgressed Neandertal ancestry segments in more than 300 genomes spanning the last 50,000 years. We study how Neandertal ancestry is shared among individuals to infer the time and duration of the Neandertal gene flow. We find the correlation of Neandertal segment locations across individuals and their divergence to sequenced Neandertals, both support a model of single major Neandertal gene flow. Our catalog of introgressed segments through time confirms that most natural selection-positive and negative-on Neandertal ancestry variants occurred immediately after the gene flow, and provides new insights into how the contact with Neandertals shaped human origins and adaptation.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression

INTRODUCTION

Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear.

MATERIALS AND METHODS

Population genetics analysis and functional genomics work was performed for this locus.

RESULTS

Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human.

CONCLUSION

This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

Pharmacogenetic Variation in Neanderthals and Denisovans and Implications for Human Health and Response to Medications

Modern humans carry both Neanderthal and Denisovan (archaic) genome elements that are part of the human gene pool and affect the life and health of living individuals. The impact of archaic DNA may be particularly evident in pharmacogenes-genes responsible for the processing of exogenous substances such as food, pollutants, and medications-as these can relate to changing environmental effects, and beneficial variants may have been retained as modern humans encountered new environments. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remain understudied. Here, we explore 11 key cytochrome P450 genes (CYP450) involved in 75% of all drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these 11 CYP450 genes in archaic individuals and find important predicted phenotypic differences relative to modern human variants. We identify several single nucleotide variants shared between archaic and modern humans in each gene, including some potentially function-altering mutations in archaic CYP450 genes, which may result in altered metabolism in living people carrying these variants. We also identified several variants in the archaic CYP450 genes that are novel and unique to archaic humans as well as one gene, CYP2B6, that shows evidence for a gene duplication found only in Neanderthals and modern Africans. Finally, we highlight CYP2A6, CYP2C9, and CYP2J2, genes which show evidence for archaic introgression into modern humans and posit evolutionary hypotheses that explain their allele frequencies in modern populations.

Archaic Introgression Shaped Human Circadian Traits

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Diverse African genomes reveal selection on ancient modern human introgressions in Neanderthals

Comparisons of Neanderthal genomes to anatomically modern human (AMH) genomes show a history of Neanderthal-to-AMH introgression stemming from interbreeding after the migration of AMHs from Africa to Eurasia. All non-sub-Saharan African AMHs have genomic regions genetically similar to Neanderthals that descend from this introgression. Regions of the genome with Neanderthal similarities have also been identified in sub-Saharan African populations, but their origins have been unclear. To better understand how these regions are distributed across sub-Saharan Africa, the source of their origin, and what their distribution within the genome tells us about early AMH and Neanderthal evolution, we analyzed a dataset of high-coverage, whole-genome sequences from 180 individuals from 12 diverse sub-Saharan African populations. In sub-Saharan African populations with non-sub-Saharan African ancestry, as much as 1% of their genomes can be attributed to Neanderthal sequence introduced by recent migration, and subsequent admixture, of AMH populations originating from the Levant and North Africa. However, most Neanderthal homologous regions in sub-Saharan African populations originate from migration of AMH populations from Africa to Eurasia ∼250 kya, and subsequent admixture with Neanderthals, resulting in ∼6% AMH ancestry in Neanderthals. These results indicate that there have been multiple migration events of AMHs out of Africa and that Neanderthal and AMH gene flow has been bi-directional. Observing that genomic regions where AMHs show a depletion of Neanderthal introgression are also regions where Neanderthal genomes show a depletion of AMH introgression points to deleterious interactions between introgressed variants and background genomes in both groups-a hallmark of incipient speciation.

Archaic Introgression Shaped Human Circadian Traits

Introduction

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultra-violet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology, and whether archaic introgression adaptively contributed to human chronotypes remains unknown.

Results

Here we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, RORC), and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among eQTLs for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency.

Conclusions

These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Methods for Assessing Population Relationships and History Using Genomic Data

Genetic data contain a record of our evolutionary history. The availability of large-scale datasets of human populations from various geographic areas and timescales, coupled with advances in the computational methods to analyze these data, has transformed our ability to use genetic data to learn about our evolutionary past. Here, we review some of the widely used statistical methods to explore and characterize population relationships and history using genomic data. We describe the intuition behind commonly used approaches, their interpretation, and important limitations. For illustration, we apply some of these techniques to genome-wide autosomal data from 929 individuals representing 53 worldwide populations that are part of the Human Genome Diversity Project. Finally, we discuss the new frontiers in genomic methods to learn about population history. In sum, this review highlights the power (and limitations) of DNA to infer features of human evolutionary history, complementing the knowledge gleaned from other disciplines, such as archaeology, anthropology, and linguistics.

Non-ultrametric phylogenetic trees shed new light on Neanderthal introgression

Ultrametric spaces are widely used to depict evolutionary times in phylogenetic trees since they assume that every population/species is located at the tips of bifurcating branches of the same length. The discrete branching of ultrametric trees permits the measurement of distances between pairs of individuals that are proportional to their divergence time. Here the traditional ultrametric concept of bifurcating and divergent phylogenetic tree is overturned and a new type of non-ultrametric diagram is introduced. The objective of this study is the description of gene flows in branching species/populations in terms of converging trees instead of bifurcating trees. To provide an operational example, the paleoanthropological issue of the date of Neanderthal genome's introgression in non-African humans is examined. Neanderthals and ancient humans are not anymore two species that exchange chunks of DNA, rather become a single, novel cluster of extant hominins that must be considered by itself. The novel converging, non-ultrametric phylogenetic trees permit the calibration of molecular clocks with a twofold benefit. When the date of the branching of two population/species from a common ancestor is known, the novel approach allows to calculate the time of subsequent introgressions. On the contrary, when the date of the introgression between two population/species is known, the novel approach allows to detect the time of their previous branching from a common ancestor.

Evolutionary analysis of JC polyomavirus in Misiones' population yields insight into the population dynamics of the early human dispersal in the Americas

BACKGROUND

JC polyomavirus (JCV) has an ethno-geographical distribution across human populations.

OBJECTIVE

Study the origins of the population of Misiones (Argentina) by using JCV as genetic marker.

METHODS

Viral detection and characterization was conducted by PCR amplification and evolutionary analysis of the intergenic region sequences.

RESULTS

22 out of 121 samples were positive for JCV, including 5 viral lineages: MY (n = 8), Eu-a (n = 7), B1-c (n = 4), B1-b (n = 2) and Af2 (n = 1). MY sequences clustered within a branch of Native American origin that diverged from its Asian counterpart about 21,914 years ago (HPD 95% interval 15,383-30,177), followed by a sustained demographic expansion around 5000 years ago.

CONCLUSIONS

JCV in Misiones reflects the multiethnic origin of the current population, with an important Amerindian contribution. Analysis of the MY viral lineage shows a pattern consistent with the arrival of early human migrations to the Americas and a population expansion by the pre-Columbian native societies.

Major Genetic Risk Factors for Dupuytren's Disease Are Inherited From Neandertals

Dupuytren's disease is characterized by fingers becoming permanently bent in a flexed position. Whereas people of African ancestry are rarely afflicted by Dupuytren's disease, up to ∼30% of men over 60 years suffer from this condition in northern Europe. Here, we meta-analyze 3 biobanks comprising 7,871 cases and 645,880 controls and find 61 genome-wide significant variants associated with Dupuytren's disease. We show that 3 of the 61 loci harbor alleles of Neandertal origin, including the second and third most strongly associated ones (P = 6.4 × 10-132 and P = 9.2 × 10-69, respectively). For the most strongly associated Neandertal variant, we identify EPDR1 as the causal gene. Dupuytren's disease is an example of how admixture with Neandertals has shaped regional differences in disease prevalence.

Pathogenic Variants Associated with Rare Monogenic Diseases Established in Ancient Neanderthal and Denisovan Genome-Wide Data

Ancient anatomically modern humans (AMHs) encountered other archaic human species, most notably Neanderthals and Denisovans, when they left Africa and spread across Europe and Asia ~60,000 years ago. They interbred with them, and modern human genomes retain DNA inherited from these interbreeding events. High quality (high coverage) ancient human genomes have recently been sequenced allowing for a direct estimation of individual heterozygosity, which has shown that genetic diversity in these archaic human groups was very low, indicating low population sizes. In this study, we analyze ten ancient human genome-wide data, including four sequenced with high-coverage. We screened these ancient genome-wide data for pathogenic mutations associated with monogenic diseases, and established unusual aggregation of pathogenic mutations in individual subjects, including quadruple homozygous cases of pathogenic variants in the PAH gene associated with the condition phenylketonuria in a ~120,000 years old Neanderthal. Such aggregation of pathogenic mutations is extremely rare in contemporary populations, and their existence in ancient humans could be explained by less significant clinical manifestations coupled with small community sizes, leading to higher inbreeding levels. Our results suggest that pathogenic variants associated with rare diseases might be the result of introgression from other archaic human species, and archaic admixture thus could have influenced disease risk in modern humans.

Advancements and Challenges in Ancient DNA Research: Bridging the Global North-South Divide

Ancient DNA (aDNA) research first began in 1984 and ever since has greatly expanded our understanding of evolution and migration. Today, aDNA analysis is used to solve various puzzles about the origin of mankind, migration patterns, and the spread of infectious diseases. The incredible findings ranging from identifying the new branches within the human family to studying the genomes of extinct flora and fauna have caught the world by surprise in recent times. However, a closer look at these published results points out a clear Global North and Global South divide. Therefore, through this research, we aim to emphasize encouraging better collaborative opportunities and technology transfer to support researchers in the Global South. Further, the present research also focuses on expanding the scope of the ongoing conversation in the field of aDNA by reporting relevant literature published around the world and discussing the advancements and challenges in the field.

Reassessment of the human mandible from Banyoles (Girona, Spain)

Since the discovery of a human mandible in 1887 near the present-day city of Banyoles, northeastern Spain, researchers have generally emphasized its archaic features, including the lack of chin structures, and suggested affinities with the Neandertals or European Middle Pleistocene (Chibanian) specimens. Uranium-series and electron spin resonance dating suggest the mandible dates to the Late Pleistocene (Tarantian), approximately ca. 45-66 ka. In this study, we reassessed the taxonomic affinities of the Banyoles mandible by comparing it to samples of Middle Pleistocene fossils from Africa and Europe, Neandertals, Early and Upper Paleolithic modern humans, and recent modern humans. We evaluated the frequencies and expressions of morphological features and performed a three-dimensional geometric morphometric analysis on a virtual reconstruction of Banyoles to capture overall mandibular shape. Our results revealed no derived Neandertal morphological features in Banyoles. While a principal component analysis based on Euclidean distances from the first two principal components clearly grouped Banyoles with both fossil and recent Homo sapiens individuals, an analysis of the Procrustes residuals demonstrated that Banyoles did not fit into any of the comparative groups. The lack of Neandertal features in Banyoles is surprising considering its Late Pleistocene age. A consideration of the Middle Pleistocene fossil record in Europe and southwest Asia suggests that Banyoles is unlikely to represent a late-surviving Middle Pleistocene population. The lack of chin structures also complicates an assignment to H. sapiens, although early fossil H. sapiens do show somewhat variable development of the chin structures. Thus, Banyoles represents a non-Neandertal Late Pleistocene European individual and highlights the continuing signal of diversity in the hominin fossil record. The present situation makes Banyoles a prime candidate for ancient DNA or proteomic analyses, which may shed additional light on its taxonomic affinities.

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.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

Eurasiaplex-2: Shifting the focus to SNPs with high population specificity increases the power of forensic ancestry marker sets

To compile a new South Asian-informative panel of forensic ancestry SNPs, we changed the strategy for selecting the most powerful markers for this purpose by targeting polymorphisms with near absolute specificity - when the South Asian-informative allele identified is absent from all other populations or present at frequencies below 0.001 (one in a thousand). More than 120 candidate SNPs were identified from 1000 Genomes datasets satisfying an allele frequency screen of ≥ 0.1 (10 % or more) allele frequency in South Asians, and ≤ 0.001 (0.1 % or less) in African, East Asian, and European populations. From the candidate pool of markers, a final panel of 36 SNPs, widely distributed across most autosomes, were selected that had allele frequencies in the five 1000 Genomes South Asian populations ranging from 0.4 to 0.15. Slightly lower average allele frequencies, but consistent patterns of informativeness were observed in gnomAD South Asian datasets used to validate the 1000 Genomes variant annotations. We named the panel of 36 South Asian-specific SNPs Eurasiaplex-2, and the informativeness of the panel was evaluated by compiling worldwide population data from 4097 samples in four genome variation databases that largely complement the global sampling of 1000 Genomes. Consistent patterns of allele frequency distribution, which were specific to South Asia, were observed in all populations in, or closely sited to, the Indian sub-continent. Pakistani populations from the HGDP-CEPH panel had markedly lower allele frequencies, highlighting the need to develop a statistical system to evaluate the ancestry inference value of counting the number of population-specific alleles present in an individual.

Midfacial Morphology and Neandertal-Modern Human Interbreeding

Ancient DNA from, Neandertal and modern human fossils, and comparative morphological analyses of them, reveal a complex history of interbreeding between these lineages and the introgression of Neandertal genes into modern human genomes. Despite substantial increases in our knowledge of these events, the timing and geographic location of hybridization events remain unclear. Six measures of facial size and shape, from regional samples of Neandertals and early modern humans, were used in a multivariate exploratory analysis to try to identify regions in which early modern human facial morphology was more similar to that of Neandertals, which might thus represent regions of greater introgression of Neandertal genes. The results of canonical variates analysis and hierarchical cluster analysis suggest important affinities in facial morphology between both Middle and Upper Paleolithic early modern humans of the Near East with Neandertals, highlighting the importance of this region for interbreeding between the two lineages.

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.

Modern human incursion into Neanderthal territories 54,000 years ago at Mandrin, France

Determining the extent of overlap between modern humans and other hominins in Eurasia, such as Neanderthals and Denisovans, is fundamental to understanding the nature of their interactions and what led to the disappearance of archaic hominins. Apart from a possible sporadic pulse recorded in Greece during the Middle Pleistocene, the first settlements of modern humans in Europe have been constrained to ~45,000 to 43,000 years ago. Here, we report hominin fossils from Grotte Mandrin in France that reveal the earliest known presence of modern humans in Europe between 56,800 and 51,700 years ago. This early modern human incursion in the Rhône Valley is associated with technologies unknown in any industry of that age outside Africa or the Levant. Mandrin documents the first alternating occupation of Neanderthals and modern humans, with a modern human fossil and associated Neronian lithic industry found stratigraphically between layers containing Neanderthal remains associated with Mousterian industries.

Multi-ancestry fine mapping implicates OAS1 splicing in risk of severe COVID-19

The OAS1/2/3 cluster has been identified as a risk locus for severe COVID-19 among individuals of European ancestry, with a protective haplotype of approximately 75 kilobases (kb) derived from Neanderthals in the chromosomal region 12q24.13. This haplotype contains a splice variant of OAS1, which occurs in people of African ancestry independently of gene flow from Neanderthals. Using trans-ancestry fine-mapping approaches in 20,779 hospitalized cases, we demonstrate that this splice variant is likely to be the SNP responsible for the association at this locus, thus strongly implicating OAS1 as an effector gene influencing COVID-19 severity.

Lower promoter activity of the ST8SIA2 gene has been favored in evolving human collective brains

ST8SIA2 is an important molecule regulating expression of the phenotype involved in schizophrenia. Lowered promoter activity of the ST8SIA2 gene is considered to be protective against schizophrenia by conferring tolerance to psychosocial stress. Here, we examined the promoter-type composition of anatomically modern humans (AMHs) and archaic humans (AHs; Neanderthals and Denisovans), and compared the promoter activity at the population level (population promoter activity; PPA) between them. In AMHs, the TCT-type, showing the second lowest promoter activity, was most prevalent in the ancestral population of non-Africans. However, the detection of only the CGT-type from AH samples and recombination tracts in AH sequences showed that the CGT- and TGT-types, exhibiting the two highest promoter activities, were common in AH populations. Furthermore, interspecies gene flow occurred into AMHs from AHs and into Denisovans from Neanderthals, influencing promoter-type compositions independently in both AMHs and AHs. The difference of promoter-type composition makes PPA unique in each population. East and Southeast Asian populations show the lowest PPA. This results from the selective increase of the CGC-type, showing the lowest promoter activity, in these populations. Every non-African population shows significantly lower PPA than African populations, resulting from the TCT-type having the highest prevalence in the ancestral population of non-Africans. In addition, PPA reduction is also found among subpopulations within Africa via a slight increase of the TCT-type. These findings indicate a trend toward lower PPA in the spread of AMHs, interpreted as a continuous adaptation to psychosocial stress arising in migration. This trend is considered as genetic tuning for the evolution of collective brains. The inferred promoter-type composition of AHs differed markedly from that of AMHs, resulting in higher PPA in AHs than in AMHs. This suggests that the trend toward lower PPA is a unique feature in AMH spread.

Geometric morphometric variability in the supraorbital and orbital region of Middle Pleistocene hominins: Implications for the taxonomy and evolution of later Homo

This study assessed variation in the supraorbital and orbital region of the Middle Pleistocene hominins (MPHs), sometimes called Homo heidelbergensis s.l., to test whether it matched the expectations of intraspecific variation. The morphological distinctiveness and relative variation of this region, which is relatively well represented in the hominin fossil record, was analyzed quantitatively in a comparative taxonomic framework. Coordinates of 230 3D landmarks (20) and sliding semilandmarks (210) were collected from 704 specimens from species of Homo, Australopithecus, Paranthropus, Gorilla, Pan, Papio, and Macaca. Results showed that the MPHs had expected levels of morphological distinctiveness and intragroup and intergroup variation in supraorbital and orbital morphology, relative to commonly recognized non-hominin catarrhine species. However, the Procrustes distances between this group and H. sapiens were significantly higher than expected for two closely related catarrhine species. Furthermore, this study showed that variation within the MPH could be similarly well contained within existing hypodigms of H. sapiens, H. neanderthalensis, and H. erectus s.l. Although quantitative assessment of supraorbital and orbital morphology did not allow differentiation between taxonomic hypotheses in later Homo, it could be used to test individual taxonomic affiliation and identify potentially anomalous individuals. This study confirmed a complicated pattern of supraorbital and orbital morphology in the MPH fossil record and raises further questions over our understanding of the speciation of H. sapiens and H. neanderthalensis and taxonomic diversity in later Homo.

An Extended Admixture Pulse Model Reveals the Limitations to Human-Neandertal Introgression Dating

Neandertal DNA makes up 2-3% of the genomes of all non-African individuals. The patterns of Neandertal ancestry in modern humans have been used to estimate that this is the result of gene flow that occurred during the expansion of modern humans into Eurasia, but the precise dates of this event remain largely unknown. Here, we introduce an extended admixture pulse model that allows joint estimation of the timing and duration of gene flow. This model leads to simple expressions for both the admixture segment distribution and the decay curve of ancestry linkage disequilibrium, and we show that these two statistics are closely related. In simulations, we find that estimates of the mean time of admixture are largely robust to details in gene flow models, but that the duration of the gene flow can only be recovered if gene flow is very recent and the exact recombination map is known. These results imply that gene flow from Neandertals into modern humans could have happened over hundreds of generations. Ancient genomes from the time around the admixture event are thus likely required to resolve the question when, where, and for how long humans and Neandertals interacted.

An etiology of human modernity

Following the refutation of the replacement hypothesis, which had proposed that a ‘superior’ hominin species arose in Africa and replaced all other humans existing at the time, the auto-domestication hypothesis remains the only viable explanation for the relatively abrupt change from robust to gracile humans in the Late Pleistocene. It invokes the incidental institution of the domestication syndrome in humans, most probably by newly introduced cultural practices. It also postulates that the induction of exograms compensated for the atrophy of the brain caused by domestication. This new explanation of the origins of modernity in humans elucidates practically all its many aspects, in stark contrast to the superseded replacement hypothesis, which explained virtually nothing. The first results of the domestication syndrome’s genetic exploration have become available in recent years, and they endorse the human self-domestication hypothesis.

Local adaptation and archaic introgression shape global diversity at human structural variant loci

Large genomic insertions and deletions are a potent source of functional variation, but are challenging to resolve with short-read sequencing, limiting knowledge of the role of such structural variants (SVs) in human evolution. Here, we used a graph-based method to genotype long-read-discovered SVs in short-read data from diverse human genomes. We then applied an admixture-aware method to identify 220 SVs exhibiting extreme patterns of frequency differentiation - a signature of local adaptation. The top two variants traced to the immunoglobulin heavy chain locus, tagging a haplotype that swept to near fixation in certain southeast Asian populations, but is rare in other global populations. Further investigation revealed evidence that the haplotype traces to gene flow from Neanderthals, corroborating the role of immune-related genes as prominent targets of adaptive introgression. Our study demonstrates how recent technical advances can help resolve signatures of key evolutionary events that remained obscured within technically challenging regions of the genome.

Effects of hybridization on pelvic morphology: A macaque model

Ancient DNA analyses have shown that interbreeding between hominin taxa occurred multiple times. Although admixture is often reflected in skeletal phenotype, the relationship between the two remains poorly understood, hampering interpretation of the hominin fossil record. Direct study of this relationship is often impossible due to the paucity of hominin fossils and difficulties retrieving ancient genetic material. Here, we use a sample of known ancestry hybrids between two closely related nonhuman primate taxa (Indian and Chinese Macaca mulatta) to investigate the effect of admixture on skeletal morphology. We focus on pelvic shape, which has potential fitness implications in hybrids, as mismatches between maternal pelvic and fetal cranial morphology are often fatal to mother and offspring. As the pelvis is also one of the skeletal regions that differs most between Homo sapiens and Neanderthals, investigating the pelvic consequences of interbreeding could be informative regarding the viability of their hybrids. We find that the effect of admixture in M. mulatta is small and proportional to the relatively small morphological difference between the parent taxa. Sexual dimorphism appears to be the main determinant of pelvic shape in M. mulatta. The lack of difference in pelvic shape between Chinese and Indian M. mulatta is in contrast to that between Neanderthals and H. sapiens, despite a similar split time (in generations) between the hybridizing pairs. Greater phenotypic divergence between hominins may relate to adaptations to disparate environments but may also highlight how the unique degree of cultural buffering in hominins allowed for greater neutral divergence. In contrast to some previous work identifying extreme morphologies in first- and second-generation hybrids, here the relationship between pelvic shape and admixture is linear. This linearity may be because most sampled animals have a multigenerational admixture history or because of relatively high constraints on the pelvis compared with other skeletal regions.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.

Ancient DNA and multimethod dating confirm the late arrival of anatomically modern humans in southern China

The expansion of anatomically modern humans (AMHs) from Africa around 65,000 to 45,000 y ago (ca. 65 to 45 ka) led to the establishment of present-day non-African populations. Some paleoanthropologists have argued that fossil discoveries from Huanglong, Zhiren, Luna, and Fuyan caves in southern China indicate one or more prior dispersals, perhaps as early as ca. 120 ka. We investigated the age of the human remains from three of these localities and two additional early AMH sites (Yangjiapo and Sanyou caves, Hubei) by combining ancient DNA (aDNA) analysis with a multimethod geological dating strategy. Although U-Th dating of capping flowstones suggested they lie within the range ca. 168 to 70 ka, analyses of aDNA and direct AMS ¹⁴C dating on human teeth from Fuyan and Yangjiapo caves showed they derive from the Holocene. OSL dating of sediments and AMS ¹⁴C analysis of mammal teeth and charcoal also demonstrated major discrepancies from the flowstone ages; the difference between them being an order of magnitude or more at most of these localities. Our work highlights the surprisingly complex depositional history recorded at these subtropical caves which involved one or more episodes of erosion and redeposition or intrusion as recently as the late Holocene. In light of our findings, the first appearance datum for AMHs in southern China should probably lie within the timeframe set by molecular data of ca. 50 to 45 ka.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

Neanderthal-derived genetic variation in living humans relates to schizophrenia diagnosis, to psychotic symptom severity, and to dopamine synthesis

Schizophrenia has been hypothesized to be a human-specific condition, but experimental approaches to testing this idea have been limited. Because Neanderthals, our closest evolutionary relatives, interbred with modern humans prior to their disappearance from the fossil record, leaving a residual echo that survives in our DNA today, we leveraged new discoveries about ancient hominid DNA to explore this hypothesis in living people in three converging ways. First, in four independent case-control datasets totaling 9,362 individuals, individuals with schizophrenia had less Neanderthal-derived genetic variation than controls (p = .044). Second, in 49 unmedicated inpatients with schizophrenia, having more Neanderthal admixture predicted less severe positive symptoms (p = .046). Finally, using 18 F-fluorodopa PET scanning in 172 healthy individuals, having greater Neanderthal introgression was significantly associated with lower dopamine synthesis capacity in the striatum and pons (p's < 2 × 10-5 ), which is fundamentally important in the pathophysiology and treatment of psychosis. These results may help to elucidate the evolutionary history of a devastating neuropsychiatric disease by supporting the notion of schizophrenia as a human-specific condition. Additionally, the relationship between Neanderthal admixture and dopamine function suggests a potential mechanism whereby Neanderthal admixture may have affected our gene pool to alter schizophrenia risk and/or course.

Signatures of Introgression across the Allele Frequency Spectrum

The detection of introgression from genomic data is transforming our view of species and the origins of adaptive variation. Among the most widely used approaches to detect introgression is the so-called ABBA-BABA test or D-statistic, which identifies excess allele sharing between nonsister taxa. Part of the appeal of D is its simplicity, but this also limits its informativeness, particularly about the timing and direction of introgression. Here we present a simple extension, D frequency spectrum or DFS, in which D is partitioned according to the frequencies of derived alleles. We use simulations over a large parameter space to show how DFS carries information about various factors. In particular, recent introgression reliably leads to a peak in DFS among low-frequency derived alleles, whereas violation of model assumptions can lead to a lack of signal at low frequencies. We also reanalyze published empirical data from six different animal and plant taxa, and interpret the results in the light of our simulations, showing how DFS provides novel insights. We currently see DFS as a descriptive tool that will augment both simple and sophisticated tests for introgression, but in the future it may be usefully incorporated into probabilistic inference frameworks.