The Genomic Legacy of Human Management and sex-biased Aurochs hybridization in Iberian Cattle

Cattle have been a valuable economic resource and cultural icon since prehistory. From the initial expansion of domestic cattle into Europe during the Neolithic period, taurine cattle (Bos taurus) and their wild ancestor, the aurochs (B. primigenius), had overlapping ranges leading to ample opportunities for intentional and unintentional hybridization. We performed a bioarchaeological analysis of 24 Bos remains from Iberia dating from the Mesolithic to the Roman period. The archaeogenomic dataset allows us to investigate the extent of domestic-wild hybridization over time, providing insight into the species’ behavior and human management by aligning changes with cultural and genomic transitions in the archaeological record. Our results show frequent hybridization during the Neolithic and Chalcolithic, likely reflecting a mix of hunting and herding or relatively unmanaged herds, with mostly male aurochs and female domestic cattle involved in hybridization. This is supported by isotopic evidence of ecological niche sharing, with only a few domestic cattle possibly being managed. The proportion of aurochs ancestry in domestic cattle remains relatively constant from about 4000 years ago, probably due to herd management and selection against hybrids, coinciding with other cultural transitions. The constant level of wild ancestry (∼20%) continues into modern western European breeds including the Spanish Lidia breed which is bred for its aggressiveness and fighting ability, but does not display elevated levels of aurochs ancestry. This study takes a genomic glance at the impact of human actions and wild introgression in the establishment of cattle as one of the most important domestic species today.

[1]  T. Günther,et al.  Ancient Sheep Genomes reveal four Millennia of North European Short-Tailed Sheep in the Baltic Sea region , 2023, bioRxiv.

[2]  R. D. da Fonseca,et al.  Iron age genomic data from Althiburos – Tunisia renew the debate on the origins of African taurine cattle , 2023, iScience.

[3]  D. Reich,et al.  On the limits of fitting complex models of population history to f-statistics , 2023, eLife.

[4]  J. Akey,et al.  The contribution of Neanderthal introgression to modern human traits , 2022, Current Biology.

[5]  J. Cañón,et al.  A novel missense variant in endothelin‐2 (EDN2) causes a growth and respiratory lethal syndrome in bovine , 2022, Animal genetics.

[6]  Andy M Jones,et al.  Dairying, diseases and the evolution of lactase persistence in Europe , 2022, Nature.

[7]  B. Guldbrandtsen,et al.  Large‐scale mitogenome sequencing reveals consecutive expansions of domestic taurine cattle and supports sporadic aurochs introgression , 2021, Evolutionary applications.

[8]  M. A. Anisimov,et al.  The origins and spread of domestic horses from the Western Eurasian steppes , 2021, Nature.

[9]  T. Sonstegard,et al.  Genome-wide local ancestry and direct evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus) , 2021, bioRxiv.

[10]  J. Arsuaga,et al.  Neolithic to Bronze Age economy and animal management revealed using analyses lipid residues of pottery vessels and faunal remains at El Portalón de Cueva Mayor (Sierra de Atapuerca, Spain) , 2021, Journal of Archaeological Science.

[11]  M. Clauss,et al.  Intensive human contact correlates with smaller brains: differential brain size reduction in cattle types , 2021, Proceedings of the Royal Society B.

[12]  D. Bradley,et al.  Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains , 2021, Proceedings of the National Academy of Sciences.

[13]  Ryan W. Schmidt,et al.  Origins and genetic legacy of prehistoric dogs , 2020, Science.

[14]  D. Hume,et al.  Whole genome analysis of water buffalo and global cattle breeds highlights convergent signatures of domestication , 2020, Nature Communications.

[15]  E. Scott,et al.  The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0–55 cal kBP) , 2020, Radiocarbon.

[16]  J. Lenstra,et al.  Evolution and domestication of the Bovini species. , 2020, Animal genetics.

[17]  J. Arsuaga,et al.  Metrical analysis of bovine bone remains from the Neolithic to the Bronze Age at the El Portalón site (Atapuerca, Burgos) in the Iberian context , 2020 .

[18]  J. Lenstra,et al.  Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication , 2020, Communications Biology.

