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

Helen Fewlass, Elena I Zavala, Yoann Fagault, Thibaut Tuna, Edouard Bard, Jean-Jacques Hublin, Mateja Hajdinjak, Jaroslaw Wilczyński

► To cite this version:

Helen Fewlass, Elena I Zavala, Yoann Fagault, Thibaut Tuna, Edouard Bard, et al.. Chronological and genetic analysis of an Upper Palaeolithic female infant burial from Borsuka Cave, Poland. iScience, 2023, 26 (12), pp.108283. 10.1016/j.isci.2023.108283. hal-05425211

HAL Id: hal-05425211

https://college-de-france.hal.science/hal-05425211v1

Submitted on 19 Dec 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

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

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

Borsuka cave

<figure><img src="image_4.png" /></figure>

Helen Fewlass,
Elena I. Zavala,
Yoann Fagault, ...,
Jean-Jacques
Hublin, Mateja
Hajdinjak,
Jarosław
Wilczyński

helen.fewlass@crick.ac.uk (H.F.)
ezavala9@berkeley.edu (E.I.Z.)

Highlights

Interdisciplinary analyses
associate Borsuka Cave
burial with the Upper
Palaeolithic

Human remains and
herbivore tooth pendants
are directly dated to
>30,000 years old

Genetic analyses
demonstrate this is the
earliest known female
infant burial

Minimally invasive
methods used <67 mg of
tooth dentine for multiple
analyses

Fewlass et al., iScience 26,
108283
December 15, 2023 © 2023 The
Author(s).
https://doi.org/10.1016/
j.isci.2023.108283
## Article

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

Helen Fewlass,<sup>1,2,9,10,*</sup> Elena I. Zavala,<sup>3,4,5,9,*</sup> Yoann Fagault,<sup>6</sup> Thibaut Tuna,<sup>6</sup> Edouard Bard,<sup>6</sup> Jean-Jacques Hublin,<sup>1,7</sup> Mateja Hajdinjak,<sup>3</sup> and Jarosław Wilczyński<sup>8</sup>

## SUMMARY

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

## INTRODUCTION

Human remains from the Upper Palaeolithic (UP) in Eurasia are sparse, but the period is marked by the increasing visibility of human burials in the archaeological record, typically associated with the Gravettian culture (~34 - 24 ka cal BP) and occasionally even earlier (~39 - 34 ka cal BP). Eastern Central Europe has a rich record of Palaeolithic archaeology with several sites containing human remains (Figure 1), which are associated with different technological assemblages. Directly dated human remains at Mladec Cave (Czech Republic),<sup>1</sup> Oblazowa Cave (Poland),<sup>2–6</sup> Peştera Muierii<sup>7</sup> and Peştera Cioclovina<sup>8</sup> in Romania, and Troisième Caverne of Goyet (individuals Q116-1 and Q376-3) in Belgium<sup>9,10</sup> all fall within the ~37 - 33 ka cal BP range and are tentatively associated with Aurignacian artifacts. The Kostenki-Borshchevo complex in Russia spans different periods of the UP and has yielded numerous human remains, including the burial of Kostenki 14 (~39 - 37 ka cal BP), which is also associated with the Aurignacian.<sup>11,12</sup> Isolated human remains from Buran Kaya III in Crimea<sup>13</sup> are associated with Gravettian assemblages, but this is in contrast with their early direct age of ~36.8–35.7 ka cal BP. In Central Europe, the Gravettian is split into the Early Gravettian (~34 - 30 ka cal BP), Pavlovian (~31 - 29 ka cal BP) and Late Gravettian/Willendorf Kostenkian (~29 - 24 ka cal BP).<sup>14,15</sup> The Pavlovian period is characterized by burials containing red ochre and grave goods such as seen at Dolní Věstonice-Pavlov<sup>16–19</sup> and Předmostí<sup>14,20</sup> in Czech Republic. Whilst UP child burials have been identified in the region, including the burial of infants at Krems-Wachtberg (Austria),<sup>21–23</sup> they are extremely rare. Much more elaborate UP burials have been found at Sunghir (Russia), containing multiple individuals, abundant ochre, and an exceptionally rich assemblage of grave goods.<sup>24,25</sup> Several attempts at directly dating these individuals produced vastly differing results, likely due to issues of contamination, thus contributing to uncertainty around the cultural association of the site. The most recent and rigorous compound specific dates place the burials at Sunghir between 35.5 and 32 ka cal BP.<sup>11,24–30</sup>

