Search

Includes supplementary materials for the online appendix. 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. We thank all members of the Archaeogenetics Department of the Max Planck Institute for Evolutionary Anthropology, especially the Population Genetics group, and we thank Paleogenomics Lab members from Instituto de Biología Evolutiva for the data discussion. We especially thank M. Hajdinjak and I. Olalde for valuable comments and suggestions. We thank G. Brandt, A. Wissgott, S. Clayton and K. Prüfer for sequencing and handling the raw data. We thank all the members of the general research project in Cueva de Ardales and Sima de las Palomas de Teba. We thank all colleagues, archaeologists, geologists and experts in different scientific disciplines, as well as students from the Universities of Cádiz, Cologne and Tübingen, who have participated in the field and laboratory works. We are grateful for the great collaboration of the Town Councils of Ardales and Teba, within the framework of collaboration agreements with the University of Cadiz and the Neanderthal Museum. This work was supported by the Max Planck Society (V.V.-M., A.B.R., J.K. and W.H.), Margarita Salas 2022 funded by Unión Europea-Next Generation EU (V.V.-M.) and the European Research Council under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 771234-PALEoRIDER to W.H. and no. 803147-RESOLUTION to S.T.). The information provided on recent dating and archaeological records from Cueva de Ardales is part of the General Research Project, authorized by the Consejería de Cultura y Patrimonio Histórico de la Junta de Andalucía: entitled ‘Prehistoric societies (from the Middle Palaeolithic to the Final Neolithic) in the Cueva de Ardales and Sima de las Palomas de Teba (Málaga, Spain)’. A geoarchaeological, chronological and environmental study was conducted under the direction of J.R.-M. and G.-C.W. between 2015 and 2021. Excavations and field and laboratory studies have been carried out in Cueva de Ardales within the framework of the projects entitled ‘Analysis of prehistoric societies from the Middle Palaeolithic to the Late Neolithic on the two shores of the Strait of Gibraltar’. Relaciones y contactos—HAR2017-8734P (Ministerio de Economía, Industria y Competitividad-Agencia Estatal de Investigación, cofinanced by FEDER funds) for which J.R.-M. and Salvador Domínguez-Bella are the responsible researchers. Collaborative Research Centre—CRC 806 Our Way to Europe, funded by the Deutsche Forschungsgemeinschaft (DFG-German Research Foundation). Project no. 57444011, under the responsibility of G.-C.W.

Robust and reliable estimates of how individuals are biologically related to each other are a key source of information when reconstructing pedigrees. In combination with contextual data, reconstructed pedigrees can be used to infer possible kinship practices in prehistoric populations. However, standard methods to estimate biological relatedness from genome sequence data cannot be applied to low coverage sequence data, such as are common in ancient DNA (aDNA) studies. Critically, a statistically robust method for assessing and quantifying the confidence of a classification of a specific degree of relatedness for a pair of individuals, using low coverage genome data, is lacking.In this paper we present the R-package BREADR (Biological RElatedness from Ancient DNA in R), which leverages the so-called pairwise mismatch rate, calculated on optimally-thinned genome-wide pseudo-haploid sequence data, to estimate genetic relatedness up to the second degree, assuming an underlying binomial distribution. BREADR also returns a posterior probability for each degree of relatedness, from identical twins/same individual, first-degree, second-degree or “unrelated” pairs, allowing researchers to quantify and report the uncertainty, even for particularly low-coverage data. We show that this method accurately recovers degrees of relatedness for sequence data with coverage as low as 0.04Xusing simulated data, and then compare the performance of BREADR on empirical data from Bronze Age Iberian human sequence data. The BREADR package is designed for pseudo-haploid genotype data, common in aDNA studies.

