# A new late Neanderthal from Crimea reveals long-distance connections across Eurasia

Emily M Pigott, Konstantina Cheshmedzhieva, Elke Zeller, Laura G van der Sluis, Manasij Pal Chowdhury, Maddalena Giann, Emese Végh, Thorsten Uthmeier, Victor Chabai, Marylène Patou-Mat, et al.

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

## ► To cite this version:

Emily M Pigott, Konstantina Cheshmedzhieva, Elke Zeller, Laura G van der Sluis, Manasij Pal Chowdhury, et al.. A new late Neanderthal from Crimea reveals long-distance connections across Eurasia. Proceedings of the National Academy of Sciences of the United States of America, 2025, 122 (45), 10.1073/pnas.25189741221of10. halshs-05373503

HAL Id: halshs-05373503

https://shs.hal.science/halshs-05373503v1

Submitted on 19 Nov 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# A new late Neanderthal from Crimea reveals long-distance connections across Eurasia

Emily M. Pigott<sup>a,b,1,2</sup> ID, Konstantina Cheshmedzhieva<sup>a,b,1</sup> ID, Elke Zeller<sup>c</sup> ID, Laura G. van der Sluis<sup>a,b</sup>, Manasij Pal Chowdhury<sup>a,b</sup>, Maddalena Gianni<sup>a,b</sup>, Emese Végh<sup>a,b</sup>, Thorsten Uthmeier<sup>c</sup> ID, Victor Chabai<sup>d</sup> ID, Marylène Patou-Mathis<sup>e,f</sup>, Petra G. Šimková<sup>a,b</sup>, Jana N. Voglmayr<sup>a,b</sup>, Gerhard W. Weber<sup>a,b</sup> ID, Ron Pinhasi<sup>a,b</sup> ID, Axel Timmermann<sup>g,h</sup>, Martin Kuhlwilm<sup>a,b</sup>, Katerina Douka<sup>a,b</sup> ID, and Thomas Higham<sup>a,b,2</sup> ID

Affiliations are included on p. 8.

Edited by James O'Connell, The University of Utah, Salt Lake City, UT; received July 17, 2025; accepted September 22, 2025

The Crimean Peninsula contains several important Middle and Upper Paleolithic sites, including Starosele, Kabazi II, and Siuren I. The region has been considered a potential refugium for Neanderthals before their replacement by *Homo sapiens*. However, no genetic data have been obtained from any of these late Neanderthals, some being inaccessible or badly preserved. Starosele is a notable site which has undergone excavations in recent years. We used collagen peptide mass fingerprinting (Zooarchaeology by Mass Spectrometry, ZooMS), to screen for potential human remains among thousands of fragmented bones from the site. Of the 150 bone fragments we analyzed, 97.3% had sufficient collagen preservation for taxonomic identification. Our results suggest Paleolithic humans primarily hunted horses. One ~5 cm bone fragment yielded peptide mass fingerprints matching *Hominidae*. Radiocarbon dating revealed an age range of 46 to 45,000 y old, close to the transition from the disappearance of Neanderthals to the dispersal of *H. sapiens* in western Europe. We sequenced a twofold coverage mitochondrial genome from this bone, indicating the individual belongs to the Neanderthal lineage. The mitogenome clusters with other Neanderthal mitogenomes previously generated from the Russian Altai region. Alongside this, an analysis of the lithic corpus from both regions suggests that a wider Neanderthal dispersal, linked to the Micoquian stone tool industry, occurred after ~60,000 y ago. We assessed the paleoclimate connection (temperature and precipitation) between these locations and identified a high habitat suitability corridor along 55°N, suggesting that the long-distance movement of Neanderthals would have been facilitated by periods of favorable climate.

Paleolithic archaeology | Neanderthals | paleoproteomics | radiocarbon dating | ancient DNA

The Middle to Upper Paleolithic transition marks the biocultural transformation from a Neanderthal-dominated western Eurasia to one exclusively populated by *Homo sapiens*. Documenting the precise nature and timing of this transition requires a careful analysis of the excavated remains from the numerous archaeological sites in the wider region, along with the application of the latest biomolecular and chronometric methods. Precise radiocarbon dating (1) and ancient genomics (2–4) have made significant contributions to our understanding of the disappearance of Neanderthals. However, Hominin remains are extremely rare and limit the impact of these methods.

There are currently very few known remains of late Neanderthals and early modern humans across Eurasia. A significant proportion of the faunal remains excavated at Paleolithic sites are not identifiable using traditional zooarchaeological methods (5) This is mostly due to the activity of carnivores (6), anthropogenic fragmentation (7), post-depositional movement and erosion, and taphonomic influences (8). At some sites, the proportion of unidentifiable bones is more than 95% (e.g., Denisova Cave), resulting in substantial loss of archaeological information. Collagen peptide mass fingerprinting, also known as Zooarchaeology by Mass Spectrometry (ZooMS), enables the identification of unidentifiable bone fragments to species/genus level with a high rate of throughput (9, 10). This method has allowed archaeologists to explore the full spectrum of faunal remains from archaeological contexts and in some instances to identify rare hominin bone fragments (11–14).

We used ZooMS to analyze Paleolithic bone assemblages from the Crimean site of Starosele (Fig. 1). The aims were threefold; first, to widen our understanding of the subsistence strategies of the human occupants, second, to screen bone fragments to find potential hominin bone remains which could be used for direct dating, genetic sequencing, and paleodietary analysis, and finally to determine the history of hominin occupation on the Crimean Peninsula.

## Significance

The Crimean Peninsula contains key Middle to Upper Paleolithic transitional archaeological sites, including the site of Starosele, where we identified a new Neanderthal; Star 1. This study highlights the integration of Zooarchaeology by Mass Spectrometry (ZooMS), radiocarbon dating, and ancient DNA analysis to uncover rare hominin remains and enhance our understanding of late Neanderthal populations in this region. Genetically, Star 1 is closely related to Neanderthals from the Altai via its mitochondrial DNA, suggesting long-range migrations of Neanderthal groups across Eurasia. These migrations during favorable climatic conditions likely involved the spread of the Micoquian lithic tradition, indicating both cultural continuity and regional mobility during the Late Pleistocene.

The authors declare no competing interest.

This article is a PNAS Direct Submission.

Copyright © 2025 the Author(s). Published by PNAS. This open access article is distributed under Creative Commons Attribution License 4.0 (CC BY).

PNAS policy is to publish maps as provided by the authors.

¹E.M.P. and K.C. contributed equally to this work.

²To whom correspondence may be addressed. Email: emily.pigott@univie.ac.at or thomas.higham@univie.ac.at.

This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2518974122/-/DCSupplemental.

Published October 27, 2025.
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** The location of the sites of Starosele, Siuren I, Kabazi II, Kabazi V, and Chokurcha in the Crimean Peninsula. The red to brown colors represent the first and highest mountain range, which reaches elevations over 1,000 m, the lighter and less frequent red to orange color indicates the second middle mountain range; due to low elevations, the lowermost third range in the North of the middle range is almost not visible in this map. Map made with Mapbox and Topographic map.</figcaption></figure>

**The Archaeological Site of Starosele.** The Crimean Peninsula contains numerous well-preserved, stratified Paleolithic sites, with many spanning within the transitional period in terms of bioculture and hominin occupations of approximately 47,000 to 42,000 cal BP (15). Based on previous radiocarbon dates, the Crimean Peninsula has been described as a refugium for late surviving Neanderthals just before their disappearance. However, new research, including radiocarbon dating using compound-specific techniques, has shown that this may need to be modified because some of the previous dates determined are, in fact, much younger in age in comparison to the compound-specific results (16).

Starosele is a rock-shelter positioned within a steep canyon (17, 18) and comprises four distinct cultural layers (19). The stratigraphy and site plan are shown in *SI Appendix, section A*. Level 2 and 3 are separated by a rockfall event, with evidence suggesting that Level 4 was briefly occupied (17). Previously, human remains within Level 1 were assumed to be Medieval Muslim burials after many decades of debate and uncertainty (20).