Borsuka Cave (50°9'53.94"N; 19°42'12.23"E), located in the Szklarka River valley in southern Poland (Figure 2A), adds important data to this discussion. Excavations next to the entrance of the cave in 2008–2010 (Figure 2B) uncovered six deciduous human teeth in association with 112 pierced teeth of Steppe wisent/aurochs (Bison priscus/Bos bison; n = 78) and European elk (Alces alces; n = 34), which have traces of ochre on

<footnote n="1"><sup>1</sup>Department of Human Evolution, Max Planck Institute for Evolutionary Anthropology, 04103 Leipzig, Germany</footnote>

<footnote n="2"><sup>2</sup>Ancient Genomics Lab, Francis Crick Institute, London NW1 1AT, UK</footnote>

<footnote n="3"><sup>3</sup>Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, 04103 Leipzig, Germany</footnote>

<footnote n="4"><sup>4</sup>Department of Cell and Molecular Biology, University of California, Berkeley, Berkeley, CA 94720-3200, USA</footnote>

<footnote n="5"><sup>5</sup>Department of Biology, San Francisco State University, San Francisco, CA 94132, USA</footnote>

<footnote n="6"><sup>6</sup>CEREGE, Aix Marseille Université, CNRS, IRD, INRA, Collège de France, Technopôle de l'Arbois BP 80, 13545 Aix-en-Provence Cedex 4, France</footnote>

<footnote n="7"><sup>7</sup>Chaire de Paléoanthropologie, CIRB (UMR 7241 - U1050), Collège de France, 75231 Paris, France</footnote>

<footnote n="8"><sup>8</sup>Institute of Systematics and Evolution of Animals, Polish Academy of Sciences, Sławkowska 17, 31-016 Krakow, Poland</footnote>

<footnote n="9"><sup>9</sup>These authors contributed equally</footnote>

<footnote n="10"><sup>10</sup>Lead contact</footnote>

<footnote>*Correspondence: helen.fewlass@crick.ac.uk (H.F.), ezavala9@berkeley.edu (E.I.Z.)
https://doi.org/10.1016/j.isci.2023.108283</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 1. Map of Upper Palaeolithic sites mentioned in the text from Central and Eastern Europe containing human remains and the archaeological complex with which they are currently associated**</figcaption></figure>

Sites number key: 1. Goyet Cave, 2. Krems-Wachtberg, 3. Dolní Věstonice, 4. Pavlov, 5. Mladeč, 6. Předmostí, 7. Oblazowa Cave, 8. Borsuka Cave, 9. Peştera Cioclovina, 10. Peştera Muierii, 11. Bacho Kiro Cave, 12. Buran Kaya III, 13. Kostenki, 14. Sunghir.

them.³¹,³² Seven layers were excavated and separated into a lower Pleistocene unit (Layers VII - V) and an upper Holocene unit (Layers IV - I). At the base of the sequence, overlying bedrock, Layer VII was archaeologically sterile. In Layer VI (from 150 to 250 cm depth), the human remains and pierced teeth were found spread over several m² in an NW linear distribution down the slope from the cave entrance. The pendants were concentrated by the southern wall of the trench at depths of 160–220 cm, with a few in the central and northern parts of the trench (Figure 2). The human teeth were located to the west of the concentration of pendants, in the northwest corner of the trench within squares C7, D7, and E7 (Figures 2C and 2D). These teeth were all identified as belonging to a child aged 12–18 months. Given the young age, it was not possible to determine the sex based on morphology, but the consistent phase of dentition indicated they likely belonged to one individual. Layer VI contained the human remains, perforated ungulate teeth, and one distal fragment of a flint blade but was otherwise archaeologically sterile. The overlying Layer V was also archaeologically and palaeontologically sterile but associated with Last Glacial Maximum (LGM) cooling. Layers IV-I contained Holocene assemblages spanning from the Mesolithic to modern (see³¹,³³).

In addition to the human remains, over 2,000 faunal skeletal fragments were excavated from Layer VI. The assemblage contained an unusual collection of fauna including cold steppe-tundra species, such as *Rangifer tarandus* (reindeer), *Vulpes alopex* (Arctic fox), *Equus ferus* (wild horse) and *Coelodonta antiquitatis* (woolly rhino), alongside taxa adapted to a warmer forested environment, including *Alces alces* (elk), *Bos primigenius* (aurochs), *Martes martes* (pine marten), *Lynx lynx* (Eurasian lynx), *Castor fiber* (beaver) and *Meles meles* (badger).³³ The majority of the Pleistocene skeletal remains likely represent a natural accumulation, with signs of carnivore modifications on 1.7% of the mammalian and bird bones and no signs of anthropogenic modification aside from the pierced teeth.

Two of the herbivore tooth pendants were directly radiocarbon dated to 32,890 - 30,760 cal BP at 95.4% probability (Poz-32394: 27,350 ± 450 ¹⁴C BP; ¹⁴C errors reported at 1σ) and 29,900–29,120 cal BP (Poz-38236: 25,150 ± 160 ¹⁴C BP). A reindeer metatarsus from the same layer was dated to 31,060 - 30,270 cal BP (Poz-38237: 26,430 ± 180 ¹⁴C BP).³²,³³ While the dates indicated a Mid-Upper Palaeolithic age for the Layer VI assemblage, the lack of agreement between the dates from the two pendants at the 95% range suggested that the pendants are not strictly contemporaneous. This, along with the distribution of the herbivore and human teeth, leaves open the question of the association between the pendants and the human remains, necessitating further chronological investigation. However, due to their small size, direct dating of the human teeth was initially determined to be impossible.

Although alternative theories have been discussed,³² the absence of associated “domestic” archaeology led to the interpretation that the assemblage represents an infant burial and grave goods that were disturbed post-deposition, resulting in the absence of a burial pit and their linear spread down the slope.³² The lack of any Palaeolithic archaeology necessitates relying on the dating and stylistic analysis of the pendants for determining the cultural association of the burial. The existing pendant ¹⁴C dates are roughly contemporaneous with Pavlovian burials ~300-200 km away at Dolní Věstonice-Pavlov and Předmostí (Czech Republic)¹⁴,¹⁶⁻²⁰ and infant burials ~400 km away at Krems-Wachtberg (Austria).²¹⁻²³ However, the typology of the pendants themselves are similar to those found at Mladec³² and suggest a Late Aurignacian association. In addition, there have been numerous examples where human remains assumed to be Palaeolithic in origin in fact
<figure><img src="image_6.png" /><figcaption>**Figure 2. Visuals of Borsuka Cave location, site plans, and stratigraphic profile**</figcaption></figure>

(A) Map showing the location of Borsuka Cave within Europe.
(B) Plan view of Borsuka Cave showing the 2008–2010 excavation trench (gray shaded rectangle).
(C) Stratigraphic profile with the location of ungulate tooth pendants (black circles) and human remains (red circles).
(D) Plan view of Layer VI showing the location of human remains and pendants. Note that three human teeth found *in situ* are plotted but three were found during wet-sieving so lack exact 3D coordinates. In C and D, each square is 1 m².

represent Holocene intrusions into Palaeolithic contexts.<sup>34–39</sup> Given the small nature of fragmented deciduous human teeth and presence of Holocene archaeology in the upper layers, the possibility of this being the case at Borsuka Cave could not be ruled out.

Significant advances have been made in recent years for the radiocarbon dating and genetic analysis of ancient and highly degraded skel-
etal material, including reducing sample sizes for dating,<sup>19,40–42</sup> combining multi-method extractions from a single sample,<sup>43</sup> non-destructive
pre-screening of organic preservation,<sup>44</sup> treatments for removing contaminant DNA and/or enriching for endogenous ancient DNA
(aDNA),<sup>45–49</sup> and non-destructive aDNA extraction from osseous artifacts.<sup>50</sup> Such developments enable biomolecular analyses of material
which were previously impossible. Considering the expected age and very limited amount of material available from the human and pierced
herbivore teeth from Borsuka, we devised a multi-disciplinary study to investigate (i) if the human teeth originate from a single individual, (ii)
the chronological relationship of the human teeth and pendants, and (iii) the association of the assemblage to different Palaeolithic cultures.
In order to confirm the Palaeolithic origin of the teeth and ascertain if the human remains are contemporaneous with the pendants, we un-
dertook radiocarbon dating of two deciduous human teeth (Figure 3) and six of the pierced herbivore teeth (Figure 4). From the human ma-
terial, we selected two teeth (C7/675, C7/683) with partially intact roots to preserve the crown morphology and selected herbivore pendants
from different squares to test across the scattered assemblage. Ancient DNA analysis was attempted on the same two human teeth to deter-
mine if they belonged to one individual and to see how this individual related to other Upper Palaeolithic individuals. The assemblage at
Borsuka Cave represents a challenge in terms of the extremely limited material available, and the inherent issues of poor preservation
and contamination of ancient samples for both <sup>14</sup>C and aDNA. Nevertheless, we aimed to maximise data retrieval from extremely small sam-
ples to provide additional information about the Borsuka assemblage within the broader context of the UP.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** The two human teeth from Borsuka cave analyzed in the study
(A) The deciduous molar C7/675 (uldm1) and (B) molar fragment C7/683. Scale bar is 1 cm.</figcaption></figure>

RESULTS

Radiocarbon dating

The collagen yield of the two human teeth (4.3–8.5%) was within the range of the herbivore teeth (3.0–9.1%), consistent with degraded Palae-
olithic material but all falling well above the ~1% minimum required for dating. The elemental values (C%, N%, C/N) of all extracts fall within
the commonly accepted ranges of well-preserved collagen (C%: ~30–45%; N%: ~11–16%; C/N: 2.9–3.6), indicating that the collagen extracts
were suitable for radiocarbon dating<sup>51</sup> (Table 1). However, several fall at the limit (3.5–3.6) of the accepted C/N range, indicating low levels of
contamination may be present in the samples.<sup>52,53</sup> For C7/675, C7/658, C6/494, sufficient collagen was available for additional quality analysis
using Fourier Transform Infrared Spectroscopy (FTIR). All extracts had FTIR spectra characteristic of ancient bone collagen,<sup>54–57</sup> but the higher
intensity of the peak at ~2900 cm<sup>-1</sup> in the spectra of the human tooth C7/675 may be indicative of external contamination (Figure S1).

The bulk collagen δ¹³C and δ¹⁵N values of the two human teeth fall within the range of stable isotopic values seen for other Mid-Upper Palaeolithic humans in Eurasia,¹⁹,²³,⁵⁸⁻⁶² with very high δ¹⁵N values that are well above the herbivorous signature of the Bos/Bison and elk teeth (Figure S2). Although the δ¹³C values of the two human teeth agree within instrumental error, the δ¹⁵N value of C7/675 is higher than C7/683 by 1.3‰. Given that the teeth are both deciduous molars with similar formation times, the difference is more likely to be indicative of contamination, or of separate individuals, rather than palaeodietary. The elevated δ¹⁵N values are at the highest end of the range seen in other Mid-UP humans and in the range of two trophic levels (trophic level increase: ~3‰–5‰) above the herbivores. This is consistent with a breastfeeding signal given the age of the infant, but without comparative adult data and further chemical analyses it is not possible to elucidate this further.

Given the small size of the dentine samples (23.4–67.4 mg) and extracted collagen (1–4.5 mg), we followed a previously established approach for small samples.<sup>41,42</sup> Multiple <sup>14</sup>C dates were obtained from each collagen extract using both graphitization and CO<sub>2</sub> gas ion source (GIS) methods, and subsequently combined. X<sup>2</sup> tests were carried out to determine the statistical agreement between replicate dates from the same extract. All of the <sup>14</sup>C dates obtained fall within the Mid-UP range, from 31,070 ± 770 <sup>14</sup>C BP to 24,830 ± 290 <sup>14</sup>C BP, or ~35,290 to 29,090 cal BP (Figure 5; Table 1; Tables S1 and S2). For the herbivore pendants, all six pendants were dated to 35,290 – 33,260 cal BP (95.4% probability), much older than dates obtained previously on two of the pendants, overlapping with the Late Aurignacian/Early Gravettian periods. For the two human teeth, replicate measurements agree statistically for each collagen extract but the dates from C7/683 (Aix-12058 combined age: 26,610 ± 240 <sup>14</sup>C BP; 31,160 – 30,320 cal BP at 95.4% probability) are older than C7/675 (Aix-12047 combined age: 25,100 ± 140 <sup>14</sup>C BP; 29,860 – 29,090 cal BP at 95.4% probability) by ~1500 <sup>14</sup>C years (Figure 5; Table S1). When the dates from Layer VI obtained thus far are modeled together in one phase in OxCal, the assemblage spans 6,290–4,770 years, between 36,650 and 27,620 cal BP (Figure 5; Table S3). No dates are available from the overlying and underlying layers to constrain the modeled range, due to a lack of material suitable for dating so none of the dates were identified as outliers. The radiocarbon data is included in Table S1, the background collagen measurements in Table S2 and the modeled ranges of the Borsuka data in Table S3.

Ancient DNA analysis

Two microsamples between 4.7 and 12.9 mg of dentine powder were drilled from each of the two human teeth (C7/675 and C7/683). The resulting four tooth powder sub-samples underwent DNA extraction, single-stranded library preparation and shotgun sequencing (see method details). Between 1,628,190 and 8,575,605 sequences were generated per library, of which between 2,495 and 37,436 mapped to the human reference genome (hg19⁶³) for each library. While no ancient DNA was detected in C7/683, one of the sub-samples from C7/675 was identified as containing ancient DNA based on the observation of between 12.4% and 20.6% C-to-T substitutions on the terminal ends of recovered DNA (Table S4), which is characteristic of ancient DNA.⁶⁴,⁶⁵ The duplication rates, or number of times each DNA fragment was sequenced, for these libraries were all around 1 (Table S4), indicating that not all recovered DNA fragments were sequenced. Modern human DNA contamination for libraries produced from this subsample was estimated at 54–56% (Table S4) using the method AuthentiCT⁶⁶ which calculates levels of contamination based on deamination patterns.

Despite not detecting ancient DNA in C7/683 to the limits of our resolution, libraries from both teeth were subsequently enriched for hu-
man mitochondrial (mt) DNA via hybridization capture,⁶⁷ in an attempt to determine if both likely stem from the same individual (Table S5).
Libraries from C7/683 were still included in hybridization captures as the shotgun sequencing was not exhaustive, and the content of nuclear
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Four pendants made of large ungulate teeth from Borsuka Cave dated in the study</figcaption></figure>

(A) European elk (Alces alces) left lower incisor, sample ID: C7/656; (B) European elk (Alces alces) left lower incisor, sample ID: B5/910; (C) Bos/Bison (Bos primigenius/Bos bison) right lower incisor, sample ID: B5/913; (D) Bos/Bison (Bos primigenius/Bos bison), sample ID: C7/658. Scale bar is 2 cm.

and mitochondrial DNA has been shown to defer in ancient specimens,⁶⁸ leaving the possibility that small amounts of endogenous DNA remained undetected. The enriched libraries contained between 243,078 and 938,984 sequences longer than 34 base pairs that mapped to the human mtDNA revised Cambridge Reference Sequence⁶⁹ with elevated C-to-T substitutions (21.3–35.1% on the 5' end and 14.5 to 24.6% on the 3' end). As modern human contamination was detected in each library, ranging from 2.93% to 39.94%, all subsequent analyses were restricted to putatively deaminated fragments that contained a C-to-T substitution within the first three and/or last three terminal bases. Between 15,269 and 60,728 deaminated fragments (corresponding to mtDNA coverages of 46- and 223-fold, respectively) were used to reconstruct near-complete consensus mtDNA genomes. The resulting mtDNA haplotypes from each tooth were the same, consistent with the teeth either being from the same individuals or the same maternal lineage. However, as neither mtDNA genome was complete, we cannot exclude that individually discriminating variants were not detected. The more complete mtDNA genome from C7/675 was then used for tree building and molecular branch shortening (Tables S6 and S7). The resulting molecular date was estimated at 33,533 years BP (95% highest posterior
<table><caption>Table 1. Pretreatment and radiocarbon dating of humans and herbivore pendants from Borsuka Cave</caption><thead><tr><th>Sample ID</th><th>Material</th><th>Sampled (mg)</th><th>Collagen (mg)</th><th>% Collagen</th><th>δ<sup>13</sup>C (‰)</th><th>δ<sup>15</sup>N (‰)</th><th>%C</th><th>%N</th><th>C/N</th><th>AMS lab number</th><th>combined BP</th><th>combined err</th><th>95.4% cal BP</th><th>n</th></tr></thead><tbody><tr><td>C7/675</td><td>Human tooth - uldm1</td><td>41.1</td><td>3.5</td><td>8.5</td><td>-20.3</td><td>14.9</td><td>42.9</td><td>13.7</td><td>3.6</td><td>Aix-12047</td><td>25100</td><td>140</td><td>29860-29090</td><td>2</td></tr><tr><td>C7/683</td><td>Human tooth - dm</td><td>23.4</td><td>1</td><td>4.3</td><td>-20.1</td><td>13.6</td><td>39.7</td><td>13.3</td><td>3.5</td><td>Aix-12058</td><td>26610</td><td>240</td><td>31160-30320</td><td>3</td></tr><tr><td>C7/656</td><td>Elk incisor pendant</td><td>34.3</td><td>2.6</td><td>7.6</td><td>-20.6</td><td>4.7</td><td>42.3</td><td>14.2</td><td>3.5</td><td>Aix-12041</td><td>29310</td><td>240</td><td>34350-33260</td><td>2</td></tr><tr><td>C7/658</td><td>Bos/Bison incisor pendant</td><td>48.5</td><td>4.4</td><td>9.1</td><td>-20.3</td><td>7.4</td><td>42.7</td><td>15.2</td><td>3.3</td><td>Aix-12042</td><td>30100</td><td>260</td><td>35190-34150</td><td>2</td></tr><tr><td>C6/479</td><td>Bos/Bison incisor pendant</td><td>53.9</td><td>3.6</td><td>6.7</td><td>-19.6</td><td>6.5</td><td>41.8</td><td>15</td><td>3.3</td><td>Aix-12043</td><td>30080</td><td>260</td><td>35170-34130</td><td>2</td></tr><tr><td>C6/494</td><td>Elk incisor pendant</td><td>67.4</td><td>4.5</td><td>6.7</td><td>-19.8</td><td>2.3</td><td>42.7</td><td>15.1</td><td>3.3</td><td>Aix-12044</td><td>30460</td><td>190</td><td>35290-34480</td><td>3</td></tr><tr><td>B5/910</td><td>Elk incisor pendant</td><td>49.2</td><td>1.5</td><td>3</td><td>-20.4</td><td>3.1</td><td>37.8</td><td>13.4</td><td>3.3</td><td>Aix-12045</td><td>29310</td><td>370</td><td>34610-33860</td><td>3</td></tr><tr><td>B5/913</td><td>Bos/Bison incisor pendant</td><td>49.9</td><td>2.9</td><td>5.8</td><td>-20.1</td><td>6.3</td><td>40.4</td><td>14.5</td><td>3.3</td><td>Aix-12046</td><td>30080</td><td>260</td><td>35170-34140</td><td>2</td></tr></tbody></table>

Combined dates of multiple replicates (n) are reported here with individual replicate dates and X² values reported in Table S1.
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Single phase outlier model showing the posterior probability distributions of the ¹⁴C dates from Borsuka Cave Layer VI generated in this study (Aix-) and previously (Poz-)</figcaption></figure>

The bracket beneath each distribution shows the 95% probability range. The outlier probability (O) is given in the format "[O: posterior outlier probability/prior outlier probability]." The ranges are given in Table S3.

density interval: 28,200 - 38,935 years BP) and the mtDNA genome falls within haplogroup U6 (Figure 6A). Haplogroup U6 is most commonly observed in Northern Africa in present-day humans,⁷⁰,⁷¹ and has previously only been observed in Palaeolithic Europe in the specimens from the site of Peştera Muierii in Romania (~34,000 years old).⁹,⁷²,⁷³ Similar to the Muierii 1 and 2 mtDNA genomes, which were subsequently determined to belong to the same individual,⁷³ the Borsuka mtDNA genome is basal to the U6 haplogroup, sharing the 3348G, 10517A, and 16172C positions. Importantly, among the overlapping positions (6 missing positions) no differences are observed between the Muierii 2 and Borsuka mtDNA genomes. Even if both Borsuka teeth and Muierii 2 share an mtDNA genome, this only indicates that they likely lived within ~2,500 years of each other (Figure S3, method details).

In order to determine the genetic affinity of the child buried in Borsuka Cave to other individuals from the Upper Palaeolithic, we enriched the libraries from C7/675 for ~1.2 million single nucleotide polymorphisms (SNPs) in the human genome known to be informative for studying human population history.⁷⁴,⁷⁵ While the resulting libraries covered 513,237 SNPs, they were also identified as containing approximately 56% modern human contamination (Table S8). This required all downstream analysis to be restricted to the 46,286 putatively deaminated fragments, which were used to determine the sex of the child. As the ratio of autosome, X, Y chromosome SNPs in the 1240k does not follow conventional expectations due to the respective target size, we used an adjustment of the X and Y-rate as has been performed previously with this type of data.⁹ This resulted in an "X-rate" of 0.68 and "Y-rate" of 0.04, consistent with the individual being female (Table 2). When projected onto a principal component analysis (PCA) of modern-day individuals from West Eurasia, Central Asia and East Asia, the Borsuka individual clustered with previously published Upper Palaeolithic Western Eurasians and near modern day Western Eurasians (Figure 6B). The genetic sharing between the Borsuka individual and previously published modern and ancient individuals was calculated with f3-statistics, using present-day Mbuti individuals as an outgroup population (Figure 6A). Consistent with the PCA analysis, the greatest affinity among modern day populations was to Western Eurasians. Basal Eurasian ancestry was not detected in the Borsuka individual (Tables S9 and S10, method details). Among ancient individuals, the Borsuka individual shared the most alleles with the ~35 ka cal BP Bacho Kiro Cave individual (BK1653), ~34 ka cal BP Muierii 1, ~31 ka cal BP Vêstonice 16 and ~34 ka cal BP Sunghir 3 individuals when tested with f3-statistics (Table S11). When directly comparing the genetic affinity of the Borsuka individual with other ancient individuals with D-statistics, the Borsuka female has a greater affinity to the Gravettian and Aurignacian individuals (Figures S4 and S5; Table S12). However, no significant difference was observed in the affinity of the Borsuka child to the Vêstonice vs. Sunghir vs. BK1653 vs. Muierii 1 individuals, precluding the direct association of the Borsuka girl with one of these groups.

DISCUSSION

The chronological and genetic data presented here confirms the Upper Palaeolithic origin of the infant remains from Borsuka Cave, making these
remains currently the oldest female infant burial identified to date. While it is not possible to confirm that the teeth come from the same individual
without nuclear DNA, the mitochondrial analysis is consistent with the two teeth belonging to the same individual or the same maternal lineage.

Radiocarbon dating is based on the exponential decay of ¹⁴C over time. As the concentration of ¹⁴C in modern carbon is therefore much
larger than in ancient samples, any contamination with modern carbon will make ¹⁴C dates younger than the true age of the sample, with the
effects getting progressively worse for more ancient material. Given the ubiquity of modern carbon introduced via the burial environment,
handling, storage and analysis of artifacts, when inconsistent dating results are obtained from Palaeolithic material, older ¹⁴C ages are gener-
ally considered more accurate.²⁹,⁷⁶ The direct date from human tooth C7/683 is older than the date from C7/675. If the two teeth originate
from the same individual, this implies that either 1) the older age of C7/683 is correct and the date from C7/675 is an under-estimation of the
true age due to the presence of external carbon contamination, or 2) both dates are under-estimations.
A