Acknowledgements: The authors thank G. Marciani and O. Jöris for comments on archaeology; C. Jeong, M. Spyrou and K. Prüfer for comments on genetics; M. O’Reilly for graphical support for Fig. 5 and Extended Data Fig. 9; the entire IT and laboratory teams at the Department of Archaeogenetics of MPI-SHH for technical assistance; M. Meyer and S. Nagel for support with single-stranded library preparation; K. Post, P. van Es, J. Glimmerveen, M. Medendorp, M. Sier, S. Dikstra, M. Dikstra, R. van Eerden, D. Duineveld and A. Hoekman for providing access to human specimens from the North Sea (The Netherlands); M. D. Garralda and A. Estalrrich for providing access to human specimens from La Riera (Spain); J. Górski and M. Zając for providing access to human specimens from Maszycka cave; C. Di Patti for providing access to human specimens from San Teodoro 2 (Italy); P. Blaževičius for providing access to the Donkalnis human remains and the new radiocarbon dates; the Italian Ministry of Culture and Soprintendenza Archeologia Belle Arti e Paesaggio for the Provinces of Verona, Rovigo, and Vicenza for granting access to the human remains of Tagliente 2; F. Fontana, who carries out investigations of the Riparo Tagliente site (Italy); the Friuli Venezia Giulia Superintendency for providing access to the human tooth Pradis 1; and the Soprintendenza Archeologia Belle Arti e Paesaggio for the Provinces of Barletta-Andria-Trani and Foggia for providing access to the Paglicci human remains. This project has received funding by the European Research Council under the European Union’s Horizon 2020 research and innovation programme under grant agreements no. 803147-RESOLUTION (to S.T.), no. 771234-PALEoRIDER (to W.H.), no. 864358 (to K.M.), no. 724703 and no. 101019659 (to K.H.). K.H. is also supported by the Deutsche Forschungsgemeinschaft (DFG FOR 2237). E.A. has received funding from the Van de Kamp fonds. PACEA co-authors of this research benefited from the scientific framework of the University of Bordeaux’s IdEx Investments for the Future programme/GPR Human Past. A.G.-O. is supported by a Ramón y Cajal fellowship (RYC-2017-22558). L. Sineo, M.L. and D.C. have received funding from the Italian Ministry of University and Research (MUR) PRIN 2017 grants 20177PJ9XF and 20174BTC4R_002. H. Rougier received support from the College of Social and Behavioral Sciences of CSUN and the CSUN Competition for RSCA Awards. C.L.S. and T. Saupe received support from the European Union through the European Regional Development Fund (project no. 2014-2020.4.01.16-0030) and C.L.S. received support from the Estonian Research Council grant PUT (PRG243). S. Shnaider received support from the Russian Science Foundation (no. 19-78-10053)., 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.

The rules and structure of human culture impact health and disease as much as genetics or the natural environment. To study the origin and evolution of these patterns, we take a multidisciplinary approach combining ancient DNA, skeletal metrics, paleopathology, and stable isotopes. Our analysis focuses on cultural, environmental, and genetic contributions to variation in stature in four populations of Early Neolithic Europe. In Central Europe, low female stature is likely due to male preference in resource allocation under conditions of stress. In contrast, shorter male stature in Mediterranean populations may reflect a lack of preference. Our analysis suggests that biological consequences of sex-specific inequities can be linked to culture as early as 7000 years before present. Understanding these patterns is key to interpreting the evolution of genetic and socio-cultural determinants of health, and our results show that culture, more than environment or genetics, drove height disparities in Early Neolithic Europe.

The history of the British Isles and Ireland is characterized by multiple periods of major cultural change, including the influential transformation after the end of Roman rule, which precipitated shifts in language, settlement patterns and material culture1. The extent to which migration from continental Europe mediated these transitions is a matter of long-standing debate2–4. Here we study genome-wide ancient DNA from 460 medieval northwestern Europeans—including 278 individuals from England—alongside archaeological data, to infer contemporary population dynamics. We identify a substantial increase of continental northern European ancestry in early medieval England, which is closely related to the early medieval and present-day inhabitants of Germany and Denmark, implying large-scale substantial migration across the North Sea into Britain during the Early Middle Ages. As a result, the individuals who we analysed from eastern England derived up to 76% of their ancestry from the continental North Sea zone, albeit with substantial regional variation and heterogeneity within sites. We show that women with immigrant ancestry were more often furnished with grave goods than women with local ancestry, whereas men with weapons were as likely not to be of immigrant ancestry. A comparison with present-day Britain indicates that subsequent demographic events reduced the fraction of continental northern European ancestry while introducing further ancestry components into the English gene pool, including substantial southwestern European ancestry most closely related to that seen in Iron Age France

