Paleogenomics. Genomic structure in Europeans dating back at least 36,200 years

Urbanization and genetic homogenization in the medieval Low Countries revealed through a ten-century paleogenomic study of the city of Sint-Truiden

BACKGROUND

Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes.

RESULTS

We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes.

CONCLUSIONS

This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

sedimix : A workflow for the analysis of hominin nuclear DNA sequences from sediments

Summary

Sediment DNA--the ability to extract DNA from archaeological sediments-- is an exciting new frontier in ancient DNA research, offering the potential to study individuals at a given archaeological site without destructive sampling. In recent years, several studies have demonstrated the promise of this approach by recovering hominin DNA from prehistoric sediments, including those dating back to the Middle or Late Pleistocene. However, a lack of open-source workflows for analysis of hominin sediment DNA samples poses a challenge for data processing and reproducibility of findings across studies. Here we introduce a snakemake workflow, sedimix , for processing genomic sequences from archaeological sediment DNA samples to identify hominin sequences and generate relevant summary statistics to assess the reliability of the pipeline. By performing simulations and comparing to published studies, we show that sedimix has high sensitivity and precision. sedimix offers a reliable and adaptable framework to aid in the analysis of sediment DNA datasets and improve reproducibility across studies.

Availability and Implementation

sedimix is available as an open-source software with the associated code and user manual available at https://github.com/jierui-cell/sedimix Contact: Jierui Xu ( jierui.xu@berkeley.edu ) Supplementary information: Supplementary data are available at Bioinformatics online.

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.

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.

Filtering out the noise: metagenomic classifiers optimize ancient DNA mapping

Contamination with exogenous DNA presents a significant challenge in ancient DNA (aDNA) studies of single organisms. Failure to address contamination from microbes, reagents, and present-day sources can impact the interpretation of results. Although field and laboratory protocols exist to limit contamination, there is still a need to accurately distinguish between endogenous and exogenous data computationally. Here, we propose a workflow to reduce exogenous contamination based on a metagenomic classifier. Unlike previous methods that relied exclusively on DNA sequencing reads mapping specificity to a single reference genome to remove contaminating reads, our approach uses Kraken2-based filtering before mapping to the reference genome. Using both simulated and empirical shotgun aDNA data, we show that this workflow presents a simple and efficient method that can be used in a wide range of computational environments-including personal machines. We propose strategies to build specific databases used to profile sequencing data that take into consideration available computational resources and prior knowledge about the target taxa and likely contaminants. Our workflow significantly reduces the overall computational resources required during the mapping process and reduces the total runtime by up to ~94%. The most significant impacts are observed in low endogenous samples. Importantly, contaminants that would map to the reference are filtered out using our strategy, reducing false positive alignments. We also show that our method results in a negligible loss of endogenous data with no measurable impact on downstream population genetics analyses.

The rise and transformation of Bronze Age pastoralists in the Caucasus

The Caucasus and surrounding areas, with their rich metal resources, became a crucible of the Bronze Age1 and the birthplace of the earliest steppe pastoralist societies2. Yet, despite this region having a large influence on the subsequent development of Europe and Asia, questions remain regarding its hunter-gatherer past and its formation of expansionist mobile steppe societies3-5. Here we present new genome-wide data for 131 individuals from 38 archaeological sites spanning 6,000 years. We find a strong genetic differentiation between populations north and south of the Caucasus mountains during the Mesolithic, with Eastern hunter-gatherer ancestry4,6 in the north, and a distinct Caucasus hunter-gatherer ancestry7 with increasing East Anatolian farmer admixture in the south. During the subsequent Eneolithic period, we observe the formation of the characteristic West Eurasian steppe ancestry and heightened interaction between the mountain and steppe regions, facilitated by technological developments of the Maykop cultural complex8. By contrast, the peak of pastoralist activities and territorial expansions during the Early and Middle Bronze Age is characterized by long-term genetic stability. The Late Bronze Age marks another period of gene flow from multiple distinct sources that coincides with a decline of steppe cultures, followed by a transformation and absorption of the steppe ancestry into highland populations.

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.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Investigating food production-associated DNA methylation changes in paleogenomes: Lack of consistent signals beyond technical noise