f3(X, Borsuka; Mbuti)

<figure><img src="image_11.png" /></figure>

B

<figure><img src="image_12.png" /></figure>

C

<figure><img src="image_13.png" /><figcaption>Figure 6. The genetic analysis of the Borsuka C7/675 individual in relation to other modern and ancient humans</figcaption></figure>

(A) The calculated allele sharing between the Borsuka individual and a selection of modern and ancient individuals (X) where warmer colors represent higher amounts of shared alleles. The f3 statistic f3(X, Borsuka; Mbuti) was calculated using 28,213 - 31,393 overlapping SNPs for modern populations and 6,697 - 28,642 SNPs for ancient individuals. Ancient individual key: 1- Goyet Q116, 2- EL Miron, 3- Dolní Věstonice16, 4- Bichon, 5- Oase1, 6- Bichon, 7- Mal'ta1, 8- Karelia, 9- Yana2, 10- Tianyuan, 11- Ust-Ishim, 12- Kostenki14, 13- Sunghir3.

(B) The projection of 20 ancient individuals onto a principal component analysis of 1,267 modern individuals. Ancient individuals are colored based on their association with the Aurignacian (dark blue), Gravettian (orange), or IUP (light blue) archaeological assemblages or in green if their technological association is uncertain.

(C) A subset of the phylogenetic tree generated with BEAST2. Ancient humans are colored in blue (Borsuka in red). Nodes are labeled with their posterior probability. The x axis indicates calibrated years before present.
<table><caption>Table 2. Sex determination based on coverage of nuclear SNPs in the '1240k' array using putatively deaminated reads</caption><thead><tr><th>Sample ID</th><th>Material</th><th>Nyo (# bases covered on Y chromosome)</th><th>Nxo (# bases covered on X chromosome)</th><th>Nautoo (# bases covered on autosomal chromosomes)</th><th>Ny expected (targets on Y chromosome)</th><th>Nx expected (targets on X chromosome)</th><th>Nauto expected (targets on autosomal chromosomes)</th><th>X-ratio (Nxo/Nautoo)/(Nxe/Nautoe)</th><th>Y-ratio (Nyo/ Nautoo)/(Nye/ Nautoe)</th></tr></thead><tbody><tr><td>C7/675</td><td>Human tooth - left udm1</td><td>49</td><td>1,312</td><td>45,032</td><td>32,670</td><td>49,704</td><td>1,150,639</td><td>0.674</td><td>0.038</td></tr></tbody><tfoot><tr><td colspan="10">A Y-ratio &lt;0.05 indicates a female and &gt;0.2 indicates a male.</td></tr></tfoot></table>
Alternatively, the teeth could originate from two infants and both dates could be correct. In addition to the discrepancy in age between the direct human dates and the range of ages from the pendant, this may also be indicated by the difference in δ¹⁵N values between the teeth and the scatter of the assemblage over several square meters. Given the consistent form of pendant manufacture and that the mtDNA analysis indicates the teeth are at least from maternal relatives (if not from the same individual), a group from the same maternal lineage would need to have used the site for 4,000–6,000 years to account for the dating results in this scenario. The elevated C/N value and additional FTIR peak at ~2900 cm⁻¹ in the C7/675 extract do suggest contamination may have affected the date, and potentially the δ¹⁵N value. Thus, the most parsimonious interpretation of the results is that the two teeth originate from the same individual and are affected by contamination. If we consider the teeth as originating from one individual, the older date from C7/683 is more likely to be accurate than the younger date from C7/675 as contamination typically decreases the age. Given the small amount of material available for sampling and the C/N values from both teeth of 3.5–3.6, we consider the ~31.2–30.3 ka cal BP date from C7/683 to be a minimum age for the human remains.

The older date from C7/683 is contemporaneous with the dates previously obtained from two Borsuka pendants and the ~31 ka cal BP Pavlovian burials at Dolní Věstonice-Pavlov and Krems-Wachtberg, which would support the association of Borsuka with the Pavlovian culture in Central Europe (Figure 7). However, the new ¹⁴C dates on the herbivore tooth pendants range in age from 35,290 to 33,260 cal BP (95.4% probability). These new dates for the pendants are older than the ages obtained from the human infant teeth, but also much older than the dates previously obtained on two pendants and a reindeer bone from the same layer (Figure 5).

It is possible that the burial assemblage dates to ~31 - 30 ka cal BP, and the necklace was constructed from a collection of herbivore teeth of varying age, including some that were pierced and used for the construction of an ornament long after the death of the animal. The solitary life-style of European elk (Alces alces) implies that the elk teeth were unlikely to have been collected during one event. A recent study of the Sr isotope compositions and trace element analysis of the enamel from four of the elk tooth pendants (at least 3 individuals) demonstrated a non-local origin of the teeth, indicating that the pendants were transported ~250 km to Borsuka Cave from an area near the Austria/Slovakia border or northern Hungary, on the southern side of the Western Carpathians.⁷⁷ This evidence for exchange or human regional mobility is supported by Borsuka Cave being the only site in the region containing evidence of the presence of elk ~27 ka BP (beginning of MIS2) despite favorable environmental conditions.³³,⁷⁸⁻⁸¹ The presence of Trans-Carpathian raw stone materials in Central Europe also shows high human mobility around the Western Carpathian mountain region.⁸⁰,⁸² This evidence suggests that the collection of the herbivore teeth over a period of time is feasible, however, it seems unlikely that the pendants represent animals varying in age by over 5,000 years. Alternatively, low levels of contamination may be affecting some of the ¹⁴C results with the ~32 - 29 ka cal BP ages under-estimating the true age of the pendants and human teeth.

The older ~35.3–33.3 ka cal BP range of dates from the pendants overlap in time with the ~36.8–34.7 ka cal BP range of directly dated human remains from Mladeč Cave (Czech Republic)¹ (Figure 7). It has been suggested that the relatively shallow human remains at Mladeč could be the result of intentional deposition,³⁵,⁸³,⁸⁴ but the lack of contextual information means this hypothesis cannot be tested. Despite their uncertain context, based on the dates, the Mladeč humans have been associated with Aurignacian artifacts found elsewhere in the cave, including Mladeč-type bone points and pendants made from the teeth of large ungulates (e.g., elk, horse, bison/aurochs).⁸⁵ In comparison to other UP perforated teeth, the Borsuka Cave pendants show close similarity to perforated teeth discovered at Mladeč and generally from Aurignacian sites.⁸⁵⁻⁸⁷ The entire roots of the perforated teeth were heavily scraped on both surfaces, remarkably reducing the root thickness and the perforations were made very close to the root apex by drilling from both sides. The only difference between pendants discovered at Borsuka and Mladeč can be observed in the larger range of species from which teeth were used at Mladeč. Thereby the assemblage from Borsuka cave represents the largest collection of pendants made of large ungulate teeth potentially associated with the Aurignacian culture. Currently, no aDNA is available from the Mladeč individuals to investigate this link further.

While the origin of the Gravettian has been widely debated,¹⁴ the relatively large collection of burials across Eurasia has provided a significant body of morphological and behavioral data on Gravettian humans across a wide geographical and chronological range. Recent aDNA studies¹⁰,⁸⁸ have demonstrated a genetic distinction between the Gravettian cultures in Central Eastern Europe and Southwestern Europe. The oldest humans discovered to date with genetic continuity to Gravettians and modern West Eurasians are represented by three different genetic ancestry components from the ~35,000-year-old GoyetQ116-1 individual,⁹ the ~35,000-year-old individual from Bacho Kiro Cave (BK1653)⁸⁹ and the ~38,000-year-old Kostenki 14⁹.⁹⁰ The BK1653 individual is related, but not identical, to the genetic ancestry of GoyetQ116-1 which contributes to the genetic ancestry of the Gravettian individuals in Southwestern Europe, while the former contributes ancestry to the Central European Gravettian individuals. Genetic ancestry from Kostenki 14 has been found in all individuals associated with the Gravettian culture. We aimed to clarify the genetic relationship of the Borsuka individual to individuals of different Palaeolithic cultures, however these associations were limited due to the low amounts of data recovered. The Borsuka individual shares a basal U6 mtDNA haplotype with Peştera Muierii 2, but there is not sufficient nuclear data to determine the full genetic relationship between these individuals in context with other ancient humans. The available data show that the Borsuka individual has the most nuclear genetic affinity to the Bacho Kiro 1653 individual (presumed to be associated with the Aurignacian), as well as individuals from Dolní Věstonice (Gravettian) and Sunghir (early UP). Interestingly, the Borsuka individual is significantly closer to the Bacho Kiro individual than they are to Kostenki 14. This is consistent with Borsuka being an UP Central European, but does not determine a closer affinity to Gravettian or Aurignacian individuals within Central Europe. The ¹⁴C dates provide a minimum age of 31.3–30.3 ka cal BP for the human remains, which agrees with the molecular date estimation of 33.5 ka that falls roughly at the border between the Late Aurignacian and Early Gravettian periods.

The Borsuka Cave assemblage represents an important data point in the catalog of Palaeolithic human remains for investigating the geographical and chronological range of UP burial practices. Of the very limited number of Upper Palaeolithic infant burials discovered and successfully analyzed with aDNA, this is the first female identified. The previously oldest confirmed female infant burial was found in
<figure><img src="image_14.png" /><figcaption>**Figure 7.** Calibrated radiocarbon dates (95.4% range) from Borsuka Cave infant teeth and pendants compared to UP human burials and human remains in Central and Eastern Eurasia mentioned in the text</figcaption></figure>

Human remains tentatively associated with Aurignacian assemblages are shown in blue and Gravettian assemblages in orange while uncertain cultural associations are shown in gray. The Borsuka infant dates are conservatively considered minimum ages.

Mesolithic Italy,⁹¹ at least 20,000 years after the burial of the Borsuka girl. This study therefore provides the first evidence that child burials incorporating grave goods were not limited to males during the UP. This study further highlights the challenges of working with limited and poorly preserved ancient materials, and emphasizes that employing an interdisciplinary, multi-method approach can still provide valuable insights to allow in-depth investigations of Palaeolithic contexts.

## Limitations of the study

Palaeolithic human remains are rare, and the remains of UP infants even more so. The small size and partial fragmentation of the infant teeth presented an exceptionally limited amount of suitable material for analysis. We aimed to take as little material as possible and therefore limited destructive sampling to two of the six teeth to preserve the morphology of the crowns and ensure their preservation for future studies.
Modern carbon contamination in a ¹⁴C sample produces an under-estimation of the true age of the material. As sample size decreases, the risk of contamination increases, which is clearly a consideration for the data presented. This is highlighted by the lack of agreement in dates from the two human teeth from the (potentially) same individual and their slightly elevated C/N values, indicating that low levels of C contamination may be present either due to incomplete removal of environmental contaminants during pretreatment or introduced in the laboratory. Given the ubiquity of carbon contaminants in the environment and laboratory, older ages produced from the same sample material are generally considered to be more accurate and reliable. Rigorous compound specific dating approaches can be employed to isolate endogenous hydroxyproline in bone collagen to overcome contamination issues.²⁷,⁹²,⁹³ However, the small proportion of hydroxyproline in collagen means this approach requires large starting sample sizes, thus making the approach not feasible for this study. As contamination of the human dates is indicated we therefore consider them to be minimum ages. The molecular dating technique (calibrated using other directly radio-carbon dated human remains) indicates an age of ~33.5 ka BP (albeit with a very large confidence interval), which supports the ~31 ka cal BP being a minimum age. The new data presented nevertheless confirm the UP origin of the human remains. Given the challenges and the limited data available, we cannot conclude a definitive age for the burial assemblage or a more definitive association with a specific Palaeolithic culture. However, even using very small sample sizes for dating (in the order of aDNA sample sizes), we confirm the placement of the assemblage around the border of the Late Aurignacian/Early Gravettian, supported by the aDNA analysis and the typology of the pendants.

The high levels of modern human contamination and limited amount of endogenous DNA present in both infant teeth limited the power of determining different genetic affinities beyond broad population genetic relationships. This resulted in reduced coverage of both the mtDNA and nuclear genomes. As shown in Figures S6-S9, decreased amounts of data result in missing some genetic affinities. It is likely that with more data it would have been possible to deconvolute the association of the Borsuka girl to Bacho Kiro 1653, or the Vêstonice and Sunghir clusters, but unfortunately the limited data prevented this.

The limitations arising due to contamination emphasize the importance of taking efforts to minimize the contamination of archaeological materials. This includes wearing gloves and face masks during excavation and handling of human remains and artifacts, appropriate storage and consideration of sampling prior to consolidants or glues being applied.

## STAR★METHODS

Detailed methods are provided in the online version of this paper and include the following:

* **KEY RESOURCES TABLE**

* **RESOURCE AVAILABILITY**

* Lead contact

* Materials availability

* Data and code availability

* **EXPERIMENTAL MODEL AND SUBJECT DETAILS**

* Archaeological information

* **METHOD DETAILS**

* Radiocarbon dating

* Ancient DNA analysis

* **QUANTIFICATION AND STATISTICAL ANALYSIS**

* Radiocarbon data analysis and reporting

* Mitochondrial DNA analysis

* Estimating TMRCA for two mtDNA genomes with no differences

* Nuclear DNA analysis

## SUPPLEMENTAL INFORMATION

Supplemental information can be found online at https://doi.org/10.1016/j.isci.2023.108283.

## ACKNOWLEDGMENTS

This research was funded by the Max Planck Society. The AixMICADAS was acquired and operated in the framework of the EQUIPEX project ASTER-CEREGE (P.I.: E. Bard) with additional matching funds from the Collège de France. HF was partially funded by an EMBO post-doctoral fellowship (EMBO ALTF 590-2021). Elena I. Zavala's work was partially funded by the Miller Institute for Basic Research in Science, University of California Berkeley. Genetic data was produced by the Ancient DNA Core Unit of the Max Planck Institute for Evolutionary Anthropology which is funded by the Max Planck Society. The authors would like to thank Laurits Skov for input in the TMRCA analysis.

## AUTHOR CONTRIBUTIONS

This study was conceived by H.F., E.I.Z., and J.W. Data generation and analysis were performed by H.F., E.I.Z., Y.F., T.T., E.B., J-J.H. and M.H. The article was written by H.F., E.I.Z., and J.W. with input by all authors.
## DECLARATION OF INTERESTS

The authors declare no competing interests.

Received: April 1, 2023

Revised: October 16, 2023

Accepted: October 18, 2023

Published: October 23, 2023

## REFERENCES

1. Wild, E.M., Teschler-Nicola, M., Kutschera, W., Steier, P., Trinkaus, E., and Wanek, W. (2005). Direct dating of Early Upper Palaeolithic human remains from Mladeč. Nature 435, 332–335.

2. Valde-Nowak, P., Nadachowski, A., and Wolsan, M. (1987). Upper Palaeolithic boomerang made of a mammoth tusk in south Poland. Nature 329, 436–438.

3. Hedges, R.E.M., Housley, R.A., Pettitt, P.B., Bronk-Ramsey, C., and van Klinken, G.J. (1996). Radiocarbon dates from the oxford AMS system: Archaeometry datelist 21. Archaeometry 38, 181–207.

4. Trinkaus, E., Haduch, E., Valde-Nowak, P.W., and Wojtal, P. (2014). The Oblazowa 1 early modern human pollical phalanx and Late Pleistocene distal thumb proportions. Homo 65, 1–12.

5. Wilczyński, J. (2015). The Gravettian and Epigravettian settlement of Poland. In Forgotten Times and Spaces, S. Sázelová, M. Novák, and A. Mizerovápp, eds. Institute of Archaeology of the CAS (Brno & Masaryk University), pp. 191–213.

6. Wilczyński, J. (2016). A new beginning: Modern humans in Poland. In The past societies Volume 1, J. Kabaciński, ed. (Institute of Archaeology and Ethnology, Polish Academy of Sciences).

7. Soficaru, A., Doboş, A., and Trinkaus, E. (2006). Early modern humans from the Peştera Muierii, Baia de Fier, Romania. Proc. Natl. Acad. Sci. USA 103, 17196–17201.

8. Soficaru, A., Petrea, C., Doboş, A., and Trinkaus, E. (2007). The Human Cranium from the Peştera Cioclovina Uscată, Romania: Context, Age, Taphonomy, Morphology, and Paleopathology. Curr. Anthropol. 48, 611–619.

9. Fu, Q., Posth, C., Hajdinjak, M., Petr, M., Mallick, S., Fernandes, D., Furtwängler, A., Haak, W., Meyer, M., Mittnik, A., et al. (2016). The genetic history of Ice Age Europe. Nature 534, 200–205.

10. Posth, C., Yu, H., Ghalichi, A., Rougier, H., Crevecoeur, I., Huang, Y., Ringbauer, H., Rohrlach, A.B., Nägele, K., Villalba-Mouco, V., et al. (2023). Palaeogenomics of Upper Palaeolithic to Neolithic European hunter-gatherers. Nature 615, 117–126.

11. Marom, A., McCullagh, J.S.O., Higham, T.F.G., Sinitsyn, A.A., and Hedges, R.E.M. (2012). Single amino acid radiocarbon dating of Upper Paleolithic modern humans. Proc. Natl. Acad. Sci. USA 109, 6878–6881.

12. Dinnis, R., Bessudnov, A., Reynolds, N., Devièse, T., Pate, A., Sablin, M., Sinitsyn, A., and Higham, T. (2019). New data for the Early Upper Paleolithic of Kostenki (Russia). J. Hum. Evol. 127, 21–40.

13. Prat, S., Péan, S., Crépin, L., Puaud, S., Drucker, D.G., Lázničková-Galetová, M., van der Plicht, J., Valladas, H., Verna, C., Patou-

Mathis, M., et al. (2018). The first anatomically modern humans from South-Eastern Europe. Contributions from the Buran-Kaya III site (Crimea). Bull. Mem. Soc. Anthropol. Paris 30, 169–179.

14. Svoboda, J. (2007). The Gravettian on the middle Danube. Paléo, 203–220.

15. Wilczyński, J., Goslar, T., Wojtal, P., Oliva, M., Göhlich, U.B., Antl-Weiser, W., Sída, P., Verpoorte, A., and Lengyel, G. (2020). New Radiocarbon Dates for the Late Gravettian in Eastern Central Europe. Radiocarbon 62, 243–259.

16. Sládek, V., Trinkaus, E., Hillson, S.W., Holliday, T.W., and Svoboda, J. (2000). The People of the Pavlovian: Skeletal Catalogue and Osteometrics of the Gravettian Fossil Hominids from Dolní Věstonice and Pavlov. Academy of Sciences of the Czech Republic (Institute of Archaeology).

17. Svoboda, J. (2016). Dolní Věstonice II: Chronostratigraphy, Paleoanthropology, Academy of Sciences of the Czech Republic (Institute of Archaeology).

18. Svoboda, J., Novák, M., Sázelová, S., and Demek, J. (2016). Pavlov I: A large Gravettian site in space and time. Quat. Int. 406, 95–105.

19. Fewlass, H., Talamo, S., Kromer, B., Bard, E., Tuna, T., Fagault, Y., Sponheimer, M., Ryder, C., Hublin, J.-J., Perri, A., et al. (2019). Direct radiocarbon dates of mid Upper Palaeolithic human remains from Dolní Věstonice II and Pavlov I, Czech Republic. J. Archaeol. Sci.: Reports 27, 102000.

20. A Svoboda, J. (2008). The Upper Paleolithic burial area at Predmosti: ritual and taphonomy. J. Hum. Evol. 54, 15–33.

21. Einwögerer, T., Friesinger, H., Händel, M., Neugebauer-Maresch, C., Simon, U., and Teschler-Nicola, M. (2006). Upper Palaeolithic infant burials. Nature 444, 285.

22. Einwögerer, T., Händel, M., Neugebauer-Maresch, C., Simon, U., Steier, P., Teschler-Nicola, M., and Wild, E.M. (2009). 14C Dating of the Upper Paleolithic Site at Krems-Wachtberg, Austria. Radiocarbon 51, 847–855.

23. Teschler-Nicola, M., Fernandes, D., Händel, M., Einwögerer, T., Simon, U., Neugebauer-Maresch, C., Tangl, S., Heimel, P., Dobsak, T., Retzmann, A., et al. (2020). Ancient DNA reveals monozygotic newborn twins from the Upper Palaeolithic. Commun. Biol. 3, 650.

24. Pettitt, P.B., and Bader, N.O. (2000). Direct AMS Radiocarbon dates for the Sungir mid Upper Palaeolithic burials. Antiquity 74, 269–270.

25. Dobrovolskaya, M., Richards, M.P., and Trinkaus, E. (2012). Direct radiocarbon dates for the Mid Upper Paleolithic (eastern Gravettian) burials from Sunghir, Russia. Bull. Mém. Soc. Anthropol. 24, 96–102.

26. Kuzmin, Y., Burr, G.S., Jull, A.J.T., and Sulerzhitsky, L.D. (2004). AMS 14C age of the Upper Palaeolithic skeletons from Sungir site, Central Russian Plain. Nucl. Instrum. Methods Phys. Res. B 223–224, 731–734.

27. Nalawade-Chavan, S., McCullagh, J., and Hedges, R. (2014). New hydroxyproline radiocarbon dates from Sungir, Russia, confirm early Mid Upper Palaeolithic burials in Eurasia. PLoS One 9, e76896.

28. Reynolds, N., Dinnis, R., Bessudnov, A.A., Devièse, T., and Higham, T. (2017). The Kostënki 18 child burial and the cultural and funerary landscape of Mid Upper Palaeolithic European Russia. Antiquity 91, 1435–1450.

29. Higham, T.F. (2019). Removing contaminants: a restatement of the value of isolating single compounds for AMS dating. Antiquity 93, 1072–1075.

30. Soldatova, T. (2019). Spatial Distribution and Problems in the Interpretation of Radiocarbon Dates of the Sungir Site, Russia. Radiocarbon 61, e1–e15.

31. Wilczyński, J., Wojenka, M., Wojtal, P., Szczepanek, A., and Sobieraj, D. (2012). Human occupation of the Borsuka Cave (Southern Poland) from Upper Paleolithic to the Post-Mediaeval Period. Eurasian Prehistory 9, 77–91.

32. Wilczyński, J., Szczepanek, A., Wojtal, P., Diakowski, M., Wojenka, M., and Sobieraj, D. (2016). A mid upper Palaeolithic child burial from Borsuka Cave (southern Poland). Int. J. Osteoarchaeol. 26, 151–162.

33. Wilczyński, J., Miękina, B., Lipecki, G., Ługas, L., Marciszak, A., Rzebik-Kowalska, B., Stworzewicz, E., Szyndlar, Z., and Wertz, K. (2012). Faunal remains from Borsuka Cave – an example of local climate variability during Late Pleistocene in southern Poland. Acta Zool. Cracov. 55, 131–155.

34. Tillier, A.-M., Mester, Z., Bocherens, H., Henry-Gambier, D., and Pap, I. (2009). Direct dating of the “Gravettian” Balla child’s skeleton from Bükk Mountains (Hungary): unexpected results. J. Hum. Evol. 56, 209–212.

35. Svoboda, J.A., van der Plicht, J., and Kuželka, V. (2002). Upper Palaeolithic and Mesolithic human fossils from Moravia and Bohemia (Czech Republic): some new 14C dates. Antiquity 76, 957–962.

36. Trinkaus, E., and Pettitt, P.B. (2000). The Krems-Hundssteig “Gravettian” human remains are Holocene. Homo 51, 258–260.

37. Douka, K., Slon, V., Stringer, C., Potts, R., Hübner, A., Meyer, M., Spoor, F., Pääbo, S., and Higham, T. (2017). Direct radiocarbon dating and DNA analysis of the Darra-i-Kur (Afghanistan) human temporal bone. J. Hum. Evol. 107, 86–93.

38. Talamo, S., Hajdinjak, M., Mannino, M.A., Fasani, L., Welker, F., Martini, F., Romagnoli,
F., Zorzin, R., Meyer, M., and Hublin, J.-J. (2016). Direct radiocarbon dating and genetic analyses on the purported Neanderthal mandible from the Monti Lessini (Italy). Sci. Rep. 6, 29144.

39. Benazzi, S., Peresani, M., Talamo, S., Fu, Q., Mannino, M.A., Richards, M.P., and Hublin, J.-J. (2014). A reassessment of the presumed Neandertal remains from San Bernardino Cave, Italy. J. Hum. Evol. 66, 89-94.

40. Cersoy, S., Zazzo, A., Rofes, J., Tresset, A., Zirah, S., Gauthier, C., Kaltnecker, E., Thil, F., and Tisnerat-Laborde, N. (2017). Radiocarbon dating minute amounts of bone (3-60 mg) with ECHoMICADAS. Sci. Rep. 7, 7141.

41. Fewlass, H., Tuna, T., Fagault, Y., Hublin, J.-J., Kromer, B., Bard, E., and Talamo, S. (2019). Pretreatment and gaseous radiocarbon dating of 40-100 mg archaeological bone. Sci. Rep. 9, 5342.

42. Fewlass, H., Talamo, S., Tuna, T., Fagault, Y., Kromer, B., Hoffmann, H., Pangrazzi, C., Hublin, J.-J., and Bard, E. (2017). Size Matters: Radiocarbon Dates of <200 µg Ancient Collagen Samples with AixMICADAS and Its Gas Ion Source. Radiocarbon 60, 425-439.

43. Korlević, P., Talamo, S., and Meyer, M. (2018). A combined method for DNA analysis and radiocarbon dating from a single sample. Sci. Rep. 8, 4127.

44. Sponheimer, M., Ryder, C.M., Fewlass, H., Smith, E.K., Pestle, W.J., and Talamo, S. (2019). Saving Old Bones: a non-destructive method for bone collagen prescreening. Sci. Rep. 9, 13928.

45. Korlević, P., and Meyer, M. (2019). Pretreatment: Removing DNA Contamination from Ancient Bones and Teeth Using Sodium Hypochlorite and Phosphate. Methods Mol. Biol. 15-19.

46. Damgaard, P.B., Margaryan, A., Schroeder, H., Orlando, L., Willerslev, E., and Allentoft, M.E. (2015). Improving access to endogenous DNA in ancient bones and teeth. Sci. Rep. 5, 11184.

47. Boessenkool, S., Hanghej, K., Nistelberger, H.M., Der Sarkissian, C., Gondek, A.T., Orlando, L., Barrett, J.H., and Star, B. (2017). Combining bleach and mild predigitation improves ancient DNA recovery from bones. Mol. Ecol. Resour. 17, 742-751.

48. Schroeder, H., de Barros Damgaard, P., and Allentoft, M.E. (2019). Pretreatment: Improving Endogenous Ancient DNA Yields Using a Simple Enzymatic Predigitation Step. Methods Mol. Biol. 1963, 21-24.

49. Essel, E., Korlević, P., and Meyer, M. (2021). A method for the temperature-controlled extraction of DNA from ancient bones. Biotechniques 71, 382-386.

50. Essel, E., Zavala, E.I., Schulz-Kornas, E., Kozlikin, M.B., Fewlass, H., Vernot, B., Shunkov, M.V., Derevianko, A.P., Douka, K., Barnes, I., et al. (2023). Ancient human DNA recovered from a Palaeolithic pendant. Nature 618, 328-332.

51. van Klinken, G.J. (1999). Bone Collagen Quality Indicators for Palaeodietary and Radiocarbon Measurements. J. Archaeol. Sci. 26, 687-695.

52. Guiry, E.J., Szpak, P. (2020). Quality control for modern bone collagen stable carbon and nitrogen isotope measurements. Methods Ecol. Evol. 11, 1049-1060.

53. Schwarz, H.P., and Nahal, H. (2021). Theoretical and observed C/N ratios in

human bone collagen. J. Archaeol. Sci. 131, 105396.

54. DeNiro, M.J., and Weiner, S. (1988). Chemical, enzymatic and spectroscopic characterization of "collagen" and other organic fractions from prehistoric bones. Geochem. Cosmochim. Acta 52, 2197-2206.

55. Yizhaq, M., Mintz, G., Cohen, I., Khalaily, H., Weiner, S., and Boaretto, E. (2005). Quality Controlled Radiocarbon Dating of Bones and Charcoal from the Early Pre-Pottery Neolithic B (PPNB) of Motza (Israel). Radiocarbon 47, 193-206.

56. D'Elia, M., Gianfrate, G., Quarta, G., Giotta, L., Giancane, G., and Calcagnile, L. (2007). Evaluation of Possible Contamination Sources in the 14C Analysis of Bone Samples by FTIR Spectroscopy. Radiocarbon 49, 201-210.

57. Talamo, S., Fewlass, H., Maria, R., and Jaouen, K. (2021). "Here we go again": the inspection of collagen extraction protocols for 14C dating and palaeodietary analysis. Sci. Technol. Archaeol. Res. 7, 62-77.

58. Richards, M.P., Pettitt, P.B., Stiner, M.C., and Trinkaus, E. (2001). Stable isotope evidence for increasing dietary breadth in the European mid-Upper Paleolithic. Proc. Natl. Acad. Sci. USA 98, 6528-6532.

59. Pettitt, P.B., Richards, M., Maggi, R., and Formicola, V. (2003). The Gravettian burial known as the Prince ("Il Principe"): new evidence for his age and diet. Antiquity 77, 15-19.

60. Richards, M.P., and Trinkaus, E. (2009). Isotopic evidence for the diets of European Neanderthals and early modern humans. Proc. Natl. Acad. Sci. USA 106, 16034-16039.

61. Bocherens, H., Drucker, D.G., Germonpré, M., Lázničková-Galetová, M., Naito, Y.I., Wissing, C., Brúzek, J., and Oliva, M. (2015). Reconstruction of the Gravettian food-web at Predmostí I using multi-isotopic tracking (13C, 15N, 34S) of bone collagen. Quat. Int. 359-360, 211-228.

62. Wißing, C., Rougier, H., Baumann, C., Comeyne, A., Crevecoeur, I., Drucker, D.G., Gaudzinski-Windheuser, S., Germonpré, M., Gómez-Olivencia, A., Krause, J., et al. (2019). Stable isotopes reveal patterns of diet and mobility in the last Neandertals and first modern humans in Europe. Sci. Rep. 9, 4433. https://doi.org/10.1038/s41598-019-41033-3.

63. Church, D.M., Schneider, V.A., Graves, T., Auger, K., Cunningham, F., Bouk, N., Chen, H.-C., Agarwala, R., McLaren, W.M., Ritchie, G.R.S., et al. (2011). Modernizing reference genome assemblies. PLoS Biol. 9, e1001091.

64. Lindahl, T. (1993). Instability and decay of the primary structure of DNA. Nature 362, 709-715.

65. Briggs, A.W., Stenzel, U., Johnson, P.L.F., Green, R.E., Kelso, J., Prüfer, K., Meyer, M., Krause, J., Ronan, M.T., Lachmann, M., and Pääbo, S. (2007). Patterns of damage in genomic DNA sequences from a Neandertal. Proc. Natl. Acad. Sci. USA 104, 14616-14621.

66. Peyrégne, S., and Peter, B.M. (2020). AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination. Genome Biol. 21, 246.

67. Zavala, E.I., Aximu-Petri, A., Richter, J., Nickel, B., Vernot, B., and Meyer, M. (2022). Quantifying and reducing cross-contamination in single- and multiplex

hybridization capture of ancient DNA. Mol. Ecol. Resour. 22, 2196-2207. https://doi.org/10.1111/1755-0998.13607.

68. Furtwängler, A., Reiter, E., Neumann, G.U., Siebke, I., Steuri, N., Hafner, A., Lösch, S., Anthes, N., Schuenemann, V.J., and Krause, J. (2018). Ratio of mitochondrial to nuclear DNA affects contamination estimates in ancient DNA analysis. Sci. Rep. 8, 14075.

69. Andrews, R.M., Kubacka, I., Chinnery, P.F., Lightowlers, R.N., Turnbull, D.M., and Howell, N. (1999). Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet. 23, 147.

70. Secher, B., Fregel, R., Larruga, J.M., Cabrera, V.M., Endicott, P., Pestano, J.J., and González, A.M. (2014). The history of the North African mitochondrial DNA haplogroup U6 gene flow into the African, Eurasian and American continents. BMC Evol. Biol. 14, 109.

71. Pereira, L., Silva, N.M., Franco-Duarte, R., Fernandes, V., Pereira, J.B., Costa, M.D., Martins, H., Soares, P., Behar, D.M., Richards, M.B., and Macaulay, V. (2010). Population expansion in the North African late Pleistocene signaled by mitochondrial DNA haplogroup U6. BMC Evol. Biol. 10, 390.

72. Hervella, M., Svensson, E.M., Alberdi, A., Günther, T., Izagirre, N., Munters, A.R., Alonso, S., Ioana, M., Ridiche, F., Soficaru, A., et al. (2016). The mitogenome of a 35,000-year-old Homo sapiens from Europe supports a Palaeolithic back-migration to Africa. Sci. Rep. 6, 25501.

73. Svensson, E., Günther, T., Hoischen, A., Hervella, M., Munters, A.R., Ioana, M., Ridiche, F., Edlund, H., van Deuren, R.C., Soficaru, A., et al. (2021). Genome of Peştera Muierii skull shows high diversity and low mutational load in pre-glacial Europe. Curr. Biol. 31, 2973-2983.e9.

74. Haak, W., Lazaridis, I., Patterson, N., Rohland, N., Mallick, S., Llamas, B., Brandt, G., Nordenfelt, S., Harney, E., Stewardson, K., et al. (2015). Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522, 207-211.

75. Fu, Q., Hajdinjak, M., Moldovan, O.T., Constantin, S., Mallick, S., Skoglund, P., Patterson, N., Rohland, N., Lazaridis, I., Nickel, B., et al. (2015). An early modern human from Romania with a recent Neanderthal ancestor. Nature 524, 216-219.

76. Higham, T. (2011). European Middle and Upper Palaeolithic radiocarbon dates are often older than they look: problems with previous dates and some remedies. Antiquity 85, 235-249.

77. Kowalik, N., Anczkiewicz, R., Wilczyński, J., Wojtal, P., Müller, W., Bondioli, L., Nava, A., and Gasparik, M. (2020). Tracing human mobility in central Europe during the Upper Paleolithic using sub-seasonally resolved Sr isotope records in ornaments. Sci. Rep. 10, 10386.

78. Willis, K., and Vanandel, T. (2004). Trees or no trees? The environments of central and eastern Europe during the Last Glaciation. Quat. Sci. Rev. 23, 2369-2387.

79. Sommer, R.S., and Nadachowski, A. (2006). Glacial refugia of mammals in Europe: evidence from fossil records. Mamm Rev. 36, 251-265.

80. Prichystal, A. (2014). Lithic Raw Materials in Prehistoric Times of Eastern Central Europe (Masaryk University Press).
<figure><img src="image_15.png" /></figure>

81. Stefaniak, K., Pawłowska, K., Ratajczak, U., Roblícková, M., Gumírišk, W., and Wojtal, P. (2014). Middle and Late Pleistocene elks (Cervales Scott, 1855 and Alces Gray, 1821) from Poland: palaeoenvironmental and palaeogeographic implications. In *Annales Societatis Geologorum Poloniae* (Polskie Towarzystwo Geologiczne), pp. 341–362.

82. Kozłowski, J.K. (2013). Raw materials procurement in the Late Gravettian of the Carpathian Basin. In *The Lithic Raw Material Sources and the Interregional Human Contacts in the Northern Carpathian Regions* (Polska Akademia Umiejętności). 63–85.

83. Svoboda, J. (2000). The depositional context of the Early Upper Paleolithic human fossils from the Koněpryso (Zlatý kůn) and Mladeč Caves, Czech Republic. *J. Hum. Evol.* 38, 523–536.

84. Svoboda, J.A. (2006). The Structure of the Cave, Stratigraphy, and Depositional Context. In *Early Modern Humans at the Moravian Gate: The Mladeč Caves and their Remains*, M. Teschler-Nicola, ed. (Springer Vienna), pp. 27–40.

85. Oliva, M. (2006). The Upper Paleolithic Finds from the Mladeč Cave. In *Early Modern Humans at the Moravian Gate: The Mladeč Caves and their Remains*, M. Teschler-Nicola, ed. (Springer Vienna), pp. 41–74.

86. Vanhaeren, M., and d'Errico, F. (2006). Aurignacian ethno-linguistic geography of Europe revealed by personal ornaments. *J. Archaeol. Sci.* 33, 1105–1128.

87. White, R., and Normand, C. (2015). Early and archaic Aurignacian personal ornaments from Isturitz cave: Technological and regional perspectives. *Palethnologie* 7. https://doi.org/10.4000/palethnologie.789.

88. Villalba-Mouco, V., van de Loosdrecht, M.S., Rohrlach, A.B., Fewlass, H., Talamo, S., Yu, H., Aron, F., Lalueza-Fox, C., Cabello, L., Cantalejo Duarte, P., et al. (2023). A 23,000-year-old southern Iberian individual links human groups that lived in Western Europe before and after the Last Glacial Maximum. *Nat. Ecol. Evol.* 7, 597–609. https://doi.org/10.1038/s41559-023-01987-0.

89. Hajdinjak, M., Mafessoni, F., Skov, L., Vernot, B., Hübner, A., Fu, Q., Essel, E., Nagel, S., Nickel, B., Richter, J., et al. (2021). Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry. *Nature* 592, 253–257.

90. Seguin-Orlando, A., Korneliussen, T.S., Sikora, M., Malaspinas, A.-S., Manica, A., Moltke, I., Albrechtsen, A., Ko, A., Margaryan, A., Moiseyev, V., et al. (2014). Paleogenomics. *Genomic structure in Europeans dating back at least 36,200 years*. *Science* 346, 1113–1118.

91. Hodgkins, J., Orr, C.M., Gravel-Miguel, C., Riel-Salvatore, J., Miller, C.E., Bondioli, L., Nava, A., Lugli, F., Talamo, S., Hajdinjak, M., et al. (2021). An infant burial from Arma Veirana in northwestern Italy provides insights into funerary practices and female personhood in early Mesolithic Europe. *Sci. Rep.* 11, 23735.

92. Devièse, T., Massilani, D., Yi, S., Comeskey, D., Nagel, S., Nickel, B., Ribechini, E., Lee, J., Tseveendorj, D., Gunchinsuren, B., et al. (2019). Compound-specific radiocarbon dating and mitochondrial DNA analysis of the Pleistocene hominin from Salkhit Mongolia. *Nat. Commun.* 10, 274.

93. Devièse, T., Abrams, G., Hajdinjak, M., Pirson, S., De Groote, I., Di Modica, K.,

Toussaint, M., Fischer, V., Comeskey, D., Spindler, L., et al. (2021). Reevaluating the timing of Neanderthal disappearance in Northwest Europe. *Proc. Natl. Acad. Sci. USA* 118, e2022466118. https://doi.org/10.1073/pnas.2022466118.

94. Stuiver, M., and Polach, H.A. (1977). Discussion Reporting of 14C Data. *Radiocarbon* 19, 355–363.

95. Millard, A.R. (2014). Conventions for Reporting Radiocarbon Determinations. *Radiocarbon* 56, 555–559.

96. Katoh, K., and Standley, D.M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* 30, 772–780.

97. Bouckaert, R., Vaughan, T.G., Barido-Sottani, J., Duchéne, S., Fourment, M., Gavryushkina, A., Heled, J., Jones, G., Kühnert, D., De Maio, N., et al. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. *PLoS Comput. Biol.* 15, e1006650.

98. Leaché, A.D., Fujita, M.K., Minin, V.N., and Bouckaert, R.R. (2014). Species delimitation using genome-wide SNP data. *Syst. Biol.* 63, 534–542.

99. Rambaut, A., Drummond, A.J., Xie, D., Baele, G., and Souchard, M.A. (2018). Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. *Syst. Biol.* 67, 901–904.

100. Skov, L., Peyrégne, S., Popli, D., Iasi, L.N.M., Devièse, T., Slon, V., Zavala, E.I., Hajdinjak, M., Sümér, A.P., Grote, S., et al. (2022). Genetic insights into the social organization of Neanderthals. *Nature* 610, 519–525.

101. Tavaré, S., Balding, D.J., Griffiths, R.C., and Donnelly, P. (1997). Inferring coalescence times from DNA sequence data. *Genetics* 145, 505–518.

102. Fenner, J.N. (2005). Cross-cultural estimation of the human generation interval for use in genetics-based population divergence studies. *Am. J. Phys. Anthropol.* 128, 415–423.

103. Patterson, N., Price, A.L., and Reich, D. (2006). Population structure and eigenanalysis. *PLoS Genet.* 2, e190.

104. Renaud, G., Stenzel, U., and Kelso, J. (2014). *leeHom: adaptor trimming and merging for Illumina sequencing reads*. *Nucleic Acids Res.* 42, e141.

105. Wacker, L., Christl, M., and Synal, H.-A. (2010). Bats: A new tool for AMS data reduction. *Nucl. Instrum. Methods Phys. Res. B* 268, 976–979.

106. Gansauge, M.-T., Aximu-Petri, A., Nagel, S., and Meyer, M. (2020). Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA. *Nat. Protoc.* 15, 2279–2300.

107. Bronk-Ramsey, C., Higham, T., Bowles, A., and Hedges, R. (2004). Improvements to the Pretreatment of Bone at Oxford. *Radiocarbon* 46, 155–163.

108. Bard, E., Tuna, T., Fagault, Y., Bonvalot, L., Wacker, L., Fahri, S., and Synal, H.-A. (2015). AixMICADAS, the accelerator mass spectrometer dedicated to 14C recently installed in Aix-en-Provence, France. *Nucl. Instrum. Methods Phys. Res. B* 361, 80–86.

109. Wacker, L., Nèmec, M., and Bourquin, J. (2010). A revolutionary graphitisation system: Fully automated, compact and simple. *Nucl. Instrum. Methods Phys. Res. B* 268, 931–934.

110. Tuna, T., Fagault, Y., Bonvalot, L., Capano, M., and Bard, E. (2018). Development of small CO2 gas measurements with AixMICADAS. *Nucl. Instrum. Methods Phys. Res. B* 437, 93–97.

111. Ruff, M., Szidat, S., Gåggeler, H., Suter, M., Synal, H.-A., and Wacker, L. (2010). Gaseous radiocarbon measurements of small samples. *Nucl. Instrum. Methods Phys. Res. B* 268, 790–794.

112. Wacker, L., Fahrni, S.M., Hajdas, I., Molnar, M., Synal, H.-A., Szidat, S., and Zhang, Y.L. (2013). A versatile gas interface for routine radiocarbon analysis with a gas ion source. *Nucl. Instrum. Methods Phys. Res. B* 294, 315–319.

113. Rohland, N., Glocke, I., Aximu-Petri, A., and Meyer, M. (2018). Extraction of highly degraded DNA from ancient bones, teeth and sediments for high-throughput sequencing. *Nat. Protoc.* 13, 2447–2461.

114. Bronk-Ramsey, C. (2009). Bayesian Analysis of Radiocarbon Dates. *Radiocarbon* 51, 337–360.

115. Reimer, P.J., Austin, W.E.N., Bard, E., Bayliss, A., Blackwell, P.G., Bronk-Ramsey, C., Butzin, M., Cheng, H., Edwards, R.L., Friedrich, M., et al. (2020). The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0–55 cal kBP). *Radiocarbon* 62, 725–757.