The recent discovery of an exceptionally rich grave at La Almoloya in south-eastern Spain illuminates the political context of Early Bronze Age El Argar society. The quantity, variety and opulence of the grave goods emphasise the technological, economic and social dimensions of this unique culture. The assemblage includes politically and ideologically emblematic objects, among which a silver diadem stands out. Of equally exceptional character is the building under which the grave was found—possibly one of the first Bronze Age palaces identified in Western Europe. The architecture and artefacts from La Almoloya provide new insight into emblematic individuals and the exercise of power in societies of marked economic asymmetry. Introduction Argaric diadems as embodied emblems La Almoloya grave 38 Spaces of power Conclusions

The Middle Neolithic in western Europe is characterized by monumental funerary structures, known as megaliths, along the Atlantic façade. The first manifestations of this phenomenon occurred in modern-day France with the long mounds of the Cerny culture. Here, we present genome-wide data from the fifth-millennium BCE site of Fleury-sur-Orne in Normandy (France), famous for its impressively long monuments built for selected individuals. The site encompasses 32 monuments of variable sizes, containing the burials of 19 individuals from the Neolithic period. To address who was buried at the site, we generated genome-wide data for 14 individuals, of whom 13 are males, completing previously published data [M. Rivollat et al., Sci. Adv. 6, eaaz5344 (2020)]. Population genetic and Y chromosome analyses show that the Fleury-sur-Orne group fits within western European Neolithic genetic diversity and that the arrival of a new group is detected after 4,000 calibrated BCE. The results of analyzing uniparentally inherited markers and an overall low number of long runs of homozygosity suggest a patrilineal group practicing female exogamy. We find two pairs of individuals to be father and son, buried together in the same monument/grave. No other biological relationship can link monuments together, suggesting that each monument was dedicated to a genetically independent lineage. The combined data and documented father–son line of descent suggest a male-mediated transmission of sociopolitical authority. However, a single female buried with an arrowhead, otherwise considered a symbol of power of the male elite of the Cerny culture, questions a strictly biological sex bias in the burial rites of this otherwise “masculine” monumental cemetery. Data Overview Genetic Identity, Demographic History Genetic Relatedness, Spatial Organization, and Social Inferences Concluding Remarks Materials and Methods

The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague’s formative years in terms of its early evolution and ecology.

It has been hypothesized that the Neolithic transition towards an agricultural and pastoralist economy facilitated the emergence of human adapted pathogens. Here, we recovered eight Salmonella enterica subsp. enterica genomes from human skeletons of transitional foragers, pastoralists, and agro-pastoralists in western Eurasia that were up to 6,500 years old. Despite the high genetic diversity of S. enterica all ancient bacterial genomes clustered in a single previously uncharacterized branch that contains S. enterica adapted to multiple mammalian species. All ancient bacterial genomes from prehistoric (agro-)pastoralists fall within a part of this branch that also includes the human-specific S. enterica Paratyphi C, illustrating the evolution of a human pathogen over a period of five thousand years. Bacterial genomic comparisons suggest that the earlier ancient strains were not host specific, differed in pathogenic potential, and experienced convergent pseudogenization that accompanied their downstream host adaptation. These observations support the concept that the emergence of human adapted S. enterica is linked to human cultural transformations.

L’exploration des structures de parenté dans les sociétés du passé a été au centre des études en archéologie et en anthropologie. Cependant, la reconstruction des liens de parenté biologique en contextes archéologiques a rarement été accessible. Avec le développement des méthodes d’analyses en ADN ancien, il est désormais possible d’obtenir des données génomiques pour de nombreux individus d’un même groupe, même dans le cas d’une faible conservation de l’ADN. Nous présentons ici de nouvelles ...

