The Great Hungarian Plain was a crossroads of cultural transformations that have shaped European prehistory. Here we analyse a 5,000-year transect of human genomes, sampled from petrous bones giving consistently excellent endogenous DNA yields, from 13 Hungarian Neolithic, Copper, Bronze and Iron Age burials including two to high (~22 × ) and seven to ~1 × coverage, to investigate the impact of these on Europe’s genetic landscape. These data suggest genomic shifts with the advent of the Neolithic, Bronze and Iron Ages, with interleaved periods of genome stability. The earliest Neolithic context genome shows a European hunter-gatherer genetic signature and a restricted ancestral population size, suggesting direct contact between cultures after the arrival of the first farmers into Europe. The latest, Iron Age, sample reveals an eastern genomic influence concordant with introduced Steppe burial rites. We observe transition towards lighter pigmentation and surprisingly, no Neolithic presence of lactase persistence.

Ancient DNA studies have established that Neolithic European populations were descended from Anatolian migrants1–8 who received a limited amount of admixture from resident hunter-gatherers3–5,9. Many open questions remain, however, about the spatial and temporal dynamics of population interactions and admixture during the Neolithic period. Using the highest-resolution genome-wide ancient DNA data set assembled to date—a total of 180 samples, 130 newly reported here, from the Neolithic and Chalcolithic of Hungary (6000–2900 BCE, n = 100), Germany (5500–3000 BCE, n = 42), and Spain (5500–2200 BCE, n = 38)—we investigate the population dynamics of Neolithization across Europe. We find that genetic diversity was shaped predominantly by local processes, with varied sources and proportions of hunter-gatherer ancestry among the three regions and through time. Admixture between groups with different ancestry profiles was pervasive and resulted in observable population transformation across almost all cultural transitions. Our results shed new light on the ways that gene flow reshaped European populations throughout the Neolithic period and demonstrate the potential of time-series-based sampling and modeling approaches to elucidate multiple dimensions of historical population interactions.

The invention and development of next or second generation sequencing methods has resulted in a dramatic transformation of ancient DNA research and allowed shotgun sequencing of entire genomes from fossil specimens. However, although there are exceptions, most fossil specimens contain only low (~ 1% or less) percentages of endogenous DNA. The only skeletal element for which a systematically higher endogenous DNA content compared to other skeletal elements has been shown is the petrous part of the temporal bone. In this study we investigate whether (a) different parts of the petrous bone of archaeological human specimens give different percentages of endogenous DNA yields, (b) there are significant differences in average DNA read lengths, damage patterns and total DNA concentration, and (c) it is possible to obtain endogenous ancient DNA from petrous bones from hot environments. We carried out intra-petrous comparisons for ten petrous bones from specimens from Holocene archaeological contexts across Eurasia dated between 10,000-1,800 calibrated years before present (cal. BP). We obtained shotgun DNA sequences from three distinct areas within the petrous: a spongy part of trabecular bone (part A), the dense part of cortical bone encircling the osseous inner ear, or otic capsule (part B), and the dense part within the otic capsule (part C). Our results confirm that dense bone parts of the petrous bone can provide high endogenous aDNA yields and indicate that endogenous DNA fractions for part C can exceed those obtained for part B by up to 65-fold and those from part A by up to 177-fold, while total endogenous DNA concentrations are up to 126-fold and 109-fold higher for these comparisons. Our results also show that while endogenous yields from part C were lower than 1% for samples from hot (both arid and humid) parts, the DNA damage patterns indicate that at least some of the reads originate from ancient DNA molecules, potentially enabling ancient DNA analyses of samples from hot regions that are otherwise not amenable to ancient DNA analyses.

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 last 10,000 years 1 . While LP selection and prehistoric milk consumption must be linked, considerable uncertainty remains concerning their spatiotemporal configuration and specific interactions 2,3 . We provide detailed distributions of milk exploitation across Europe over the last 9k years using c. 7,000 pottery fat residues from >550 archaeological sites. European milk use was widespread from the Neolithic period onwards but varied spatially and temporally in intensity. Surprisingly, comparison of model likelihoods indicates that LP selection varying with levels of prehistoric milk exploitation provides no better explanation of LP allele frequency trajectories than uniform selection since the Neolithic. In the UK Biobank 4,5 cohort of ~500K 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 other hypotheses on 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 that under particular conditions and microbiological milieux this was disadvantageous, driving LP selection in prehistoric Europe. Comparison of model likelihoods indicates that population fluctuations, settlement density and wild animal exploitationproxies for these driversprovide better explanations of LP selection than the extent of milk exploitation. These findings offer new perspectives on prehistoric milk exploitation and LP evolution.

The chief objective of this paper is to improve our understanding of the Neolithic in Eastern Hungary using absolute chronological data. To accomplish this we calibrated new measurements as well as previously published dates. The up-to-date, standardized evaluation of 261 calibrated measurements showed temporal overlaps between archaeological cultures defined on the basis of ceramic styles. The increasing number of dates suggest that the Neolithic period began at the turn of the 6th and 7th millennia bc and lasted for ca. 1500 yr in the present area of the Great Hungarian Plain (Alföld). Further research should be aimed at complementing the current data set with dates from western Hungary and establishing additional correlations among stratigraphic, typological and radiocarbon dates.