The Neolithic transition introduced major diet and lifestyle changes to human populations across continents. Beyond well-documented bioarcheological and genetic effects, whether these changes also had molecular-level epigenetic repercussions in past human populations has been an open question. In fact, methylation signatures can be inferred from UDG-treated ancient DNA through postmortem damage patterns, but with low signal-to-noise ratios; it is thus unclear whether published paleogenomes would provide the necessary resolution to discover systematic effects of lifestyle and diet shifts. To address this we compiled UDG-treated shotgun genomes of 13 pre-Neolithic hunter-gatherers (HGs) and 21 Neolithic farmers (NFs) individuals from West and North Eurasia, published by six different laboratories and with coverage c.1×-58× (median = 9×). We used epiPALEOMIX and a Monte Carlo normalization scheme to estimate methylation levels per genome. Our paleomethylome dataset showed expected genome-wide methylation patterns such as CpG island hypomethylation. However, analyzing the data using various approaches did not yield any systematic signals for subsistence type, genetic sex, or tissue effects. Comparing the HG-NF methylation differences in our dataset with methylation differences between hunter-gatherers versus farmers in modern-day Central Africa also did not yield consistent results. Meanwhile, paleomethylome profiles did cluster strongly by their laboratories of origin. Using larger data volumes, minimizing technical noise and/or using alternative protocols may be necessary for capturing subtle environment-related biological signals from paleomethylomes.

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.

The Persian plateau served as hub for Homo sapiens after the main out of Africa dispersal

A combination of evidence, based on genetic, fossil and archaeological findings, indicates that Homo sapiens spread out of Africa between ~70-60 thousand years ago (kya). However, it appears that once outside of Africa, human populations did not expand across all of Eurasia until ~45 kya. The geographic whereabouts of these early settlers in the timeframe between ~70-60 to 45 kya has been difficult to reconcile. Here we combine genetic evidence and palaeoecological models to infer the geographic location that acted as the Hub for our species during the early phases of colonisation of Eurasia. Leveraging on available genomic evidence we show that populations from the Persian Plateau carry an ancestry component that closely matches the population that settled the Hub outside Africa. With the paleoclimatic data available to date, we built ecological models showing that the Persian Plateau was suitable for human occupation and that it could sustain a larger population compared to other West Asian regions, strengthening this claim.

Ancient genomes illuminate Eastern Arabian population history and adaptation against malaria

The harsh climate of Arabia has posed challenges in generating ancient DNA from the region, hindering the direct examination of ancient genomes for understanding the demographic processes that shaped Arabian populations. In this study, we report whole-genome sequence data obtained from four Tylos-period individuals from Bahrain. Their genetic ancestry can be modeled as a mixture of sources from ancient Anatolia, Levant, and Iran/Caucasus, with variation between individuals suggesting population heterogeneity in Bahrain before the onset of Islam. We identify the G6PD Mediterranean mutation associated with malaria resistance in three out of four ancient Bahraini samples and estimate that it rose in frequency in Eastern Arabia from 5 to 6 kya onward, around the time agriculture appeared in the region. Our study characterizes the genetic composition of ancient Arabians, shedding light on the population history of Bahrain and demonstrating the feasibility of studies of ancient DNA in the region.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

Chronological and genetic analysis of an Upper Palaeolithic female infant burial from Borsuka Cave, Poland

Six infant human teeth and 112 animal tooth pendants from Borsuka Cave were identified as the oldest burial in Poland. However, uncertainties around the dating and the association of the teeth to the pendants have precluded their association with an Upper Palaeolithic archaeological industry. Using <67 mg per tooth, we combined dating and genetic analyses of two human teeth and six herbivore tooth pendants to address these questions. Our interdisciplinary approach yielded informative results despite limited sampling material, and high levels of degradation and contamination. Our results confirm the Palaeolithic origin of the human remains and herbivore pendants, and permit us to identify the infant as female and discuss the association of the assemblage with different Palaeolithic industries. This study exemplifies the progress that has been made toward minimally destructive methods and the benefits of integrating methods to maximize data retrieval from precious but highly degraded and contaminated prehistoric material.

Genome sequences of 36,000- to 37,000-year-old modern humans at Buran-Kaya III in Crimea

Populations genetically related to present-day Europeans first appeared in Europe at some point after 38,000-40,000 years ago, following a cold period of severe climatic disruption. These new migrants would eventually replace the pre-existing modern human ancestries in Europe, but initial interactions between these groups are unclear due to the lack of genomic evidence from the earliest periods of the migration. Here we describe the genomes of two 36,000-37,000-year-old individuals from Buran-Kaya III in Crimea as belonging to this newer migration. Both genomes share the highest similarity to Gravettian-associated individuals found several thousand years later in southwestern Europe. These genomes also revealed that the population turnover in Europe after 40,000 years ago was accompanied by admixture with pre-existing modern human populations. European ancestry before 40,000 years ago persisted not only at Buran-Kaya III but is also found in later Gravettian-associated populations of western Europe and Mesolithic Caucasus populations.

Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes

