Ancient DNA from Chalcolithic Israel reveals the role of population mixture in cultural transformation

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Reconstructing the Genetic Relationship between Ancient and Present-Day Siberian Populations

Human populations across a vast area in northern Eurasia, from Fennoscandia to Chukotka, share a distinct genetic component often referred to as the Siberian ancestry. Most enriched in present-day Samoyedic-speaking populations such as Nganasans, its origins and history still remain elusive despite the growing list of ancient and present-day genomes from Siberia. Here, we reanalyze published ancient and present-day Siberian genomes focusing on the Baikal and Yakutia, resolving key questions regarding their genetic history. First, we show a long-term presence of a unique genetic profile in southern Siberia, up to 6,000 yr ago, which distinctly shares a deep ancestral connection with Native Americans. Second, we provide plausible historical models tracing genetic changes in West Baikal and Yakutia in fine resolution. Third, the Middle Neolithic individual from Yakutia, belonging to the Belkachi culture, serves as the best source so far available for the spread of the Siberian ancestry into Fennoscandia and Greenland. These findings shed light on the genetic legacy of the Siberian ancestry and provide insights into the complex interplay between different populations in northern Eurasia throughout history.

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.

Benchmarking software tools for trimming adapters and merging next-generation sequencing data for ancient DNA

Ancient DNA is highly degraded, resulting in very short sequences. Reads generated with modern high-throughput sequencing machines are generally longer than ancient DNA molecules, therefore the reads often contain some portion of the sequencing adaptors. It is crucial to remove those adaptors, as they can interfere with downstream analysis. Furthermore, overlapping portions when DNA has been read forward and backward (paired-end) can be merged to correct sequencing errors and improve read quality. Several tools have been developed for adapter trimming and read merging, however, no one has attempted to evaluate their accuracy and evaluate their potential impact on downstream analyses. Through the simulation of sequencing data, seven commonly used tools were analyzed in their ability to reconstruct ancient DNA sequences through read merging. The analyzed tools exhibit notable differences in their abilities to correct sequence errors and identify the correct read overlap, but the most substantial difference is observed in their ability to calculate quality scores for merged bases. Selecting the most appropriate tool for a given project depends on several factors, although some tools such as fastp have some shortcomings, whereas others like leeHom outperform the other tools in most aspects. While the choice of tool did not result in a measurable difference when analyzing population genetics using principal component analysis, it is important to note that downstream analyses that are sensitive to wrongly merged reads or that rely on quality scores can be significantly impacted by the choice of tool.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f 3 -statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P -values of qpAdm models, and f 3 -statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f 2 -statistics, demonstrating the improved demographic inference that could be achieved with advancements in f -statistic estimations.

False discovery rates of qpAdm-based screens for genetic admixture

Although a broad range of methods exists for reconstructing population history from genome-wide single nucleotide polymorphism data, just a few methods gained popularity in archaeogenetics: principal component analysis (PCA); ADMIXTURE, an algorithm that models individuals as mixtures of multiple ancestral sources represented by actual or inferred populations; formal tests for admixture such as f3-statistics and D/f4-statistics; and qpAdm, a tool for fitting two-component and more complex admixture models to groups or individuals. Despite their popularity in archaeogenetics, which is explained by modest computational requirements and ability to analyze data of various types and qualities, protocols relying on qpAdm that screen numerous alternative models of varying complexity and find "fitting" models (often considering both estimated admixture proportions and p-values as a composite criterion of model fit) remain untested on complex simulated population histories in the form of admixture graphs of random topology. We analyzed genotype data extracted from such simulations and tested various types of high-throughput qpAdm protocols ("rotating" and "non-rotating", with or without temporal stratification of target groups and proxy ancestry sources, and with or without a "model competition" step). We caution that high-throughput qpAdm protocols may be inappropriate for exploratory analyses in poorly studied regions/periods since their false discovery rates varied between 12% and 68% depending on the details of the protocol and on the amount and quality of simulated data (i.e., >12% of fitting two-way admixture models imply gene flows that were not simulated). We demonstrate that for reducing false discovery rates of qpAdm protocols to nearly 0% it is advisable to use large SNP sets with low missing data rates, the rotating qpAdm protocol with a strictly enforced rule that target groups do not pre-date their proxy sources, and an unsupervised ADMIXTURE analysis as a way to verify feasible qpAdm models. Our study has a number of limitations: for instance, these recommendations depend on the assumption that the underlying genetic history is a complex admixture graph and not a stepping-stone model.

Mandibular morphology and the Mesolithic-Neolithic transition in Westernmost Iberia