116. Bronk-Ramsey, C. (2009). Dealing with Outliers and Offsets in Radiocarbon Dating. *Radiocarbon* 51, 1023–1045.

117. Weissensteiner, H., Pacher, D., Kloss-Brandstätter, A., Forer, L., Specht, G., Bandelt, H.-J., Kronenberg, F., Salas, A., and Schönherr, S. (2016). *HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing*. *Nucleic Acids Res.* 44, W58–W63.

118. Dür, A., Huber, N., and Parson, W. (2021). Fine-Tuning Phylogenetic Alignment and Haplogrouping of mtDNA Sequences. *Int. J. Mol. Sci.* 22, 5747. https://doi.org/10.3390/ijms22115747.

119. Ingman, M., Kaessmann, H., Pääbo, S., and Gyllensten, U. (2000). Mitochondrial genome variation and the origin of modern humans. *Nature* 408, 708–713.

120. Fu, Q., Li, H., Moorjani, P., Jay, F., Slepchenko, S.M., Bondarev, A.A., Johnson, P.L.F., Aximu-Petri, A., Prüfer, K., de Filippo, C., et al. (2014). Genome sequence of a 45,000-year-old modern human from western Siberia. *Nature* 514, 445–449.

121. Fu, Q., Mittnik, A., Johnson, P.L.F., Bos, K., Lari, M., Bollongino, R., Sun, C., Giemsch, L., Schmitz, R., Burger, J., et al. (2013). A revised timescale for human evolution based on ancient mitochondrial genomes. *Curr. Biol.* 23, 553–559.