f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.

Genetic continuity, isolation, and gene flow in Stone Age Central and Eastern Europe

The genomic landscape of Stone Age Europe was shaped by multiple migratory waves and population replacements, but different regions do not all show similar patterns. To refine our understanding of the population dynamics before and after the dawn of the Neolithic, we generated and analyzed genomic sequence data from human remains of 56 individuals from the Mesolithic, Neolithic, and Eneolithic across Central and Eastern Europe. We found that Mesolithic European populations formed a geographically widespread isolation-by-distance zone ranging from Central Europe to Siberia, which was already established 10,000 years ago. We found contrasting patterns of population continuity during the Neolithic transition: people around the lower Dnipro Valley region, Ukraine, showed continuity over 4000 years, from the Mesolithic to the end of the Neolithic, in contrast to almost all other parts of Europe where population turnover drove this cultural change, including vast areas of Central Europe and around the Danube River.

Early presence of Homo sapiens in Southeast Asia by 86-68 kyr at Tam Pà Ling, Northern Laos

The timing of the first arrival of Homo sapiens in East Asia from Africa and the degree to which they interbred with or replaced local archaic populations is controversial. Previous discoveries from Tam Pà Ling cave (Laos) identified H. sapiens in Southeast Asia by at least 46 kyr. We report on a recently discovered frontal bone (TPL 6) and tibial fragment (TPL 7) found in the deepest layers of TPL. Bayesian modeling of luminescence dating of sediments and U-series and combined U-series-ESR dating of mammalian teeth reveals a depositional sequence spanning ~86 kyr. TPL 6 confirms the presence of H. sapiens by 70 ± 3 kyr, and TPL 7 extends this range to 77 ± 9 kyr, supporting an early dispersal of H. sapiens into Southeast Asia. Geometric morphometric analyses of TPL 6 suggest descent from a gracile immigrant population rather than evolution from or admixture with local archaic populations.

A 23,000-year-old southern Iberian individual links human groups that lived in Western Europe before and after the Last Glacial Maximum

Human populations underwent range contractions during the Last Glacial Maximum (LGM) which had lasting and dramatic effects on their genetic variation. The genetic ancestry of individuals associated with the post-LGM Magdalenian technocomplex has been interpreted as being derived from groups associated with the pre-LGM Aurignacian. However, both these ancestries differ from that of central European individuals associated with the chronologically intermediate Gravettian. Thus, the genomic transition from pre- to post-LGM remains unclear also in western Europe, where we lack genomic data associated with the intermediate Solutrean, which spans the height of the LGM. Here we present genome-wide data from sites in Andalusia in southern Spain, including from a Solutrean-associated individual from Cueva del Malalmuerzo, directly dated to ~23,000 cal yr BP. The Malalmuerzo individual carried genetic ancestry that directly connects earlier Aurignacian-associated individuals with post-LGM Magdalenian-associated ancestry in western Europe. This scenario differs from Italy, where individuals associated with the transition from pre- and post-LGM carry different genetic ancestries. This suggests different dynamics in the proposed southern refugia of Ice Age Europe and posits Iberia as a potential refugium for western European pre-LGM ancestry. More, individuals from Cueva Ardales, which were thought to be of Palaeolithic origin, date younger than expected and, together with individuals from the Andalusian sites Caserones and Aguilillas, fall within the genetic variation of the Neolithic, Chalcolithic and Bronze Age individuals from southern Iberia.

Palaeogenomics of Upper Palaeolithic to Neolithic European hunter-gatherers

Modern humans have populated Europe for more than 45,000 years1,2. Our knowledge of the genetic relatedness and structure of ancient hunter-gatherers is however limited, owing to the scarceness and poor molecular preservation of human remains from that period3. Here we analyse 356 ancient hunter-gatherer genomes, including new genomic data for 116 individuals from 14 countries in western and central Eurasia, spanning between 35,000 and 5,000 years ago. We identify a genetic ancestry profile in individuals associated with Upper Palaeolithic Gravettian assemblages from western Europe that is distinct from contemporaneous groups related to this archaeological culture in central and southern Europe4, but resembles that of preceding individuals associated with the Aurignacian culture. This ancestry profile survived during the Last Glacial Maximum (25,000 to 19,000 years ago) in human populations from southwestern Europe associated with the Solutrean culture, and with the following Magdalenian culture that re-expanded northeastward after the Last Glacial Maximum. Conversely, we reveal a genetic turnover in southern Europe suggesting a local replacement of human groups around the time of the Last Glacial Maximum, accompanied by a north-to-south dispersal of populations associated with the Epigravettian culture. From at least 14,000 years ago, an ancestry related to this culture spread from the south across the rest of Europe, largely replacing the Magdalenian-associated gene pool. After a period of limited admixture that spanned the beginning of the Mesolithic, we find genetic interactions between western and eastern European hunter-gatherers, who were also characterized by marked differences in phenotypically relevant variants.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