Neolithic farming and animal husbandry were first developed in the Near East ~ 10,000 BCE and expanded westwards, reaching westernmost Iberia no later than 5500 BCE. It resulted in major social, cultural, economic and dietary changes. Yet, the impact of this change on human mandibular morphology in Iberia is yet to be assessed, which is regrettable because mandible form is impacted by population history and diet. In this study we used Mesolithic to Chalcolithic Iberian samples to examine the impact of this transition on mandibular morphology. We also compared these samples with a Southern Levantine Chalcolithic population to assess their relationship. Lastly, we assessed dental wear to determine if the morphological differences identified were related to the material properties of the diet. We found differences between samples in mandibular shape but not size, which we attribute to contrasting population histories between Mesolithic and later populations. Some differences in the severity of dental wear were also found between Mesolithic and later Iberian samples, and smaller between the Mesolithic Iberians and southern Levantines. Little relationship was found between wear magnitude and mandibular shape. Altogether, our results show that the Mesolithic-Neolithic Iberian transition resulted in a meaningful change in mandibular morphology, which was likely driven more by population history than by dietary change.

Northwest African Neolithic initiated by migrants from Iberia and Levant

In northwestern Africa, lifestyle transitioned from foraging to food production around 7,400 years ago but what sparked that change remains unclear. Archaeological data support conflicting views: (1) that migrant European Neolithic farmers brought the new way of life to North Africa^1-3 or (2) that local hunter-gatherers adopted technological innovations^4,5. The latter view is also supported by archaeogenetic data⁶. Here we fill key chronological and archaeogenetic gaps for the Maghreb, from Epipalaeolithic to Middle Neolithic, by sequencing the genomes of nine individuals (to between 45.8- and 0.2-fold genome coverage). Notably, we trace 8,000 years of population continuity and isolation from the Upper Palaeolithic, via the Epipaleolithic, to some Maghrebi Neolithic farming groups. However, remains from the earliest Neolithic contexts showed mostly European Neolithic ancestry. We suggest that farming was introduced by European migrants and was then rapidly adopted by local groups. During the Middle Neolithic a new ancestry from the Levant appears in the Maghreb, coinciding with the arrival of pastoralism in the region, and all three ancestries blend together during the Late Neolithic. Our results show ancestry shifts in the Neolithization of northwestern Africa that probably mirrored a heterogeneous economic and cultural landscape, in a more multifaceted process than observed in other regions.

Genetic population structure of the Xiongnu Empire at imperial and local scales

The Xiongnu established the first nomadic imperial power, controlling the Eastern Eurasian steppe from ca. 200 BCE to 100 CE. Recent archaeogenetic studies identified extreme levels of genetic diversity across the empire, corroborating historical records of the Xiongnu Empire being multiethnic. However, it has remained unknown how this diversity was structured at the local community level or by sociopolitical status. To address this, we investigated aristocratic and local elite cemeteries at the western frontier of the empire. Analyzing genome-wide data from 18 individuals, we show that genetic diversity within these communities was comparable to the empire as a whole, and that high diversity was also observed within extended families. Genetic heterogeneity was highest among the lowest-status individuals, implying diverse origins, while higher-status individuals harbored less genetic diversity, suggesting that elite status and power was concentrated within specific subsets of the broader Xiongnu population.

Genomic analysis of indigenous goats in Southwest Asia reveals evidence of ancient adaptive introgression related to desert climate

Understanding how evolutionary pressures related to climate change have shaped the current genetic background of domestic animals is a fundamental pursuit of biology. Here, we generated whole-genome sequencing data from native goat populations in Iraq and Pakistan. Combined with previously published data on modern, ancient (Late Neolithic to Medieval periods), and wild Capra species worldwide, we explored the genetic population structure, ancestry components, and signatures of natural positive selection in native goat populations in Southwest Asia (SWA). Results revealed that the genetic structure of SWA goats was deeply influenced by gene flow from the eastern Mediterranean during the Chalcolithic period, which may reflect adaptation to gradual warming and aridity in the region. Furthermore, comparative genomic analysis revealed adaptive introgression of the KITLG locus from the Nubian ibex ( C. nubiana) into African and SWA goats. The frequency of the selected allele at this locus was significantly higher among goat populations located near northeastern Africa. These results provide new insights into the genetic composition and history of goat populations in the SWA region.

Methodological Changes in the Field of Paleogenetics

