Genomic Diversity of Aurochs From a Mediterranean Ice‐Age Refugium

The aurochs, the wild ancestor of domestic cattle, was a keystone herbivore in Late Pleistocene Eurasian ecosystems and a major prey species for Palaeolithic hunter-gatherers. Despite its significance, the genetic structure of aurochs populations that survived the Last Glacial Maximum (LGM) remains poorly understood, especially in southeastern Europe. Here, we present the first directly dated ancient genomes of aurochs from the eastern Adriatic region, recovered from the Upper Palaeolithic site of Šandalja (Istria, Croatia). Two female individuals, dated to approximately 14,800-14,200 and 11,800-11,400 calibrated years before present, were sequenced for low-coverage whole genomes and near-complete mitochondrial genomes. Bayesian phylogenetic analyses and median-joining network reconstruction place both specimens within mitochondrial haplogroup P, the dominant European aurochs lineage. However, they do not cluster within the main P sub-haplogroup observed in most ancient aurochs samples and in modern cattle carrying P-type mitochondrial lineages. Instead, one specimen is placed within an 'alternative' P sub-haplogroup, whereas the position of the other appears more isolated and should be interpreted cautiously, as it may be influenced by limited sequence coverage and the resulting uncertainty in phylogenetic placement. At the nuclear genomic level, the two Šandalja aurochs show affinity to Late Pleistocene and Early Holocene aurochs from Italy. Although based on a limited number of specimens, this pattern is consistent with possible genetic connectivity across the Adriatic region, potentially associated with the now-submerged Great Adriatic Plain (GAP). Overall, our results suggest regional structure among Late Pleistocene aurochs, potentially associated with the exposed Adriatic Plain as a refugium and dispersal corridor between the Apennine and Balkan Peninsulas. By filling a major geographic and temporal gap in the aurochs genomic record, this study highlights the Adriatic Basin as a potentially overlooked centre of Pleistocene megafaunal diversity and refines models of postglacial recolonization and cattle evolutionary history.

Authors

Institutions

Publication Details

Journal
Molecular Ecology
Published
2026-06-29
DOI
https://doi.org/10.1111/mec.70449
Primary Topic
Forensic and Genetic Research
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Genomic Diversity of Aurochs From a Mediterranean Ice‐Age Refugium

Ivana Keršić, Preston Miracle, Ino Čurik, Ana Prohaska et al.
Molecular Ecology
Forensic and Genetic Research
article

Genomic Diversity of Aurochs From a Mediterranean Ice‐Age Refugium

Ivana Keršić, Preston Miracle, Ino Čurik, Ana Prohaska, Dinko Novosel, Maja Ferenčaković, Vlatka Čubrić-Ćurik, Zvonimir Mijadžiković, Aurélien Capitan, Siniša Radović, Mekki Boussaha, Eske Willerslev, Jesper Stenderup
article en

Abstract

The aurochs, the wild ancestor of domestic cattle, was a keystone herbivore in Late Pleistocene Eurasian ecosystems and a major prey species for Palaeolithic hunter-gatherers. Despite its significance, the genetic structure of aurochs populations that survived the Last Glacial Maximum (LGM) remains poorly understood, especially in southeastern Europe. Here, we present the first directly dated ancient genomes of aurochs from the eastern Adriatic region, recovered from the Upper Palaeolithic site of Šandalja (Istria, Croatia). Two female individuals, dated to approximately 14,800-14,200 and 11,800-11,400 calibrated years before present, were sequenced for low-coverage whole genomes and near-complete mitochondrial genomes. Bayesian phylogenetic analyses and median-joining network reconstruction place both specimens within mitochondrial haplogroup P, the dominant European aurochs lineage. However, they do not cluster within the main P sub-haplogroup observed in most ancient aurochs samples and in modern cattle carrying P-type mitochondrial lineages. Instead, one specimen is placed within an 'alternative' P sub-haplogroup, whereas the position of the other appears more isolated and should be interpreted cautiously, as it may be influenced by limited sequence coverage and the resulting uncertainty in phylogenetic placement. At the nuclear genomic level, the two Šandalja aurochs show affinity to Late Pleistocene and Early Holocene aurochs from Italy. Although based on a limited number of specimens, this pattern is consistent with possible genetic connectivity across the Adriatic region, potentially associated with the now-submerged Great Adriatic Plain (GAP). Overall, our results suggest regional structure among Late Pleistocene aurochs, potentially associated with the exposed Adriatic Plain as a refugium and dispersal corridor between the Apennine and Balkan Peninsulas. By filling a major geographic and temporal gap in the aurochs genomic record, this study highlights the Adriatic Basin as a potentially overlooked centre of Pleistocene megafaunal diversity and refines models of postglacial recolonization and cattle evolutionary history.

Molecular EcologyVol. 35(13)
University of Copenhagen (DK), Croatian Academy of Sciences and Arts (HR), University of Zagreb (HR), AgroParisTech (FR), University of Cambridge (GB), Microbiologie de l’alimentation au service de la santé (FR), Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement (FR), Sustainability Institute (ZA), Croatian Veterinary Institute (HR), University of Kaposvár (HU)
Life below water
Openalex Percentile: Top 7%
Forensic and Genetic Research
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.