Admixture has obscured signals of historical hard sweeps in humans

The role of natural selection in shaping biological diversity is an area of intense interest in modern biology. To date, studies of positive selection have primarily relied on genomic datasets from contemporary populations, which are susceptible to confounding factors associated with complex and often unknown aspects of population history. In particular, admixture between diverged populations can distort or hide prior selection events in modern genomes, though this process is not explicitly accounted for in most selection studies despite its apparent ubiquity in humans and other species. Through analyses of ancient and modern human genomes, we show that previously reported Holocene-era admixture has masked more than 50 historic hard sweeps in modern European genomes. Our results imply that this canonical mode of selection has probably been underappreciated in the evolutionary history of humans and suggest that our current understanding of the tempo and mode of selection in natural populations may be inaccurate.

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.

Ancient Yersinia pestis and Salmonella enterica genomes from Bronze Age Crete

During the late 3^rd millennium BCE, the Eastern Mediterranean and Near East witnessed societal changes in many regions, which are usually explained with a combination of social and climatic factors.^1-4 However, recent archaeogenetic research forces us to rethink models regarding the role of infectious diseases in past societal trajectories.⁵ The plague bacterium Yersinia pestis, which was involved in some of the most destructive historical pandemics,^5-8 circulated across Eurasia at least from the onset of the 3^rd millennium BCE,^9-13 but the challenging preservation of ancient DNA in warmer climates has restricted the identification of Y.pestis from this period to temperate climatic regions. As such, evidence from culturally prominent regions such as the Eastern Mediterranean is currently lacking. Here, we present genetic evidence for the presence of Y. pestis and Salmonella enterica, the causative agent of typhoid/enteric fever, from this period of transformation in Crete, detected at the cave site Hagios Charalambos. We reconstructed one Y. pestis genome that forms part of a now-extinct lineage of Y. pestis strains from the Late Neolithic and Bronze Age that were likely not yet adapted for transmission via fleas. Furthermore, we reconstructed two ancient S. enterica genomes from the Para C lineage, which cluster with contemporary strains that were likely not yet fully host adapted to humans. The occurrence of these two virulent pathogens at the end of the Early Minoan period in Crete emphasizes the necessity to re-introduce infectious diseases as an additional factor possibly contributing to the transformation of early complex societies in the Aegean and beyond.

Ancient Maltese genomes and the genetic geography of Neolithic Europe

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

An integrative skeletal and paleogenomic analysis of stature variation suggests relatively reduced health for early European farmers

Human culture, biology, and health were shaped dramatically by the onset of agriculture ∼12,000 y B.P. This shift is hypothesized to have resulted in increased individual fitness and population growth as evidenced by archaeological and population genomic data alongside a decline in physiological health as inferred from skeletal remains. Here, we consider osteological and ancient DNA data from the same prehistoric individuals to study human stature variation as a proxy for health across a transition to agriculture. Specifically, we compared “predicted” genetic contributions to height from paleogenomic data and “achieved” adult osteological height estimated from long bone measurements for 167 individuals across Europe spanning the Upper Paleolithic to Iron Age (∼38,000 to 2,400 B.P.). We found that individuals from the Neolithic were shorter than expected (given their individual polygenic height scores) by an average of −3.82 cm relative to individuals from the Upper Paleolithic and Mesolithic (P = 0.040) and −2.21 cm shorter relative to post-Neolithic individuals (P = 0.068), with osteological vs. expected stature steadily increasing across the Copper (+1.95 cm relative to the Neolithic), Bronze (+2.70 cm), and Iron (+3.27 cm) Ages. These results were attenuated when we additionally accounted for genome-wide genetic ancestry variation: for example, with Neolithic individuals −2.82 cm shorter than expected on average relative to pre-Neolithic individuals (P = 0.120). We also incorporated observations of paleopathological indicators of nonspecific stress that can persist from childhood to adulthood in skeletal remains into our model. Overall, our work highlights the potential of integrating disparate datasets to explore proxies of health in prehistory.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

Human BRCA pathogenic variants were originated during recent human history