Paleogenetics has significantly changed since its inception almost forty years ago. Initially, molecular techniques available to the researchers offered minimal possibilities for ancient DNA analysis. The subsequent expansion of the scientific tool cabinet allowed for more remarkable achievements, combined has with the newfound popularity of this budding field of science. Finally, a breakthrough was made with the development of next-generation sequencing (NGS) technologies and the update of DNA isolation protocols, through which even very fragmented aDNA samples could be used to sequence whole genomes. In this paper, we review the achievements made thus far and compare the methodologies utilized in this field of science, discussing their benefits and challenges.

Spatial and temporal heterogeneity in human mobility patterns in Holocene Southwest Asia and the East Mediterranean

We present a spatiotemporal picture of human genetic diversity in Anatolia, Iran, Levant, South Caucasus, and the Aegean, a broad region that experienced the earliest Neolithic transition and the emergence of complex hierarchical societies. Combining 35 new ancient shotgun genomes with 382 ancient and 23 present-day published genomes, we found that genetic diversity within each region steadily increased through the Holocene. We further observed that the inferred sources of gene flow shifted in time. In the first half of the Holocene, Southwest Asian and the East Mediterranean populations homogenized among themselves. Starting with the Bronze Age, however, regional populations diverged from each other, most likely driven by gene flow from external sources, which we term "the expanding mobility model." Interestingly, this increase in inter-regional divergence can be captured by outgroup-f₃-based genetic distances, but not by the commonly used F_ST statistic, due to the sensitivity of F_ST, but not outgroup-f₃, to within-population diversity. Finally, we report a temporal trend of increasing male bias in admixture events through the Holocene.

4000-year-old hair from the Middle Nile highlights unusual ancient DNA degradation pattern and a potential source of early eastern Africa pastoralists

Petrous bones and teeth are the skeletal elements most often targeted by researchers for ancient DNA (aDNA) extraction, and the sources of the majority of previously published ancient African genomes. However, the high temperature environments that characterise much of Africa often lead to poor preservation of skeletal remains. Here, we successfully reconstruct and analyse genome-wide data from the naturally mummified hair of a 4000-year-old individual from Sudan in northeastern Africa, after failed attempts at DNA extraction from teeth, petrous, and cranium of this and other individuals from the Kadruka cemeteries. We find that hair DNA extracted with an established single-stranded library protocol is unusually enriched in ultra-short DNA molecules and exhibits substantial interior molecular damage. The aDNA was nonetheless amenable to genetic analyses, which revealed that the genome is genetically indistinguishable from that of early Neolithic eastern African pastoralists located 2500 kms away. Our findings are consistent with established models for the southward dispersal of Middle Nile Valley pastoral populations to the Rift Valley of eastern Africa, and provide a possible genetic source population for this dispersal. Our study highlights the value of mummified hair as an alternate source of aDNA from regions with poor bone preservation.

Ancient DNA from Mesopotamia suggests distinct Pre-Pottery and Pottery Neolithic migrations into Anatolia

We present the first ancient DNA data from the Pre-Pottery Neolithic of Mesopotamia (Southeastern Turkey and Northern Iraq), Cyprus, and the Northwestern Zagros, along with the first data from Neolithic Armenia. We show that these and neighboring populations were formed through admixture of pre-Neolithic sources related to Anatolian, Caucasus, and Levantine hunter-gatherers, forming a Neolithic continuum of ancestry mirroring the geography of West Asia. By analyzing Pre-Pottery and Pottery Neolithic populations of Anatolia, we show that the former were derived from admixture between Mesopotamian-related and local Epipaleolithic-related sources, but the latter experienced additional Levantine-related gene flow, thus documenting at least two pulses of migration from the Fertile Crescent heartland to the early farmers of Anatolia.

The genetic history of the Southern Arc: A bridge between West Asia and Europe

By sequencing 727 ancient individuals from the Southern Arc (Anatolia and its neighbors in Southeastern Europe and West Asia) over 10,000 years, we contextualize its Chalcolithic period and Bronze Age (about 5000 to 1000 BCE), when extensive gene flow entangled it with the Eurasian steppe. Two streams of migration transmitted Caucasus and Anatolian/Levantine ancestry northward, and the Yamnaya pastoralists, formed on the steppe, then spread southward into the Balkans and across the Caucasus into Armenia, where they left numerous patrilineal descendants. Anatolia was transformed by intra-West Asian gene flow, with negligible impact of the later Yamnaya migrations. This contrasts with all other regions where Indo-European languages were spoken, suggesting that the homeland of the Indo-Anatolian language family was in West Asia, with only secondary dispersals of non-Anatolian Indo-Europeans from the steppe.

A genetic probe into the ancient and medieval history of Southern Europe and West Asia