122. Tamura, K., and Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. *Mol. Biol. Evol.* 10, 512–526.

123. Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., and Durbin, R.; 1000 Genome Project Data Processing Subgroup (2009). Genome Project Data Processing Subgroup (2009). The Sequence Alignment/Map format and SAMtools. *Bioinformatics* 25, 2078–2079.

124. Price, A.L., Patterson, N.J., Plenge, R.M., Weinblatt, M.E., Shadick, N.A., and Reich, D. (2006). Principal components analysis
corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909.

125. Patterson, N., Moorjani, P., Luo, Y., Mallick, S., Rohland, N., Zhan, Y., Genschoreck, T., Webster, T., and Reich, D. (2012). Ancient admixture in human history. Genetics 192, 1065–1093.

126. Petr, M., Vernot, B., and Kelso, J. (2019). admixr—R package for reproducible analyses using ADMIXTOOLS. Bioinformatics 35, 3194–3195.

127. Mallick, S., Li, H., Lipson, M., Mathieson, I., Gymrek, M., Racimo, F., Zhao, M., Chennagiri, N., Nordenfelt, S., Tandon, A., et al. (2016). The Simons Genome Diversity Project: 300 genomes from 142 diverse populations. Nature 538, 201–206.
## STAR★METHODS

### KEY RESOURCES TABLE

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td colspan="3">Biological samples</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C7/675</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C7/683</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C7/656</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C7/658</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C6/479</td></tr><tr><td>Osteological remain</td><td>This study</td><td>C6/494</td></tr><tr><td>Osteological remain</td><td>This study</td><td>B5/910</td></tr><tr><td>Osteological remain</td><td>This study</td><td>B5/913</td></tr><tr><td colspan="3">Chemicals, peptides, and recombinant proteins</td></tr><tr><td>0.5M EDTA pH 8.0</td><td>AppliChem</td><td>Cat.no. A4892</td></tr><tr><td>100% Tween 20</td><td>Merck</td><td>Cat.no. P5927-100ML</td></tr><tr><td>10 mg/mL Proteinase K</td><td>Sigma Aldrich</td><td>Cat.no. P6556</td></tr><tr><td>H20, HPLC grade</td><td>Merck</td><td>Cat.no. 270733</td></tr><tr><td>Guanidine hydrochloride</td><td>Sigma Aldrich</td><td>Cat.no. G3272</td></tr><tr><td>Isopropanol</td><td>Merck</td><td>Cat.no. 109634</td></tr><tr><td>100% Ethanol</td><td>Merck</td><td>Cat.no. 100983</td></tr><tr><td>PE Buffer</td><td>Qiagen</td><td>Cat.no. 19065</td></tr><tr><td>1 M Tris-HCl pH 8.0</td><td>AppliChem</td><td>Cat.no. A4577</td></tr><tr><td>3M Sodium Acetate, pH 5.2</td><td>Merck</td><td>Cat.no. S7899-500ML</td></tr><tr><td>5 M NaCl</td><td>Merck</td><td>Cat.no. S5150-1L</td></tr><tr><td>20% SDS</td><td>Thermo Fisher Scientific</td><td>Cat.no. AM9820</td></tr><tr><td>10x AmpliTaq Gold buffer (without Mg)</td><td>Life Technologies</td><td>Cat.no. 4379874</td></tr><tr><td>1M NaOH</td><td>Roth</td><td>Cat.no. 9062.3</td></tr><tr><td>20x SSC</td><td>Thermo Fisher Scientific</td><td>Cat.no. AM9763</td></tr><tr><td>Silica magnetic beads</td><td>G-Biosciences, VWR</td><td>Cat.no. 786-915</td></tr><tr><td>T4 RNA ligase reaction buffer including 50% (wt/vol) PEG 8000</td><td>NEB</td><td>Cat.no. B0216L</td></tr><tr><td>ATP solution (100 mM)</td><td>Thermo Fisher Scientific</td><td>Cat.no. R0441</td></tr><tr><td>T4 polynucleotide kinase 10 U uL<sup>-1</sup></td><td>Thermo Fisher Scientific</td><td>Cat.no. EK0031</td></tr><tr><td>T4 DNA ligase high concentrated 30 U uL<sup>-1</sup></td><td>Thermo Fisher Scientific</td><td>Cat.no. EL0013</td></tr><tr><td>Fast AP (1U/uL)</td><td>Thermo Fisher Scientific</td><td>Cat.no. EF0651</td></tr><tr><td>Klenow fragment (10 U/uL), including 10X reaction buffer</td><td>Thermo Fisher Scientific</td><td>Cat.no. EP0052</td></tr><tr><td>T4 DNA Ligase, 5 U/uL, including 10X reaction buffer and 50% (wt/vol) PEG-4000</td><td>Thermo Fisher Scientific</td><td>Cat.no.EL0012</td></tr><tr><td>Dynabeads MyOne Streptavidin C1 beads</td><td>Life Technologies</td><td>Cat.no. 65001</td></tr><tr><td>PEG-8000 powder</td><td>Promega</td><td>Cat.no. V3011</td></tr><tr><td>Sera-Mag SpeedBeads</td><td>Sigma Aldrich</td><td>Cat.no. GE6515-2105-050250</td></tr><tr><td>UltraPure DNase/RNase-Free Distilled Water</td><td>Invitrogen</td><td>Cat.no. 10977049</td></tr><tr><td>dNTPs 25mM each</td><td>Thermo Fisher Scientific</td><td>Cat.no. R1121</td></tr><tr><td>Herculase II Fusion DNA polymerase, including 5X Herculase II reaction buffer</td><td>Agilent</td><td>Cat.no. 600679</td></tr></tbody></table>