BRCA1 and BRCA2 (BRCA) play essential roles in maintaining genome stability. BRCA germline pathogenic variants increase cancer risk. However, the evolutionary origin of human BRCA pathogenic variants remains largely elusive. We tested the 2,972 human BRCA1 and 3,652 human BRCA2 pathogenic variants from ClinVar database in 100 vertebrates across eight clades, but failed to find evidence to show cross-species evolution conservation as the origin; we searched the variants in 2,792 ancient human genome data, and identified 28 BRCA1 and 22 BRCA2 pathogenic variants in 44 cases dated from 45,000 to 300 yr ago; we analyzed the haplotype-dated human BRCA pathogenic founder variants, and observed that they were mostly arisen within the past 3,000 yr; we traced ethnic distribution of human BRCA pathogenic variants, and found that the majority were present in single or a few ethnic populations. Based on the data, we propose that human BRCA pathogenic variants were highly likely arisen in recent human history after the latest out-of-Africa migration, and the expansion of modern human population could largely increase the variation spectrum.

Refining models of archaic admixture in Eurasia with ArchaicSeeker 2.0

We developed a method, ArchaicSeeker 2.0, to identify introgressed hominin sequences and model multiple-wave admixture. The new method enabled us to discern two waves of introgression from both Denisovan-like and Neanderthal-like hominins in present-day Eurasian populations and an ancient Siberian individual. We estimated that an early Denisovan-like introgression occurred in Eurasia around 118.8-94.0 thousand years ago (kya). In contrast, we detected only one single episode of Denisovan-like admixture in indigenous peoples eastern to the Wallace-Line. Modeling ancient admixtures suggested an early dispersal of modern humans throughout Asia before the Toba volcanic super-eruption 74 kya, predating the initial peopling of Asia as proposed by the traditional Out-of-Africa model. Survived archaic sequences are involved in various phenotypes including immune and body mass (e.g., ZNF169), cardiovascular and lung function (e.g., HHAT), UV response and carbohydrate metabolism (e.g., HYAL1/HYAL2/HYAL3), while "archaic deserts" are enriched with genes associated with skin development and keratinization.

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.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Insights into human history from the first decade of ancient human genomics

[Figure: see text].

Landscape genetics and the genetic legacy of Upper Paleolithic and Mesolithic hunter-gatherers in the modern Caucasus

This study clarifies the role of refugia and landscape permeability in the formation of the current genetic structure of peoples of the Caucasus. We report novel genome-wide data for modern individuals from the Caucasus, and analyze them together with available Paleolithic and Mesolithic individuals from Eurasia and Africa in order (1) to link the current and ancient genetic structures via landscape permeability, and (2) thus to identify movement paths between the ancient refugial populations and the Caucasus. The ancient genetic ancestry is best explained by landscape permeability implying that human movement is impeded by terrain ruggedness, swamps, glaciers and desert. Major refugial source populations for the modern Caucasus are those of the Caucasus, Anatolia, the Balkans and Siberia. In Rugged areas new genetic signatures take a long time to form, but once they do so, they remain for a long time. These areas act as time capsules harboring genetic signatures of ancient source populations and making it possible to help reconstruct human history based on patterns of variation today.

Different historical generation intervals in human populations inferred from Neanderthal fragment lengths and mutation signatures

After the main Out-of-Africa event, humans interbred with Neanderthals leaving 1–2% of Neanderthal DNA scattered in small fragments in all non-African genomes today. Here we investigate what can be learned about human demographic processes from the size distribution of these fragments. We observe differences in fragment length across Eurasia with 12% longer fragments in East Asians than West Eurasians. Comparisons between extant populations with ancient samples show that these differences are caused by different rates of decay in length by recombination since the Neanderthal admixture. In concordance, we observe a strong correlation between the average fragment length and the mutation accumulation, similar to what is expected by changing the ages at reproduction as estimated from trio studies. Altogether, our results suggest differences in the generation interval across Eurasia, by up 10–20%, over the past 40,000 years. We use sex-specific mutation signatures to infer whether these changes were driven by shifts in either male or female age at reproduction, or both. We also find that previously reported variation in the mutational spectrum may be largely explained by changes to the generation interval. We conclude that Neanderthal fragment lengths provide unique insight into differences among human populations over recent history.

Historical interbreeding between Neanderthals and humans should leave signatures of historical demographics in modern human genomes. Analysing the size distribution of Neanderthal fragments in non-African genomes suggests consistent differences in the generation interval across Eurasia, and that this could explain mutational spectrum variation.

Where Asia meets Europe - recent insights from ancient human genomics

CONTEXT

The peopling of Europe by modern humans is a widely debated topic in the field of modern and ancient genomics. While several recent syntheses have focussed on this topic, little has been discussed about the genetic history of populations in the continent's surrounding regions.