Literary and archaeological sources have preserved a rich history of Southern Europe and West Asia since the Bronze Age that can be complemented by genetics. Mycenaean period elites in Greece did not differ from the general population and included both people with some steppe ancestry and others, like the Griffin Warrior, without it. Similarly, people in the central area of the Urartian Kingdom around Lake Van lacked the steppe ancestry characteristic of the kingdom's northern provinces. Anatolia exhibited extraordinary continuity down to the Roman and Byzantine periods, with its people serving as the demographic core of much of the Roman Empire, including the city of Rome itself. During medieval times, migrations associated with Slavic and Turkic speakers profoundly affected the region.

Dairying, diseases and the evolution of lactase persistence in Europe

In European and many African, Middle Eastern and southern Asian populations, lactase persistence (LP) is the most strongly selected monogenic trait to have evolved over the past 10,000 years1. Although the selection of LP and the consumption of prehistoric milk must be linked, considerable uncertainty remains concerning their spatiotemporal configuration and specific interactions2,3. Here we provide detailed distributions of milk exploitation across Europe over the past 9,000 years using around 7,000 pottery fat residues from more than 550 archaeological sites. European milk use was widespread from the Neolithic period onwards but varied spatially and temporally in intensity. Notably, LP selection varying with levels of prehistoric milk exploitation is no better at explaining LP allele frequency trajectories than uniform selection since the Neolithic period. In the UK Biobank4,5 cohort of 500,000 contemporary Europeans, LP genotype was only weakly associated with milk consumption and did not show consistent associations with improved fitness or health indicators. This suggests that other reasons for the beneficial effects of LP should be considered for its rapid frequency increase. We propose that lactase non-persistent individuals consumed milk when it became available but, under conditions of famine and/or increased pathogen exposure, this was disadvantageous, driving LP selection in prehistoric Europe. Comparison of model likelihoods indicates that population fluctuations, settlement density and wild animal exploitation-proxies for these drivers-provide better explanations of LP selection than the extent of milk exploitation. These findings offer new perspectives on prehistoric milk exploitation and LP evolution.

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.

Hunter-gatherer genomes reveal diverse demographic trajectories during the rise of farming in Eastern Africa

The fate of hunting and gathering populations following the rise of agriculture and pastoralism remains a topic of debate in the study of human prehistory. Studies of ancient and modern genomes have found that autochthonous groups were largely replaced by expanding farmer populations with varying levels of gene flow, a characterization that is influenced by the almost universal focus on the European Neolithic.^1-5 We sought to understand the demographic impact of an ongoing cultural transition to farming in Southwest Ethiopia, one of the last regions in Africa to experience such shifts.⁶ Importantly, Southwest Ethiopia is home to several of the world's remaining hunter-gatherer groups, including the Chabu people, who are currently transitioning away from their traditional mode of subsistence.⁷ We generated genome-wide data from the Chabu and four neighboring populations, the Majang, Shekkacho, Bench, and Sheko, to characterize their genetic ancestry and estimate their effective population sizes over the last 60 generations. We show that the Chabu are a distinct population closely related to ancient people who occupied Southwest Ethiopia >4,500 years ago. Furthermore, the Chabu are undergoing a severe population bottleneck, which began approximately 1,400 years ago. By analyzing eleven Eastern African populations, we find evidence for divergent demographic trajectories among hunter-gatherer-descendant groups. Our results illustrate that although foragers respond to encroaching agriculture and pastoralism with multiple strategies, including cultural adoption of agropastoralism, gene flow, and economic specialization, they often face population decline.

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.

Social stratification without genetic differentiation at the site of Kulubnarti in Christian Period Nubia

Relatively little is known about Nubia's genetic landscape prior to the influence of the Islamic migrations that began in the late 1st millennium CE. Here, we increase the number of ancient individuals with genome-level data from the Nile Valley from three to 69, reporting data for 66 individuals from two cemeteries at the Christian Period (~650-1000 CE) site of Kulubnarti, where multiple lines of evidence suggest social stratification. The Kulubnarti Nubians had ~43% Nilotic-related ancestry (individual variation between ~36-54%) with the remaining ancestry consistent with being  introduced through Egypt and ultimately deriving from an ancestry pool like that found in the Bronze and Iron Age Levant. The Kulubnarti gene pool - shaped over a millennium - harbors disproportionately female-associated West Eurasian-related ancestry. Genetic similarity among individuals from the two cemeteries supports a hypothesis of social division without genetic distinction. Seven pairs of inter-cemetery relatives suggest fluidity between cemetery groups. Present-day Nubians are not directly descended from the Kulubnarti Nubians, attesting to additional genetic input since the Christian Period.