(Continued on next page)
Continued

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td>Human Cot-1 DNA</td><td>Thermo Fisher Scientific</td><td>Cat.no. 15279011</td></tr><tr><td>2x HI-RPM hybridization buffer, including Agilent blocking agent (dissolve in 1250uL water)</td><td>Agilent</td><td>Cat.no. 5188-5380_25 mL</td></tr><tr><td>Maxima Probe qPCR Master Mix (2x)</td><td>Thermo Fisher Scientific</td><td>Cat.no. K0261</td></tr><tr><td>AccuPrime Pfx DNA polymerase, including 10x AccuPrime reaction mix</td><td>Thermo Fisher Scientific</td><td>Cat.no. 12344-024</td></tr><tr><td>Maxima SYBR Green qPCR Master Mix (2X)</td><td>Thermo Fisher Scientific</td><td>Cat.no. K0253</td></tr><tr><td>Ultrapure H₂O (resistivity value: 18.2 MΩ·cm @ 25 °C; ≤ 5 ppb)</td><td>MilliporeSigma</td><td>MilliQ Element A10 Water Purification system</td></tr><tr><td>HCl 37%</td><td>Roth</td><td>Nr. 9277.1</td></tr><tr><td>HCl (6N)</td><td>Roth</td><td>Nr. 0281.1</td></tr><tr><td>NaOH (2N)</td><td>Roth</td><td>Nr. T135.1</td></tr><tr><td>Ezee Filters 9ml</td><td>Elkay Labs</td><td>Cat.no. 127-3193-000</td></tr><tr><td>Vivaspin Turbo 15 Ultrafilters PES (30 kDa MWCO)</td><td>Sartorius</td><td>Cat.no. VS15T01</td></tr><tr><td>AMS standard reference material Oxalic acid</td><td>NIST</td><td>SRM 4990C</td></tr><tr><td>IAEA-CH-6 (EA-IRMS standard)</td><td>IAEA</td><td>https://nucleus.iaea.org/sites/ReferenceMaterials/Pages/IAEA-CH-6.aspx</td></tr><tr><td>IAEA-CH-7 (EA-IRMS standard)</td><td>IAEA</td><td>https://nucleus.iaea.org/sites/ReferenceMaterials/Pages/IAEA-CH-7.aspx</td></tr><tr><td>IAEA-N-1 (EA-IRMS standard)</td><td>IAEA</td><td>https://nucleus.iaea.org/sites/ReferenceMaterials/Pages/IAEA-N-1.aspx</td></tr><tr><td>IAEA-N-2 (EA-IRMS standard)</td><td>IAEA</td><td>https://nucleus.iaea.org/sites/ReferenceMaterials/Pages/IAEA-N-2.aspx</td></tr><tr><td>Methionine</td><td>Elemental Microanalysis, Okehampton, UK</td><td>B2100</td></tr></tbody></table>