Human expansion in the course of the Neolithic transition in western Eurasia has been one of the major topics in ancient DNA research in the last 10 years. Multiple studies have shown that the spread of agriculture and animal husbandry from the Near East across Europe was accompanied by large-scale human expansions. Moreover, changes in subsistence and migration associated with the Neolithic transition have been hypothesized to involve genetic adaptation. Here, we present high quality genome-wide data from the Linear Pottery Culture site Derenburg-Meerenstieg II (DER) (N = 32 individuals) in Central Germany. Population genetic analyses show that the DER individuals carried predominantly Anatolian Neolithic-like ancestry and a very limited degree of local hunter-gatherer admixture, similar to other early European farmers. Increasing the Linear Pottery culture cohort size to ∼100 individuals allowed us to perform various frequency- and haplotype-based analyses to investigate signatures of selection associated with changes following the adoption of the Neolithic lifestyle. In addition, we developed a new method called Admixture-informed Maximum-likelihood Estimation for Selection Scans that allowed us test for selection signatures in an admixture-aware fashion. Focusing on the intersection of results from these selection scans, we identified various loci associated with immune function (JAK1, HLA-DQB1) and metabolism (LMF1, LEPR, SORBS1), as well as skin color (SLC24A5, CD82) and folate synthesis (MTHFR, NBPF3). Our findings shed light on the evolutionary pressures, such as infectious disease and changing diet, that were faced by the early farmers of Western Eurasia. Introduction New Approaches Results and discussion - Data generation and authentication. - Genome-wide data analysis. - Uniparentally inherited markers. - Biological Relatedness. - LBK inter-site comparisons. - Phenotypic analysis. - HLA class-I and II allele analysis. - Selection Scans Conclusions Material and methods - Archaeological background/site description. - Sample preparation. - Genotyping. - F-statistics - Biological Relatedness and ROH Analyses. - Mitochondrial and Y haplotype analysis. - Shotgun data. - HLA typing. - Selection scans. - AIMLESS Admixture-Informed Maximum-Likelihood Estimation for Selection Scans. - Haplotype-based scan.

[Significance] The bacterium Yersinia pestis has caused numerous historically documented outbreaks of plague and research using ancient DNA could demonstrate that it already affected human populations during the Neolithic. However, the pathogen’s genetic diversity, geographic spread, and transmission dynamics during this early period of Y. pestis evolution are largely unexplored. Here, we describe a set of ancient plague genomes up to 5,000 y old from across Eurasia. Our data demonstrate that two genetically distinct forms of Y. pestis evolved in parallel and were both distributed across vast geographic distances, potentially occupying different ecological niches. Interpreted within the archeological context, our results suggest that the spread of plague during this period was linked to increased human mobility and intensification of animal husbandry., The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague’s formative years in terms of its early evolution and ecology., This study was funded by the Max Planck Society, Max Planck Harvard Research Center for the Archaeoscience of the Ancient Mediterranean and the European Research Council under the European Union’s Horizon 2020 research and innovation program under Grant Agreement 771234 – PALEoRIDER (to W.H.), 856453 – HistoGenes (to J.K.), and 834616 – ARCHCAUCASUS (to S.H.). The Heidelberg Academy of Science financed the genetic and archeological research on human individuals from the Augsburg region within the project WIN Kolleg: “Times of Upheaval: Changes of Society and Landscape at the Beginning of the Bronze Age. M.E. was supported by the award “Praemium Academiae” of the Czech Academy of Sciences. M.D. was supported by the project RVO 67985912 of the Institute of Archaeology of the Czech Academy of Sciences, Prague. I.O. was supported by the Ramón y Cajal grant from Ministerio de Ciencia e Innovación, Spanish Government (RYC2019-027909-I). A. H€ubner was supported by the Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy (EXC 2051 – Project-ID 390713860). J.F.-E. and J.A.M.-A. were supported by the Diputación Foral de Alava, IT 1223-19, Gobierno Vasco. A. Buzhilova was supported by the Center of Information Technologies and Systems (CITIS), Moscow, Russia 121041500329-0. L. M., L.B.D., and E. Khussainova were supported by the Grant AP08856654, Ministry of Education and Science of the Republic of Kazakhstan. A. Beisenov was supported by the Grant AP08857177, Ministry of Education and Science of the Republic of Kazakhstan.