Parental relatedness through time revealed by runs of homozygosity in ancient DNA

Parental relatedness of present-day humans varies substantially across the globe, but little is known about the past. Here we analyze ancient DNA, leveraging that parental relatedness leaves genomic traces in the form of runs of homozygosity. We present an approach to identify such runs in low-coverage ancient DNA data aided by haplotype information from a modern phased reference panel. Simulation and experiments show that this method robustly detects runs of homozygosity longer than 4 centimorgan for ancient individuals with at least 0.3 × coverage. Analyzing genomic data from 1,785 ancient humans who lived in the last 45,000 years, we detect low rates of first cousin or closer unions across most ancient populations. Moreover, we find a marked decay in background parental relatedness co-occurring with or shortly after the advent of sedentary agriculture. We observe this signal, likely linked to increasing local population sizes, across several geographic transects worldwide.

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.

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 2^nd 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.

Ancient genomes reveal structural shifts after the arrival of Steppe-related ancestry in the Italian Peninsula

Across Europe, the genetics of the Chalcolithic/Bronze Age transition is increasingly characterized in terms of an influx of Steppe-related ancestry. The effect of this major shift on the genetic structure of populations in the Italian Peninsula remains underexplored. Here, genome-wide shotgun data for 22 individuals from commingled cave and single burials in Northeastern and Central Italy dated between 3200 and 1500 BCE provide the first genomic characterization of Bronze Age individuals (n = 8; 0.001-1.2× coverage) from the central Italian Peninsula, filling a gap in the literature between 1950 and 1500 BCE. Our study confirms a diversity of ancestry components during the Chalcolithic and the arrival of Steppe-related ancestry in the central Italian Peninsula as early as 1600 BCE, with this ancestry component increasing through time. We detect close patrilineal kinship in the burial patterns of Chalcolithic commingled cave burials and a shift away from this in the Bronze Age (2200-900 BCE) along with lowered runs of homozygosity, which may reflect larger changes in population structure. Finally, we find no evidence that the arrival of Steppe-related ancestry in Central Italy directly led to changes in frequency of 115 phenotypes present in the dataset, rather that the post-Roman Imperial period had a stronger influence, particularly on the frequency of variants associated with protection against Hansen's disease (leprosy). Our study provides a closer look at local dynamics of demography and phenotypic shifts as they occurred as part of a broader phenomenon of widespread admixture during the Chalcolithic/Bronze Age transition.

Y-chromosome variability and genetic history of Commons from Northern Italy

Objectives

Genetic drift and admixture are driving forces in human evolution, but their concerted impact to population evolution in historical times and at a micro‐geographic scale is poorly assessed. In this study we test a demographic model encompassing both admixture and drift to the case of social‐cultural isolates such as the so‐called “Commons.”

Materials and methods

Commons are peculiar institutions of medieval origins whose key feature is the tight relationship between population and territory, mediated by the collective property of shared resources. Here, we analyze the Y‐chromosomal genetic structure of four Commons (for a total of 366 samples) from the Central and Eastern Padana plain in Northern Italy.

Results

Our results reveal that all these groups exhibit patterns of significant diversity reduction, peripheral/outlier position within the Italian/European genetic space and high frequency of Common‐specific haplogroups. By explicitly testing different drift‐admixture models, we show that a drift‐only model is more probable for Central Padana Commons, while additional admixture (~20%) from external population around the same time of their foundation cannot be excluded for the Eastern ones.

Discussion

Building on these results, we suggest central Middle Ages as the most probable age of foundation for three of the considered Commons, the remaining one pointing to late antiquity. We conclude that an admixture‐drift model is particularly useful for interpreting the genetic structure and recent demographic history of small‐scale populations in which social‐cultural features play a significant role.

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 genetic history of the pre-contact Caribbean