The Crimean Micoquian lasted from the Last Interglacial until approximately 40 kyr cal BP (21, 22). This technocomplex is characterized by soft hammer bifacial thinning and retouching, to produce bifacial tools (23). Micoquian lithic assemblages are exclusively linked with Neanderthals across Eurasia, as evidenced by the lithic assemblages co-occurring with Neanderthal remains at various sites, such as Kiik Koba (24), Zaskalnaya VI (25), Stajnia Cave (26), Mezmaiskaya (27) Chagyrskaya and Okladnikov (28, 29). In Starosele, Levels 1, 2, and 4 are assigned to the Crimean Micoquian, while Level 3 remains unassigned to a previously known lithic industry (Starosele-Level.3-Industry). Bifacial tools are present at all levels, except Level 3, which exhibits blank production methods (30). This may indicate that another group of Neanderthals was present at the site at this time. However, further analysis is needed to characterize the lithic assemblage fully.

Zooarchaeological analyses have previously been conducted at Starosele. A significant percentage of the faunal spectrum could not be identified to species through morphological inspection alone. Of the identified fauna, *Equus hydruntinus* dominates the assemblages. Level 1 comprises ~75% Equidae (19, 31). Cervus sp. (deer), various carnivores, and several bird species have also been identified in all levels (31, 32). Many of the faunal remains exhibit anthropogenic cut marks, indicating that butchering activities occurred at the site. Level 1 has a large assemblage of faunal material, with extensive evidence of animal butchery. However, there is also a high degree of fragmentation among the identified remains, which include *Sus* sp., *Vulpes* sp., *Ursus* sp., *Crocuta* sp., *Rupicapra rupicapra*, Cervid, Elephantidae, *Bos/Bison*, and *Saiga* sp. (17).

Here we investigate the potential of the fragmented faunal assemblage for hominin bones and apply paleoproteomic analyses to unidentified fauna from Starosele. Following this, we apply DNA and radiocarbon analyses to a small number of bones.

## Results

**Paleoproteomics.** On the basis of morphological and preservation characteristics, i.e., unburnt, sizeable bones over 2 cm long, we applied ZooMS to 150 bone fragments from all layers of Starosele. For these, we obtained 146 ZooMS identifications (97.3% success rate) (SI Appendix, section C and Dataset S1). The vast majority of these were identified as animal bones. Overall, the ZooMS results broadly confirm what was known from formal zooarchaeological analyses (17, 32) with six taxa forming the bulk of the identifications and *Equus* sp. dominating (93%). However, within Level 3, Square F21, a higher percentage of Bovidae or Cervidae/Bovidae is observed in comparison to the other stratigraphic contexts. Two faunal species were found only using ZooMS; Rhinocerotidae and *Canis lupus*. The taxonomic identifications are consistent with a classic MOIS3 paleoenvironment (33), characterized by
fluctuating warmer and colder environments during stadial and interstadial episodes (34).

**Starosele Human Remains.** Among the ZooMS-analyzed bones, we identified a Hominidae bone fragment (Fig. 2 A and B), which originated from Level 1 (Micoquian), Square I23 (SI Appendix, section A), close to previously identified human burials (see below).

We refer to this new specimen as Star 1. See SI Appendix, section B for stable isotope information. The ZooMS spectra for Star 1 (coded Z00113) shows the six peptide markers (A-D & F and G) associated with this Family, including 1,235.71, 1,478.74, 1,580.86, 2,115.19, 2,869.56, 2,957.68 (Fig. 3C) (35, 36). Peptide marker E is missing, which may be due to collagen degradation and its scarce presence in humans (12). The human bone has no

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** (A) Photograph of the hominin bone, Star 1, from Level 1, Square I23 (Z00113) shortly after it was identified using biomolecular methods. (B) µCT scans of the hominin bone. (C) MALDI mass spectrometry data from Z00113, Star 1. Previously published peptide markers are labeled A-D and F and G. Peptide marker E is missing. (D) Summary of the ZooMS results from Starosele plotted by the stratigraphic context and percentage. The exact percentage for each taxon can be found in SI Appendix, section C.</figcaption></figure>
clear recognizable morphological characteristics or human-induced modifications, which is why it was overlooked during excavation and zooarchaeological analysis (Fig. 2A).

The fragment was microcomputed tomography (μCT) scanned and virtually analyzed prior to any further invasive examinations. It measures 49.8 mm in maximum length and 18.8 mm in maximum width. Two relatively flat surfaces meet at an approximately orthogonal angle (Fig. 2 A and B). The fragment exhibits minimal signs of fossilization and is composed primarily of cortical bone, with residual trabecular bone lining the entire inner surface. Cross-sectional measurements of the cortical bone

yielded a maximum thickness of 6.21 mm at the angled midsection. The thickness gradually decreases along both surfaces extending from the corner, with the thinnest part measuring 2.78 mm. According to Gosman et al. (2013) (37), the individual to whom this fragment belonged was likely not younger than an adolescent. The surface data can be found at the following link https://www.virtual-anthropology.com/3d-data/free-data/.

Based on the size, shape, and curvature of the fragment, as well as its ZooMS classification, we suggest it could only originate from a hominin femur or humerus. According to the comparative analysis we performed, the fragment most likely derives from the

<figure><img src="image_6.png" /><figcaption>**Fig. 3.** (A) Phylogenetic tree of the mitochondrial genome featuring Star 1, together with other ancient hominins and modern humans. The Starosele, Star 1 branch, is highlighted in orange, and the bootstrap values are marked in a color gradient from red (minimum) to green (maximum). For the actual bootstrap values, see *SI Appendix, section D*. (B) Number of pairwise differences between Star 1 and other ancient and modern humans, sorted from most to least dissimilar.</figcaption></figure>
anterior distal shaft of a femur, since it closely matches our comparable specimens in cross sectional angle, continuity of the surfaces, and cortical thickness, while it does not align with humeri with respect to these characteristics. However, it should be noted that these results are based on a small comparative sample.

Radiocarbon dating was undertaken on the bone using the standard ultrafiltration approach at the University of Vienna's Higham Lab. The result was 39,858 ± 736 BP (VIE-1203). The calibrated age ranges between 43,860 to 42,690 cal BP (at 68.3% prob.) and 44,390 to 42,450 cal BP (at 95.4% prob.), a result consistent with previous radiocarbon dates from Level 1 (SI Appendix, section A). We also treated the collagen from the bone sample following ultrafiltration with acid hydrolysis and purified the amino acids using XAD-2 resin chromatography. This resulted in a radiocarbon date of 43,212 ± 295 BP (VIE-1541). The calibrated age ranges between 45,910 to 45,340 cal BP (at 68.3% prob.) and 46,130 to 44,990 cal BP (at 95.4% prob.). We conclude that the XAD radiocarbon date is the most accurate of the two, as it has been shown in the past to remove environmental and museum-derived contaminants successfully and thereby yield more accurate results compared with bulk collagen determinations (38, 39).

The bone was processed for DNA extraction and sequencing and two libraries were prepared, using 40 mg (STS1.1) and 10 mg (STS1.2) of bone powder, respectively (40). One of the subsamples (STS1.2) was subjected to a bleach pretreatment (41). Two single-stranded DNA libraries were constructed using an optimal method for ancient DNA (42). Diagnostic positions for derived Neanderthal alleles compared to modern humans (43) indicated that the sample contains DNA matching the Neanderthal allele (diagnostic positions with Neanderthal reads: 10/22 for STS1.1 and 11/52 for STS1.2). Due to the extremely low endogenous DNA content in the constructed libraries (0.038% STS1.1; 0.043% STS1.2), we performed mitochondrial target enrichment, resulting in an average coverage of 2.28-fold for the mitochondrial genome.

This mitochondrial genome confirms that the bone is most consistent with belonging to the Neanderthal lineage, not modern human or Denisovan-like lineages (44–46). Although due to relatively low coverage the placement in a phylogenetic tree cannot be fully resolved, Star 1 appears basal to European Neanderthals (Goyet Q305-4/Feldhofer2), and derived to Siberian Neanderthals (Denisova 11, Okladnikov A, 2). When comparing the mtDNA sequence to complete Neanderthal mtDNA data (Fig. 3B), Star 1 appears most closely related to five hominins from three sites; Denisova 11, Chagyrskaya E, and Okladnikov 2, A and B (12, 47, 48). Notably, all three are in the Altai region of Russian Siberia despite being ~3,000 km distant.

## Discussion