[19]  D. Reich,et al.  Assessing the performance of qpAdm: a statistical tool for studying population admixture , 2020, bioRxiv.

[20]  L. Orlando,et al.  Animal domestication in the era of ancient genomics , 2020, Nature Reviews Genetics.

[21]  Cristina E. Valdiosera,et al.  Ancient pigs reveal a near-complete genomic turnover following their introduction to Europe , 2019, Proceedings of the National Academy of Sciences.

[22]  Benjamin S. Arbuckle,et al.  Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent , 2019, Science.

[23]  Arne Ludwig,et al.  Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series , 2019, Cell.

[24]  A. B. Marín‐Arroyo,et al.  Investigating prehistoric diet and lifeways of early farmers in central northern Spain (3000–1500 CAL BC) using stable isotope techniques , 2019, Archaeological and Anthropological Sciences.

[25]  L. Chikhi,et al.  Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle , 2018, BMC Genomics.

[26]  K. Gregersen,et al.  Ancient DNA analysis of Scandinavian medieval drinking horns and the horn of the last aurochs bull , 2018, Journal of Archaeological Science.

[27]  Alexandra Jamieson,et al.  Paleogenomics of Animal Domestication , 2018 .

[28]  Domingo C. Salazar-García,et al.  Reconstruction of human subsistence and husbandry strategies from the Iberian Early Neolithic: A stable isotope approach. , 2018, American journal of physical anthropology.

[29]  Benjamin S. Arbuckle,et al.  Ancient goat genomes reveal mosaic domestication in the Fertile Crescent , 2018, Science.

[30]  Hong Chen,et al.  Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia , 2018, Nature Communications.

[31]  Mattias Jakobsson,et al.  Four millennia of Iberian biomolecular prehistory illustrate the impact of prehistoric migrations at the far end of Eurasia , 2018, Proceedings of the National Academy of Sciences.

[32]  Marie Besse,et al.  The Beaker Phenomenon and the Genomic Transformation of Northwest Europe , 2018, Nature.

[33]  R. Schulting,et al.  Living different lives: Early social differentiation identified through linking mortuary and isotopic variability in Late Neolithic/ Early Chalcolithic north-central Spain , 2017, PloS one.

[34]  A. Rosell‐Melé,et al.  Feeding management strategies among the early Neolithic pigs in the NE of the Iberian Peninsula , 2017 .

[35]  Domingo C. Salazar-García,et al.  Late Neolithic-Chalcolithic socio-economical dynamics in Northern Iberia. A multi-isotope study on diet and provenance from Santimamiñe and Pico Ramos archaeological sites (Basque Country, Spain) , 2017, Quaternary International.

[36]  T. Stafford,et al.  Genome data on the extinct Bison schoetensacki establish it as a sister species of the extant European bison (Bison bonasus) , 2017, BMC Evolutionary Biology.

[37]  J. Krause,et al.  A Molecular Approach to the Sexing of the Triple Burial at the Upper Paleolithic Site of Dolní Věstonice , 2016, PloS one.

[38]  D. A. Magee,et al.  Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle , 2016, Heredity.

[39]  M. Sinding,et al.  The Draft Genome of Extinct European Aurochs and its Implications for De-Extinction , 2016 .

[40]  Matthew Oetjens,et al.  Ancient European dog genomes reveal continuity since the Early Neolithic , 2016, Nature Communications.

[41]  Ulrike H. Taron,et al.  Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild , 2016, bioRxiv.

[42]  M. Cubas,et al.  Re-evaluating the Neolithic: The Impact and the Consolidation of Farming Practices in the Cantabrian Region (Northern Spain) , 2016 .

[43]  M. Cubas,et al.  Erratum to: Re-evaluating the Neolithic: The Impact and the Consolidation of Farming Practices in the Cantabrian Region (Northern Spain) , 2016, Journal of World Prehistory.

[44]  Stinus Lindgreen,et al.  AdapterRemoval v2: rapid adapter trimming, identification, and read merging , 2016, BMC Research Notes.

[45]  Urmila Kulkarni-Kale,et al.  Analysis of Next-generation Sequencing Data in Virology - Opportunities and Challenges , 2016 .