OBJECTIVE

We explore genetic transformations in three key areas that played an essential role in the formation of the European genetic landscape through time, focussing on the periods spanning from the Epipalaeolithic/Mesolithic and up until the Iron Age.

METHODS

We review published ancient genomic studies and integrate the associated data to provide a quantification and visualisation of major trends in the population histories of the Near East, the western Eurasian Steppe and North East Europe.

RESULTS

We describe cross-regional as well as localised prehistoric demographic shifts and discuss potential research directions while highlighting geo-temporal gaps in the data.

CONCLUSION

In recent years, archaeogenetic studies have contributed to the understanding of human genetic diversity through time in regions located at the doorstep of Europe. Further studies focussing on these areas will allow for a better characterisation of genetic shifts and regionally-specific patterns of admixture across western Eurasia.

Paleogenomics of the prehistory of Europe: human migrations, domestication and disease

A substantial portion of ancient DNA research has been centred on understanding European populations' origin and evolution. A rchaeological evidence has already shown that the peopling of Europe involved an intricate pattern of demic and/or cultural diffusion since the Upper Palaeolithic, which became more evident during the Neolithic and Bronze Age periods. However, ancient DNA data has been crucial in determining if cultural changes occurred due to the movement of ideas or people. With the advent of next-generation sequencing and population-based paleogenomic research, ancient DNA studies have been directed not only at the study of continental human migrations, but also to the detailed analysis of particular archaeological sites, the processes of domestication, or the spread of disease during prehistoric times. With this vast paleogenomic effort added to a proper archaeological contextualisation of results, a deeper understanding of Europe's peopling is starting to emanate.

Inferring Population Histories for Ancient Genomes Using Genome-Wide Genealogies

Ancient genomes anchor genealogies in directly observed historical genetic variation and contextualize ancestral lineages with archaeological insights into their geography and cultural associations. However, the majority of ancient genomes are of lower coverage and cannot be directly built into genealogies. Here, we present a fast and scalable method, Colate, the first approach for inferring ancestral relationships through time between low-coverage genomes without requiring phasing or imputation. Our approach leverages sharing patterns of mutations dated using a genealogy to infer coalescence rates. For deeply sequenced ancient genomes, we additionally introduce an extension of the Relate algorithm for joint inference of genealogies incorporating such genomes. Application to 278 present-day and 430 ancient DNA samples of >0.5x mean coverage allows us to identify dynamic population structure and directional gene flow between early farmer and European hunter-gatherer groups. We further show that the previously reported, but still unexplained, increase in the TCC/TTC mutation rate, which is strongest in West Eurasia today, was already present at similar strength and widespread in the Late Glacial Period ~10k-15k years ago, but is not observed in samples >30k years old. It is strongest in Neolithic farmers, and highly correlated with recent coalescence rates between other genomes and a 10,000-year-old Anatolian hunter-gatherer. This suggests gene-flow among ancient peoples postdating the last glacial maximum as widespread and localizes the driver of this mutational signal in both time and geography in that region. Our approach should be widely applicable in future for addressing other evolutionary questions, and in other species.

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.

Genetic insights into the paternal admixture history of Chinese Mongolians via high-resolution customized Y-SNP SNaPshot panels

The Mongolian people, one of the Mongolic-speaking populations, are native to the Mongolian Plateau in North China and southern Siberia. Many ancient DNA studies recently reported extensive population transformations during the Paleolithic to historic periods in this region, while little is known about the paternal genetic legacy of modern geographically different Mongolians. Here, we genotyped 215 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) and 37 Y-chromosomal short tandem repeats (Y-STRs) among 679 Mongolian individuals from Hohhot, Hulunbuir, and Ordos in North China using the AGCU Y37 kit and our developed eight Y-SNP SNaPshot panels (including two panels first reported herein). The C-M130 Y-SNP SNaPshot panel defines 28 subhaplogroups, and the N/O/Q complementary Y-SNP SNaPshot panel defines 30 subhaplogroups of N1b-F2930, N1a1a1a1a3-B197, Q-M242, and O2a2b1a1a1a4a-CTS4658, which improved the resolution our developed Y-SNP SNaPshot panel set and could be applied for dissecting the finer-scale paternal lineages of Mongolic speakers. We found a strong association between Mongolian-prevailing haplogroups and some observed microvariants among the newly generated Y-STR haplotype data, suggesting the possibility of haplogroup prediction based on the distribution of Y-STR haplotypes. We identified three main ancestral sources of the observed Mongolian-dominant haplogroups, including the local lineage of C2*-M217 and incoming lineages from other regions of southern East Asia (O2*-M122, O1b*-P31, and N1*-CTS3750) and western Eurasia (R1*-M173). We also observed DE-M145, D1*-M174, C1*-F3393, G*-M201, I-M170, J*-M304, L-M20, O1a*-M119, and Q*-M242 at relatively low frequencies (< 5.00%), suggesting a complex admixture history between Mongolians and other incoming Eurasians from surrounding regions. Genetic clustering analyses indicated that the studied Mongolians showed close genetic affinities with other Altaic-speaking populations and Sinitic-speaking Hui people. The Y-SNP haplotype/haplogroup-based genetic legacy not only revealed that the stratification among geographically/linguistically/ethnically different Chinese populations was highly consistent with the geographical division and language classification, but also demonstrated that patrilineal genetic materials could provide fine-scale genetic structures among geographically different Mongolian people, suggesting that our developed high-resolution Y-SNP SNaPshot panels have the potential for forensic pedigree searches and biogeographical ancestry inference.