Humans settled the Caribbean about 6,000 years ago, and ceramic use and intensified agriculture mark a shift from the Archaic to the Ceramic Age at around 2,500 years ago1-3. Here we report genome-wide data from 174 ancient individuals from The Bahamas, Haiti and the Dominican Republic (collectively, Hispaniola), Puerto Rico, Curaçao and Venezuela, which we co-analysed with 89 previously published ancient individuals. Stone-tool-using Caribbean people, who first entered the Caribbean during the Archaic Age, derive from a deeply divergent population that is closest to Central and northern South American individuals; contrary to previous work4, we find no support for ancestry contributed by a population related to North American individuals. Archaic-related lineages were >98% replaced by a genetically homogeneous ceramic-using population related to speakers of languages in the Arawak family from northeast South America; these people moved through the Lesser Antilles and into the Greater Antilles at least 1,700 years ago, introducing ancestry that is still present. Ancient Caribbean people avoided close kin unions despite limited mate pools that reflect small effective population sizes, which we estimate to be a minimum of 500-1,500 and a maximum of 1,530-8,150 individuals on the combined islands of Puerto Rico and Hispaniola in the dozens of generations before the individuals who we analysed lived. Census sizes are unlikely to be more than tenfold larger than effective population sizes, so previous pan-Caribbean estimates of hundreds of thousands of people are too large5,6. Confirming a small and interconnected Ceramic Age population7, we detect 19 pairs of cross-island cousins, close relatives buried around 75 km apart in Hispaniola and low genetic differentiation across islands. Genetic continuity across transitions in pottery styles reveals that cultural changes during the Ceramic Age were not driven by migration of genetically differentiated groups from the mainland, but instead reflected interactions within an interconnected Caribbean world1,8.

Assessing the performance of qpAdm: a statistical tool for studying population admixture

qpAdm is a statistical tool for studying the ancestry of populations with histories that involve admixture between two or more source populations. Using qpAdm, it is possible to identify plausible models of admixture that fit the population history of a group of interest and to calculate the relative proportion of ancestry that can be ascribed to each source population in the model. Although qpAdm is widely used in studies of population history of human (and nonhuman) groups, relatively little has been done to assess its performance. We performed a simulation study to assess the behavior of qpAdm under various scenarios in order to identify areas of potential weakness and establish recommended best practices for use. We find that qpAdm is a robust tool that yields accurate results in many cases, including when data coverage is low, there are high rates of missing data or ancient DNA damage, or when diploid calls cannot be made. However, we caution against co-analyzing ancient and present-day data, the inclusion of an extremely large number of reference populations in a single model, and analyzing population histories involving extended periods of gene flow. We provide a user guide suggesting best practices for the use of qpAdm.

Genomic History of Neolithic to Bronze Age Anatolia, Northern Levant, and Southern Caucasus

Here, we report genome-wide data analyses from 110 ancient Near Eastern individuals spanning the Late Neolithic to Late Bronze Age, a period characterized by intense interregional interactions for the Near East. We find that 6^th millennium BCE populations of North/Central Anatolia and the Southern Caucasus shared mixed ancestry on a genetic cline that formed during the Neolithic between Western Anatolia and regions in today's Southern Caucasus/Zagros. During the Late Chalcolithic and/or the Early Bronze Age, more than half of the Northern Levantine gene pool was replaced, while in the rest of Anatolia and the Southern Caucasus, we document genetic continuity with only transient gene flow. Additionally, we reveal a genetically distinct individual within the Late Bronze Age Northern Levant. Overall, our study uncovers multiple scales of population dynamics through time, from extensive admixture during the Neolithic period to long-distance mobility within the globalized societies of the Late Bronze Age. VIDEO ABSTRACT.

A Dynamic 6,000-Year Genetic History of Eurasia's Eastern Steppe

The Eastern Eurasian Steppe was home to historic empires of nomadic pastoralists, including the Xiongnu and the Mongols. However, little is known about the region’s population history. Here, we reveal its dynamic genetic history by analyzing new genome-wide data for 214 ancient individuals spanning 6,000 years. We identify a pastoralist expansion into Mongolia ca. 3000 BCE, and by the Late Bronze Age, Mongolian populations were biogeographically structured into three distinct groups, all practicing dairy pastoralism regardless of ancestry. The Xiongnu emerged from the mixing of these populations and those from surrounding regions. By comparison, the Mongols exhibit much higher eastern Eurasian ancestry, resembling present-day Mongolic-speaking populations. Our results illuminate the complex interplay between genetic, sociopolitical, and cultural changes on the Eastern Steppe.

•

Genome-wide analysis of 214 ancient individuals from Mongolia and the Baikal region

•

Three genetically distinct dairy pastoralist groups in Late Bronze Age Mongolia

•

Xiongnu nomadic empire formed through mixing of distinct local and distant groups

•

No selection on the lactase persistence alleles despite 5,000 years of dairy culture

A systematic investigation of human DNA preservation in medieval skeletons