The Neanderthal bone fragment identified in this work is a key finding in the light of previous research conducted at Starosele. The site was known for many decades for the discovery of the widely cited “Starosele child” in the 1950s, originally excavated by Formozov (32). The remains, estimated to be of an 18 to 19-mo-old child, have been controversial and strongly debated. They were initially identified as a “transitional Neanderthal” or *H. sapiens* individual by Formozov and others, who saw similarities with the remains found in the Levantine site of Skhul (32). Howell agreed, comparing the Starosele child to Skhul, Qafzeh, and modern humans in Russia (49–51). Later, in the 1990s, Soffer suggested that the child was an anatomically modern human (19). Confusion regarding biological identification was mirrored by uncertainty in the age of the remains and their context, owing to excavation methods that were not of a robust nature and a lack of

photographs or detailed field notes that could confirm the exact positioning of the remains (19). New excavations by a joint Ukrainian/American team between 1993 to 1995, however, showed that it was most likely to be a much later burial, perhaps late Medieval (19). Marks and Chabai found two other burials at the site in 1993 to 1994: one infant and one adult, in squares H25 (adult) and I22 (infant). They suggested that these were probably 17/18th century Muslim burials, based on the stratigraphy of the burial pits, layout of the bodies and their direction and the typology of the ceramic pieces from the burial pits. Hence, it was inferred that the Starosele child burial was probably from a similar period (SI Appendix, section A) (19). No dating was performed on the new human remains, however, and while it would be interesting to demonstrate their late age, it has not been possible to do this in the current study. The most likely explanation is that the burials are late and were interred into Middle Paleolithic cultural horizons (Level 1), where the Star 1 human bone originates.

Level 1 has been dated using various methods and protocols (SI Appendix, section A). The calibrated ages of the latest three radiocarbon dates using the ultrafiltration protocol from Level 1 range between 47,710 to 42,490 cal BP (at 68.3% prob.) and 51,010 to 44,280 cal BP (at 95.4% prob.) (SI Appendix, section A). New compound specific protocols are currently being applied to date other bones.

Despite the extremely low endogenous DNA content and the lack of nuclear DNA, the 2.28-fold coverage after targeted capture shows that Star 1 did not belong to the divergent lineage of Hohlenstein Stadel (52). The sequence instead places it as being most similar to Denisova 11, Chagyrskaya E, and Okladnikov 2, A and B, which are all in the Altai region (Fig. 3). Because the genome data are incomplete, especially in parts of the D-loop and the hypervariable regions, we cannot confirm whether Star 1 can be assigned to this lineage of Siberian Neanderthals or represents a sister lineage. However, in the currently available record of Neanderthal diversity, the Siberian Neanderthals are the closest related individuals. Okladnikov and Chagyrskaya have been described as short-term, seasonal hunting camps, occupied by Neanderthals (47, 53). One suggestion, based on the close affinity of these Neanderthal mitochondrial genomes, is that they might represent an expansive movement of a Neanderthal group, most likely in an easterly direction, at some point during MIS4 ranging from ~71 to 57 ka (28).

To demonstrate further potential dispersal corridors and identify paleoclimatic and paleoenvironmental preferences, we constructed two Mahalanobis distance habitat suitability models covering the window between MIS6-MIS3, ~150 to 30,000 BP (Methods). This method has previously been used to explore habitat overlap between Neanderthals and Denisovans (54). Our habitat suitability models suggest that climatic conditions linking central Eurasia (Altai) and eastern Europe were most favorable during interglacial periods, supporting MIS 5e as a potentially attractive time for dispersal along this route [Fig. 4 and SI Appendix, section F, (55)]. We identified the most likely route between the Crimean Peninsula and Altai region along longitude 55°N (Fig. 4 A and E). This dating estimate appears to fit the evidence in the Altai for the presence of Neanderthals at Denisova Cave, where genetic and chronometric evidence places them between ~130,000 to 100,000 y ago. There are hints that Neanderthals might have been present prior to this, based on mtDNA in ~170,000 y ago sediment at Denisova Cave, but no human remains of this antiquity have been found yet (56). Later, Neanderthals appear at sites such as Okladnikov Cave. The initial radiocarbon dating of Neanderthal remains there disclosed a wide range of results, including ages younger than 40,000 cal BP. Okladnikov 2, 11, and 15 have been used to date carbon dated using robust single compound hydroxyproline dating to older than 42,000 cal BP (57, 58). At other Altai sites, such as Chagyrskaya,
<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Average habitat suitability of Star 1 and genetically related hominins from Croatia, Poland, the Crimean Peninsula, and Altai region (A). With time series of global mean average temperature from the 3 My transient CESM1.2 simulation (61, 62) (B), and results of the habitat suitability model over the area in Box 1, 2, 3, and 4 (C–F) where the blue line is the average suitability and the gray lines are the results of the individual models of the Monte Carlo analysis (Methods). In box 1,2, and 4 each sample included in the primary habitat suitability model was added with the entire time uncertainty (black bars) and at the top of the timeline plots the marine isotope stages (MIS) are indicated in black and white, representing warmer and colder periods, respectively.</figcaption></figure>

the most likely age of the Neanderthals is more precise and dates to between 59 and 51,000 cal BP (47). This suggests two, perhaps three, dispersals of Neanderthals from the west to the Altai at 170, 120-100, and ~60,000 cal BP. This hypothesis is supported by genetic data which suggest that the Chagyrskaya and Okladnikov Neanderthals were more closely related to the Vindija Cave 33.19 individual (59), than to the earlier Altai Neanderthal. Similarly, Denisova 11's Neanderthal mother is also genetically more closely related to European populations of Neanderthals, represented by Vindija 33.19, having 19.6% alleles present compared with 12.4% in Altai Neanderthals (60). It seems most parsimonious to conclude, therefore, that movement occurred from west to east and during times of favorable climate, particularly during MIS5e. Potential later movement, however, accounting for the Chagyrskaya and Okladnikov Neanderthals, seems more likely to have occurred during less benign climate phases, perhaps during MIS4.

In addition to the genetic links between Starosele, Chagyrskaya, and Okladnikov, these Neanderthal groups appear to have had

similar subsistence strategies, focused upon exploiting mainly horses, as well as bison and small ungulates, characteristic of steppe and piedmont environments. Sites dating to these periods have been interpreted as hunting and butchering sites focusing on large herbivores and ungulates within steppe environments (23, 63). Lithic assemblages excavated from Okladnikov, Chagyrskaya, and Starosele are also very similar, and characterized by artifacts of the Micoquian technocomplex (22, 28, 64, 65). A statistical study examining the lithic components of Micoquian assemblages from Chagyrskaya and comparing them to other Micoquian assemblages, showed that the Crimean Micoquian was significantly closer to Chagyrskaya than Denisova (64). The lithic assemblages from Chagyrskaya, in turn, compare closely to Okladnikov assemblages, being a unique variant in the Altai region, possibly representing small Neanderthal groups (63). The Neanderthals at Denisova Cave had a completely different lithic assemblage, indicating a different population of Neanderthals were present, and at a much earlier time (29, 66).

We may conclude, from the archaeological and genetic evidence, therefore, that it is likely that Neanderthals were in the Altai at two, perhaps even three, different periods. Our results show that Starosele was part of a wide-ranging dispersal network of Neanderthals, linked by genes, lithics and behavior.

## Methods

**Sample Selection.** All faunal remains were excavated between 1993 to 1994 and sampled by EP at the University of Erlangen-Nuremberg, Germany. The samples were collected for ZooMS analysis depending on the availability and layers present at the institution. Approximately 20 to 25 bones were selected from each layer/level within different squares, when possible. Some bags containing the faunal remains had a large number of bones (50 to 100), while others contained approximately 20. A further 50 bones were sampled from Level 1, square I23, after the hominin bone was found (Table 1). Star 1 originally weighed 6.66 g, measuring a maximum 52.3 mm in length and 16.5 mm in width. In total, 26 mg was taken for ZooMS analysis, with ~900 mg taken for radiocarbon dating. Stable isotope analysis was undertaken on collagen from the radiocarbon dating preparation.

**Sampling.** The bone samples were manually cleaned with air abrasion mixed with aluminum oxide powder at the Douka/Higham labs, Department of Evolutionary Anthropology at the University of Vienna, Austria. This enabled the removal of superficial contaminants before sampling, with around 30 mg of bone being clipped off with clean dental pliers for ZooMS analysis. For the human bone, 26 mg was clipped for ZooMS and ~900 mg sampled for radiocarbon dating.