Human mobility at Tell Atchana (Alalakh), Hatay, Turkey during the 2nd millennium BC: Integration of isotopic and genomic evidence

The Middle and Late Bronze Age, a period roughly spanning the 2nd millennium BC (ca. 2000-1200 BC) in the Near East, is frequently referred to as the first 'international age', characterized by intense and far-reaching contacts between different entities from the eastern Mediterranean to the Near East and beyond. In a large-scale tandem study of stable isotopes and ancient DNA of individuals excavated at Tell Atchana (Alalakh, located in Hatay, Turkey), we explored the role of mobility at the capital of a regional kingdom, named Mukish during the Late Bronze Age, which spanned the Amuq Valley and some areas beyond. We generated strontium and oxygen isotope data from dental enamel for 53 individuals and 77 individuals, respectively, and added ancient DNA data of 10 newly sequenced individuals to a dataset of 27 individuals published in 2020. Additionally, we improved the DNA coverage of one individual from this 2020 dataset. The DNA data revealed a very homogeneous gene pool. This picture of an overwhelmingly local ancestry was consistent with the evidence of local upbringing in most of the individuals indicated by the isotopic data, where only five were found to be non-local. High levels of contact, trade, and exchange of ideas and goods in the Middle and Late Bronze Ages, therefore, seem not to have translated into high levels of individual mobility detectable at Tell Atchana.

Human population history at the crossroads of East and Southeast Asia since 11,000 years ago

Past human genetic diversity and migration between southern China and Southeast Asia have not been well characterized, in part due to poor preservation of ancient DNA in hot and humid regions. We sequenced 31 ancient genomes from southern China (Guangxi and Fujian), including two ∼12,000- to 10,000-year-old individuals representing the oldest humans sequenced from southern China. We discovered a deeply diverged East Asian ancestry in the Guangxi region that persisted until at least 6,000 years ago. We found that ∼9,000- to 6,000-year-old Guangxi populations were a mixture of local ancestry, southern ancestry previously sampled in Fujian, and deep Asian ancestry related to Southeast Asian Hòabìnhian hunter-gatherers, showing broad admixture in the region predating the appearance of farming. Historical Guangxi populations dating to ∼1,500 to 500 years ago are closely related to Tai-Kadai and Hmong-Mien speakers. Our results show heavy interactions among three distinct ancestries at the crossroads of East and Southeast Asia.

The deep population history of northern East Asia from the Late Pleistocene to the Holocene

Northern East Asia was inhabited by modern humans as early as 40 thousand years ago (ka), as demonstrated by the Tianyuan individual. Using genome-wide data obtained from 25 individuals dated to 33.6-3.4 ka from the Amur region, we show that Tianyuan-related ancestry was widespread in northern East Asia before the Last Glacial Maximum (LGM). At the close of the LGM stadial, the earliest northern East Asian appeared in the Amur region, and this population is basal to ancient northern East Asians. Human populations in the Amur region have maintained genetic continuity from 14 ka, and these early inhabitants represent the closest East Asian source known for Ancient Paleo-Siberians. We also observed that EDAR V370A was likely to have been elevated to high frequency after the LGM, suggesting the possible timing for its selection. This study provides a deep look into the population dynamics of northern East Asia.

Genome of Peştera Muierii skull shows high diversity and low mutational load in pre-glacial Europe