The methods presented here seek to maximize the chances for the recovery of human DNA from ancient archaeological remains while limiting input sample material. This was done by targeting anatomical sampling locations previously determined to yield the highest amounts of ancient DNA (aDNA) in a comparative analysis of DNA recovery across the skeleton. Prior research has suggested that these protocols maximize the chances for the successful recovery of ancient human and pathogen DNA from archaeological remains. DNA yields were previously assessed by Parker et al. 2020 in a broad survey of aDNA preservation across multiple skeletal elements from 11 individuals recovered from the medieval (radiocarbon dated to a period of circa (ca.) 1040-1400 CE, calibrated 2-sigma range) graveyard at Krakauer Berg, an abandoned medieval settlement near Peißen Germany. These eight sampling spots, which span five skeletal elements (pars petrosa, permanent molars, thoracic vertebra, distal phalanx, and talus) successfully yielded high-quality ancient human DNA, where yields were significantly greater than the overall average across all elements and individuals. Yields were adequate for use in most common downstream population genetic analyses. Our results support the preferential use of these anatomical sampling locations for most studies involving the analyses of ancient human DNA from archaeological remains. Implementation of these methods will help to minimize the destruction of precious archaeological specimens. Introduction Protocol 1. Considerations before sample processing 2. Pretreatment 3. Bone powder generation Representative Results Discussion

Description, Steppe ancestry among the non–Indo-European–speaking Etruscans challenges previous hypotheses on their recent Anatolian origin., The origin, development, and legacy of the enigmatic Etruscan civilization from the central region of the Italian peninsula known as Etruria have been debated for centuries. Here we report a genomic time transect of 82 individuals spanning almost two millennia (800 BCE to 1000 CE) across Etruria and southern Italy. During the Iron Age, we detect a component of Indo-European–associated steppe ancestry and the lack of recent Anatolian-related admixture among the putative non–Indo-European–speaking Etruscans. Despite comprising diverse individuals of central European, northern African, and Near Eastern ancestry, the local gene pool is largely maintained across the first millennium BCE. This drastically changes during the Roman Imperial period where we report an abrupt population-wide shift to ~50% admixture with eastern Mediterranean ancestry. Last, we identify northern European components appearing in central Italy during the Early Middle Ages, which thus formed the genetic landscape of present-day Italian populations.

Uniparentally-inherited markers on mitochondrial DNA (mtDNA) and the non-recombining regions of the Y chromosome (NRY), have been used for the past 30 years to investigate the history of humans from a maternal and paternal perspective. Researchers have preferred mtDNA due to its abundance in the cells, and comparatively high substitution rate. Conversely, the NRY is less susceptible to back mutations and saturation, and is potentially more informative than mtDNA owing to its longer sequence length. However, due to comparatively poor NRY coverage via shotgun sequencing, and the relatively low and biased representation of Y-chromosome variants on capture assays such as the 1240 k, ancient DNA studies often fail to utilize the unique perspective that the NRY can yield. Here we introduce a new DNA enrichment assay, coined YMCA (Y-mappable capture assay), that targets the "mappable" regions of the NRY. We show that compared to low-coverage shotgun sequencing and 1240 k capture, YMCA significantly improves the mean coverage and number of sites covered on the NRY, increasing the number of Y-haplogroup informative SNPs, and allowing for the identification of previously undiscovered variants. To illustrate the power of YMCA, we show that the analysis of ancient Y-chromosome lineages can help to resolve Y-chromosomal haplogroups. As a case study, we focus on H2, a haplogroup associated with a critical event in European human history: the Neolithic transition. By disentangling the evolutionary history of this haplogroup, we further elucidate the two separate paths by which early farmers expanded from Anatolia and the Near East to western Europe., Open Access funding enabled and organized by Projekt DEAL. This study was funded by the Max Planck Society, the French (ANR) and German (DFG) Research Foundations under the INTERACT project (ANR-17-FRAL-0010, DFG-HA-5407/4-1, 2018-2021) to M.R. and W.H., the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under Grant agreement no. 771234-PALEoRIDER to W.H., the award Praemium Academiae of the Czech Academy of Sciences to M.E. and the project RVO 67985912 of the Institute of Archaeology of the Czech Academy of Sciences, Prague to M.S.