**CT Scanning.** The Star 1 hominin bone was µCT scanned at the Vienna Micro-CT Lab with a VISCOM X8060 Scanner (39 µm, 110 kV, 400 µA, 1,400 ms, 0.75 mm copper filter), before further invasive sampling for dating and paleoproteomics were undertaken. For further comparative analysis a femur and humerus of one Sub-Saharan African and one European individual respectively were scanned at the same facility (63 µm, 130 kV, 320 µA, 1,400 ms, 0.75 mm copper filter). Additionally, scans of the Neanderthal 1 specimen from Feldhofer Grotto were provided from the Max Planck Institute for Evolutionary Anthropology Leipzig with a resolution of 91 µm. All data were virtually analyzed and measured in the Amira 3D 2024.1 software from Thermo Fisher.

**ZooMS.** The ZooMS analysis was conducted at the Douka Lab, Department of Evolutionary Anthropology, University of Vienna, Austria. Bone samples were analyzed using an acid insoluble protocol (9, 67). Following wet chemistry, the spectra were analyzed using flexAnalysis 3.4 (Brucker Daltonics) and mMass software (68). Subsequently, the spectra were compared against known peptide markers in a reference library (9). Some species/family faunal remains had already been identified from previous zoological analyses.

ZooMS analysis was undertaken using the acid insoluble collagen protocol (9). The bone specimens were first sampled using pliers, chipping off 15 to 25 mg of bone. The samples were demineralized up to 48 h with 500 µL of cold (4 °C) 0.6 M HCl (covering the sample), until somewhat spongy. After demineralization, samples were vortexed for approximately 10 s and centrifuged for 1 min at 12 RPM. The samples then went
<table><caption>**Table 1.** Number of bones analyzed for ZooMS analysis from five different square numbers from Levels 1 and 3</caption><thead><tr><th>Site</th><th>Typology</th><th>Level</th><th>Square</th><th>No. of bones</th></tr></thead><tbody><tr><td>Starosele</td><td>Starosele-Level-3-Industry</td><td>3</td><td>F21</td><td>20</td></tr><tr><td>Starosele</td><td>Micoquian</td><td>1</td><td>G21</td><td>20</td></tr><tr><td>Starosele</td><td>Micoquian</td><td>1</td><td>H23</td><td>20</td></tr><tr><td>Starosele</td><td>Micoquian</td><td>1</td><td>I23</td><td>70</td></tr><tr><td>Starosele</td><td>Micoquian</td><td>1</td><td>I22</td><td>20</td></tr></tbody></table>

through 2 to 3 washes with 200 µL of 50 mM AmBic, followed by centrifuging and discarding the AmBic solution. Following this, 100 µL of 50 mM AmBic (ammonium bicarbonate buffer, pH 8.0) was inserted into each sample and then incubated for 1 h at 65 °C for gelatinization. Half (50 µL) of the sample was removed in case a backup was needed in the future. Trypsin digestion was undertaken by inserting 1 µL of trypsin solution 0.4 µg (Thermo Scientific Pierce™ Trypsin Protease), which was added to each supernatant and digested for 18 h at 37 °C. The digested samples were treated with 1 µL of 5% Trifluoroacetic Acid (TFA) to stop further enzymatic digestion. ZipTipping was performed for purification and peptide extraction, using C18 ZipTips (Thermo Scientific Pierce™ C18 Tips), using 50 mL of 50% ACN in 0.5% TFA. A calibrant at 1:2 ratio of 0.5 µL was added to the plate, with 0.5 µL matrix solution. The samples (0.5 µL) were spotted with 0.5 µL of matrix solution. Thereafter, the processed samples were analyzed using a Bruker Autoflex Speed Matrix-Assisted Laser Desorption/Ionization Time-of-Flight (LRF MALDI ToF/ToF) mass spectrometer (69). The peptide marker (A-G) was manually inspected in mMass 5.5.0 (70).

After manual taxon identification, SpecieScan software was used to confirm taxon identifications (71). SpecieScan is a semiautomated process using R and Python to output correlated identifications from reference databases.

**Radiocarbon Dating.** Collagen extraction and pretreatment of the bone was conducted at the Higham Lab, Department of Evolutionary Anthropology, University of Vienna, Austria. The human bone underwent ABA (Acid-Base-Acid) and ultrafiltration protocols (72, 73). The bone was first demineralized with 0.5 M HCl twice for 2 h and then for ~24 h. Then the sample was washed three times with ultrapure Milli-Q water. Following this, a base wash was performed with 0.1 M NaOH for 30 min. Three more washes with Milli-Q water were then undertaken, with a 15-min treatment of 0.5 M HCl, followed by three Milli-Q water washes. The collagen was then gelatinized in pH3 water for 20 h at 75 °C. The gelatin was then filtered with a precleaned Eezee-filter, before being ultrafiltered using a precleaned Vivaspin® Turbo 15 30 kD MWCO and centrifuged until 1.0 mL of the >30 kD gelatin fraction remained (~20 min at 2,700 RPM). The gelatin was then placed into a clean glass tube with Milli-Q water and frozen before being freeze-dried for 48 h.

After freeze-drying the collagen yielded 24 mg of collagen. A subsample of the collagen was weighed into precleaned tin capsules then combusted and graphitized using an AGE3 system, before being accelerator mass spectrometry (AMS) dated at VERA, University of Vienna.

XAD-2 resin treatment was also undertaken on a subsample of the leftover collagen (8.33 mg) after freeze-drying. This method involves hydrolysis of the collagen into amino acids and their recovery using ion-exchange columns (74). The sample was hydrolyzed for 24 h at 110 °C, in 6 M HCl. A precleaned and preconditioned XAD column was used to allow the passage of the hydrolyzed amino acids. The column had a pre-fitted filter frit at the bottom, and a filter frit was placed above ~100 µL of XAD-2 resin. To clean the resin, 20 washes of 1 M HCl and 10 washes of 6 M HCl were passed through the column. Then the sample passed through the column, with one more wash of 1 mL of 6 M HCl, to collect any leftover amino acids. The sample was then placed into a Genevac system for 12 h, to remove the HCl acid through evaporation. ~100 µL of Milli-Q water was placed into the tube and then placed back in the Genevac for 1.5 h until full evaporation to remove trace remnants of HCl that might otherwise interfere with the graphitization process later. Following this, ~100 µL of Milli-Q was placed again into the sample tube to collect the sample, before transferal of the sample into a precleaned tin boat together with precleaned chromosorb for combustion, graphitization, and AMS dating.

OxCal v. 4.4.3 software (75) and the INTCAL20 calibration curve (76) was used to calibrate the result.

**Stable Isotopes.** After the collagen was weighed out for radiocarbon dating, 300 µg was weighed into a tin capsule for isotopic analysis of carbon and nitrogen by elemental analyzer-isotope ratio mass spectrometry (Thermo Scientific EA-Isolink with a Flash 2000 coupled to a Delta V Advantage isotope ratio mass spectrometer) in the Silver Laboratory (Large-Instrument Facility for Advanced Isotope Research) at the Center of Microbiology and Environmental Systems Science of the University of Vienna. Stable isotope values of all samples are measured relative to the laboratory standard alanine, which has a reproducibility of 0.07 wt% for N and 0.1 wt% for δ¹³C and δ¹⁵N values are reported relative to the Vienna Pee Dee belemnite and Atmospheric air standard, respectively, and are measured with an analytical precision of ± 0.1‰ for both δ¹³C and δ¹⁵N values.

**DNA Extraction, Library Construction, and Target Enrichment.** Ancient DNA extraction was carried out in the Pinhasi Lab, with dedicated clean-room facilities at the University of Vienna as per the protocol first described by Dabney et al. (2013) (40). Two bone powder samples were obtained using a dental drill and incubated overnight (18 h) at 37 °C on a rotation plate (800 rpm) in 1 mL of extraction buffer (final concentrations: 0.45 M Ethylenediaminetetraacetic acid, 0.25 mg/mL proteinase K). The supernatant was mixed with 13 mL of binding buffer (final concentrations: 5 M guanidine hydrochloride, 40% [v/v] isopropanol, 0.05% Tween-20, and 90 mM sodium acetate) in preassembled silica spin columns (Roche) following Korlić et al. (2015), (41). Following a 4-min centrifugation at 1,500 rpm, the columns were disassembled. The remaining spin columns were washed twice with 650 µL of PE purification wash buffer (Qiagen; in-house) with an intermediate centrifugation step at 6,000 rpm, followed by a dry spin. DNA was eluted in 50 µL EBT (Qiagen with 0.05% Tween-20 added) following a 5-min incubation at 37 °C and centrifugation at maximum speed for 1 min. The DNA concentrations were measured via the QuBit dsDNA DNA High Sensitivity kit. Blank samples at the beginning and the end of the batches were included to control potential cross-contamination.