Few complete human genomes from the European Early Upper Palaeolithic (EUP) have been sequenced. Using novel sampling and DNA extraction approaches, we sequenced the genome of a woman from "Peştera Muierii," Romania who lived ∼34,000 years ago to 13.5× coverage. The genome shows similarities to modern-day Europeans, but she is not a direct ancestor. Although her cranium exhibits both modern human and Neanderthal features, the genome shows similar levels of Neanderthal admixture (∼3.1%) to most EUP humans but only half compared to the ∼40,000-year-old Peştera Oase 1. All EUP European hunter-gatherers display high genetic diversity, demonstrating that the severe loss of diversity occurred during and after the Last Glacial Maximum (LGM) rather than just during the out-of-Africa migration. The prevalence of genetic diseases is expected to increase with low diversity; however, pathogenic variant load was relatively constant from EUP to modern times, despite post-LGM hunter-gatherers having the lowest diversity ever observed among Europeans.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

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.

Widespread Denisovan ancestry in Island Southeast Asia but no evidence of substantial super-archaic hominin admixture

The hominin fossil record of Island Southeast Asia (ISEA) indicates that at least two endemic 'super-archaic' species-Homo luzonensis and H. floresiensis-were present around the time anatomically modern humans arrived in the region >50,000 years ago. Intriguingly, contemporary human populations across ISEA carry distinct genomic traces of ancient interbreeding events with Denisovans-a separate hominin lineage that currently lacks a fossil record in ISEA. To query this apparent disparity between fossil and genetic evidence, we performed a comprehensive search for super-archaic introgression in >400 modern human genomes, including >200 from ISEA. Our results corroborate widespread Denisovan ancestry in ISEA populations, but fail to detect any substantial super-archaic admixture signals compatible with the endemic fossil record of ISEA. We discuss the implications of our findings for the understanding of hominin history in ISEA, including future research directions that might help to unlock more details about the prehistory of the enigmatic Denisovans.

Systematic benchmark of ancient DNA read mapping

The current standard practice for assembling individual genomes involves mapping millions of short DNA sequences (also known as DNA 'reads') against a pre-constructed reference genome. Mapping vast amounts of short reads in a timely manner is a computationally challenging task that inevitably produces artefacts, including biases against alleles not found in the reference genome. This reference bias and other mapping artefacts are expected to be exacerbated in ancient DNA (aDNA) studies, which rely on the analysis of low quantities of damaged and very short DNA fragments (~30-80 bp). Nevertheless, the current gold-standard mapping strategies for aDNA studies have effectively remained unchanged for nearly a decade, during which time new software has emerged. In this study, we used simulated aDNA reads from three different human populations to benchmark the performance of 30 distinct mapping strategies implemented across four different read mapping software-BWA-aln, BWA-mem, NovoAlign and Bowtie2-and quantified the impact of reference bias in downstream population genetic analyses. We show that specific NovoAlign, BWA-aln and BWA-mem parameterizations achieve high mapping precision with low levels of reference bias, particularly after filtering out reads with low mapping qualities. However, unbiased NovoAlign results required the use of an IUPAC reference genome. While relevant only to aDNA projects where reference population data are available, the benefit of using an IUPAC reference demonstrates the value of incorporating population genetic information into the aDNA mapping process, echoing recent results based on graph genome representations.

A genome sequence from a modern human skull over 45,000 years old from Zlatý kůň in Czechia

Modern humans expanded into Eurasia more than 40,000 years ago following their dispersal out of Africa. These Eurasians carried ~2–3% Neanderthal ancestry in their genomes, originating from admixture with Neanderthals that took place sometime between 50,000 and 60,000 years ago, probably in the Middle East. In Europe, the modern human expansion preceded the disappearance of Neanderthals from the fossil record by 3,000–5,000 years. The genetic makeup of the first Europeans who colonized the continent more than 40,000 years ago remains poorly understood since few specimens have been studied. Here, we analyse a genome generated from the skull of a female individual from Zlatý kůň, Czechia. We found that she belonged to a population that appears to have contributed genetically neither to later Europeans nor to Asians. Her genome carries ~3% Neanderthal ancestry, similar to those of other Upper Palaeolithic hunter-gatherers. However, the lengths of the Neanderthal segments are longer than those observed in the currently oldest modern human genome of the ~45,000-year-old Ust’-Ishim individual from Siberia, suggesting that this individual from Zlatý kůň is one of the earliest Eurasian inhabitants following the expansion out of Africa.

The authors present the genome sequence of a >45,000-year-old female Homo sapiens individual from the site of Zlatý kůň, Czechia. Although radiometric dating of the human remains was inconclusive, the authors were able to use molecular methods to demonstrate that she was probably among the earliest Eurasian inhabitants following expansion out of Africa.

Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.