Critical commercial assays

<table><tr><td>Custom capture probes</td><td>Agilent</td><td>N/A</td></tr><tr><td>P5 8bp primer plate (384)</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>P7 8bp primer plate (384)</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>MinElute PCR Purification kit</td><td>Qiagen</td><td>Cat.no. 28006</td></tr><tr><td>QIAquick Nucleotide Removal kit</td><td>Qiagen</td><td>Cat.no. 28306</td></tr></table>

Deposited data

<table><tr><td>European Nucleotide Archive (nuclear DNA data)</td><td>This study</td><td>PRJEB66365</td></tr><tr><td>Dryad (mtDNA genome)</td><td>This study</td><td>https://doi.org/10.5061/dryad.47d7wm3m4</td></tr></table>

Oligonucleotides

<table><tr><td>(ddC, dideoxycytidine; TEG, triethylene glycol spacer; *, phosphothioate linkage; [N], 2'-O-methyl-RNA; {N}, locked nucleic acid (LNA); SpacerC12, 12-carbon spacer; FAM, fluorescein amidite; BHQ1, Black Hole Quencher-1)</td><td></td><td></td></tr><tr><td>TL 181 100uM (desalted), 1st adapter Phosphate-AGATCGGAAGAAA[A][A][A][A][A][A][A]-TEG-Biotin</td><td>IDT</td><td>N/A</td></tr><tr><td>TL 159 100uM (desalted), splinter SpacerC12-[A][A][A]CTTCCGATCTNNNNNNNN[A]-AminoC6</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>CL 128 100uM (HPLC), extension primer GTGACTGGAGTTCAGACGTGTGCTCTTCC*G*A*T*C*T</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>CL 53 500uM (HPLC), 2nd adapter strand 1 CGACGCTCTTC-ddC</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr></table>

(Continued on next page)
Continued

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td>TL 178 500uM (desalted), 2<sup>nd</sup> adapter strand 2<br>Phosphate-GGAAGAGCGTCGTGTAGGGAAAGAGTGTA</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>CL 304 100uM (HPLC), control DNA<br>Phosphate-ATTCAGCTCCGGTTCCCAACGATCAAGGC<br>GAGTTACATGA-Phosphate</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS5 100uM (desalted), forward primer<br>AATGATACGGCGACCACCGA</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS5 biotinylated 100uM (desalted), forward primer<br>Biotin-AATGATACGGCGACCACCGA</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS6 100uM (desalted), reverse primer<br>CAAGCAGAAGACGGCATACGA</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>CL105 100uM (desalted), forward primer<br>ACACTCTTTCCCTACACGACGCTCTTCCTCGTCGTTT<br>GGTATGGCTTC</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>CL106 100uM (desalted), reverse primer<br>GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTTC<br>ATGTAACTCGCCTTGATCGT</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS7 100uM (HPLC), forward primer A&B<br>ACACTCTTTCCCTACACGAC</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS8 100uM (HPLC), reverse primer A<br>GTGACTGGAGTTCAGACGTGT</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>IS10 100uM (reverse phase - HPLC), probe A<br>FAM-A(G)A(T)C(G)GAAGAGC(A)CAC-BHQ1</td><td>Eurogentec</td><td>N/A</td></tr><tr><td>CL107 100uM (HPLC), reverse primer B<br>TCATGTAACTCGCCTTGATCGT</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr><tr><td>CL118 100uM (HPLC), probe B<br>FAM-TTCAGCTCCGGTTCCCAACGAT-BHQ1</td><td>Sigma-Adrich (Merck)</td><td>N/A</td></tr></tbody></table>

Software and algorithms

<table><tr><td>leeHom</td><td>Renaud et al.<sup>94</sup></td><td>https://bioinf.eva.mpg.de/leehom/</td></tr><tr><td>BWA version 0.5</td><td>N/A</td><td>https://github.com/mpieva/network-aware-bwa</td></tr><tr><td>SAMtools version 1.3.1</td><td>Li, H. et al.<sup>95</sup></td><td>http://www.htslib.org/doc/1.3.1/samtools.html</td></tr><tr><td>Bam-rmdup version 0.6.1</td><td>NA</td><td>https://github.com/mpieva/biohazard-tools</td></tr><tr><td>MAFFT v.7</td><td>Katoh, K. et al.<sup>96</sup></td><td>https://mafft.cbrc.jp/alignment/server/index.html</td></tr><tr><td>BEAST V2.6.6</td><td>Bouckaert, R. et al.<sup>97</sup></td><td>http://www.beast2.org/</td></tr><tr><td>MODEL SELECTION</td><td>Leaché, A.D. et al.<sup>98</sup></td><td>https://beast community/workshop_model_selection</td></tr><tr><td>Tracer v1.7</td><td>Rambaut, A. et al.<sup>99</sup></td><td>https://beast community/tracer</td></tr><tr><td>Bam-caller v0.1</td><td></td><td>https://github.com/bodkan/bam-caller</td></tr><tr><td>SmartPCA</td><td>Patterson, N. et al., 2006 and<br>Price, A.L. et al., 2006<sup>100,101</sup></td><td>https://reich.hms.harvard.edu/software</td></tr><tr><td>ADMIXTOOLS</td><td>Patterson, N. et al.<sup>102</sup></td><td>https://reich.hms.harvard.edu/software</td></tr><tr><td>Admixr v0.7.1</td><td>Petr, M. et al.<sup>103</sup></td><td>https://github.com/bodkan/admixr</td></tr><tr><td>OxCal v4.4</td><td>Bronk-Ramsey,<sup>104</sup></td><td>https://c14.arch.ox.ac.uk/oxcal.html</td></tr><tr><td>BATS v3.64</td><td>Wacker et al.<sup>105</sup></td><td></td></tr><tr><td>AGE3 software v4.666</td><td>IonPlus AG</td><td>www.ionplus.ch</td></tr><tr><td>GIS software v4.15</td><td>IonPlus AG</td><td>www.ionplus.ch</td></tr></table>

(Continued on next page)
<table><thead><tr><th colspan="3">Continued</th></tr><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td>Other</td><td></td><td></td></tr><tr><td>Allen ancient DNA resource (aadr) for previously published ancient and modern DNA data</td><td>N/A</td><td>https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data</td></tr><tr><td>IntCal20 Calibration curve</td><td>Reimer et al.<sup>106</sup></td><td></td></tr></tbody></table>

# RESOURCE AVAILABILITY

## Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Helen Fewlass (helen.fewlass@crick.ac.uk).

## Materials availability

This study did not generate new unique reagents.

## Data and code availability

* Nuclear DNA data have been deposited at the ENA (European Nucleotide Archive) repository: PRJEB66365 and mitochondrial DNA data at the Dryad repository: https://doi.org/10.5061/dryad.47d7wm3m4. Both are publicly available as of the date of publication. The accession number and DOI are also listed in the key resources table. This paper also analyses existing, publicly available data from the Allen Ancient DNA Resource. The link to the dataset is listed in the key resources table. The radiocarbon data is available in this paper's supplemental information.

* All original code is available in this paper's supplemental information.

* Any additional information required to reanalyse the data reported in this paper is available from the lead contact upon request.

# EXPERIMENTAL MODEL AND SUBJECT DETAILS

## Archaeological information

Three Elk and three Bos/Bison incisor tooth pendants were selected for dating from squares B5, C6 and C7 to represent the spread of the 112 pendants across the trench, each originating from a different animal based on tooth morphology. Five of the six teeth selected had pre-existing broken roots with only part of the drilled hole preserved. The two human teeth C7/675 (uldm1) and C7/683 (fragment of a cusp of a deciduous molar) were selected for analysis as they had the largest amount of dentine preserved and were microCT-scanned prior to the study to preserve their 3D morphology. Following excavation (2009-2010), the material was washed and wet-sieved in a local stream and handled without gloves prior to laboratory analysis, as per standard archaeological procedure. No consolidants were applied to the surface of the human or herbivore teeth.

# METHOD DETAILS

## Radiocarbon dating

### Pretreatment for radiocarbon dating

Small samples of the dentine (human: 23.4 - 41.1 mg; herbivore: 34.4 - 67.4 mg) were removed using a dentist's rotary drill. Dentine samples were pretreated in the Human Evolution department at the Max Planck Institute for Evolutionary Anthropology, Leipzig, using a previously published protocol for small samples.⁴¹ Briefly, sample chunks were demineralised in HCl 0.5M at 4°C until mechanically soft and CO₂ effervescence stopped. Samples were treated with NaOH 0.1M to remove humic acid contamination and re-acidified in HCl 0.5M. The samples were gelatinised in HCl pH3 for 3-10 h at 70°C. The gelatin extracts were then filtered (Ezee-filters, Elkay labs, UK) to remove >80 µm particles and ultrafiltered to concentrate the >30 kDa weight fraction. Filters were pre-cleaned according to Bronk-Ramsey et al.¹⁰⁷ Samples were rinsed to a neutral pH between each step. The extracts were freeze-dried for >48 hours and weighed immediately to determine the collagen yield. Small aliquots of bones dating beyond the ¹⁴C method (background) were pretreated and measured alongside the samples.

### Collagen quality assessment

To assess the quality of the extracts, ~0.4 mg collagen was weighed into a tin capsule and measured on a ThermoFinnigan Flash 2000 elemental analyser (EA) coupled to a Thermo Delta plus XP isotope ratio mass spectrometer (IRMS) via a Conflo III interface (Thermo Fisher Scientific, Germany) to determine the elemental (C%, N%, C/N) and stable isotopic values (δ¹³C, δ¹⁵N). Stable carbon isotope ratios were expressed relative to Vienna PeeDee Belemnite (VPDB) and stable nitrogen isotope ratios were measured relative to atmospheric N₂ (AIR), using the delta notation (δ) in parts per thousand (‰). Stable isotope delta values were two-point scale normalised using international
reference material IAEA-CH-6 (sucrose, δ¹³C = -10.449 ± 0.033 ‰), IAEA-CH-7 (polyethylene, δ¹³C = -32.151 ± 0.050 ‰), IAEA-N-1 (ammonium sulphate, δ¹⁵N = 0.4 ± 0.2 ‰) and IAEA-N-2 (ammonium sulphate, δ¹⁵N = 20.3 ± 0.2 ‰). Repeated analysis of internal and international standards indicates an analytical error of ± 0.2‰ (1σ). Two in-house quality control standards were used to quality check the scale normalisation and evaluate analytical precision: 1) EVA-0012 methionine (Elemental Microanalysis, Okehampton, UK), n = 12, δ¹³C = -28.0 ± 0.14 ‰ (1 s.d.), δ¹⁵N = -6.5 ± 0.12 ‰ (1 s.d.); and EVA MRG pig gelatin, n = 11, δ¹³C = -19.6 ± 0.24‰ (1 s.d.) and δ¹⁵N = 4.9 ± 0.11 ‰ (1 s.d.). These values compare well with the long-term average values of δ¹³C = -28.0 ± 0.1‰ (1 s.d.) and δ¹⁵N = -6.4 ± 0.1‰ (1 s.d.) for EVA-0012 and δ¹³C = -19.7 ± 0.3‰ (1 s.d.) and δ¹⁵N = 5.0 ± 0.1‰ (1 s.d.) for EVA MRG.