Ancient DNA (aDNA) analyses necessitate the destructive sampling of archaeological material. Currently, the cochlea, part of the osseous inner ear located inside the petrous pyramid, is the most sought after skeletal element for molecular analyses of ancient humans as it has been shown to yield high amounts of endogenous DNA. However, destructive sampling of the petrous pyramid may not always be possible, particularly in cases where preservation of skeletal morphology is of top priority. To investigate alternatives, we present a survey of human aDNA preservation for each of ten skeletal elements in a skeletal collection from Medieval Germany. Through comparison of human DNA content and quality we confirm best performance of the petrous pyramid and identify seven additional sampling locations across four skeletal elements that yield adequate aDNA for most applications in human palaeogenetics. Our study provides a better perspective on DNA preservation across the human skeleton and takes a further step toward the more responsible use of ancient materials in human aDNA studies.

A Genetic History of the Near East from an aDNA Time Course Sampling Eight Points in the Past 4,000 Years

The Iron and Classical Ages in the Near East were marked by population expansions carrying cultural transformations that shaped human history, but the genetic impact of these events on the people who lived through them is little-known. Here, we sequenced the whole genomes of 19 individuals who each lived during one of four time periods between 800 BCE and 200 CE in Beirut on the Eastern Mediterranean coast at the center of the ancient world’s great civilizations. We combined these data with published data to traverse eight archaeological periods and observed any genetic changes as they arose. During the Iron Age (∼1000 BCE), people with Anatolian and South-East European ancestry admixed with people in the Near East. The region was then conquered by the Persians (539 BCE), who facilitated movement exemplified in Beirut by an ancient family with Egyptian-Lebanese admixed members. But the genetic impact at a population level does not appear until the time of Alexander the Great (beginning 330 BCE), when a fusion of Asian and Near Easterner ancestry can be seen, paralleling the cultural fusion that appears in the archaeological records from this period. The Romans then conquered the region (31 BCE) but had little genetic impact over their 600 years of rule. Finally, during the Ottoman rule (beginning 1516 CE), Caucasus-related ancestry penetrated the Near East. Thus, in the past 4,000 years, three limited admixture events detectably impacted the population, complementing the historical records of this culturally complex region dominated by the elite with genetic insights from the general population.

Ancient genomes reveal complex patterns of population movement, interaction, and replacement in sub-Saharan Africa

Africa hosts the greatest human genetic diversity globally, but legacies of ancient population interactions and dispersals across the continent remain understudied. Here, we report genome-wide data from 20 ancient sub-Saharan African individuals, including the first reported ancient DNA from the DRC, Uganda, and Botswana. These data demonstrate the contraction of diverse, once contiguous hunter-gatherer populations, and suggest the resistance to interaction with incoming pastoralists of delayed-return foragers in aquatic environments. We refine models for the spread of food producers into eastern and southern Africa, demonstrating more complex trajectories of admixture than previously suggested. In Botswana, we show that Bantu ancestry post-dates admixture between pastoralists and foragers, suggesting an earlier spread of pastoralism than farming to southern Africa. Our findings demonstrate how processes of migration and admixture have markedly reshaped the genetic map of sub-Saharan Africa in the past few millennia and highlight the utility of combined archaeological and archaeogenetic approaches.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

West Asian sources of the Eurasian component in Ethiopians: a reassessment

The presence of genomic signatures of Eurasian origin in contemporary Ethiopians has been reported by several authors and estimated to have arrived in the area from 3000 years ago. Several studies reported plausible source populations for such a signature, using haplotype based methods on modern data or single-site methods on modern or ancient data. These studies did not reach a consensus and suggested an Anatolian or Sardinia-like proxy, broadly Levantine or Neolithic Levantine as possible sources. We demonstrate, however, that the deeply divergent, autochthonous African component which accounts for ~50% of most contemporary Ethiopian genomes, affects the overall allele frequency spectrum to an extent that makes it hard to control for it and, at once, to discern between subtly different, yet important, Eurasian sources (such as Anatolian or Levant Neolithic ones). Here we re-assess pattern of allele sharing between the Eurasian component of Ethiopians (here called “NAF” for Non African) and ancient and modern proxies. Our results unveil a genomic legacy that may connect the Eurasian genetic component of contemporary Ethiopians with Sea People and with population movements that affected the Mediterranean area and the Levant after the fall of the Minoan civilization.

Ancient DNA from the skeletons of Roopkund Lake reveals Mediterranean migrants in India

Situated at over 5,000 meters above sea level in the Himalayan Mountains, Roopkund Lake is home to the scattered skeletal remains of several hundred individuals of unknown origin. We report genome-wide ancient DNA for 38 skeletons from Roopkund Lake, and find that they cluster into three distinct groups. A group of 23 individuals have ancestry that falls within the range of variation of present-day South Asians. A further 14 have ancestry typical of the eastern Mediterranean. We also identify one individual with Southeast Asian-related ancestry. Radiocarbon dating indicates that these remains were not deposited simultaneously. Instead, all of the individuals with South Asian-related ancestry date to ~800 CE (but with evidence of being deposited in more than one event), while all other individuals date to ~1800 CE. These differences are also reflected in stable isotope measurements, which reveal a distinct dietary profile for the two main groups.