Single-stranded libraries were built as per Kapp et al. (2021) (42). Based on initial DNA concentration, appropriate dilutions of single-stranded DNA binding protein (SSB), P5 and P7 adapter dilutions were prepared. 20 µL of DNA extract were mixed with 2 µL of the appropriate SSB dilution. The mix was incubated at 95 °C for 3 min and transferred onto an ice bath. While on ice, × 1 µL of the appropriate dilutions of P5 and P7 adapter dilutions were added together with 26 µL of the Santa Cruz Reaction mix. An incubation at 37 °C for 45 min followed. MinElute columns and buffers from the Qiagen MinElute PCR Purification kits were used for intermediate clean-up steps. qPCR was performed on 1:40 dilution of the resulting libraries in order to estimate the number of amplification cycles to take place post unique molecular indexing (UMI). UMIs from the Ancient Meyer collection were used to double-index the samples, following an amplification with the NebNext Q5U Master Mix DNA Polymerase (NEB). Indexed libraries were pooled with other samples and submitted to the Vienna Biocenter Core Facility for Next Generation Sequencing (VBCF NGS) where the samples were sequenced on a single lane of an Illumina NovaSeq X Plus single-read mode 100 bp. For mitochondrial enrichment we used the myBaits Expert Mito kit (Arbor Biosciences), following the protocols recommended by the provider.

**Bioinformatic Analysis.** Standard Illumina adapters were cut with leeHom v.1.2.18 in -ancient-dna mode and filtering for sequences longer than 30 bp was performed with seqtk (77). Mapping to hg19 and Revised Cambridge reference sequences was performed with bwa aln v. 0.7.18 with ancient DNA settings (SI Appendix, section E). Unmapped reads were removed with samtools v.1.20 and PCR duplicates were filtered out with GATK MarkDuplicates v.3.1.1. Typical Ancient DNA deamination patterns were confirmed with MapDamage v.2.2.2, contamination was assessed with schmutzi v.1.5.7 and filtering out of nondeaminated reads was performed with PMDtools python3 updated version (78-80). Mitochondrial genome sequence and quality was confirmed with samtools, mosdepth, and qualimap, and visualized with BamPlotter (81-83). A consensus sequence was called with ANGSD -doFasta (SI Appendix, section E) and multiple sequence alignment including a collection of present-day humans, Neanderthals, and Pan troglodytes was performed via the MEGA11 in-built MUSCLE algorithm (84, 85). The phylogenetic tree was constructed with IQ-tree (SI Appendix, section E) and visualized with Interactive Tree of Life (86, 87). Pairwise nucleotide differences between Star 1 and other genomes in the alignment were calculated using a custom Python script available on the GitHub page of the project. All code used is available at https://github.com/KonstantinaChe/Star1_Neanderthal.
<table><caption>**Table 2.** Specimens and dated materials including their time uncertainty used for constructing the primary habitat suitability model</caption><thead><tr><th>Specimen/dated material</th><th>latitude</th><th>longitude</th><th>min_age</th><th>max_age</th><th>Location</th><th>Age estimate method</th></tr></thead><tbody><tr><td>Vi-33-19</td><td>46.301562</td><td>16.079687</td><td>50,901</td><td>45,561</td><td>Box 1</td><td>Age at 68.3% probability</td></tr><tr><td>Vi-208</td><td>46.301562</td><td>16.079687</td><td>47,256</td><td>44,126</td><td>Box 1</td><td>Age at 68.3% probability</td></tr><tr><td>Vi-207</td><td>46.301562</td><td>16.079687</td><td>48,904</td><td>44,579</td><td>Box 1</td><td>Age at 68.3% probability</td></tr><tr><td>Mammoth from layer D - same layer as the discovered Neanderthals OxA-24944</td><td>50.366</td><td>19.294</td><td>52,257</td><td>46,646</td><td>Box 1</td><td>Age at 68.3% probability</td></tr><tr><td>Star 1</td><td>44.75</td><td>33.92</td><td>43,858</td><td>42,695</td><td>Box 2</td><td>Age at 68.3% probability</td></tr><tr><td>Mezmaiskaya 2</td><td>44.5540</td><td>39.5534</td><td>43,986</td><td>42,526</td><td>Box 2</td><td>Age at 68.3% probability</td></tr><tr><td>Chagyrskaya Level 6a where Neanderthal remains were found</td><td>51.263299</td><td>83.091628</td><td>53,000</td><td>48,200</td><td>Box 4</td><td>Weighted arithmetic mean age from four samples</td></tr><tr><td>Chagyrskaya Level 6b where Neanderthal remains were found</td><td>51.263299</td><td>83.091628</td><td>55,000</td><td>49,200</td><td>Box 4</td><td>Weighted arithmetic mean age from two samples</td></tr><tr><td>Chagyrskaya Level 6C where Neanderthal remains were found</td><td>51.263299</td><td>83.091628</td><td>57,400</td><td>52,400</td><td>Box 4</td><td>Weighted arithmetic mean age from nine samples</td></tr><tr><td>Layer 7 sediment where a Neanderthal was found</td><td>51.40</td><td>84.20</td><td>51,338</td><td>44,716</td><td>Box 4</td><td>Age at 68.3% probability</td></tr><tr><td>Layer 7 sediment where a Neanderthal was found</td><td>51.40</td><td>84.20</td><td>51,903</td><td>44,801</td><td>Box 4</td><td>Age at 68.3% probability</td></tr><tr><td>Denisova 11 age estimates</td><td>51.235</td><td>84.403</td><td>118,100</td><td>79,300</td><td>Box 4</td><td></td></tr></tbody></table>

Raw untrimmed sequencing data are available at the European Nucleotide Archive (ENA) under the accession number: PRJEB93921. All code used for data processing is available at https://github.com/KonstantinaChe/Star1_Neanderthal.

**Mahalanobis Distance Habitat Suitability Modeling.** Two Mahalanobis distance habitat suitability models were constructed to assess the ecological connections between three geographic regions: Croatia/Poland, the Crimean Peninsula, and the Altai region.

The primary model utilized only dated fossil data, comprising a limited dataset of 13 samples (Table 2). To assess the validity of the primary model, a second, more comprehensive model was developed. This validation model incorporated data from an extensive fossil and archaeological database (47, 48, 57, 66, 88) focused on the regions surrounding the sample locations of the primary model. The following coordinate boundaries were used for the geographical extent of these regions:

Croatia/Poland: longitude: 12.5°E-21°E, latitude: 44°N-53°N

Crimean Peninsula: longitude: 31°E-43.5°E, latitude: 43.5°N-49°N

Altai region: longitude: 80°E-89°E, latitude: 47.5°N-53°N

For both models, we employed a Monte Carlo sampling approach to address two key challenges: 1) the uneven distribution of data points at the three locations and 2) the variable temporal uncertainties associated with the data.

For each ensemble member, we randomly sampled an equal number of samples from each location to ensure balanced spatial representation. The primary model utilized two samples per location (totaling six samples per ensemble member), while the validation model used 10 samples per location (totaling 30 samples per ensemble member). Additionally, for each sample, a random age was selected from within the age uncertainty.

Sampling the data using this approach ensured that we account for the uneven sample size and varying length in age uncertainty without weighting specific locations or specimens more heavily. This sampling process was repeated 5,000 times.

For each ensemble member, we constructed habitat suitability models by extracting 1,000-y averages of annual mean temperature (T) and precipitation (P) from a 3 My quasi-transient climate simulation at a 5.75 × 3.75 resolution (61, 62). The extracted temperature and precipitation correspond to the nearest spatial-temporal point for each sample in the ensemble member constructed using the Monte Carlo sampling approach. Environmental conditions across the selected points were concatenated into a single matrix $C' = [T'(Z'), P(Z')]$, and used to compute the mean vector $-C'$ and inverse covariance matrix $\Sigma^{-1}$ of the sampled

conditions. Mahalanobis distance was then calculated across the entire climate grid using these statistics $d_M = \sqrt{(C' - C'')^T \Sigma^{-1} (C' - C'')}$). Suitability was derived as the exponential decay of Mahalanobis distance, i.e., $\exp(-d_M)$, where lower distances result in higher suitability.