Sufficient collagen was available from C7/675, C7/658, C6/494 for additional quality checks using Fourier Transform Infrared Spectroscopy (FTIR).⁵⁴⁻⁵⁶ Roughly 0.3 mg collagen was homogenised and mixed with ~40 mg of IR grade KBr powder in an agate mortar and pestle. The powder was pressed into a pellet using a manual hydraulic press (Wasserman) and analysed with an Agilent Technologies Cary FTIR with a DTGS detector. Spectra were recorded in transmission mode at 4 cm⁻¹ resolution with averaging of 34 scans between 4000 and 400 cm⁻¹ using Resolution Pro software (Agilent Technologies) and were compared to library spectra of well preserved collagen.

## AMS dating

The extracts were dated on the AixMICADAS¹⁰⁸ installed at CEREGE (Centre de Recherche et d'Enseignement de Geosciences de l'Environnement) in Aix-en-Provence, France. Each collagen extract was split and measured in several aliquots as per the protocol for small collagen extracts described in Fewlass et al.⁴¹,⁴² (Table S1). Where sufficient material was available, collagen (combustion weight ~700 ug C) was weighed into aluminium cups and graphitized using the AGE 3¹⁰⁹ (Automated Graphitisation Equipment) prior to dating on the AixMICADAS. Oxalic acid II standards and background collagen samples (Table S2) were measured in the same session and used in the age calculation of the archaeological samples in BATS.¹⁰⁵ A relative error of 30% was applied to the blank value and an additional uncertainty of 1.6‰ was propagated in the F¹⁴C error calculation as per standard practice.

Gas measurements were performed using the protocol described in Tuna et al.¹¹⁰ and Fewlass et al.⁴¹ Small aliquots of collagen (combustion weight ~60 ug C) were introduced into cleaned silver cups (800°C, 2 h) and combusted in an Elementar Vario MICRO cube EA (Elementar Analysensysteme GmbH, Germany) which was directly coupled to the gas ion source (GIS) of the AixMICADAS.¹¹¹,¹¹² Oxalic acid II NIST standards (from a gas canister) were measured to normalize and correct samples for fractionation. The age of each sample was corrected with background measurements from the same batch (Table S2). A relative error of 30% was applied to the blank value and an additional uncertainty of 4‰ was propagated in the F¹⁴C error.

## Ancient DNA analysis

### Data generation and shotgun sequencing

Sampling of the teeth for radiocarbon dating and genetic analysis was performed in a dedicated clean room at the Max Planck Institute for Evolutionary Anthropology. A microsampling approach was used for sample collection from the two of the infant teeth (C7/675 and C7/683) for downstream genetic analyses. Between 7.2 and 12.9 mg of dentine was collected from the two roots of C7/675 and two microsamples of 4.7 and 5.2 mg of dentine were drilled from the crown of C7/683 for a total of four samples. To extract DNA from the tooth powder, 500 uL of lysis buffer (for 5 mL lysis buffer: 4.5 mL 0.5M EDTA pH 8.0, 2.5 uL of Tween 20, 125 uL of 10mg/mL proteinase K, 373 uL water) was added to each tooth powder aliquot for overnight lysis at 37°C. The resulting lysate was purified following the automated extraction protocol described,¹¹³ using binder buffer 'D'. Each DNA extract was then converted into dual-indexed single-stranded libraries¹⁰⁶ and sequenced on Illumina MiSeq or HiSeq2500 platforms in pools of 63-96 libraries that included libraries from other projects not discussed in this paper. Negative controls were included starting at both the DNA extraction and library preparation steps. Base calling was performed using Bustard (Illumina) and *leeHom*¹⁰⁴ was used to trim adapters and merge overlapping forward and reverse sequences. The resulting sequences were mapped to the human reference genome (hg19⁶³) using BWA version 0.5 (https://github.com/mpieva/network-aware-bwa) with ancient DNA parameters (-n 0.001 -o 2 -l 16500). Reads that were not a perfect match to expected index combinations, shorter than 35 base pairs, or had a mapping quality less than 25 were removed using SAMtools version 1.3.1¹⁰⁴. *Bam-rndup* version 0.6.1 (https://github.com/mpieva/biohazard-tools) was used to remove PCR duplicates. Each library (including controls) was evaluated for the presence of ancient human DNA based on the presence of elevated C-to-T substitutions on the terminal ends of DNA fragments (Table S4). Potential cross-contamination due to index swapping between libraries was evaluated using a previously published method⁶⁷ that evaluates all sequenced index pairs in a sequencing run to estimate the number of reads assigned to a library that may originate from another library. All libraries were estimated to have less than 6 cross-contaminated reads. The negative controls were also examined for evidence of ancient DNA based on elevated C-to-T substitutions and no evidence of aDNA was identified in these libraries (Table S4).

## QUANTIFICATION AND STATISTICAL ANALYSIS

### Radiocarbon data analysis and reporting

The radiocarbon dates were calibrated and analysed using OxCal 4.4¹¹⁴ and the IntCal20 calibration curve.¹¹⁵ The dates were input into a one phase General t-type Outlier Model with a 5% prior likelihood of being an outlier, except C7/675 which was given an outlier likelihood of 100% due to the C/N value and FTIR analysis¹¹⁶ (Table S3). As per convention (e.g. Stuiver & Polach,⁹⁴ Millard⁹⁵), throughout the text uncalibrated ¹⁴C ages are reported with the abbreviation 'BP', meaning 'radiocarbon years before AD 1950' and are reported with 1σ errors, whereas
calibrated ranges are denoted as 'cal BP' through the text and are given at the 95.4% (2σ) probability range. The agreement of the replicate dates from a single extract was tested using a X² test using the R_Combine function. All previously reported dates have been re-calibrated with the IntCal20 curve.¹¹⁵

## Mitochondrial DNA analysis

Six libraries (three from the samples and three negative controls) were enriched for human mitochondrial DNA via hybridization capture on an automated platform as described in⁶⁷ and sequenced on Illumina's MiSeq platform in pools of 32 to 76 libraries. Analysis of the resulting sequences was performed as described for the shotgun sequencing with the exception that libraries were mapped to the rCRS (Table S5). The reconstruction of mitochondrial genomes was performed on one library from each C7/675 and C7/683 where a coverage of at least 5-fold and support of at least 80% was required for calling a consensus base using only putatively deaminated reads.

No evidence of cross-contamination from index swapping between libraries (< 1 contaminating read per library) was identified in the sequencing runs. Reads with elevated C-to-T substitutions were also absent from the negative controls (Table S5). While the library negative controls had relatively few unique mtDNA reads (5 to 36) there were 680 unique mtDNA reads in the extraction negative control. Support for the variants observed in C7/675 and C7/683 were examined in the extraction negative control. The only shared variant was a G at position 73, which is a common variant across human populations,¹¹⁷,¹¹⁸ supporting that laboratory contamination did not impact downstream analyses.

The resulting mtDNA genome from C7/675 was aligned to the rCRS with mtDNA genomes from one Neanderthal (Vindija33.19), 54 present-day¹¹⁹ and 18 ancient humans (Table S6) with MAFFT⁹⁶ for tree building and molecular dating using the program BEAST V2.6.6⁹⁷. Present-day individuals were set to 0 and the Borsuka individual was constrained to an age range of 0 to 60,000. A model using a strict clock and constant population was selected after using the path sampling approach with 40 path steps from the MODEL SELECTION⁹⁸ package in BEAST2 to evaluate different clock and tree models (Table S7). These models were selected because there was no significant difference identified between the models and therefore the simplest were selected. In this testing a chain length of 25 million iterations, 0.3 alpha parameter, pre-burn of 75,000 iterations, and an 80% whole chain burn-in were used. For both testing and the final analysis, clock models were set with a normal distribution around a mean mutation rate of 2.67E-08 per base per year¹²⁰ with a sigma of 1.00E-10.¹²¹ The Tamura-Nei 193 (TN93)¹²² substitution model was used for all models. For the final tree building and molecular dating, 30 million iterations (1 million pre-burn iterations) of three Markov chain Monte Carlo runs with sampling every 2,000 trees were combined using Logcombiner2 from BEAST2. The tree annotator program from BEAST2 was used to identify a single tree from the combined tree file, which was visualised with Figtree v1.4.4. The Tracer v1.7¹¹⁹ program was used to examine the resulting molecular dates.

## Estimating TMRCA for two mtDNA genomes with no differences

In order to estimate the potential time span between two mtDNA genomes that share an mtDNA haplotype, we followed a previously published approach¹⁰⁰ to estimate the time to the most recent common ancestor (TMRCA) in this scenario. Following Skov et al., 2022,¹⁰⁰ we used a probability density (see below), which combines equations 7, 8, and 9 from Tavare et al., 1996¹⁰¹ to determine the mean coalescence time for two sequences with 0 differences. To calculate Θ we used a mutation rate (μ) of 2.67e-8 mutations per base pair per year, a generation time of 29 years,¹⁰² length of 15,569, and set the effective population size to 1,000 individuals (Θ=2*N*μ*L). As the resulting mean time (Figure S3) represents one branch length, the total branch length is then 2 * 1,206 years, or 2,412 years.

$$ \text{Probability density function: } f_{\text{coalescence time}}(t|k) = \frac{(1+\Theta)^{1+k}}{k!} t^{k} e^{-(1+\Theta)t} $$

## Nuclear DNA analysis

Twelve libraries (six samples and six negative controls) were also enriched for ~1.2 million SNPs ('1240k array') previously described in⁷⁴ and⁷⁵ via hybridization capture. They were sequenced on Illumina's MiSeq and NextSeq platforms and analysed as described in the shotgun sequencing section (Table S8). For these libraries less than 0.001% of total reads (<15 reads) were estimated as originating from index swapping. The data from the two libraries from C7/675 identified as containing ancient DNA were merged for downstream analysis using samtools.¹²³ Pseudo-haplotype genotyping was performed using random read sampling with bam-caller (https://github.com/bodkan/bam-caller; v0.1) after trimming the three termini bases of each fragment to minimise the effects of deamination on downstream analyses. These data were converted into EIGENSTRAT format and merged with previously published genotypes of ancient and modern humans previously published and compiled (https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data). The principal component analysis was performed using smartpca¹⁰³,¹²⁴ on the 597,573 SNPs in the Affymetrix Human Origins array¹²⁵ on 1,395 modern-day humans from West Eurasia, the Americas, Central Asia and Siberia, and East Asia. The Borsuka individual was then projected onto the PCA with 22 ancient individuals. The genetic relationship between the Borsuka individual and modern and ancient individuals was explored further with SNPs from the '1240k' array using ADMIXTOOLS¹²⁵ with the R package admixr v0.7.1.¹²⁶ The modern populations used were from Simons Genome Diversity Project.¹²⁷

The high amount of modern human contamination in the genetic data recovered from C7/675 required all downstream analysis to be restricted to putatively deaminated fragments only. This limited the data to only 46,286 SNPs in the '1240k' panel. As this only represents 3.8% of the SNPs within the panel, this raises concern as to what conclusions are possible to draw with this limited amount of data. Therefore, comparisons between other ancient humans and the Borsuka individual were restricted only to individuals with at least 500,000 SNPs with one
exception. The Vêstonice complex includes multiple individuals with varying qualities of genetic data. Four of these individuals were included for analysis to serve as a proxy for Borsuka for the input of limited data on each analysis.

Basal Eurasian ancestry evaluation

D-statistics was used to explore how the Borsuka individual's genetic affinities compared to other previously published ancient humans. It has previously been shown that the D statistic is ~ 0 for both D(Kostenki14, X; Han/Dai, Mbuti) and D(W, Stuttgart; Han, Mbuti), when W/X is a European from ~39,000 to 14,000 years ago without Basal Eurasian ancestry[9]. Calculating these D-statistics with Borsuka in the X or W position is consistent with an individual that also does not contain basal Eurasian ancestry (Tables S9 and S10).