Remains of several hundred humans are scattered around Roopkund Lake, situated over 5,000 meters above sea level in the Himalayan Mountains. Here the authors analyze genome-wide data from 38 skeletons and find 3 clusters with different ancestries and dates, showing the people were desposited in multiple catastrophic events.

Ancient DNA sheds light on the genetic origins of early Iron Age Philistines

The ancient Mediterranean port city of Ashkelon, identified as "Philistine" during the Iron Age, underwent a marked cultural change between the Late Bronze and the early Iron Age. It has been long debated whether this change was driven by a substantial movement of people, possibly linked to a larger migration of the so-called "Sea Peoples." Here, we report genome-wide data of 10 Bronze and Iron Age individuals from Ashkelon. We find that the early Iron Age population was genetically distinct due to a European-related admixture. This genetic signal is no longer detectible in the later Iron Age population. Our results support that a migration event occurred during the Bronze to Iron Age transition in Ashkelon but did not leave a long-lasting genetic signature.

Changes in human mandibular shape during the Terminal Pleistocene-Holocene Levant

The transition to food production, exploitation of ‘secondary’ products (e.g., milk), and advances in cookware technology have affected all aspects of human life. The aim of the present study was to follow changes in mandibular form and shape throughout the terminal Pleistocene-Holocene Levant. The hemimandibles of four populations were included in this study: Natufian hunter-gatherers (n = 10), Pre-pottery Neolithic early farmers (n = 6), Chalcolithic farmers (n = 9), Roman-Byzantine (n = 16), and modern (n = 63) populations. A surface mesh of each mandible was reconstructed from CT or surface scans. Changes in mandibular form and shape were studied using the Procrustes-based geometric morphometrics method. Univariate and multivariate analyses were carried out to examine differences in size and shape between the studied populations. Our results reveal considerable temporal changes in mandibular shape throughout the Holocene Levant, mainly between the pre-agricultural population (the Natufian) and the succeeding ones, and between the post-industrial (the Modern) and the pre-industrial populations. A tendency for a reduction in mandibular size was identified between the pre-agricultural population and the farmers. Most regions of the mandible underwent shape changes. In conclusion, substantial changes in mandibular shape occurred throughout the Holocene Levant, especially following the agricultural revolution. These changes can be explained by the “masticatory-functional hypothesis”.

Ancient DNA reveals a multistep spread of the first herders into sub-Saharan Africa

How food production first entered eastern Africa ~5000 years ago and the extent to which people moved with livestock is unclear. We present genome-wide data from 41 individuals associated with Later Stone Age, Pastoral Neolithic (PN), and Iron Age contexts in what are now Kenya and Tanzania to examine the genetic impacts of the spreads of herding and farming. Our results support a multiphase model in which admixture between northeastern African-related peoples and eastern African foragers formed multiple pastoralist groups, including a genetically homogeneous PN cluster. Additional admixture with northeastern and western African-related groups occurred by the Iron Age. These findings support several movements of food producers while rejecting models of minimal admixture with foragers and of genetic differentiation between makers of distinct PN artifacts.

Ancient human mitochondrial genomes from Bronze Age Bulgaria: new insights into the genetic history of Thracians

One of the best documented Indo-European civilizations that inhabited Bulgaria is the Thracians, who lasted for more than five millennia and whose origin and relationships with other past and present-day populations are debated among researchers. Here we report 25 new complete mitochondrial genomes of ancient individuals coming from three necropolises located in different regions of Bulgaria – Shekerdja mogila, Gabrova mogila and Bereketska mogila – dated to II-III millennium BC. The identified mtDNA haplogroup composition reflects the mitochondrial variability of Western Eurasia. In particular, within the ancient Eurasian genetic landscape, Thracians locate in an intermediate position between Early Neolithic farmers and Late Neolithic-Bronze Age steppe pastoralists, supporting the scenario that the Balkan region has been a link between Eastern Europe and the Mediterranean since the prehistoric time. Spatial Principal Component Analysis (sPCA) performed on Thracian and modern mtDNA sequences, confirms the pattern highlighted on ancient populations, overall indicating that the maternal gene pool of Thracians reflects their central geographical position at the gateway of Europe.