[46]  Steven G. Schroeder,et al.  Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle , 2015, Genome Biology.

[47]  R. Murphy,et al.  DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals , 2015, Molecular ecology resources.

[48]  M. Bonilla,et al.  Investigating palaeodietary and social differences between two differentiated sectors of a Neolithic community, La Bòbila Madurell-Can Gambús (north-east Iberian Peninsula) , 2015 .

[49]  G. Barker,et al.  Detecting the T1 cattle haplogroup in the Iberian Peninsula from Neolithic to medieval times: new clues to continuous cattle migration through time , 2015 .

[50]  Swapan Mallick,et al.  Massive migration from the steppe was a source for Indo-European languages in Europe , 2015, Nature.

[51]  Anders Albrechtsen,et al.  ANGSD: Analysis of Next Generation Sequencing Data , 2014, BMC Bioinformatics.

[52]  J. Schibler,et al.  Incorporation of aurochs into a cattle herd in Neolithic Europe: single event or breeding? , 2014, Scientific Reports.

[53]  Domingo C. Salazar-García,et al.  Isotope evidence for the use of marine resources in the Eastern Iberian Mesolithic , 2014 .

[54]  Cristina E. Valdiosera,et al.  Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments , 2013, Proceedings of the National Academy of Sciences.

[55]  Jianquan Liu,et al.  The Yak genome database: an integrative database for studying yak biology and high-altitude adaption , 2012, BMC Genomics.

[56]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[57]  Matthias Meyer,et al.  Illumina sequencing library preparation for highly multiplexed target capture and sequencing. , 2010, Cold Spring Harbor protocols.

[58]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[59]  David R. Kelley,et al.  A whole-genome assembly of the domestic cow, Bos taurus , 2009, Genome Biology.

[60]  R. Hedges,et al.  Where the wild things are: aurochs and cattle in England , 2008, Antiquity.

[61]  J. Cañón,et al.  Genetic variation within the Lidia bovine breed. , 2008, Animal genetics.

[62]  Hans-Jürgen Bandelt,et al.  Mitochondrial genomes of extinct aurochs survive in domestic cattle , 2008, Current Biology.

[63]  M. Richards,et al.  Palaeodiets of Humans and Fauna at the Spanish Mesolithic Site of El Collado , 2006, Current Anthropology.

[64]  P. Taberlet,et al.  The origin of European cattle: evidence from modern and ancient DNA. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[65]  H. Ellegren,et al.  Prehistoric contacts over the Straits of Gibraltar indicated by genetic analysis of Iberian Bronze Age cattle. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[66]  T. Price,et al.  Diet of aurochs and early cattle in southern Scandinavia: evidence from 15N and 13C stable isotopes , 2005 .

[67]  R. Hedges,et al.  Stable isotope evidence for similarities in the types of marine foods used by Late Mesolithic humans at sites along the Atlantic coast of Europe , 1999 .

[68]  J. Arsuaga,et al.  Upper Pleistocene and Holocene palaeoenvironmental records in Cueva Mayor karst (Atapuerca, Spain) from different proxies: speleothem crystal fabrics, palynology and archaeology. , 2014 .

[69]  Martin Kircher,et al.  Analysis of high-throughput ancient DNA sequencing data. , 2012, Methods in molecular biology.

[70]  D. Carlos,et al.  Aproximación a la dieta de la población calcolítica de La Vital a través del análisis de isótopos estables del carbono y del nitrógeno sobre restos óseos , 2011 .

[71]  S. Salzberg,et al.  FLASH: fast length adjustment of short reads to improve genome assemblies , 2011, Bioinform..

[72]  Daniel Helmer,et al.  Identifying early domestic cattle from Pre-Pottery Neolithic sites on the Middle Euphrates using sexual dimorphism , 2005 .

[73]  N. Noe-Nygaarda,et al.  Diet of aurochs and early cattle in southern Scandinavia : evidence from 15 N and 13 C stable isotopes , 2005 .

[74]  A. V. D. Driesch,et al.  Early Animal Husbandry in the Northern Levant , 1999 .