**Data, Materials, and Software Availability.** Raw untrimmed sequencing data have been deposited in Star_1_Neanderthal (PRJEB93921-https://github.com/KonstantinaChe/Star1_Neanderthal)(89). Anonymized surface scan data have been deposited in Star 1 (https://www.virtual-anthropology.com/3d-data/free-data/)(90).

**ACKNOWLEDGMENTS.** This research was funded by the European Union (ERC grant number 101142352 "DISPERSE" awarded to T.H.). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union. Neither the European Union nor the granting authority can be held responsible for them. We thank Martin Dockner (Vienna Micro-CT Lab, University of Vienna) for scanning specimens. We are grateful to Ralf W. Schmitz (LVR-Landesmuseum Bonn) and Tracy Kivell, Philipp Gunz, and Mykolas Imbrasas (Max Planck Institute for Evolutionary Anthropology, Leipzig) for providing access to CT-scans of the Neanderthal 1 from Feldhofer Grotto and Bärbel Auffermann for providing surface scans of the same specimen. This project has been funded by the Vienna Science and Technology Fund (WWTF) [10.47379/VRG20001] to M.K. The computational results of this work have been achieved using the Life Science Compute Cluster of the University of Vienna. A.T. was supported by the Institute for Basic Science (IBS) Center for Climate Physics grant IBS-R028-D1.

Author affiliations: <sup>a</sup>Department of Evolutionary Anthropology, University of Vienna, Vienna 1030, Austria; <sup>b</sup>Human Evolution and Archaeological Science Network, Vienna 1030, Austria; <sup>c</sup>Department of Geosciences, University of Arizona, Tucson, AZ 85721; <sup>d</sup>Institut für Ur- und Frühgeschichte, Friedrich-Alexander-Universität (FAU) Erlangen-Nürnberg, Erlangen D-91054, Germany; <sup>e</sup>Institute of Archaeology, National Ukrainian Academy of Science, Kyiv 04210, Ukraine; <sup>f</sup>Muséum National d'Histoire Naturelle, Institut de Paléontologie Humaine, Paris 75013, France; <sup>g</sup>Institute for Basic Science (IBS) Center for Climate Physics, Busan 46241, South Korea; and <sup>h</sup>Department of Integrated Climate System Science, Pusan National University, Busan 46241, South Korea

Author contributions: E.M.P., K.C., E.Z., A.T., M.K., K.D., and T.H. designed research; E.M.P., K.C., E.Z., L.G.v.d.S., P.G.Š., J.N.V., G.W.W., A.T., and M.K. performed research; E.M.P., E.Z., M.P.C., M.G., E.V., and M.P.-M. contributed new reagents/analytic tools; E.M.P., K.C., E.Z., M.P.C., M.G., E.V., T.U., V.C., P.G.Š., J.N.V., G.W.W., and M.K. analyzed data; T.U., V.C., G.W.W., and A.T. edited paper; R.P. DNA lab provider and edited paper; and E.M.P., K.C., E.Z., M.K., K.D., and T.H. wrote the paper.
1. T. F. G. Higham *et al.*, The timing and spatio-temporal patterning of Neanderthal disappearance. *Nature* **512**, 306–309 (2014), 10.1038/nature13621.

2. R. E. Green *et al.*, A draft sequence of the neanderthal genome. *Science* **328**, 710–722 (2010), 10.1126/science.1188021.

3. K. Prüfer *et al.*, The complete genome sequence of a Neanderthal from the Altai Mountains. *Nature* **505**, 43–49 (2014), 10.1038/nature12886.

4. M. Hajdinjak *et al.*, Reconstructing the genetic history of late Neanderthals. *Nature* **555**, 652–656 (2018), 10.1038/nature26151.

5. E. Morin, E. Ready, A. Boilleau, C. Beauval, M.-P. Coumont, Problems of identification and quantification in archaeozoological analysis, Part II: Presentation of an alternative counting method. *J. Archaeol. Method. Theory* **24**, 938–973 (2017).

6. G. M. Smith *et al.*, The ecology, subsistence and diet of -45, 000-year-old Homo sapiens at Ilsenhöhle in Ranis, Germany. *Nat. Ecol. Evol.* **8**, 564–577 (2024), 10.1038/s41559-023-02303-6.

7. Z. Sin-Metviot *et al.*, Combining ZooMS and zoarchaeology to study Late Pleistocene hominin behaviour at Fumane (Italy). *Sci. Rep.* **9**, 12350 (2019), 10.1038/s41598-019-48706-z.

8. A. Pineda, P. Saladié, Beyond the problem of bone surface preservation in taphonomic studies of early and middle Pleistocene open-air sites. *J. Archaeol. Method Theory* **29**, 1090–1130 (2022), 10.1007/s10816-022-09550-0.

9. F. Welker, M. Soressi, W. Rendu, J.-J. Hublin, M. Collins, Using ZooMS to identify fragmentary bone from the late middle/early Upper Palaeolithic sequence of Les Cottés. *France. J. Archaeol. Sci.* **54**, 279–286 (2015), 10.1016/j.jas.2014.12.010.

10. K. K. Richter, M. C. Codlin, M. Seabrook, C. Warinner, A primer for ZooMS applications in archaeology. *Proc. Natl. Acad. Sci. U.S.A.* **119**, 1–10 (2022), 10.1073/pnas.2109323119.

11. F. Welker *et al.*, Palaeoproteomic evidence identifies archaic hominins associated with the Châtelperronian at the Grotte du Renne. *Proc. Natl. Acad. Sci. U.S.A.* **113**, 11162–11167 (2016), 10.1073/pnas.1605834113.

12. S. Brown *et al.*, Identification of a new hominin bone from Denisova Cave, Siberia using collagen fingerprinting and mitochondrial DNA analysis. *Sci. Rep.* **6**, 23559 (2016), 10.1007/srep23559.

13. M. Buckley, "Zooarchaeology by Mass Spectrometry (ZooMS) collagen fingerprinting for the species identification of archaeological bone fragments" in *Zooarchaeology in Practice*, C. Giovas, M. LeFebvre, Eds. (Springer, 2018), pp. 227–247, 10.1007/978-3-319-64763-0–12.

14. S. Brown *et al.*, Zooarchaeology through the lends of collagen fingerprinting at Denisova Cave. *Sci. Rep.* **11**, 15457 (2021), 10.1038/s41598-021-94731-2.

15. M. Fernández-García, M. Vidal-Cordasco, J. R. Jones, A. B. Marín-Arroyo, Reassessing palaeoenvironmental conditions during the Middle to Upper Palaeolithic transition in the Cantabrian region (Southwestern Europe). *Quat. Sci. Rev.* **301**, 107928 (2023), 10.1016/j.quascirev.2022.107928.

16. L. Spindler *et al.*, Dating the last Middle Palaeolithic of the Crimean Peninsula: New hydroxyproline AMS dates from the site of Kabazi II. *J. Hum. Evol.* **156**, 1–10 (2021), 10.1016/j.jhevol.2021.102996.

17. A. Burke, The view from Starosel'e: Faunal exploitation at a middle palaeolithic site in Western Crimea. *Int. J. Ostearchaeol.* **11**, 325–355 (2000).

18. A. Marks, Yu. Demidenko, K. Monigal, V. Usik, C. Ferring, "Starosel'e: The 1993–1996 excavations" in *The Middle Palaeolithic of the Western Crimea*, A. Marks, V. Chabai, Eds. (ERAUL 84, Liège, 1998), vol. 1, pp. 67–99.

19. A. E. Marks *et al.*, Starosel'e and the Starosel'e child: New excavations, new results. *Curr. Anthropol.* **38**, 112–123 (1997), 10.1086/204591.

20. K. Monigal *et al.*, Nouvelles découvertes de restes humains au site Paléolithique moyen de Starosel'e, Crimée (Ukraine). *Préhist. Eur.* **11**, 11–31 (1997).

21. V. P. Chabai, "The Late Middle and Early Upper Palaeolithic in Crimea (Ukraine)" in *Les premiers harnes modernes de la péninsule ibérienne Vol 17* "Actes du colloque de la commission VII de l'UISSP, Vila nova de Foz Cóa, 22-24 October 1998., J. Zilhão, T. Aubry, A. F. Carvalho, Eds. (Trabalhos de Arqueologia, 2005).

22. E. M. Pigott, T. Uthmeier, V. Chabai, T. F. G. Higham, The late middle and early Upper Palaeolithic in Crimea (Ukraine). Revier of the Neanderthal refugium hypothesis. *J. Paleolit. Archaeol.* **7**, 27 (2024), 10.1007/s41982-024-00194-y.

23. V. P. Chabai, A. E. Marks, K. Monigal, "Crimea in the context of the eastern European Middle Palaeolithic and early Upper Palaeolithic" in *The Middle Palaeolithic and early Upper Palaeolithic of Eastern Europe: Micoquian layers IIIa and III of the sites of Zaskalnaya VI (Kolosovskaya), anthropological records and context*. *Quat. Int.* **428**, 132–150 (2017), 10.1016/j.quaint.

24. G. A. Bonch-Osmolowski, The hand of fossil man from Kiik-KOBA rock-shelter, Vol. 2 of Palaeolithic of the Crimea series (USSR Acad of Science, Moscow-Leningrad: Vypusk 1941), Vol. 1 (in Russian).

25. V. N. Stepanchuk, S. V. Silyev, N. I. Khaldeeva, N. V. Kharlamova, S. B. Borutskaya, The last Neanderthals of Eastern Europe: Micoquian layers IIIa and III of the sites of Zaskalnaya VI (Kolosovskaya), anthropological records and context. *Quat. Int.* **428**, 132–150 (2017), 10.1016/j.quaint.

26. A. Picin *et al.*, New perspectives on Neanderthal dispersal and turnover from Stajnia Cave (Poland). *Sci. Rep.* **10**, 14778 (2020), 10.1038/s41598-020-7154-9.

27. T. V. Andreeva *et al.*, Genomic analysis of a novel Neanderthal from Mezmaiskaya Cave provides insights into the genetic relationships of Middle Palaeolithic population. *Sci. Rep.* **12**, 13106 (2012), 10.1038/s41598-022-16164-9.

28. K. A. Kolobova *et al.*, New data on Neanderthal behaviour in the Altai region, Russia. *Archaeol. Res. Asia* **3**, 100489 (2023), 10.1016/j.ara.2023.100489.

29. K. Kolobova *et al.*, Archaeological evidence for two separate dispersals of Neanderthals into Southern Siberia. *Proc. Natl. Acad. Sci. U.S.A.* **117**, 2879–2885 (2020).

30. A. E. Marks, K. Monigal, "Starosel'e 1993–1995: The lithic artifacts" in *The Middle Palaeolithic of Western Crimea*, A. E. Marks, V. P. Chabai, Eds. (ERAUL 84, Liège, 1998), vol. 1, pp. 117–166.

31. A. Burke, "Butchering and scavenging at the Middle Palaeolithic site of Starosel'e" in *The Middle Palaeolithic of Western Crimea*, V. P. Chabai, K. Monigal, Eds. (ERAUL 87, Liège, 1999), vol. 2, pp. 1–27.

32. B. L. Hardy, M. Kay, A. E. Marks, K. Monigal, Stone tool function at the palaeolithic sites of Starosel'e and Buran Kaya III, Crimea: Behavioural implications. *Proc. Natl. Acad. Sci. U.S.A.* **98**, 10972–10977 (2001), 10.1073/pnas.191384498.

33. A. Shpansky, Y. V. Kuzmin, Chronology of the MIS 3 megafauna in southeastern west Siberia and the possibility of late survival of the Khozarian steppe mammoth (Mammuthus trogonotherii chosaricus). *Radiocarbon* **63**, 575–584 (2021), 10.1017/RDC.2021.202.

34. B. Albouy, S. Paquin, M. Hinz, C. D. Wren, A. Burke, "The last of them: Investigating the Palaeogeography of the Last Neanderthals in Europe (Marine Isotopic Stage 3)" in *Modeling Human-Environment Interactions in and beyond Prehistoric Europe*, S. Seuru, B. Albouy, Eds. (Springer, 2023), pp. 27–45.

35. M. Buckley, S. W. Kansa, Collagen fingerprinting of bone and teeth remains from Domuztepe, South Eastern Turkey. *Archaeol. Anthropol. Sci.* **3**, 271–280 (2011), 10.1007/s12520-011-0066-z.

36. K. Douka *et al.*, FINDER project: Collagen fingerprinting (ZooMS) for the identification of new human fossils. *Antiquity* **93**, 1–6 (2019), 10.15184/agt.2019.3.

37. J. H. Gosman, Z. R. Hubbell, C. N. Shaw, R. M. Ryan, Development of cortical bone geometry in the human femoral and tibial diaphysis. *Anat. Rec. (Hoboken)* **296**, 774–787 (2013), 10.1002/ar.22688.

38. T. Devièse *et al.*, Increasing accuracy for the radiocarbon dating of sites occupied by the first Americans. *Quart. Sci. Rev.* **198**, 171–180 (2018), 10.1016/j.quascirev.2018.08.023.

39. L. G. van der Sluis, A. Zazzo, O. Tombret, F. Thil, J.-M. Pétillon, Testing the use of XAD resin to remove synthetic contamination from archaeological bone prior to radiocarbon dating. *Radiocarbon* **65**, 1160–1175 (2023), 10.1017/RDC.2023.100.

40. J. Nabday *et al.*, Complete mitochondrial genome sequence of a Middle Pleistocene bone bear reconstructed from ultraash DNA fragments. *Proc. Natl. Acad. Sci. U.S.A.* **110**, 15758–15763 (2013), 10.1073/pnas.1314445110.

41. P. Korlević *et al.*, Reducing microbial and human contamination in DNA extractions from ancient bones and teeth. *Biotechniques* **59**, 57–93 (2015), 10.2144/000114320.

42. J. D. Kapp, R. E. Green, B. Shapiro, A fast and efficient single-stranded genomic library preparation method optimized for ancient DNA. *J. Heredity* **112**, 241–249 (2021), 10.1093/jhered/esab012.

43. L. F. K. Kuderna *et al.*, A global catalog of whole-genome diversity from 233 primate species. *Science* **380**, 906–913 (2023), 10.1126/science.abn7829.

44. J. Krause *et al.*, The complete mitochondrial DNA genome of an unknown hominin from southern Siberia. *Nature* **464**, 894–897 (2010), 10.1038/nature08976.

45. M. Meyer *et al.*, A mitochondrial genome sequence of a hominin from Sima de los Huesos. *Nature* **505**, 403–406 (2013), 10.1038/nature12788.

46. Q. Fu *et al.*, Denisovan mitochondrial DNA from dental calculus of the >146,000-year-old Harbin cyanin. *Cell* **188**, 3919–3926 (2025), 10.1016/j.cell.2025.05.040.

47. J. Krause *et al.*, Neanderthals in Central Asia and Siberia. *Nature* **448**, 902–904 (2007), 10.1038/nature06193.

48. A. A. Formozov, The cave site of Starosel'e and its place in the Palaeolithic (Materials and investigations of the Archaeology of the USSR, Moscow, 1958), vol. 71.

49. V. P. Alexeyev, "Fossil man on the territory of the USSR and related problems" in *Les processus de l'hominiation*, D. Ferrembach, Ed. (CNRS, Paris, 1981), pp. 183–188.

50. V. P. Alexeyev, "Fossil man on the territory of the USSR and related problems" in *Les processus de l'hominiation*, D. Ferrembach, Ed. (CNRS, Paris, 1981), pp. 183–188.

51. V. P. Alexeyev, Position of the Staroselye find in the hominin sysem. *J. Hum. Evol.* **5**, 413–421 (1976), 10.1016/0047-2484(76)90083-X.

52. C. Posth *et al.*, Deeply divergent archaic mitochondrial genome provides lower time boundary for African gene flow into Neanderthals. *Nat. Commun.* **8**, 1 (2017), 10.1038/ncomms16046.

53. K. Kolobova *et al.*, Exploitation of the natural environment by Neanderthals from Chagyrskaya Cave (Altai). *Quat.* **66**, 7–21 (2021), 10.7485/QU66_1.

54. J. Ruan *et al.*, Climate shift orchestrated hominin interbreeding events across Eurasia. *Science* **381**, 699–704 (2023), 10.1126/science.add4459.

55. E. Coco, R. Lovita, Agent-based simulations reveal the possibility of multiple rapid northern routes for the second Neanderthal dispersal from Western to Eastern Crimea. *PLoS One* **20**, e0325693 (2025), 10.1371/journal.pone.0325693.

56. E. I. Zavala *et al.*, Pleistocene sediment DNA reveals hominin and faunal turnovers at Denisova Cave. *Nature* **595**, 399–403 (2021), https://doi.org/1038/s41586-021-03675-0.

57. T. Devièse, D. Comeskey, J. McCullagh, C. Bronk Ramsey, T. Higham, New protocol for compound-specific radiocarbon analysis of archaeological bones. *Rapid Commun. Mass Spectrom.* **32**, 373–379 (2017), 10.1002/rcom.8047.

58. A. Marom, J. O. S. McCullagh, T. F. G. Higham, A. A. Sinitsyn, R. E. M. Hedges, Single amino acid radiocarbon dating of Upper Palaeolithic modern humans. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 6787–6881 (2012), 10.1073/pnas.116328109.

59. F. Maffessoni *et al.*, A high-coverage Neandertal genome from Chagyrskaya Cave. *Proc. Natl. Acad. Sci. U.S.A.* **117**, 15132–15136 (2020), 10.1073/pnas.2004944117.

60. V. Slon *et al.*, The genome of the offspring of a Neandertal mother and a Denisovan father. *Nature* **561**, 113–116 (2018), 10.1038/s41586-018-0455-x.

61. A. Timmermann *et al.*, Climate effects on archaic human habitats and species successions. *Nature* **604**, 495–501 (2022), 10.1038/s41586-022-04600-9.

62. K.-S. Yun *et al.*, A transient coupled general circulation model (CGCM) simulation of the past 3 million years. * Clim. Past. *19*, 1951–1974 (2023), 10.5194/clp-19-191-2023.

63. A. P. Derevianko, S. V. Markin, N. A. Kulik, K. A. Kolobova, Lithic raw material exploitation in the Sibiriyackhia facies, the middle Palaeolithic of Altai. *Archaeol. Ethnol. Anthropol. Eurasia.* **41**, 89–103 (2013), 10.1016/j.aeae.2013.07.008.

64. K. A. Kolobova *et al.*, Exploitation of the natural environment by Neanderthals from Chagyrskaya Cave. *Quat. Int.* **66**, 7–32 (2019), 10.7485/QU66_1.

65. R. Lovita *et al.*, A newly Dated Late Pleistocene and holocene archaeological assemblage from Bukhtarma Cave in the Southern Altai Piedmont, East Kazakhstan. *J. Paleolit. Archaeol.* **7**, 21 (2024), 10.1007/s41982-024-00187-x.

66. K. Douka *et al.*, Age estimates for hominin fossils and the onset of the upper Palaeolithic at Denisova Cave. *Nature* **565**, 640–644 (2019), 10.1038/s41586-018-0870-z.

67. S. Brown *et al.*, Zoarchaeology by Mass Spectrometry (ZooMS) for Bone Material–AmBic Protocol v1 (protocols.io.bffddij16). (2020). https://doi.org/10.17504/protocols.io.bffddij16. Accessed 10 October 2025.

68. M. Strohalm, M. Hassman, B. Kosata, M. Kodicek, mMass data miner: An open source alternative for mass spectrometric data analysis. *Rapid Commun. Mass Spectrom.* **22**, 905–908 (2008), 10.1002/rcom.3444.

69. M. Bruckley, M. Collins, J. Thomas-Oates, J. C. Wilson, Species identification by analysis of bone collagen using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. *Rapid Commun. Mass Spectrom.* **23**, 3843–3854 (2009), 10.1002/rcom.4316.

70. M. Strohalm, D. Kavan, P. Novak, M. Volny, V. Havlicek, mMass 3: A crossplatform software environment for precise analysis of mass spectrometric data. *Anal. Chem.* **82** **11**, 4648–4651 (2010), 10.1021/ac100818g.

71. E. Végh, K. Douka, SpecieScan: Semi-automated taxonomic identification of bone collagen peptides from MALDI-ToF-MS. *Bioinformatics* **40**, btae054 (2024), 10.1093/bioinformatics/btae054.
72. F. Brock, T. Higham, P. Ditchfield, C. Bronk Ramsey, Current pretreatment methods for AMD radiocarbon dating at Oxford Radiocarbon Accelerator Unit (ORAU). *Radiocarbon* **52**, 103–112 (2010), 10.1017/S0033822200045069.

73. T. Higham, R. M. Jacobi, C. Bronk Ramsey, AMS radiocarbon dating of ancient bone using ultrafiltration. *Radiocarbon* **48**, 179–195 (2016), 10.1017/S0033822200066388.

74. T. W. Stafford, P. E. Hare, L. Currie, A. J. T. Jull, D. J. Donahue, Accelerator radiocarbon dating at the molecular level. *J. Archaeol. Sci.* **18**, 35–72 (1991), 10.1016/0305-4403(91)90078-4.

75. C. Bronk Ramsey, Bayesian analysis of radiocarbon dates. *Radiocarbon* **51**, 337–360 (2009), 10.1017/S0033822200033865.

76. P. Reimer et al., "The IntCal20 Northern manufacture?" in *Flakes Not Blades: The Role of Flake Production at the Onset of the Upper Palaeolithic*, A. Pastoors, M. Peresani, Eds. (Wissenschaftliche Schriften des Neanderthal Museums 5. Neanderthal Museum, Mettmann, Germany, 2020), pp. 239–260.

77. G. Renaud, U. Stenzel, J. Kelso, IeeHom: Adaptor trimming and merging for Illumina sequencing reads. *Nucleic Acids Res.* **42**, e141 (2014), 10.1093/nar/gku699.

78. H. Jónsson, A. Ginolhac, M. Schubert, P. L. F. Johnson, L. Orlando, MapDamage2.0: Fast approximate Bayesian estimates of ancient DNA damage parameters. *Bioinformatics* **29**, 1682–1684 (2013), 10.1093/bioinformatics/btt193.

79. G. Renaud, V. Slon, A. T. Duggan, J. Kelso, Schmutzi: Estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA. *Genome Biol.* **16**, 1 (2015), 10.1186/s13059-015-0776-0.

80. P. Skoglund et al., Separating endogenous ancient DNA from modern day contamination in a Siberian Neandertal. *Proc. Natl. Acad. Sci. U.S.A.* **111**, 2229–2234 (2014), 10.1073/pnas.1318934111.

81. B. S. Pedersen, A. R. Quinlan, Mosdepth: Quick coverage calculation for genomes and exomes. *Bioinformatics* **34**, 867–868 (2017), 10.1093/bioinformatics/btx699.

82. K. Okonechnikov, A. Conesa, F. García-Alcalde, Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. *Bioinformatics* **32**, 292–294 (2015), 10.1093/bioinformatics/btv566.

83. G. Meriam, MeriamGuellil/aDNA-BAMPlotter: aDNA-BAMPlotter. Zenodo (2021), 10.5281/zenodo.5676092.

84. T. S. Korneliussen, A. Albrechtsen, R. Nielsen, ANGSD: Analysis of next generation sequencing data. *BMC Bioinformatics* **15**, 1 (2014), 10.1186/s12859-014-0356-4.

85. K. Tamura, G. Stecher, S. Kumar, MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. **38**, 3022–3027 (2021), 10.1093/molbev/msab120.

86. J. Trifinopoulos, L. Nguyen, A. Von Haeseler, B. Q. Minh, W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. *Nucleic Acids Res.* **44**, W232-W235 (2016), 10.1093/nar/gkw256.

87. I. Letunic, P. Bork, Interactive Tree of Life (iTOL) v6: Recent updates to the phylogenetic tree display and annotation tool. *Nucleic Acids Res.* **52**, W78-W82 (2024), 10.1093/nar/gkae268.

88. M. Żarski, Stratigraphy and palaeoenvironment of Stajnia Cave (Southern Poland) with regard to habitation of the site by Neanderthals. *Geo. Quart.* **61**, 350–369 (2017), 10.7306/gq.1355.

89. K. Cheshmedzhieva, Raw untrimmed sequencing of Star 1. ENA. https://www.ebi.ac.uk/ena/browse/view/PRJEB93921. Accessed 10 October 2025.

90. G. Weber, Anonymized surface scan data of Star 1. Virtual Anthropology, University of Vienna. https://www.virtual-anthropology.com/3d-data/free-data/. Accessed 10 October 2025.