# Stable population structure in Europe since the Iron Age, despite high mobility

Margaret L Antonio, Clemens L Weiss, Ziyue Gao, Susanna Sawyer, Victoria Oberreiter, Hannah M Moots, Jeffrey P Spence, Olivia Cheronet, Brina Zagorc, Elisa Praxmarer, et al.

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

## To cite this version:

Margaret L Antonio, Clemens L Weiss, Ziyue Gao, Susanna Sawyer, Victoria Oberreiter, et al. Stable population structure in Europe since the Iron Age, despite high mobility. eLife, 2024, 12, pp.1-41. 10.7554/eLife.79714. halshs-04378683

HAL Id: halshs-04378683

https://shs.hal.science/halshs-04378683

Submitted on 8 Jan 2024

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_2.png" /></figure>
# Stable population structure in Europe since the Iron Age, despite high mobility

Margaret L Antonio, Clemens L Weiss, Ziyue Gao, Susanna Sawyer, Victoria Oberreiter, Hannah M Moots, Jeffrey P Spence, Olivia Cheronet, Brina Zagorc, Elisa Praxmarer, et al.

► To cite this version:

Margaret L Antonio, Clemens L Weiss, Ziyue Gao, Susanna Sawyer, Victoria Oberreiter, et al. Stable population structure in Europe since the Iron Age, despite high mobility. eLife, 2024, 12, pp.1-41.
10.7554/eLife.79714. halshs-04378683

HAL Id: halshs-04378683

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

Submitted on 8 Jan 2024

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
<figure><img src="image_3.png" /></figure>

# Stable population structure in Europe since the Iron Age, despite high mobility

Margaret L Antonio¹†, Clemens L Weiss²†, Ziyue Gao³†, Susanna Sawyer⁴†, Victoria Oberreiter⁴, Hannah M Moots⁵,⁶, Jeffrey P Spence², Olivia Cheronet⁴, Brina Zagorc⁴, Elisa Praxmarer⁴, Kadir Toykan Özdogan⁷, Lea Demetz⁴, Pere Gelabert⁴, Daniel Fernandes⁴, Michaela Lucci⁸, Timka Alihodžić⁹, Selma Amrani¹⁰, Pavel Avetisyan¹¹, Christèle Baillif-Ducros¹², Željka Bedić¹³, Audrey Bertrand¹⁴, Maja Bilić¹⁵, Luca Bondioli¹⁶, Paulina Borówka¹⁷, Emmanuel Botte¹⁸, Josip Burmaz¹⁹, Domagoj Bužanić²⁰, Francesca Candilio²¹, Mirna Cvetko²⁰, Daniela De Angelis²², Ivan Drnić²³, Kristián Elschek²⁴, Mounir Fantar²⁵, Andrej Gaspari²⁶, Gabriella Gasperetti²⁷, Francesco Genchi²⁸, Snežana Golubović²⁹, Zuzana Hukeľová²⁴, Rimantas Jankauskas³⁰, Kristina Jelinčić Vučković³¹, Gordana Jeremić²⁹, Iva Kaić²⁰, Kevin Kazek³², Hamazasp Khachatryan³³, Anahit Khudaverdyan³⁴, Sylvia Kirchengast⁴, Miomir Korać²⁹, Valérie Kozlowski³⁵, Mária Krošláková²⁴, Dora Kušan Špalj³⁶, Francesco La Pastina³⁷, Marie Laguardia³⁸, Sandra Legrand³⁵, Tino Lelekovic³⁹, Tamara Leskovar²⁶, Wiesław Lorkiewicz¹⁷, Dženi Los¹⁹, Ana Maria Silva⁴⁰,⁴¹,⁴², Rene Masaryk⁴³, Vinka Matijević²⁰, Yahia Mehdi Seddik Cherifi⁴,⁴⁴,⁴⁵, Nicholas Meyer¹², Ilija Mikić²⁹, Nataša Miladinović-Radmilović²⁹, Branka Milošević Zakić⁴⁶, Lina Nacouzi⁴⁷, Magdalena Natuniewicz-Sekuła⁴⁸, Alessia Nava⁴⁹, Christine Neugebauer-Maresch⁵⁰,⁵¹, Jan Nováček⁵², Anna Osterholtz⁵³, Julianne Paige⁵⁴, Lujana Paraman⁵⁵, Dominique Pieri⁴⁷, Karol Pieta²⁴, Stefan Pop-Lazić²⁹, Matej Ruttkay²⁴, Mirjana Sanader²⁰, Arkadiusz Sołtysiak⁵⁶, Alessandra Sperduti²¹,⁵⁷, Tijana Stankovic Pesterac⁵⁸, Maria Teschler-Nicola⁴,⁵⁹, Iwona Teul⁶⁰, Domagoj Tončinić²⁰, Julien Trapp⁶¹, Dragana Vulović²⁹, Tomasz Waliszewski⁵⁶, Diethard Walter⁵², Milos Zivanovic⁶², Mohamed el Mostefa Filah⁶³, Morana Čaušević-Bully⁶⁴, Mario Šlaus¹³, Dušan Borić³⁷,⁶⁵, Mario Novak⁶⁶,⁶⁷, Alfredo Coppa⁴,³⁷,⁶⁶, Ron Pinhasi⁴,⁶⁸, Jonathan K Pritchard²,⁶⁹*

¹Biomedical Informatics Program, Stanford University, Stanford, United States;

²Department of Genetics, Stanford University, Stanford, United States; ³Department

of Genetics, University of Pennsylvania, Perelman School of Medicine, Philadelphia,

United States; ⁴Department of Evolutionary Anthropology, University of Vienna,

Vienna, Austria; ⁵Stanford Archaeology Center, Stanford University, Stanford,

United States; ⁶University of Chicago, Department of Human Genetics, Chicago,

United States; ⁷Department of History and Art History, Utrecht University, Utrecht,

Netherlands; ⁸Dipartimento di Storia Antropologia Religioni Arte Spettacolo,

Sapienza University, Rome, Italy; ⁹Archaeological Museum Zadar, Zadar, Croatia;

¹⁰LBEIG, Population Genetics & Conservation Unit, Department of Cellular and

Molecular Biology - Faculty of Biological Sciences, University of Sciences and

Technology Houari Boumediene, Algiers, Algeria; ¹¹National Academy of Sciences

of Armenia, Institute of Archaeology and Ethnography, Yerevan, Armenia; ¹²French

National Institute for Preventive Archaeological Research (INRAP), Zagreb, Croatia;

¹³Anthropological Centre, Croatian Academy of Sciences and Arts, Zagreb, Croatia;

¹⁴Université Gustave Eiffel - Laboratoire ACP, Paris, France; ¹⁵Palisada Ltd, Split,

**Competing interest:** See page 20

**Funding:** See page 20

**Received:** 23 April 2022

**Preprinted:** 16 May 2022

**Accepted:** 12 December 2023

**Reviewing Editor:** Christian R Landry, Université Laval, Canada

© Copyright Antonio, Weiss, Gao et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

<footnote>*For correspondence: pritch@stanford.edu</footnote>

<footnote>†These authors contributed
equally to this work</footnote>
<figure><img src="image_2.png" /></figure>

# Stable population structure in Europe since the Iron Age, despite high mobility

Margaret L Antonio¹†, Clemens L Weiss²†, Ziyue Gao³†, Susanna Sawyer⁴†, Victoria Oberreiter⁴, Hannah M Moots⁵,⁶, Jeffrey P Spence², Olivia Cheronet⁴, Brina Zagorc⁴, Elisa Praxmarer⁴, Kadir Toykan Özdogan⁷, Lea Demetz⁴, Pere Gelabert⁴, Daniel Fernandes⁴, Michaela Lucci⁸, Timka Alihodžić⁹, Selma Amrani¹⁰, Pavel Avetisyan¹¹, Christèle Baillif-Ducros¹², Željka Bedić¹³, Audrey Bertrand¹⁴, Maja Bilić¹⁵, Luca Bondioli¹⁶, Paulina Borówka¹⁷, Emmanuel Botte¹⁸, Josip Burmaz¹⁹, Domagoj Bužanić²⁰, Francesca Candilio²¹, Mirna Cvetko²⁰, Daniela De Angelis²², Ivan Drnić²³, Kristián Elschek²⁴, Mounir Fantar²⁵, Andrej Gaspari²⁶, Gabriella Gasperetti²⁷, Francesco Genchi²⁸, Snežana Golubović²⁹, Zuzana Hukeľová²⁴, Rimantas Jankauskas³⁰, Kristina Jelinčić Vučković³¹, Gordana Jeremić²⁹, Iva Kaić²⁰, Kevin Kazek³², Hamazasp Khachatryan³³, Anahit Khudaverdyan³⁴, Sylvia Kirchengast⁴, Miomir Korać²⁹, Valérie Kozlowski³⁵, Mária Krošláková²⁴, Dora Kušan Špalj³⁶, Francesco La Pastina³⁷, Marie Laguardia³⁸, Sandra Legrand³⁵, Tino Lelekovic³⁹, Tamara Leskovar²⁶, Wiesław Lorkiewicz¹⁷, Dženi Los¹⁹, Ana Maria Silva⁴⁰,⁴¹,⁴², Rene Masaryk⁴³, Vinka Matijević²⁰, Yahia Mehdi Seddik Cherifi⁴,⁴⁴,⁴⁵, Nicholas Meyer¹², Ilija Mikić²⁹, Nataša Miladinović-Radmilović²⁹, Branka Milošević Zakić⁴⁶, Lina Nacouzi⁴⁷, Magdalena Natuniewicz-Sekuła⁴⁸, Alessia Nava⁴⁹, Christine Neugebauer-Maresch⁵⁰,⁵¹, Jan Nováček⁵², Anna Osterholtz⁵³, Julianne Paige⁵⁴, Lujana Paraman⁵⁵, Dominique Pieri⁴⁷, Karol Pieta²⁴, Stefan Pop-Lazić²⁹, Matej Ruttkay²⁴, Mirjana Sanader²⁰, Arkadiusz Sołtysiak⁵⁶, Alessandra Sperduti²¹,⁵⁷, Tijana Stankovic Pesterac⁵⁸, Maria Teschler-Nicola⁴,⁵⁹, Iwona Teul⁶⁰, Domagoj Tončinić²⁰, Julien Trapp⁶¹, Dragana Vulović²⁹, Tomasz Waliszewski⁵⁶, Diethard Walter⁵², Milos Zivanovic⁶², Mohamed el Mostefa Filah⁶³, Morana Čaušević-Bully⁶⁴, Mario Šlaus¹³, Dušan Borić³⁷,⁶⁵, Mario Novak⁶⁶,⁶⁷, Alfredo Coppa⁴,³⁷,⁶⁶, Ron Pinhasi⁴,⁶⁸, Jonathan K Pritchard²,⁶⁹*

¹Biomedical Informatics Program, Stanford University, Stanford, United States;

²Department of Genetics, Stanford University, Stanford, United States; ³Department

of Genetics, University of Pennsylvania, Perelman School of Medicine, Philadelphia, United States; ⁴Department of Evolutionary Anthropology, University of Vienna,

Vienna, Austria; ⁵Stanford Archaeology Center, Stanford University, Stanford, United States; ⁶University of Chicago, Department of Human Genetics, Chicago,

United States; ⁷Department of History and Art History, Utrecht University, Utrecht, Netherlands; ⁸Dipartimento di Storia Antropologia Religioni Arte Spettacolo,

Sapienza University, Rome, Italy; ⁹Archaeological Museum Zadar, Zadar, Croatia;

¹⁰LBEIG, Population Genetics & Conservation Unit, Department of Cellular and Molecular Biology - Faculty of Biological Sciences, University of Sciences and

Technology Houari Boumediene, Algiers, Algeria; ¹¹National Academy of Sciences of Armenia, Institute of Archaeology and Ethnography, Yerevan, Armenia; ¹²French

National Institute for Preventive Archaeological Research (INRAP), Zagreb, Croatia;

¹³Anthropological Centre, Croatian Academy of Sciences and Arts, Zagreb, Croatia;

¹⁴Université Gustave Eiffel - Laboratoire ACP, Paris, France; ¹⁵Palisada Ltd, Split,

¹⁶French National Institute for Preventive Archaeological Research (INRAP), Zagreb, Croatia;

¹⁷Anthropological Centre, Croatian Academy of Sciences and Arts, Zagreb, Croatia;

¹⁸Université Gustave Eiffel - Laboratoire ACP, Paris, France; ¹⁹Palisada Ltd, Split,

**Competing interest:** See page 20

**Funding:** See page 20

**Received:** 23 April 2022

**Preprinted:** 16 May 2022

**Accepted:** 12 December 2023

**Reviewing Editor:** Christian R Landry, Université Laval, Canada

© Copyright Antonio, Weiss, Gao et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

<footnote>*For correspondence: pritch@stanford.edu</footnote>

<footnote>†These authors contributed equally to this work</footnote>
Croatia; ¹⁶Dipartimento dei Beni Culturali, Archeologia, Storia dell'arte, del Cinema e della Musica, Università di Padova, Padova, Italy; ¹⁷Department of Anthropology, Faculty of Biology and Environmental Protection, University of Lodz, Łódź, Poland; ¹⁸Aix Marseille Université, CNRS, Centre Camille Jullian, Aix-en-Provence, France; ¹⁹Kaducej Ltd, Split, Croatia; ²⁰Faculty of Humanities and Social Sciences, University of Zagreb, Zagreb, Croatia; ²¹Bioarchaeology Service, Museum of Civilizations, Rome, Italy; ²²Museo Nazionale Etrusco di Tarquinia, Direzione Generale Musei Lazio, Rome, Italy; ²³Archaeological Museum, Zagreb, Croatia; ²⁴Institute of Archaeology, Slovak Academy of Sciences, Bratislava, Slovakia; ²⁵Département des Monuments et des Sites Antiques - Institut National du Patrimoine INP, Tunis, Tunisia; ²⁶University of Ljubljana, Faculty of Arts, Department for Archaeology, Ljubljana, Slovenia; ²⁷Soprintendenza Archeologia, belle arti e paesaggio per le province di Sassari e Nuoro, Sassari, Italy; ²⁸Department of Oriental Studies, Sapienza University of Rome, Rome, Italy; ²⁹Institute of Archaeology Belgrade, Belgrade, Serbia; ³⁰Department of Anatomy, Histology and Anthropology, Vilnius University, Vilnius, Lithuania; ³¹Institute of Archaeology, Zagreb, Croatia; ³²Université de Lorraine, Centre de Recherche Universitaire Lorrain d' Histoire (CRULH), Nancy, France; ³³Department of Archaeologi, Shirak Centere of Armenological Studies National Academy of Sciences Republic of Armenia, Gyumri, Armenia; ³⁴Institute of Archaeology and Ethnography of the National Academy of Sciences of the Republic of Armenia, Yerevan, Armenia; ³⁵Musée Archéologique de l'Oise, Vendeuil- Caply, France; ³⁶Archaeological Museum in Zagreb, Zagreb, Croatia; ³⁷Department of Environmental Biology, Sapienza University of Rome, Rome, Italy; ³⁸Panthéon- Sorbonne University / UMR 7041 ArScAn, Paris, France; ³⁹Archaeology Division, Croatian Academy of Sciences and Arts, Zagreb, Croatia; ⁴⁰CIAS, Department of Life Sciences, University of Coimbra, Coimbra, Portugal; ⁴¹CEF - University of Coimbra, Coimbra, Portugal; ⁴²UNIARQ - University of Lisbon, Lisbon, Portugal; ⁴³Skupina STIK Zavod za preučevanje povezovalnih področij preteklosti in sedanjosti, Ljubljana, Slovenia; ⁴⁴Cardiolo-Oncology Research Collaborative Group (CORCG), Faculty of Medicine, Benyoucef Benkhedda University, Algiers, Algeria; ⁴⁵Molecular Pathology, University Paul Sabatier Toulouse III, Toulouse, France; ⁴⁶Museum of Croatian Archaeological Monuments, Split, Croatia; ⁴⁷L'Institut français du Proche- Orient, Beirut, Lebanon; ⁴⁸Institute of Archaeology and Ethnology Polish Academy of Sciences, Centre of Interdisciplinary Archaeological Research, Warsaw, Poland; ⁴⁹Department of Odontostomatological and Maxillofacial Sciences, Sapienza University of Rome, Rome, Italy; ⁵⁰Austrian Archaeological Institute, Austrian Academy of Sciences, Vienna, Austria; ⁵¹Institute of Prehistory and Early History, University of Vienna, Vienna, Austria; ⁵²Thuringia State Service for Cultural Heritage and Archaeology Weimar, Thuringia, Germany; ⁵³Mississippi State University, Starkville, United States; ⁵⁴University of Nevada, Las Vegas, United States; ⁵⁵Trogir Town Museum, Trogir, Croatia; ⁵⁶Faculty of Archaeology, University of Warsaw, Warsaw, Poland; ⁵⁷Dipartimento Asia, Africa e Mediterraneo, Università degli Studi di Napoli "L'Orientale", Naples, Italy; ⁵⁸Museum of Vojvodina, Novi Sad, Serbia; ⁵⁹Department of Anthropology, Natural History Museum Vienna, Vienna, Austria; ⁶⁰Department of General and Clinical Anatomy, Faculty of Medicine and Dentistry, Pomeranian Medical University in Szczecin, Szczecin, Poland; ⁶¹Musée de La Cour d'Or, Eurométropole de Metz, Metz, France; ⁶²Department of Archeology, Center for Conservation and Archeology of Montenegro, Cetinje, Montenegro; ⁶³Institut d'Archeologie, University Algiers 2, Algiers, Algeria; ⁶⁴Université de Franche Comté / UMR Chrono-Environnement, Besançon, France; ⁶⁵Department of Anthropology, New York University, New York, United States; ⁶⁶Department of Genetics, Harvard Medical
Croatia; ¹⁶Dipartimento dei Beni Culturali, Archeologia, Storia dell'arte, del Cinema e della Musica, Università di Padova, Padova, Italy; ¹⁷Department of Anthropology, Faculty of Biology and Environmental Protection, University of Lodz, Łódź, Poland; ¹⁸Aix Marseille Université, CNRS, Centre Camille Jullian, Aix-en-Provence, France; ¹⁹Kaducej Ltd, Split, Croatia; ²⁰Faculty of Humanities and Social Sciences, University of Zagreb, Zagreb, Croatia; ²¹Bioarchaeology Service, Museum of Civilizations, Rome, Italy; ²²Museo Nazionale Etrusco di Tarquinia, Direzione Generale Musei Lazio, Rome, Italy; ²³Archaeological Museum, Zagreb, Croatia; ²⁴Institute of Archaeology, Slovak Academy of Sciences, Bratislava, Slovakia; ²⁵Département des Monuments et des Sites Antiques - Institut National du Patrimoine INP, Tunis, Tunisia; ²⁶University of Ljubljana, Faculty of Arts, Department for Archaeology, Ljubljana, Slovenia; ²⁷Soprintendenza Archeologia, belle arti e paesaggio per le province di Sassari e Nuoro, Sassari, Italy; ²⁸Department of Oriental Studies, Sapienza University of Rome, Rome, Italy; ²⁹Institute of Archaeology Belgrade, Belgrade, Serbia; ³⁰Department of Anatomy, Histology and Anthropology, Vilnius University, Vilnius, Lithuania; ³¹Institute of Archaeology, Zagreb, Croatia; ³²Université de Lorraine, Centre de Recherche Universitaire Lorrain d' Histoire (CRULH), Nancy, France; ³³Department of Archaeologi, Shirak Centere of Armenological Studies National Academy of Sciences Republic of Armenia, Gyumri, Armenia; ³⁴Institute of Archaeology and Ethnography of the National Academy of Sciences of the Republic of Armenia, Yerevan, Armenia; ³⁵Musée Archéologique de l'Oise, Vendeuil- Caply, France; ³⁶Archaeological Museum in Zagreb, Zagreb, Croatia; ³⁷Department of Environmental Biology, Sapienza University of Rome, Rome, Italy; ³⁸Panthéon- Sorbonne University / UMR 7041 ArScAn, Paris, France; ³⁹Archaeology Division, Croatian Academy of Sciences and Arts, Zagreb, Croatia; ⁴⁰CIAS, Department of Life Sciences, University of Coimbra, Coimbra, Portugal; ⁴¹CEF - University of Coimbra, Coimbra, Portugal; ⁴²UNIARQ - University of Lisbon, Lisbon, Portugal; ⁴³Skupina STIK Zavod za preučevanje povezovalnih področij preteklosti in sedanjosti, Ljubljana, Slovenia; ⁴⁴Cardiolo-Oncology Research Collaborative Group (CORCG), Faculty of Medicine, Benyoucef Benkhedda University, Algiers, Algeria; ⁴⁵Molecular Pathology, University Paul Sabatier Toulouse III, Toulouse, France; ⁴⁶Museum of Croatian Archaeological Monuments, Split, Croatia; ⁴⁷L'Institut français du Proche- Orient, Beirut, Lebanon; ⁴⁸Institute of Archaeology and Ethnology Polish Academy of Sciences, Centre of Interdisciplinary Archaeological Research, Warsaw, Poland; ⁴⁹Department of Odontostomatological and Maxillofacial Sciences, Sapienza University of Rome, Rome, Italy; ⁵⁰Austrian Archaeological Institute, Austrian Academy of Sciences, Vienna, Austria; ⁵¹Institute of Prehistory and Early History, University of Vienna, Vienna, Austria; ⁵²Thuringia State Service for Cultural Heritage and Archaeology Weimar, Thuringia, Germany; ⁵³Mississippi State University, Starkville, United States; ⁵⁴University of Nevada, Las Vegas, United States; ⁵⁵Trogir Town Museum, Trogir, Croatia; ⁵⁶Faculty of Archaeology, University of Warsaw, Warsaw, Poland; ⁵⁷Dipartimento Asia, Africa e Mediterraneo, Università degli Studi di Napoli "L'Orientale", Naples, Italy; ⁵⁸Museum of Vojvodina, Novi Sad, Serbia; ⁵⁹Department of Anthropology, Natural History Museum Vienna, Vienna, Austria; ⁶⁰Department of General and Clinical Anatomy, Faculty of Medicine and Dentistry, Pomeranian Medical University in Szczecin, Szczecin, Poland; ⁶¹Musée de La Cour d'Or, Eurométropole de Metz, Metz, France; ⁶²Department of Archeology, Center for Conservation and Archeology of Montenegro, Cetinje, Montenegro; ⁶³Institut d'Archeologie, University Algiers 2, Algiers, Algeria; ⁶⁴Université de Franche Comté / UMR Chrono-Environnement, Besançon, France; ⁶⁵Department of Anthropology, New York University, New York, United States; ⁶⁶Department of Genetics, Harvard Medical
School, Boston, United States; ⁶⁷Centre for Applied Bioanthropology, Institute for Anthropological Research, Zagreb, Croatia; ⁶⁸Human Evolution and Archaeological Sciences, University of Vienna, Vienna, Austria; ⁶⁹Department of Biology, Stanford University, Stanford, United States

**Abstract** Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000–3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

**Editor's evaluation**

This important study provides an impressive dataset containing more than 200 novel ancient human genome sequences and a creative, robust, novel approach for studying human migration across time. The authors' conclusions are well supported by the data, and the methods used are convincing and solid. This paper will be of great interest to population geneticists and other scholars in the field of paleogenomics.

**Introduction**

Ancient DNA (aDNA) sequencing has provided immense insight into previously unanswered questions about human population history. Initially, sequencing efforts were focused on identifying the main ancestry groups and transitions during prehistoric times, for which there is no written record. Recently, aDNA sampling has expanded to more recent times, allowing the study of movements of people using genetic data alongside the well-studied historical record. However, we lack a comprehensive assessment of historical genetic structure, including characterizing genetic heterogeneity and interactions across regions. Integrating historical period genetics will be instrumental to better understanding the development of European and Mediterranean population structure from prehistoric to present-day.

Prehistoric ancient genomes have allowed disentangling the movements of people and technologies across two major demographic transitions in prehistoric western Eurasia: first the farming transition ~7500 BCE (Lazaridis et al., 2014; Skoglund et al., 2012), and later the Bronze Age Steppe migrations ~3500 BCE (Haak et al., 2015). Over the course of generations, genetically differentiated peoples across western Eurasia came together and admixed. As a result, most present-day European genomes can be modeled as a three-way mixture of these prehistoric groups: Western Hunter-Gatherers, Neolithic farmers, and Bronze Age Herders from the Steppe (Haak et al., 2015; Lazaridis et al., 2014) with minor contributions from other groups (Antonio et al., 2019; Fernandes et al., 2020; Lazaridis et al., 2016; Morrison et al., 2020; Mathieson et al., 2018). These ancestry components are present at different proportions across western Eurasia, leading to a pattern where the genetic structure of Europe mirrors its geography (Novembre et al., 2008).
School, Boston, United States; ⁶⁷Centre for Applied Bioanthropology, Institute for Anthropological Research, Zagreb, Croatia; ⁶⁸Human Evolution and Archaeological Sciences, University of Vienna, Vienna, Austria; ⁶⁹Department of Biology, Stanford University, Stanford, United States

**Abstract** Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000–3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

**Editor's evaluation**

This important study provides an impressive dataset containing more than 200 novel ancient human genome sequences and a creative, robust, novel approach for studying human migration across time. The authors' conclusions are well supported by the data, and the methods used are convincing and solid. This paper will be of great interest to population geneticists and other scholars in the field of paleogenomics.

**Introduction**

Ancient DNA (aDNA) sequencing has provided immense insight into previously unanswered questions about human population history. Initially, sequencing efforts were focused on identifying the main ancestry groups and transitions during prehistoric times, for which there is no written record. Recently, aDNA sampling has expanded to more recent times, allowing the study of movements of people using genetic data alongside the well-studied historical record. However, we lack a comprehensive assessment of historical genetic structure, including characterizing genetic heterogeneity and interactions across regions. Integrating historical period genetics will be instrumental to better understanding the development of European and Mediterranean population structure from prehistoric to present-day.

Prehistoric ancient genomes have allowed disentangling the movements of people and technologies across two major demographic transitions in prehistoric western Eurasia: first the farming transition ~7500 BCE (Lazaridis et al., 2014; Skoglund et al., 2012), and later the Bronze Age Steppe migrations ~3500 BCE (Haak et al., 2015). Over the course of generations, genetically differentiated peoples across western Eurasia came together and admixed. As a result, most present-day European genomes can be modeled as a three-way mixture of these prehistoric groups: Western Hunter-Gatherers, Neolithic farmers, and Bronze Age Herders from the Steppe (Haak et al., 2015; Lazaridis et al., 2014) with minor contributions from other groups (Antonio et al., 2019; Fernandes et al., 2020; Lazaridis et al., 2016; Morrison et al., 2020; Mathieson et al., 2018). These ancestry components are present at different proportions across western Eurasia, leading to a pattern where the genetic structure of Europe mirrors its geography (Novembre et al., 2008).
Given that the major ancestry components of present-day west Eurasians were largely established by the end of the Bronze Age, it is unclear how and what types of demographic processes impacted the genetic make-up of western Eurasia over the last ~3000 years, from the end of the Bronze Age to present-day. Recent studies of historical period genomes from individual regions shed light on this question; they paint a picture of heterogeneity and mobility, rather than of stable population structure. In the city of Rome alone, the population was dynamic and harbored a large diversity of ancestries from across Europe and the Mediterranean from the Iron Age (~1000 BCE) through the Imperial Roman period (27 BCE-300 CE; Antonio et al., 2019). Historical genomes from the Iberian Peninsula also highlight gene flow from across the Mediterranean (Olalde et al., 2019).

These regional reports fit well with archaeological and historical records. By the Iron Age, sea travel was already common, enabling peoples from across the Mediterranean to come into contact for trade (Abulafia, 2011; Broodbank, 2013). Subsequently, the Roman Empire leveraged its organization, labor force, and military prowess to build upon existing waterways and roads throughout Europe and create a united Mediterranean for the only time in history (Beard, 2015; Harper, 2017; Symonds, 2017). Not only did the Empire provide a means for movement, it also provided a reason for individuals to move. Empire building activities, broadly categorized into military, labor, and trade, pulled in people and resources from inside and outside the Empire (Scheidel, 2019).

We sequenced 204 new historical period genomes from across Europe and the Mediterranean to more comprehensively investigate the Roman Empire's impact on the genetic landscape suggested by these regional reports of heterogeneous, mobile populations. By analyzing genetic similarities between individuals across historical Eurasia, we were able to quantify individual movements during this time. Based on population genetic simulations, we explore potential explanations of how population structure may be maintained in the face of frequent individual dispersal.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Timeline of new and published genomes. (A) 204 newly reported genomes (black circles) are shown alongside published genomes (gray circles), ordered by time and region (colored the same way as in B). (B) Sampling locations of newly reported (black) and published (gray) genomes are indicated by diamonds, sized according to the number of genomes at each location.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 1:

**Figure supplement 1.** Detailed map of locations for newly reported samples.
Given that the major ancestry components of present-day west Eurasians were largely established by the end of the Bronze Age, it is unclear how and what types of demographic processes impacted the genetic make-up of western Eurasia over the last ~3000 years, from the end of the Bronze Age to present-day. Recent studies of historical period genomes from individual regions shed light on this question; they paint a picture of heterogeneity and mobility, rather than of stable population structure. In the city of Rome alone, the population was dynamic and harbored a large diversity of ancestries from across Europe and the Mediterranean from the Iron Age (~1000 BCE) through the Imperial Roman period (27 BCE-300 CE; Antonio et al., 2019). Historical genomes from the Iberian Peninsula also highlight gene flow from across the Mediterranean (Olalde et al., 2019).

These regional reports fit well with archaeological and historical records. By the Iron Age, sea travel was already common, enabling peoples from across the Mediterranean to come into contact for trade (Abulafia, 2011; Broodbank, 2013). Subsequently, the Roman Empire leveraged its organization, labor force, and military prowess to build upon existing waterways and roads throughout Europe and create a united Mediterranean for the only time in history (Beard, 2015; Harper, 2017; Symonds, 2017). Not only did the Empire provide a means for movement, it also provided a reason for individuals to move. Empire building activities, broadly categorized into military, labor, and trade, pulled in people and resources from inside and outside the Empire (Scheidel, 2019).

We sequenced 204 new historical period genomes from across Europe and the Mediterranean to more comprehensively investigate the Roman Empire's impact on the genetic landscape suggested by these regional reports of heterogeneous, mobile populations. By analyzing genetic similarities between individuals across historical Eurasia, we were able to quantify individual movements during this time. Based on population genetic simulations, we explore potential explanations of how population structure may be maintained in the face of frequent individual dispersal.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Timeline of new and published genomes. (A) 204 newly reported genomes (black circles) are shown alongside published genomes (gray circles), ordered by time and region (colored the same way as in B). (B) Sampling locations of newly reported (black) and published (gray) genomes are indicated by diamonds, sized according to the number of genomes at each location.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 1:

**Figure supplement 1.** Detailed map of locations for newly reported samples.
# Results

## 204 new historical genomes from Europe and the Mediterranean

We collected whole genomes from 204 individuals across 53 archaeological sites in 18 countries spanning Europe and the Mediterranean (Figure 1—figure supplement 1), 26 of these individuals were recently reported (Moots et al., 2022). This collection includes the first historical genomes (Iron Age and later, i.e. after 1000 BCE) from present-day Armenia, Algeria, Austria, and France. Dates for 126 samples were directly determined through radiocarbon dating, and were used alongside archaeological contexts to infer dates for the remaining samples.

DNA was extracted from either the powdered cochlear portion of the petrous bone (n=203) or from teeth (n=1). Libraries were partially treated with uracil-DNA glycosylase (UDG) and screened for ancient DNA damage patterns, high endogenous DNA content, and low contamination. We performed whole genome sequencing to a median depth of 0.92 x (0.16x to 2.38x).

For downstream integration with published data, pseudohaploid genotypes were called for the 1240 k SNP panel (Mathieson et al., 2015), resulting in a median of 685,058 SNPs (167,000–1,029,345) per sample. We analyzed newly reported genomes in conjunction with 2033 present-day genomes, 1998 prehistoric genomes, and 764 published historical period genomes (Clemente et al., 2021; Kovacevic et al., 2014, Mallick et al., 2023; Pagani et al., 2016; Saupe et al., 2021; Žegarac et al., 2021, primary AADR sources cited in Materials and methods). Genomes were grouped by regions and time periods (Figure 1) and analyzed using principal component analysis (PCA) and *qpAdm* modeling (Haak et al., 2015).

## Local historical population structure varies across regions

To investigate historical population structure, we categorized the data into 14 geographical regions, split into three sub-periods of the historical period: Iron Age (1000–1 BCE), Imperial Rome & Late

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Armenia: two homogeneous genetic clusters distinguished by a temporal shift. (A) Sampling locations of ancient genomes (open circles) colored by their genetic cluster identified using *qpAdm* modeling. (B) Date ranges for the genomes: each line represents the 95% confidence interval for the radiocarbon date or the upper and lower limit of the inferred date, and the point represents the midpoint of that range. (C) Projections of the genomes onto a PCA of present-day genomes (gray points labeled by their population). Present-day genomes from Armenia are shown with dark gray open circles.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 2:

**Figure supplement 1.** Principal component analysis of present-day genomes from Europe and the Mediterranean.

**Figure supplement 2.** Ancestry clusters identified within regions.

**Figure supplement 3.** SNP coverage comparison across cluster sizes and downstream outlier status.
# Results

## 204 new historical genomes from Europe and the Mediterranean

We collected whole genomes from 204 individuals across 53 archaeological sites in 18 countries spanning Europe and the Mediterranean (Figure 1—figure supplement 1), 26 of these individuals were recently reported (Moots et al., 2022). This collection includes the first historical genomes (Iron Age and later, i.e. after 1000 BCE) from present-day Armenia, Algeria, Austria, and France. Dates for 126 samples were directly determined through radiocarbon dating, and were used alongside archaeological contexts to infer dates for the remaining samples.

DNA was extracted from either the powdered cochlear portion of the petrous bone (n=203) or from teeth (n=1). Libraries were partially treated with uracil-DNA glycosylase (UDG) and screened for ancient DNA damage patterns, high endogenous DNA content, and low contamination. We performed whole genome sequencing to a median depth of 0.92 x (0.16x to 2.38x).

For downstream integration with published data, pseudohaploid genotypes were called for the 1240 k SNP panel (Mathieson et al., 2015), resulting in a median of 685,058 SNPs (167,000–1,029,345) per sample. We analyzed newly reported genomes in conjunction with 2033 present-day genomes, 1998 prehistoric genomes, and 764 published historical period genomes (Clemente et al., 2021; Kovacevic et al., 2014, Mallick et al., 2023; Pagani et al., 2016; Saupe et al., 2021; Žegarac et al., 2021, primary AADR sources cited in Materials and methods). Genomes were grouped by regions and time periods (Figure 1) and analyzed using principal component analysis (PCA) and *qpAdm* modeling (Haak et al., 2015).

## Local historical population structure varies across regions

To investigate historical population structure, we categorized the data into 14 geographical regions, split into three sub-periods of the historical period: Iron Age (1000–1 BCE), Imperial Rome & Late

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Armenia: two homogeneous genetic clusters distinguished by a temporal shift. (A) Sampling locations of ancient genomes (open circles) colored by their genetic cluster identified using *qpAdm* modeling. (B) Date ranges for the genomes: each line represents the 95% confidence interval for the radiocarbon date or the upper and lower limit of the inferred date, and the point represents the midpoint of that range. (C) Projections of the genomes onto a PCA of present-day genomes (gray points labeled by their population). Present-day genomes from Armenia are shown with dark gray open circles.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 2:

**Figure supplement 1.** Principal component analysis of present-day genomes from Europe and the Mediterranean.

**Figure supplement 2.** Ancestry clusters identified within regions.

**Figure supplement 3.** SNP coverage comparison across cluster sizes and downstream outlier status.
Antiquity (1–700 CE), and Medieval Ages & Early Modern (700–1950 CE). We then characterized inter-individual heterogeneity within these spatio-temporal groups by examining (1) variation of projections onto a PCA space of present-day genomes (Figure 2—figure supplement 1), (2) genetic groups identified by *qpAdm* and clustering across time within a region, and (3) admixture modeling of genetic groups.

A majority of regions have highly heterogeneous populations in at least one historical time period (Figure 2—figure supplement 2). This is illustrated by both the visual spread in PCA and the genetically distinct clusters of individuals based on pairwise modeling with *qpAdm* (Haak et al., 2015; Harney et al., 2021). On average, we identified 10 genetic clusters within each region present during the historical period, with a minimum of two and a maximum of 23. With genetically similar samples grouped together, we have more power relative to individual-level analyses when performing admixture modeling on clusters of interest using *qpAdm*.

Regional vignettes reveal various patterns of historical population structure. In Armenia, for example, the population is highly homogeneous at any given time (Figure 2). After the Copper Age, there are two distinct genetic clusters, separated by a temporal split around 772–403 BCE (Figure 2BC). The earlier cluster (C1) includes newly reported samples (n=5) from Beniamin and published ones (n=6) from five other sites. This cluster cannot be modeled by any single source of ancestry using existing data. The later cluster (C3), which contains newly reported samples (n=12) from Beniamin dating between 403 BCE-500 CE, is genetically similar to present-day Armenians (excluding two Kurdish individuals; Figure 2C). Despite the split, there is evidence of partial continuity between the earlier and later clusters: the later (C3) can be modeled using around 50% of the earlier cluster (C1) and an additional source of Steppe ancestry. Historical genomes from Northern Europe, particularly newly reported genomes from Lithuania and Poland, exhibit a similar level of homogeneity (Figure 2—figure supplement 2).

In contrast to the homogeneity of the Armenian population, most of the regions, including Italy, Southeastern Europe, and Western Europe, had strikingly heterogeneous populations. Newly collected samples reinforce previous findings of high heterogeneity in Rome, including a large portion of the population having affinities for present-day Near Eastern populations (Antonio et al., 2019; Posth et al., 2021; Figure 3—figure supplement 1). Interestingly, Southeastern European and Western

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Southeastern Europe: highly heterogeneous Imperial Roman and Late Antiquity period population. (A) Sampling locations of genetic clusters are represented by a single point per location. Outlier ancestries are black stars, all others are open circles colored by genetic cluster. (B) Colored bars span the minimum and maximum of the date ranges of samples (95% confidence interval from radiocarbon dating or archaeological range). Points are the mean of an individual's date range. (C) Projections of the ancient genomes onto a PCA of present-day genomes (gray points). Population labels for the PCA reference space are shown in Figure 2C. Present-day genomes from Southeastern Europe are shown with dark gray open circles.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 3:

**Figure supplement 1.** Population structure of Italy during the Imperial Roman and Late Antiquity period.
Antiquity (1–700 CE), and Medieval Ages & Early Modern (700–1950 CE). We then characterized inter-individual heterogeneity within these spatio-temporal groups by examining (1) variation of projections onto a PCA space of present-day genomes (Figure 2—figure supplement 1), (2) genetic groups identified by *qpAdm* and clustering across time within a region, and (3) admixture modeling of genetic groups.

A majority of regions have highly heterogeneous populations in at least one historical time period (Figure 2—figure supplement 2). This is illustrated by both the visual spread in PCA and the genetically distinct clusters of individuals based on pairwise modeling with *qpAdm* (Haak et al., 2015; Harney et al., 2021). On average, we identified 10 genetic clusters within each region present during the historical period, with a minimum of two and a maximum of 23. With genetically similar samples grouped together, we have more power relative to individual-level analyses when performing admixture modeling on clusters of interest using *qpAdm*.

Regional vignettes reveal various patterns of historical population structure. In Armenia, for example, the population is highly homogeneous at any given time (Figure 2). After the Copper Age, there are two distinct genetic clusters, separated by a temporal split around 772–403 BCE (Figure 2BC). The earlier cluster (C1) includes newly reported samples (n=5) from Beniamin and published ones (n=6) from five other sites. This cluster cannot be modeled by any single source of ancestry using existing data. The later cluster (C3), which contains newly reported samples (n=12) from Beniamin dating between 403 BCE-500 CE, is genetically similar to present-day Armenians (excluding two Kurdish individuals; Figure 2C). Despite the split, there is evidence of partial continuity between the earlier and later clusters: the later (C3) can be modeled using around 50% of the earlier cluster (C1) and an additional source of Steppe ancestry. Historical genomes from Northern Europe, particularly newly reported genomes from Lithuania and Poland, exhibit a similar level of homogeneity (Figure 2—figure supplement 2).

In contrast to the homogeneity of the Armenian population, most of the regions, including Italy, Southeastern Europe, and Western Europe, had strikingly heterogeneous populations. Newly collected samples reinforce previous findings of high heterogeneity in Rome, including a large portion of the population having affinities for present-day Near Eastern populations (Antonio et al., 2019; Posth et al., 2021; Figure 3—figure supplement 1). Interestingly, Southeastern European and Western

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Southeastern Europe: highly heterogeneous Imperial Roman and Late Antiquity period population. (A) Sampling locations of genetic clusters are represented by a single point per location. Outlier ancestries are black stars, all others are open circles colored by genetic cluster. (B) Colored bars span the minimum and maximum of the date ranges of samples (95% confidence interval from radiocarbon dating or archaeological range). Points are the mean of an individual's date range. (C) Projections of the ancient genomes onto a PCA of present-day genomes (gray points). Population labels for the PCA reference space are shown in Figure 2C. Present-day genomes from Southeastern Europe are shown with dark gray open circles.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 3:

**Figure supplement 1.** Population structure of Italy during the Imperial Roman and Late Antiquity period.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Western Europe: heterogeneous Imperial Roman and Late Antiquity period population. (A) Sampling locations of genetic clusters are represented by a single point per location. Outlier ancestries are black stars, all others are open circles colored by genetic cluster. (B) Colored bars span the minimum and maximum of the date ranges of samples (95% confidence interval from radiocarbon dating or archaeological range). Points are the mean of an individual's date range. (C) Projections of the ancient genomes onto a PCA of present-day genomes (gray points). Population labels for the PCA reference space are shown in **Figure 2C**. Present-day genomes from Southeastern Europe are shown with dark gray open circles.</figcaption></figure>

European individuals during the Imperial Roman & Late Antiquity period also exhibit high heterogeneity, on par with that of contemporaneous Italy (Figures 3 and 4).

Furthermore, these ancestries are often shared across regions. In Southeastern Europe, a core group of individuals have ancestry similar to that of present-day and contemporaneous Central Europeans (C6), while other clusters have ancestry similar to that of Northern Europeans (C7) and Eastern Mediterraneanans (C8) (Figure 3C). These ancestry groups are found in contemporaneous Italy and Western Europe as well (Figure 4C, Figure 3—figure supplement 1). We also observe individuals of eastern nomadic ancestry, similar to that of Sarmatian individuals previously reported, in both Western Europe (C8, n=2) and Southeastern Europe (C2, n=2).

Overall, we see remarkable local genetic heterogeneity as well as cross-regional similarities which point to common ancestry sources and, on a broader scale, demographic events affecting different regions in similar ways.

## At least 7-11% of historical individuals are ancestry outliers

The high regional genetic heterogeneity with long range, cross-regional similarities suggests historical populations were highly mobile. We therefore sought to quantify the amount of movement during the historical period by estimating the proportion of individuals who are ancestry outliers with respect to all individuals found in the same region. We considered an individual an outlier if they belonged to an ancestry cluster that is underrepresented (consisting of fewer than 5% of individuals in a region or at most two individuals) within their sampling region from the Bronze Age up to present-day. To focus on first-generation migrants as well as long-range movements, we further identified outlier individuals who can be modeled as 100% (i.e. 'one-component model') of a majority ancestry cluster found in a different region.

In total, we identified 11% of individuals as outliers, and could connect 7% of individuals to a putative source in a different region (Figure 5A). Based on the regions where these outliers and their sources originated, we created a network to illustrate their movements (Figure 5B). This network reveals the interconnectedness of Europe and the Mediterranean during the historical period. For example, as discussed above, the Armenian population is quite homogeneous (Figure 2). Unsurprisingly, no outliers were found within Armenia; however, we found outlier individuals in the Levant
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Western Europe: heterogeneous Imperial Roman and Late Antiquity period population. (A) Sampling locations of genetic clusters are represented by a single point per location. Outlier ancestries are black stars, all others are open circles colored by genetic cluster. (B) Colored bars span the minimum and maximum of the date ranges of samples (95% confidence interval from radiocarbon dating or archaeological range). Points are the mean of an individual's date range. (C) Projections of the ancient genomes onto a PCA of present-day genomes (gray points). Population labels for the PCA reference space are shown in **Figure 2C**. Present-day genomes from Southeastern Europe are shown with dark gray open circles.</figcaption></figure>

European individuals during the Imperial Roman & Late Antiquity period also exhibit high heterogeneity, on par with that of contemporaneous Italy (Figures 3 and 4).

Furthermore, these ancestries are often shared across regions. In Southeastern Europe, a core group of individuals have ancestry similar to that of present-day and contemporaneous Central Europeans (C6), while other clusters have ancestry similar to that of Northern Europeans (C7) and Eastern Mediterraneanans (C8) (Figure 3C). These ancestry groups are found in contemporaneous Italy and Western Europe as well (Figure 4C, Figure 3—figure supplement 1). We also observe individuals of eastern nomadic ancestry, similar to that of Sarmatian individuals previously reported, in both Western Europe (C8, n=2) and Southeastern Europe (C2, n=2).

Overall, we see remarkable local genetic heterogeneity as well as cross-regional similarities which point to common ancestry sources and, on a broader scale, demographic events affecting different regions in similar ways.

## At least 7-11% of historical individuals are ancestry outliers

The high regional genetic heterogeneity with long range, cross-regional similarities suggests historical populations were highly mobile. We therefore sought to quantify the amount of movement during the historical period by estimating the proportion of individuals who are ancestry outliers with respect to all individuals found in the same region. We considered an individual an outlier if they belonged to an ancestry cluster that is underrepresented (consisting of fewer than 5% of individuals in a region or at most two individuals) within their sampling region from the Bronze Age up to present-day. To focus on first-generation migrants as well as long-range movements, we further identified outlier individuals who can be modeled as 100% (i.e. 'one-component model') of a majority ancestry cluster found in a different region.

In total, we identified 11% of individuals as outliers, and could connect 7% of individuals to a putative source in a different region (Figure 5A). Based on the regions where these outliers and their sources originated, we created a network to illustrate their movements (Figure 5B). This network reveals the interconnectedness of Europe and the Mediterranean during the historical period. For example, as discussed above, the Armenian population is quite homogeneous (Figure 2). Unsurprisingly, no outliers were found within Armenia; however, we found outlier individuals in the Levant
<figure><img src="image_8.png" /><figcaption>**B. Network of sources for detected outliers**</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**C. Examples of individual regions**</figcaption></figure>

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

**Figure 5.** Ancestry outliers and their potential sources. (A) The proportions of outliers in each region were determined by individual pairwise qpAdm modeling followed by clustering. (B) Sources were inferred by one component qpAdm modeling of resulting clusters with all genetic clusters in the dataset. In the network visualizations, nodes are regions and directed edges are drawn from sources to outliers (i.e. potential migrants). The full network of source to outlier is shown. (C) Examples of individual regions are shown in greater detail.

The online version of this article includes the following figure supplement(s) for figure 5:

**Figure supplement 1.** Lack of sex-bias amongst outliers with valid qpAdm sources.

**Figure supplement 2.** Distances of outliers to their candidate sources.

**Figure supplement 3.** Example routes and travel times across the Roman Empire.

and Italy who can be putatively traced back to Armenia according to their ancestry (Figure 5C; blue outgoing arrows from Armenia). In contrast, the heterogeneous population in Italy connects it to many other regions, with bi-directional movement in most cases. In North Africa, outliers found in Iron Age Tunisia (Moots et al., 2022) indicate movements from many regions in Europe, and North African-like outliers were found in Italy and Austria (Western Europe). North African ancestry in Italy is supported by a single previously reported individual from the Imperial Roman period (R132; Antonio
A. Proportions of outliers

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

B. Network of sources for detected outliers

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

C. Examples of individual regions

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Ancestry outliers and their potential sources. (A) The proportions of outliers in each region were determined by individual pairwise qqAdm modeling followed by clustering. (B) Sources were inferred by one component qqAdm modeling of resulting clusters with all genetic clusters in the dataset. In the network visualizations, nodes are regions and directed edges are drawn from sources to outliers (i.e. potential migrants). The full network of source to outlier is shown. (C) Examples of individual regions are shown in greater detail.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 5:

**Figure supplement 1.** Lack of sex-bias amongst outliers with valid qqAdm sources.

**Figure supplement 2.** Distances of outliers to their candidate sources.

**Figure supplement 3.** Example routes and travel times across the Roman Empire.

and Italy who can be putatively traced back to Armenia according to their ancestry (Figure 5C; blue outgoing arrows from Armenia). In contrast, the heterogeneous population in Italy connects it to many other regions, with bi-directional movement in most cases. In North Africa, outliers found in Iron Age Tunisia (Moots et al., 2022) indicate movements from many regions in Europe, and North African-like outliers were found in Italy and Austria (Western Europe). North African ancestry in Italy is supported by a single previously reported individual from the Imperial Roman period (R132; Antonio
et al., 2019). Similar North African ancestry in Western Europe is supported by a single individual, R10667, from Wels, Austria, a site located on the frontier of the Roman Empire (C28 in Figure 4). This individual from Austria can be modeled using Canary Islander individuals from the Medieval Ages or an Iron Age outlier (distinguished by having more sub-Saharan ancestry) from Kerkouane, a Punic city near Carthage in modern-day Tunisia.

The 7% estimate for outliers with source should be considered conservative for the proportion of 'non-local' individuals. There are several cases where a cluster comprises more than 5% of the individuals in the region, but are clearly of a different ancestry than the majority and seem to be transient (only found in a single sub-period of the historical period). For example, in Southeastern Europe (Figure 3B), Imperial Roman & Late Antiquity individuals in C8 are (1) of distant ancestry (Near Eastern) and (2) not found in previous or subsequent time periods. However, since there are five individuals in this cluster, it does not meet our strict criteria for outlier consideration. Additionally, many clusters of underrepresented ancestry cannot be modeled as one-component models because they are recently admixed (i.e. require two or more ancestry components) or of ancestry not sampled elsewhere. Thus, we expect the actual proportion of individuals involved in long distance movements to be higher than reported here.

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Relatively stable population structure from Bronze Age to present-day. (A) Overall genetic differentiation between populations (measured by $F_{ST}$) and its relationship to geographical distance (spatial structure) is similar from Bronze Age onward. Confidence intervals were calculated through a bootstrap procedure, using 200 bootstrap replicates. (B) In PC space, each genome is represented by a point, colored based on their origin (for present-day individuals) or sampling location (for historical samples). The PC space is established by present-day samples (bottom), onto which either historical period (middle) or prehistoric genomes (top) were projected. For projections, the present-day samples are shown in gray, and their extent is visualized by a gray polygon.</figcaption></figure>
et al., 2019). Similar North African ancestry in Western Europe is supported by a single individual, R10667, from Wels, Austria, a site located on the frontier of the Roman Empire (C28 in Figure 4). This individual from Austria can be modeled using Canary Islander individuals from the Medieval Ages or an Iron Age outlier (distinguished by having more sub-Saharan ancestry) from Kerkouane, a Punic city near Carthage in modern-day Tunisia.

The 7% estimate for outliers with source should be considered conservative for the proportion of 'non-local' individuals. There are several cases where a cluster comprises more than 5% of the individuals in the region, but are clearly of a different ancestry than the majority and seem to be transient (only found in a single sub-period of the historical period). For example, in Southeastern Europe (Figure 3B), Imperial Roman & Late Antiquity individuals in C8 are (1) of distant ancestry (Near Eastern) and (2) not found in previous or subsequent time periods. However, since there are five individuals in this cluster, it does not meet our strict criteria for outlier consideration. Additionally, many clusters of underrepresented ancestry cannot be modeled as one-component models because they are recently admixed (i.e. require two or more ancestry components) or of ancestry not sampled elsewhere. Thus, we expect the actual proportion of individuals involved in long distance movements to be higher than reported here.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Relatively stable population structure from Bronze Age to present-day. (A) Overall genetic differentiation between populations (measured by $F_{ST}$) and its relationship to geographical distance (spatial structure) is similar from Bronze Age onward. Confidence intervals were calculated through a bootstrap procedure, using 200 bootstrap replicates. (B) In PC space, each genome is represented by a point, colored based on their origin (for present-day individuals) or sampling location (for historical samples). The PC space is established by present-day samples (bottom), onto which either historical period (middle) or prehistoric genomes (top) were projected. For projections, the present-day samples are shown in gray, and their extent is visualized by a gray polygon.</figcaption></figure>
## Spatial population structure is relatively stable in the last 3,000 years

The remarkable amount of heterogeneity and mobility in the historical period leads to the question of what impact this might have had on population structure over time. To investigate this, we sought to quantify the overall change in population structure across time, from prehistoric to present-day. To assess the spatial structure of population differentiation, we calculated $F_{ST}$ across groups of individuals on a sliding spatial grid in each time period and related it to their mean geographic distance. In each time period, we recovered the classical pattern of isolation-by-distance (Figure 6A), where individuals closer in geographic space are also more similar genetically. Across time periods, we see a large decrease in overall $F_{ST}$ from the Mesolithic & Neolithic periods to the Bronze Age (approximately 10,000–2300 BCE), coinciding with the major prehistoric migrations (Haak et al., 2015; Lazaridis et al., 2014). From the Bronze Age onward, however, $F_{ST}$ does not decrease further with time, indicating that the level of genetic differentiation across space is relatively stable from the Bronze Age to present-day.

To assess not only the amount, but also the structure of geographic population differentiation, we compared the 'genetic maps' of historical period and present-day genomes. To construct these 'maps', we performed principal component analysis on 829 present-day European and Mediterranean genomes sampled across geographical space (Figure 6B, bottom) and projected historical period genomes onto the same PC space. Echoing close correspondence between genetic structure and geographic space in present-day Europeans (Novembre et al., 2008), we recovered similar spatial structure for historical samples as well, although noisier due to a narrower sampling distribution and higher local genetic heterogeneity (Figure 6B, middle). The similarity in structure between present-day and historical period is especially striking in comparison to a projection of prehistoric genomes, which shows much weaker correspondence to the present-day PCA as well as to geographic space (Figure 6B, top). Together, our analyses indicate that European and Mediterranean population structure has been relatively stable over the last 3000 years.

This raises the question: is it surprising for stable population structure to be maintained in the presence of ~7–11% long-range migration? To address this, we simulated Wright-Fisher populations evolving neutrally in continuous space. In these simulations, spatial population structure is established through local mate choice and limited dispersal, which we calibrated to approximately match the spatial differentiation observed in historical-period Europe (Figure 6A, Figure 7A and Figure 7—figure supplement 1, maximum $F_{ST}$ of ~0.03). We then allowed a proportion of the population to disperse longer distances, empirically matching the migration distances we observed in the data during the historical period (Figure 7—figure supplement 2). Even with long-range dispersal as low as 4%, we observe decreasing $F_{ST}$ over 120 generations (~3000 years with a generation time of 25 years) as individuals become less differentiated genetically across space (Figure 7B). At 8%, $F_{ST}$ decreases dramatically within 120 generations as spatial structure collapses to the point that it is hardly detectable in the first two principal components (Figure 7C). These simulations indicate that under a basic spatial population genetics model we would expect structure to collapse by present-day given the levels of movement we observe.

## Discussion

In summary, we observed largely stable spatial population structure across western Eurasia and high mobility of people evidenced by local genetic heterogeneity and cross-regional connections. These two observations are seemingly incompatible with each other under standard population genetics assumptions.

A possible explanation for this apparent paradox is that our simulations did not capture some key features of human behavior and population dynamics. In the simulated populations, migration implies both movement and reproduction with local random mate choice. However, in real human populations migration can be more complex: people do not necessarily reproduce where they migrate, and reproduction is not necessarily random. We hypothesize that in the historical period there was an increasing decoupling of movement and reproduction, compared to prehistoric times. For the spread of Farmer and Steppe ancestry, we know that these prehistoric migrations would take hundreds of years to traverse the continent (Allentoft et al., 2015; Haak et al., 2015; Lazaridis et al., 2016). In contrast, in the historical period, there were dense travel networks of roads and waterways as well
## Spatial population structure is relatively stable in the last 3,000 years

The remarkable amount of heterogeneity and mobility in the historical period leads to the question of what impact this might have had on population structure over time. To investigate this, we sought to quantify the overall change in population structure across time, from prehistoric to present-day. To assess the spatial structure of population differentiation, we calculated $F_{ST}$ across groups of individuals on a sliding spatial grid in each time period and related it to their mean geographic distance. In each time period, we recovered the classical pattern of isolation-by-distance (Figure 6A), where individuals closer in geographic space are also more similar genetically. Across time periods, we see a large decrease in overall $F_{ST}$ from the Mesolithic & Neolithic periods to the Bronze Age (approximately 10,000–2300 BCE), coinciding with the major prehistoric migrations (Haak et al., 2015; Lazaridis et al., 2014). From the Bronze Age onward, however, $F_{ST}$ does not decrease further with time, indicating that the level of genetic differentiation across space is relatively stable from the Bronze Age to present-day.

To assess not only the amount, but also the structure of geographic population differentiation, we compared the 'genetic maps' of historical period and present-day genomes. To construct these 'maps', we performed principal component analysis on 829 present-day European and Mediterranean genomes sampled across geographical space (Figure 6B, bottom) and projected historical period genomes onto the same PC space. Echoing close correspondence between genetic structure and geographic space in present-day Europeans (Novembre et al., 2008), we recovered similar spatial structure for historical samples as well, although noisier due to a narrower sampling distribution and higher local genetic heterogeneity (Figure 6B, middle). The similarity in structure between present-day and historical period is especially striking in comparison to a projection of prehistoric genomes, which shows much weaker correspondence to the present-day PCA as well as to geographic space (Figure 6B, top). Together, our analyses indicate that European and Mediterranean population structure has been relatively stable over the last 3000 years.

This raises the question: is it surprising for stable population structure to be maintained in the presence of ~7–11% long-range migration? To address this, we simulated Wright-Fisher populations evolving neutrally in continuous space. In these simulations, spatial population structure is established through local mate choice and limited dispersal, which we calibrated to approximately match the spatial differentiation observed in historical-period Europe (Figure 6A, Figure 7A and Figure 7—figure supplement 1, maximum $F_{ST}$ of ~0.03). We then allowed a proportion of the population to disperse longer distances, empirically matching the migration distances we observed in the data during the historical period (Figure 7—figure supplement 2). Even with long-range dispersal as low as 4%, we observe decreasing $F_{ST}$ over 120 generations (~3000 years with a generation time of 25 years) as individuals become less differentiated genetically across space (Figure 7B). At 8%, $F_{ST}$ decreases dramatically within 120 generations as spatial structure collapses to the point that it is hardly detectable in the first two principal components (Figure 7C). These simulations indicate that under a basic spatial population genetics model we would expect structure to collapse by present-day given the levels of movement we observe.

## Discussion

In summary, we observed largely stable spatial population structure across western Eurasia and high mobility of people evidenced by local genetic heterogeneity and cross-regional connections. These two observations are seemingly incompatible with each other under standard population genetics assumptions.

A possible explanation for this apparent paradox is that our simulations did not capture some key features of human behavior and population dynamics. In the simulated populations, migration implies both movement and reproduction with local random mate choice. However, in real human populations migration can be more complex: people do not necessarily reproduce where they migrate, and reproduction is not necessarily random. We hypothesize that in the historical period there was an increasing decoupling of movement and reproduction, compared to prehistoric times. For the spread of Farmer and Steppe ancestry, we know that these prehistoric migrations would take hundreds of years to traverse the continent (Allentoft et al., 2015; Haak et al., 2015; Lazaridis et al., 2016). In contrast, in the historical period, there were dense travel networks of roads and waterways as well
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Simulation of population structure with and without long-range dispersal. (A) A base model of spatial structure is established by calibrating per-generation dispersal rate to generate a maximum $F_{ST}$ of ~0.03 across the maximal spatial distance, and visualized using PCA. In addition to this base dispersal, either 4% (B) or 8% (C) of individuals disperse longer distances, and the effect is tracked by analyzing spatial $F_{ST}$ through time, as well as PCA after 120 generations of long-range dispersal.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 7:

**Figure supplement 1.** A sigmaDisp - *N* parameter pair was chosen to closely approximate the observed $F_{ST}$max of ~0.03 using grid search across a range of parameter pairs.

**Figure supplement 2.** A sigmaDispLR parameter was chosen to qualitatively resemble long-range dispersal distances observed in the data, by comparing the distribution of distances under long-range dispersal (outliers) to randomly chosen distances given the spatial distribution of samples.

as clear incentives for cross-Mediterranean and cross-continental movement (Abulafia, 2011; Beard, 2015; Broodbank, 2013; Symonds, 2017). This enabled people to travel cross-continental distances on the order of weeks or months, well within their lifetimes (Figure 5—figure supplements 2 and 3; Scheidel, 2015).

The Roman Empire is particularly important in understanding how transient mobility could become a unique hallmark of this period. During the expansion of the Empire, existing and new cities quickly expanded as hubs for trade and labor. Urban-military complexes emerged along the frontier as military forces established themselves and drew in local communities which sought protection or economic benefit (Verhagen et al., 2019). To support these rapidly growing economic city-centers, human capital beyond the local population was necessary, thus drawing in people from far away places either
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Simulation of population structure with and without long-range dispersal. (A) A base model of spatial structure is established by calibrating per-generation dispersal rate to generate a maximum $F_{ST}$ of ~0.03 across the maximal spatial distance, and visualized using PCA. In addition to this base dispersal, either 4% (B) or 8% (C) of individuals disperse longer distances, and the effect is tracked by analyzing spatial $F_{ST}$ through time, as well as PCA after 120 generations of long-range dispersal.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 7:

**Figure supplement 1.** A sigmaDisp - *N* parameter pair was chosen to closely approximate the observed $F_{ST}$max of ~0.03 using grid search across a range of parameter pairs.

**Figure supplement 2.** A sigmaDispLR parameter was chosen to qualitatively resemble long-range dispersal distances observed in the data, by comparing the distribution of distances under long-range dispersal (outliers) to randomly chosen distances given the spatial distribution of samples.

as clear incentives for cross-Mediterranean and cross-continental movement (Abulafia, 2011; Beard, 2015; Broodbank, 2013; Symonds, 2017). This enabled people to travel cross-continental distances on the order of weeks or months, well within their lifetimes (Figure 5—figure supplements 2 and 3; Scheidel, 2015).

The Roman Empire is particularly important in understanding how transient mobility could become a unique hallmark of this period. During the expansion of the Empire, existing and new cities quickly expanded as hubs for trade and labor. Urban-military complexes emerged along the frontier as military forces established themselves and drew in local communities which sought protection or economic benefit (Verhagen et al., 2019). To support these rapidly growing economic city-centers, human capital beyond the local population was necessary, thus drawing in people from far away places either
freely or forcibly (e.g. slavery, military). According to a longstanding historical hypothesis, the Urban Graveyard Effect, the influx of migrants in city-centers disproportionately contributed to the death rate over the birth rate; a process which would contribute to observing individuals as 'transient' migrants (*Tacoma, 2016*). Long-range, transient migration, combined with the Roman Empire's highly efficient travel networks (*Scheidel et al., 2007; Oleson, 2009; Scheidel, 2015*) may explain the genetically heterogeneous populations, especially along the frontier regions (e.g. Serbia, Croatia, and Austria).

With transient mobility as the main contributor to the observed heterogeneity, it remains unclear what additional demographic processes contributed to the maintenance of spatial genetic structure. The collapse of the Empire involved a loss of urban-military complexes and depopulation of cities, followed by ruralization (*Burgess, 2007; Dey, 2015; Roymans et al., 2020*). Without the Empire incentivising trade and movement, there may have been little motive for individuals to remain in these suddenly remote regions.

If this hypothesis is true, we would expect a reduction in local genetic heterogeneity after the collapse of the Empire. Unfortunately, we do not have this period sampled densely enough to assess this comprehensively. The lack of samples is further amplified by the fact that ancient DNA comes from archaeological excavations, which tend to be enriched in urban areas; a stone mausoleum in the city center, for example, will produce more surface scatter than a wood farmhouse, making urban areas more likely for excavation (*Bowes, 2011*). This makes it difficult to comprehensively address differences in rural versus urban demography. Collecting more genetic data from both urban and rural contexts across the historical period will be a valuable future step in understanding how spatial population structure was maintained. Furthermore, it could elucidate the role of other historical events and peoples, such as the Franks, Lombards, Visigoths, and Huns, during the Migration Period.

Based on genetic analyses and the rich historical record, we hypothesize that both the loss of transient migrants which contributed to population heterogeneity, as well as repopulation by less heterogeneous, but temporally stable, local populations could have helped maintain overall stability of genetic structure from the Iron Age to present-day. This work highlights the utility of ancient DNA in revealing complex population dynamics through direct genetic observations through time and the importance of integrating historical contexts to understand these complexities.

# Materials and methods

## Sample collection and archaeological sites

The archaeological context for ancient individuals reported in this study is detailed in *Supplementary file 1*. Site descriptions were written by the contributing archaeologists. Descriptions of individual-level burials are included where possible.

Sampling was performed to maximize coverage across Europe and the Mediterranean, as opposed to detailed site level sampling. We particularly focused on regions where there was no published genomic data for the Imperial Roman and Late Antiquity Period at the time of sampling (e.g. France, Austria, the Balkans, Armenia, and North Africa).

Although we aimed to collect samples in the Imperial Roman and Late Antiquity period (approximately 1 CE-700 CE), some samples fall outside of this period due to limited sample availability and/or lack of date specificity at the time of sampling (prior to radiocarbon dating).

## Date determination for individuals and time periods

Determination of time periods and boundaries were influenced by (1) the wide geographic range represented in the newly reported and published data (including historical changes in those regions), (2) the amount of data available in the proposed time periods, (3) the amount of genetic change observed during a time interval, and (4) the types of temporal comparisons made in the analysis.

* Mesolithic and Neolithic 10000 BCE - 3500 BCE

* Copper Age 3500 BCE - 2300 BCE

* Bronze Age 2300 BCE - 1000 BCE

* Iron Age 1000 BCE - 1 CE

* Imperial Rome & Late Antiquity 1 CE - 700 CE

* Medieval Ages & Early Modern 700 CE - 1900 CE

* Present-day 1900 CE - onward
freely or forcibly (e.g. slavery, military). According to a longstanding historical hypothesis, the Urban Graveyard Effect, the influx of migrants in city-centers disproportionately contributed to the death rate over the birth rate; a process which would contribute to observing individuals as 'transient' migrants (*Tacoma, 2016*). Long-range, transient migration, combined with the Roman Empire's highly efficient travel networks (*Scheidel et al., 2007; Oleson, 2009; Scheidel, 2015*) may explain the genetically heterogeneous populations, especially along the frontier regions (e.g. Serbia, Croatia, and Austria).

With transient mobility as the main contributor to the observed heterogeneity, it remains unclear what additional demographic processes contributed to the maintenance of spatial genetic structure. The collapse of the Empire involved a loss of urban-military complexes and depopulation of cities, followed by ruralization (*Burgess, 2007; Dey, 2015; Roymans et al., 2020*). Without the Empire incentivising trade and movement, there may have been little motive for individuals to remain in these suddenly remote regions.

If this hypothesis is true, we would expect a reduction in local genetic heterogeneity after the collapse of the Empire. Unfortunately, we do not have this period sampled densely enough to assess this comprehensively. The lack of samples is further amplified by the fact that ancient DNA comes from archaeological excavations, which tend to be enriched in urban areas; a stone mausoleum in the city center, for example, will produce more surface scatter than a wood farmhouse, making urban areas more likely for excavation (*Bowes, 2011*). This makes it difficult to comprehensively address differences in rural versus urban demography. Collecting more genetic data from both urban and rural contexts across the historical period will be a valuable future step in understanding how spatial population structure was maintained. Furthermore, it could elucidate the role of other historical events and peoples, such as the Franks, Lombards, Visigoths, and Huns, during the Migration Period.

Based on genetic analyses and the rich historical record, we hypothesize that both the loss of transient migrants which contributed to population heterogeneity, as well as repopulation by less heterogeneous, but temporally stable, local populations could have helped maintain overall stability of genetic structure from the Iron Age to present-day. This work highlights the utility of ancient DNA in revealing complex population dynamics through direct genetic observations through time and the importance of integrating historical contexts to understand these complexities.

# Materials and methods

## Sample collection and archaeological sites

The archaeological context for ancient individuals reported in this study is detailed in *Supplementary file 1*. Site descriptions were written by the contributing archaeologists. Descriptions of individual-level burials are included where possible.

Sampling was performed to maximize coverage across Europe and the Mediterranean, as opposed to detailed site level sampling. We particularly focused on regions where there was no published genomic data for the Imperial Roman and Late Antiquity Period at the time of sampling (e.g. France, Austria, the Balkans, Armenia, and North Africa).

Although we aimed to collect samples in the Imperial Roman and Late Antiquity period (approximately 1 CE-700 CE), some samples fall outside of this period due to limited sample availability and/or lack of date specificity at the time of sampling (prior to radiocarbon dating).

## Date determination for individuals and time periods

Determination of time periods and boundaries were influenced by (1) the wide geographic range represented in the newly reported and published data (including historical changes in those regions), (2) the amount of data available in the proposed time periods, (3) the amount of genetic change observed during a time interval, and (4) the types of temporal comparisons made in the analysis.

* Mesolithic and Neolithic 10000 BCE - 3500 BCE

* Copper Age 3500 BCE - 2300 BCE

* Bronze Age 2300 BCE - 1000 BCE

* Iron Age 1000 BCE - 1 CE

* Imperial Rome & Late Antiquity 1 CE - 700 CE

* Medieval Ages & Early Modern 700 CE - 1900 CE

* Present-day 1900 CE - onward
Dates for newly reported study individuals were determined through a combination of radiocarbon dating and archaeologically inferred dates. Sample groups were created based on the finest grouping using a composite of criteria: site, projection in PCA, and archaeologically inferred date or period. We aimed to radiocarbon date at least one sample from each of these groups, summing to a total of 126 samples. For each sample, one gram of petrous bone was sent for radiocarbon dating by Accelerator Mass Spectrometry (AMS) at the Keck Carbon Cycle AMS Facility at the University of California, Irvine. Resulting dates were calibrated using https://c14.arch.ox.ac.uk/oxcal.html.

For samples not directly dated, we assigned the range of directly dated samples within their sample group (n=49), or the archaeological date where the former was not available (n=29). Note that one individual (R3477, R3476) had two radiocarbon dates since two samples (a tooth and petrous) were dated and sequenced before they were determined to be from the same individual based on genetic information. Another individual (R9818, R9823) had two radiocarbon dates because both left and right petrosals were sampled, but the sequence data for one turned out to be contaminated. In both cases of identical samples the radiocarbon date ranges were almost entirely overlapping.

In the study we use a single date estimate (for both newly reported and published samples) which is the midpoint of the 95% confidence interval when using AMS dates, and the average of the lower and upper bound inference dates when using archaeological context for dating. The dating approach used for each sample is included in files *Supplementary file 1* (archaeological context) and *Supplementary file 2* (sample metadata). The full AMS and calibration results are reported in *Supplementary file 3*.

## DNA extraction, library preparation, and sequencing

The 204 ancient genomes reported in this study, 26 of which were recently reported in *Moots et al., 2022*, represent a subset of samples screened from 53 archaeological sites across 18 countries. We isolated and finely ground the cochlear regions of the petrous bones in dedicated clean room facilities at the University of Vienna following the protocols described in *Pinhasi et al., 2019; Pinhasi et al., 2015*. Using 50 mg of bone powder, DNA was extracted by 18 hr incubation of the powder in a solution of Proteinase-K and EDTA. DNA was eluted in 50 µl 10 mM Tris-HCl, 1 mM EDTA, 0.05% Tween-20, pH 8.0 as in *Dabney et al., 2013; Rohland and Hofreiter, 2007*. 12.5–25 µL of DNA extract was used to prepare partial uracil–DNA–glycosylase (UDG) double stranded libraries as described in *Rohland et al., 2015*. After a partial (30 minute) UDG treatment, library preparation followed a modified version of the *Meyer and Kircher, 2010* protocol (*Meyer and Kircher, 2010*): the initial DNA fragmentation step was not required and MinElute PCR purification kits (Qiagen) were used for all library clean-up steps. Libraries were double-indexed with Accuprime Pfx Supermix. The PCR cycling conditions were as follows: an initial denaturation at 95 °C for 5 min followed by 12 cycles of 95 °C for 15 s, 60 °C for 30 s and 68 °C for 30 s with a final elongation at 68 °C for 5 min. After indexing, the libraries were purified using the MinElute system (Qiagen) and eluted in 25 µL of 1 mM EDTA, 0.05% Tween-20.

Libraries were screened based on Qubit concentration and visual validation of Bioanalyzer peaks for an initial low coverage (MiSeq or NextSeq) screening run.

## Processing sequence data and sample screening

Newly reported samples were initially sequenced to low coverage on MiSeq or NextSeq for screening. Following demultiplexing of the sequencing libraries, reads were trimmed, aligned, filtered for quality, and deduplicated. The following sequence data processing pipeline was applied to both screening and full sequencing runs for new data.

Adapters were removed from sequence reads using Cutadapt (v1.14) (*Martin, 2011*). Then, for each sample, reads were processed further (a) with the 2 base pairs at either end of the reads trimmed off and (b) without trimming. Since partial UDG treatment was performed on the libraries, a damage signature consisting of elevated C>T transitions on the 5' end and G>A transitions on the 3' end should remain at the ends of reads. Therefore, analyzing untrimmed, aligned reads would allow us to assess the amount of the ancient DNA damage signature present in a sample, and to use this as a criteria for authenticating that the sampled DNA is ancient. Other than the variable trimming parameter for the ends of the reads, all other parameters remained the same for both screening and high coverage sequencing data.
Dates for newly reported study individuals were determined through a combination of radiocarbon dating and archaeologically inferred dates. Sample groups were created based on the finest grouping using a composite of criteria: site, projection in PCA, and archaeologically inferred date or period. We aimed to radiocarbon date at least one sample from each of these groups, summing to a total of 126 samples. For each sample, one gram of petrous bone was sent for radiocarbon dating by Accelerator Mass Spectrometry (AMS) at the Keck Carbon Cycle AMS Facility at the University of California, Irvine. Resulting dates were calibrated using https://c14.arch.ox.ac.uk/oxcal.html.

For samples not directly dated, we assigned the range of directly dated samples within their sample group (n=49), or the archaeological date where the former was not available (n=29). Note that one individual (R3477, R3476) had two radiocarbon dates since two samples (a tooth and petrous) were dated and sequenced before they were determined to be from the same individual based on genetic information. Another individual (R9818, R9823) had two radiocarbon dates because both left and right petrosals were sampled, but the sequence data for one turned out to be contaminated. In both cases of identical samples the radiocarbon date ranges were almost entirely overlapping.

In the study we use a single date estimate (for both newly reported and published samples) which is the midpoint of the 95% confidence interval when using AMS dates, and the average of the lower and upper bound inference dates when using archaeological context for dating. The dating approach used for each sample is included in files *Supplementary file 1* (archaeological context) and *Supplementary file 2* (sample metadata). The full AMS and calibration results are reported in *Supplementary file 3*.

## DNA extraction, library preparation, and sequencing

The 204 ancient genomes reported in this study, 26 of which were recently reported in Moots et al., 2022, represent a subset of samples screened from 53 archaeological sites across 18 countries. We isolated and finely ground the cochlear regions of the petrous bones in dedicated clean room facilities at the University of Vienna following the protocols described in Pinhasi et al., 2019; Pinhasi et al., 2015. Using 50 mg of bone powder, DNA was extracted by 18 hr incubation of the powder in a solution of Proteinase-K and EDTA. DNA was eluted in 50 µl 10 mM Tris-HCl, 1 mM EDTA, 0.05% Tween-20, pH 8.0 as in Dabney et al., 2013; Rohland and Hofreiter, 2007. 12.5–25 µL of DNA extract was used to prepare partial uracil–DNA–glycosylase (UDG) double stranded libraries as described in Rohland et al., 2015. After a partial (30 minute) UDG treatment, library preparation followed a modified version of the Meyer and Kircher, 2010 protocol (Meyer and Kircher, 2010): the initial DNA fragmentation step was not required and MinElute PCR purification kits (Qiagen) were used for all library clean-up steps. Libraries were double-indexed with Accuprime Pfx Supermix. The PCR cycling conditions were as follows: an initial denaturation at 95 °C for 5 min followed by 12 cycles of 95 °C for 15 s, 60 °C for 30 s and 68 °C for 30 s with a final elongation at 68 °C for 5 min. After indexing, the libraries were purified using the MinElute system (Qiagen) and eluted in 25 µL of 1 mM EDTA, 0.05% Tween-20.

Libraries were screened based on Qubit concentration and visual validation of Bioanalyzer peaks for an initial low coverage (MiSeq or NextSeq) screening run.

## Processing sequence data and sample screening

Newly reported samples were initially sequenced to low coverage on MiSeq or NextSeq for screening. Following demultiplexing of the sequencing libraries, reads were trimmed, aligned, filtered for quality, and deduplicated. The following sequence data processing pipeline was applied to both screening and full sequencing runs for new data.

Adapters were removed from sequence reads using Cutadapt (v1.14) (Martin, 2011). Then, for each sample, reads were processed further (a) with the 2 base pairs at either end of the reads trimmed off and (b) without trimming. Since partial UDG treatment was performed on the libraries, a damage signature consisting of elevated C>T transitions on the 5' end and G>A transitions on the 3' end should remain at the ends of reads. Therefore, analyzing untrimmed, aligned reads would allow us to assess the amount of the ancient DNA damage signature present in a sample, and to use this as a criteria for authenticating that the sampled DNA is ancient. Other than the variable trimming parameter for the ends of the reads, all other parameters remained the same for both screening and high coverage sequencing data.
et al., 2014; Omrak et al., 2016; O'Sullivan et al., 2018; Patterson et al., 2012; Prüfer et al., 2017; Rivollat et al., 2020; Rodríguez-Varela et al., 2017; Saag et al., 2019; Saag et al., 2017; Sánchez-Quinto et al., 2019; Schiffels et al., 2016; Schroeder et al., 2019; Schuenemann et al., 2017; Sikora et al., 2017; Skoglund et al., 2014; Skourtanioti et al., 2020; Unterländer et al., 2017; Valdiosera et al., 2018; van den Brink et al., 2017; Veeramah et al., 2018; Villalba-Mouco et al., 2019; Wang et al., 2019; Zalloua et al., 2018). We also included relevant genetic data made available by authors that were not in the AADR: present-day genomes from the Balkans (Kovacevic et al., 2014), present-day genomes from 4 Poles, 3 Germans, and 2 Moldavians (Pagani et al., 2016), and Bronze Age Italian genomes (Saupe et al., 2021). Pseudohaploid genotypes for published Bronze Age Aegean genomes (Clemente et al., 2021) and Bronze Age Serbian genomes (Žegarac et al., 2021) were generated from BAM files using our pipeline. All published genomes were filtered for contamination based on reported contamination levels in the original study and SNP coverage based on the genomic data. All published samples that contributed to this study are listed in Supplementary file 4. To ensure maximum overlap with present-day and ancient samples in analyses, the merged dataset was subset to SNPs in the Human Origin Panel array, resulting in a total of 481,259 SNPs. For PCA and qpAdm modeling, SNPs that are transitions at CpG sites (n=76,678) were excluded since they may have arisen from DNA damage as opposed to true genetic variation.

## Principal component analysis (PCA)

### Setting up the principal component analysis

Principal component analysis was performed on genotypes from present-day and Mediterranean individuals using smartpca v16000 (https://github.com/chrchang/eigensoft/blob/master/POPGEN/README). The following parameters were used: 5 outlier iterations (numoutlieriter), 10 principal components along which to remove outliers (numoutlierevec), altnormstyle set to NO, with least squares projection turned on (lsqproject set to YES). To calculate principal components only using present-day individuals, a file (poplistname) was provided with the population names of present-day individuals, randomly subsampled per population. After outlier removal (which removed 55 samples), 829 individuals and 480,712 SNPs were used in the initial analysis. All individuals (non 'reference' present-day genomes, and all of the ancient individuals) whose population was not listed in the poplist-name file were projected onto the calculated principal components. In the paper, we refer to the individuals used in the calculation of principal components as belonging to the 'reference PCA space'. These 'reference' genomes were used to calculate the PCs because (1) they represent a wide range of present-day variation and (2) the genotypes tend to be of high quality.

### Visual representation of PCA

In the figures, present-day genomes used in the reference space are generally colored gray in order to illustrate the background space of genetic variation. To reduce visual clutter and emphasize the ancient genomes, these present-day 'reference' genomes are typically unlabeled. Labels for these populations are shown in **Figure 2—figure supplement 1**.

### Calculation of F<sub>ST</sub>

To assess the extent of genetic differentiation across geographic space within a time period, we calculated the Fixation index (F<sub>ST</sub>) between groups of individuals on a sliding spatial grid. Each grid cell measured 10 degrees longitude by 10 degrees latitude, and was slid by 1 degree in both directions (north and east) nine times to build a total of 10 spatial grids. For each of these grids, pairwise F<sub>ST</sub> was calculated between all populated 10-by-10 grid cells using Hudson's estimator, correcting for unequal sample size (Bhatia et al., 2013). In addition to F<sub>ST</sub>, we also calculated the average geographic distance (in kilometers) between all individuals across pairs of grid cells to assess how spatial distance relates to genetic differentiation. To visualize this relationship, we used lowess smoothing as implemented in python's statsmodels package (statsmodels api.nonparametric.lowess, v. 0.12.2). To infer confidence intervals for the lowess smoothing estimates, we devised a spatial bootstrapping procedure. Our bootstrap approach samples pairs of grid cells in a way that always samples all overlapping cells or none of them, so individuals are either fully included in a bootstrap replicate or not at all. This prevents double-counting individuals since they contribute to several comparisons across space due to the sliding grid.
et al., 2014; Omrak et al., 2016; O'Sullivan et al., 2018; Patterson et al., 2012; Prüfer et al., 2017; Rivollat et al., 2020; Rodríguez-Varela et al., 2017; Saag et al., 2019; Saag et al., 2017; Sánchez-Quinto et al., 2019; Schiffels et al., 2016; Schroeder et al., 2019; Schuenemann et al., 2017; Sikora et al., 2017; Skoglund et al., 2014; Skourtanioti et al., 2020; Unterländer et al., 2017; Valdiosera et al., 2018; van den Brink et al., 2017; Veeramah et al., 2018; Villalba-Mouco et al., 2019; Wang et al., 2019; Zalloua et al., 2018). We also included relevant genetic data made available by authors that were not in the AADR: present-day genomes from the Balkans (Kovacevic et al., 2014), present-day genomes from 4 Poles, 3 Germans, and 2 Moldavians (Pagani et al., 2016), and Bronze Age Italian genomes (Saupe et al., 2021). Pseudohaploid genotypes for published Bronze Age Aegean genomes (Clemente et al., 2021) and Bronze Age Serbian genomes (Žegarac et al., 2021) were generated from BAM files using our pipeline. All published genomes were filtered for contamination based on reported contamination levels in the original study and SNP coverage based on the genomic data. All published samples that contributed to this study are listed in Supplementary file 4. To ensure maximum overlap with present-day and ancient samples in analyses, the merged dataset was subset to SNPs in the Human Origin Panel array, resulting in a total of 481,259 SNPs. For PCA and qpAdm modeling, SNPs that are transitions at CpG sites (n=76,678) were excluded since they may have arisen from DNA damage as opposed to true genetic variation.

## Principal component analysis (PCA)

### Setting up the principal component analysis

Principal component analysis was performed on genotypes from present-day and Mediterranean individuals using smartpca v16000 (https://github.com/chrchang/eigensoft/blob/master/POPGEN/README). The following parameters were used: 5 outlier iterations (numoutlieriter), 10 principal components along which to remove outliers (numoutlierevec), altnormstyle set to NO, with least squares projection turned on (lsqproject set to YES). To calculate principal components only using present-day individuals, a file (poplistname) was provided with the population names of present-day individuals, randomly subsampled per population. After outlier removal (which removed 55 samples), 829 individuals and 480,712 SNPs were used in the initial analysis. All individuals (non 'reference' present-day genomes, and all of the ancient individuals) whose population was not listed in the poplist-name file were projected onto the calculated principal components. In the paper, we refer to the individuals used in the calculation of principal components as belonging to the 'reference PCA space'. These 'reference' genomes were used to calculate the PCs because (1) they represent a wide range of present-day variation and (2) the genotypes tend to be of high quality.

### Visual representation of PCA

In the figures, present-day genomes used in the reference space are generally colored gray in order to illustrate the background space of genetic variation. To reduce visual clutter and emphasize the ancient genomes, these present-day 'reference' genomes are typically unlabeled. Labels for these populations are shown in **Figure 2—figure supplement 1**.

### Calculation of F<sub>ST</sub>

To assess the extent of genetic differentiation across geographic space within a time period, we calculated the Fixation index (F<sub>ST</sub>) between groups of individuals on a sliding spatial grid. Each grid cell measured 10 degrees longitude by 10 degrees latitude, and was slid by 1 degree in both directions (north and east) nine times to build a total of 10 spatial grids. For each of these grids, pairwise F<sub>ST</sub> was calculated between all populated 10-by-10 grid cells using Hudson's estimator, correcting for unequal sample size (Bhatia et al., 2013). In addition to F<sub>ST</sub>, we also calculated the average geographic distance (in kilometers) between all individuals across pairs of grid cells to assess how spatial distance relates to genetic differentiation. To visualize this relationship, we used lowess smoothing as implemented in python's statsmodels package (statsmodels api.nonparametric.lowess, v. 0.12.2). To infer confidence intervals for the lowess smoothing estimates, we devised a spatial bootstrapping procedure. Our bootstrap approach samples pairs of grid cells in a way that always samples all overlapping cells or none of them, so individuals are either fully included in a bootstrap replicate or not at all. This prevents double-counting individuals since they contribute to several comparisons across space due to the sliding grid.
## Modeling ancestry and identifying outliers using *qpAdm*

We used the *qpAdm* tool of admixtools 2.0 to build a workflow that:

1. Identifies similar individuals within a region

2. Groups these into regional clusters

3. Compares these clusters both within and across regions

In this workflow, we heavily utilize what we call one-component models, where we test whether two individuals or clusters of individuals form a clade relative to a chosen set of reference populations (which we define below). To clarify what we mean by one-component model, assume that we have two focal individuals *i1* and *i2*, as well as a set of reference populations. Using the terminology of admixtools 2 as well as (Harney et al., 2021), the following four tests are equivalent with respect to the resulting p-value:

* `qpWave(left = c(i1, i2), right = ref)`

* `qpAdm(left = c(i1, i2), target = NULL, right = ref)`

* `qpAdm(left = i1, target = i2, right = ref)`

* `qpAdm(left = i2, target = i1, right = ref)`

where we use the *R*-style notation of *c(i1, i2)* to denote a vector consisting of *i1* and *i2*. In all four cases, we are testing against the null hypothesis that the two individuals *do* form a clade, with low p-values indicating a rejection of that hypothesis. We will call any implementation of this test a 'one-component *qpAdm* model', but the equivalence stated above shows that our one-component models are equivalent to what was called a 'qpWave analysis' by *Fernandes et al., 2020*.

For our reference populations, we chose a set similar to those previously used to model Eurasian historical genomes (*Fernandes et al., 2020*). However, we added two Asian populations (Laos_Hoabinhian and Onge) based on evidence of gene flow from further east in a subset of our data. The purpose of a set of reference populations in the *qpAdm* modeling setting is to represent components of ancestry which are differentially related to the focal individuals being tested (i.e. 'left', or 'sources' and 'target') in order to resolve differences in ancestry, but distally related enough to minimize the chance of recent gene flow between 'left' and 'right'. Our final set of references is:

Mbuti.DG (n=4), WHG (n=8), Russia_Ust_Ishim.DG (n=1), CHG (n=2), EHG (n=3), Iran_GanjDareh_N (n=8), Israel_Natufian_published (n=3), Jordan_PPNB (n=6), Laos_Hoabinhian (n=1), Russia_EBA_Yamnaya_Samara (n=9), Onge (n=6), Spain_EIMiron (n=1), Turkey_N_published (n=8), Russia_MA1_HG (n=1), Morocco_Iberomaurusian (n=6), Czech_Vestonice16 (n=1).

## Individual-based one-component models within regions

In the first step of our workflow, we perform one-component *qpAdm* tests between all pairs of individuals from the same region, from the Copper Age (inclusive) up to present-day (exclusive). As mentioned above, this approach tests against the null hypothesis that the two individuals *do* form a clade, with low p-values indicating a rejection of that hypothesis. Low *qpAdm* p-values thus suggest that the test individuals are not more closely related to each other than to one or more populations in the reference set. To convert the *qpAdm* p-value into a measure of dissimilarity (*d*), we calculate *d* = -log₁₀(*p*-value), where a large value of *d* indicates a low p-value, and thus a rejection of the null hypothesis of the two individuals forming a clade.

To cluster individuals into groups of genetically similar individuals, we performed hierarchical clustering on the dissimilarity matrix constructed from all pairwise values of *d* within a region. Hierarchical clustering was performed using the UPGMA algorithm as implemented in python's `scipy.cluster.hierarch` (v. 1.6.1). The hierarchical clustering was then split into flat clusters using a dissimilarity cutoff of 1.3, which corresponds to a nominal p-value cutoff of 0.05. Intermediate results from the pairwise analysis and clustering are shown for each region in Appendix 1.

## Identifying ancestry outliers and their potential sources

Once clusters within regions were identified as described above, we classified them into two groups based on size. Clusters consisting of less than 5% of the total population or no more than two individuals in the region (across time) were classified as outlier candidates, whose ancestry is under-represented in the region they were sampled. All other clusters were classified as majority clusters.
# Modeling ancestry and identifying outliers using *qpAdm*

We used the *qpAdm* tool of admixtools 2.0 to build a workflow that:

1. Identifies similar individuals within a region

2. Groups these into regional clusters

3. Compares these clusters both within and across regions

In this workflow, we heavily utilize what we call one-component models, where we test whether two individuals or clusters of individuals form a clade relative to a chosen set of reference populations (which we define below). To clarify what we mean by one-component model, assume that we have two focal individuals *i1* and *i2*, as well as a set of reference populations. Using the terminology of admixtools 2 as well as (Harney et al., 2021), the following four tests are equivalent with respect to the resulting p-value:

* `qpWave(left = c(i1, i2), right = ref)`

* `qpAdm(left = c(i1, i2), target = NULL, right = ref)`

* `qpAdm(left = i1, target = i2, right = ref)`

* `qpAdm(left = i2, target = i1, right = ref)`

where we use the *R*-style notation of *c(i1, i2)* to denote a vector consisting of *i1* and *i2*. In all four cases, we are testing against the null hypothesis that the two individuals *do* form a clade, with low p-values indicating a rejection of that hypothesis. We will call any implementation of this test a 'one-component *qpAdm* model', but the equivalence stated above shows that our one-component models are equivalent to what was called a 'qpWave analysis' by *Fernandes et al., 2020*.

For our reference populations, we chose a set similar to those previously used to model Eurasian historical genomes (*Fernandes et al., 2020*). However, we added two Asian populations (Laos_Hoabinhian and Onge) based on evidence of gene flow from further east in a subset of our data. The purpose of a set of reference populations in the *qpAdm* modeling setting is to represent components of ancestry which are differentially related to the focal individuals being tested (i.e. 'left', or 'sources' and 'target') in order to resolve differences in ancestry, but distally related enough to minimize the chance of recent gene flow between 'left' and 'right'. Our final set of references is:

Mbuti.DG (n=4), WHG (n=8), Russia_Ust_Ishim.DG (n=1), CHG (n=2), EHG (n=3), Iran_GanjDareh_N (n=8), Israel_Natufian_pUBLISHED (n=3), Jordan_PPNB (n=6), Laos_Hoabinhian (n=1), Russia_EBA_Yamnaya_Samara (n=9), Onge (n=6), Spain_EIMiron (n=1), Turkey_N_pUBLISHED (n=8), Russia_MA1_HG (n=1), Morocco_Iberomaurusian (n=6), Czech_Vestonice16 (n=1).

# Individual-based one-component models within regions

In the first step of our workflow, we perform one-component *qpAdm* tests between all pairs of individuals from the same region, from the Copper Age (inclusive) up to present-day (exclusive). As mentioned above, this approach tests against the null hypothesis that the two individuals *do* form a clade, with low p-values indicating a rejection of that hypothesis. Low *qpAdm* p-values thus suggest that the test individuals are not more closely related to each other than to one or more populations in the reference set. To convert the *qpAdm* p-value into a measure of dissimilarity (*d*), we calculate *d* = -log₁₀(*p*-value), where a large value of *d* indicates a low p-value, and thus a rejection of the null hypothesis of the two individuals forming a clade.

To cluster individuals into groups of genetically similar individuals, we performed hierarchical clustering on the dissimilarity matrix constructed from all pairwise values of *d* within a region. Hierarchical clustering was performed using the UPGMA algorithm as implemented in python's *scipy.cluster.hierarch* (v. 1.6.1). The hierarchical clustering was then split into flat clusters using a dissimilarity cutoff of 1.3, which corresponds to a nominal p-value cutoff of 0.05. Intermediate results from the pairwise analysis and clustering are shown for each region in Appendix 1.

# Identifying ancestry outliers and their potential sources

Once clusters within regions were identified as described above, we classified them into two groups based on size. Clusters consisting of less than 5% of the total population or no more than two individuals in the region (across time) were classified as outlier candidates, whose ancestry is under-represented in the region they were sampled. All other clusters were classified as majority clusters.
Following this classification of clusters, we then split each cluster by time period (Copper Age, Bronze Age, Iron Age, Imperial Rome & Late Antiquity, Middle Ages and Early Modern), to end up with a region_period_cluster sub-classification. All downstream analyses were done using these region_period clusters.

For each cluster identified as an outlier candidate, we then tested all one-component models involving the candidate and a majority cluster within the same region. This was done to ensure that outlier candidates were truly distinct from majority ancestries identified in the region. If an outlier candidate could be connected to a majority cluster within its region through a valid one-component model (above a p-value threshold of 0.01), it was removed from the outlier candidate list, as it represented a majority ancestry within the region.

For all outlier candidates that remained, we aimed to identify potential source ancestries that were majority ancestries in other regions. To do so, we tested all one-component models for a given outlier candidate and each of the majority clusters across all other regions. A majority cluster from a different region that had a valid one-component model with an outlier candidate was considered a potential source. To find the most likely source(s), we then subjected these potential source populations to model competition.

Model competition involves re-testing the model fit of each potential source after adding another potential source to the right group, first described in *Lazaridis et al., 2016*, and more recently detailed in *Harney et al., 2021*. The idea is that if a population in the right group has significantly more allele sharing with the target than the source, then the model will be rejected. (This is why right group populations are chosen to be distal, yet relevant, to target and source). We use this property to our advantage by rotating all n-1 valid sources for the target through the right group set, one at a time, for the same target and a source x (the one source that is not included in the rotation). If the previously valid source x is rejected when including another valid source y in the right group then we remove source x from the list of potential sources. Note that this does not make source y the best source, only a better one than x. Thus, this scheme only eliminates sub-optimal sources, rather than selecting a best source.

If there were still multiple valid sources following the model competition scheme described above, we prioritized a candidate source that is from the same time period as the target cluster, or from a previous period going backwards in time. If there were still multiple candidate sources, they were considered equivalent and all kept in downstream analyses.

Among the outliers identified, we did not find a significant sex bias compared to non-outliers. Overall, there are more males than females in the dataset. However, the proportions of males in non-outliers, outliers with source, and outliers without source do not differ significantly by a Chi-squared test (-value = 0.4117, df = 2; Figure 5—figure supplement 1). When outliers (with and without source) are treated as one group, there is still no significant association with outlier status and sex (p-value = 0.633, df = 1).

## Admixture modeling with *qpAdm*

For targeted analyses of other clusters beyond just outlier candidates, for example to annotate *Figures 2–4*, we also used cluster-based *qpAdm*. In addition to one-component models as described above, we also used two-component models of admixture. These models test the hypothesis that a focal target cluster can be modeled as a two-way admixture of two sources (or ‘left’ populations). As above, a p-value below the threshold rejects this hypothesis, that is the proposed admixture model is not a good fit and a different model needs to be considered to disentangle the admixture scenario in question.

## Simulations

To assess how spatial population structure would be impacted by different modes of dispersal, we set up forward simulations in continuous 2D space using SLiM v. 3.6 (*Haller and Messer, 2019*). The aim of these simulations was to approximate the extent of spatial population structure we observe by the beginning of the Iron Age in Western Eurasia, after the major prehistoric migrations had taken place. To achieve that, we decided not to attempt simulating the precise ancestry composition of populations in different regions at that time, but rather to simulate simply the extent of spatial structure as measured by the relationship of population differentiation (*F*<sub>ST</sub>) and geographic distance. We chose the SLiM simulation framework to make use of its extensive feature set to simulate individuals
Following this classification of clusters, we then split each cluster by time period (Copper Age, Bronze Age, Iron Age, Imperial Rome & Late Antiquity, Middle Ages and Early Modern), to end up with a region_period_cluster sub-classification. All downstream analyses were done using these region_period clusters.

For each cluster identified as an outlier candidate, we then tested all one-component models involving the candidate and a majority cluster within the same region. This was done to ensure that outlier candidates were truly distinct from majority ancestries identified in the region. If an outlier candidate could be connected to a majority cluster within its region through a valid one-component model (above a p-value threshold of 0.01), it was removed from the outlier candidate list, as it represented a majority ancestry within the region.

For all outlier candidates that remained, we aimed to identify potential source ancestries that were majority ancestries in other regions. To do so, we tested all one-component models for a given outlier candidate and each of the majority clusters across all other regions. A majority cluster from a different region that had a valid one-component model with an outlier candidate was considered a potential source. To find the most likely source(s), we then subjected these potential source populations to model competition.

Model competition involves re-testing the model fit of each potential source after adding another potential source to the right group, first described in Lazaridis et al., 2016, and more recently detailed in Harney et al., 2021. The idea is that if a population in the right group has significantly more allele sharing with the target than the source, then the model will be rejected. (This is why right group populations are chosen to be distal, yet relevant, to target and source). We use this property to our advantage by rotating all n-1 valid sources for the target through the right group set, one at a time, for the same target and a source x (the one source that is not included in the rotation). If the previously valid source x is rejected when including another valid source y in the right group then we remove source x from the list of potential sources. Note that this does not make source y the best source, only a better one than x. Thus, this scheme only eliminates sub-optimal sources, rather than selecting a best source.

If there were still multiple valid sources following the model competition scheme described above, we prioritized a candidate source that is from the same time period as the target cluster, or from a previous period going backwards in time. If there were still multiple candidate sources, they were considered equivalent and all kept in downstream analyses.

Among the outliers identified, we did not find a significant sex bias compared to non-outliers. Overall, there are more males than females in the dataset. However, the proportions of males in non-outliers, outliers with source, and outliers without source do not differ significantly by a Chi-squared test (-value = 0.4117, df = 2; Figure 5—figure supplement 1). When outliers (with and without source) are treated as one group, there is still no significant association with outlier status and sex (p-value = 0.633, df = 1).

## Admixture modeling with *qpAdm*

For targeted analyses of other clusters beyond just outlier candidates, for example to annotate Figures 2–4, we also used cluster-based *qpAdm*. In addition to one-component models as described above, we also used two-component models of admixture. These models test the hypothesis that a focal target cluster can be modeled as a two-way admixture of two sources (or ‘left’ populations). As above, a p-value below the threshold rejects this hypothesis, that is the proposed admixture model is not a good fit and a different model needs to be considered to disentangle the admixture scenario in question.

## Simulations

To assess how spatial population structure would be impacted by different modes of dispersal, we set up forward simulations in continuous 2D space using SLiM v. 3.6 (Haller and Messer, 2019). The aim of these simulations was to approximate the extent of spatial population structure we observe by the beginning of the Iron Age in Western Eurasia, after the major prehistoric migrations had taken place. To achieve that, we decided not to attempt simulating the precise ancestry composition of populations in different regions at that time, but rather to simulate simply the extent of spatial structure as measured by the relationship of population differentiation ($F_{ST}$) and geographic distance. We chose the SLiM simulation framework to make use of its extensive feature set to simulate individuals
in continuous space. We simulated diploid genomes made up of a single, 10⁸ bp long chromosome, with recombination rate and mutation rate set to 10⁻⁸. We used the default Wright-Fisher simulation mode, where a single population of constant size *N* is simulated with non-overlapping populations, that is each generation is made up of offspring generated from the previous generation. Spatial structure is established by associating each individual with a continuous 2D coordinate (i.e. latitude and longitude), and by using these coordinates to govern three demographic processes: mate choice, competition, and dispersal. An overview of how these processes can be set up to interact in SLiM can be found in Recipe 15.4 of SLiM v. 3.6 (see e.g. here: https://github.com/MesserLab/SLiM/tree/v3.6/SLiMgui/Recipes). Briefly, for mate choice, a Gaussian interaction function with *maxDistance* = 0.1, *maxStrength* = 1.0, *sigma* = 0.02 is used to govern a *mateChoice* callback using the *strength* of that interaction function. For competition, another Gaussian interaction function with *maxDistance* = 0.3, *maxStrength* = 3.0, *sigma* = 0.1 is used to calculate competition using the *totalNeighbor-Strength* vector of that interaction function to scale an individual's relative fitness as 1.1 - *competition* / *N*. Finally, we establish local dispersal through a *modifyChild* callback, where a newly generated offspring's position is drawn from a Gaussian centered at the location of the maternal individual with standard deviation *sigmaDisp*.

We let this population evolve forward in time for 2000+120 generations during which the processes outlined above lead to spatial population structure, where individuals sampled closely together in 2D space are also more closely related genetically. We do not simulate mutations in SLiM, as this poses a major computational burden. Instead, we use tree sequence recording (Haller et al., 2019; Kelleher et al., 2018) to track the full genealogy of all individuals in the simulation which are either alive at the end of the simulation, or explicitly sampled through time using the *treeSeqRememberIndividuals* function of SLiM. While 2000 generations are enough to establish spatial structure under the parameters we consider, it is by far insufficient for all sampled individuals to fully coalesce. To accurately assess neutral variation however, we need all sampled individuals to have a common ancestor at some point in the past, as mutations may have arisen at any point leading back to this ultimate coalescence event. Therefore, we approximate the deep history of our population with a panmictic population simulated backwards in time using the coalescent with recombination as implemented in msprime (Kelleher et al., 2016). This process has been referred to as 'recapitation' (Haller et al., 2019), where an incomplete genealogy with multiple roots (from SLiM) is 'recapitated' using coalescent simulation backwards in time. This is made possible by using the tree sequence data structure to record and simulate genealogies in both SLiM and msprime.

Since our simulation is only concerned with how processes such as dispersal affect neutral variation across space and through time, we can use the 'recapitated' tree sequence to overlay mutations onto the full genealogy of all sampled individuals, also using msprime. The rationale here is that under neutrality, mutations will not affect the structure of the genealogy, so we can simulate the genealogy without mutations first, and overlay neutral mutations second, thereby greatly reducing computational burden.

We then extracted the resulting genotypes of all individuals from the tree sequence for downstream analysis. We only kept sites segregating in individuals at the end of the simulation ('present day'), and filtered for minor allele frequency of at least 0.01 across the entire dataset, to make downstream analysis of simulated genomes comparable to how the empirical data was ascertained and analyzed.

To assess the relationship between $F_{ST}$ and spatial distance, we split geographic space into a 10-by-10 grid and calculated all pairwise $F_{ST}$ between inhabited grid cells using Hudson's estimator with unequal sample size correction (Bhatia et al., 2013), as well as the average geographic (euclidean) distance between individuals across grid cells. We used this $F_{ST}$ analysis to calibrate the base dispersal *sigmaDisp* as well as the population size *N*, so that $F_{ST}$ at maximum distance ($F_{STmax}$) would approximately match the $F_{STmax}$ we observed at the start of the historical period (~0.03). We used grid search with a range of *sigmaDisp* and *N* values, and found the parameter pair $N=50,000$ & *sigmaDisp* = 0.02 to qualitatively produce the closest match (Figure 7—figure supplement 1). We use this parameter set as our base model of population structure without long-range dispersal, where we allow spatial structure to establish over 2000 generations, and then observe the $F_{ST}$ - Distance relationship over the following 120 generations for a total of 2120 generations simulated in SLiM (Figure 7A).
in continuous space. We simulated diploid genomes made up of a single, 10⁸ bp long chromosome, with recombination rate and mutation rate set to 10⁻⁸. We used the default Wright-Fisher simulation mode, where a single population of constant size *N* is simulated with non-overlapping populations, that is each generation is made up of offspring generated from the previous generation. Spatial structure is established by associating each individual with a continuous 2D coordinate (i.e. latitude and longitude), and by using these coordinates to govern three demographic processes: mate choice, competition, and dispersal. An overview of how these processes can be set up to interact in SLiM can be found in Recipe 15.4 of SLiM v. 3.6 (see e.g. here: https://github.com/MesserLab/SLiM/tree/v3.6/SLiMgui/Recipes). Briefly, for mate choice, a Gaussian interaction function with *maxDistance* = 0.1, *maxStrength* = 1.0, *sigma* = 0.02 is used to govern a *mateChoice* callback using the *strength* of that interaction function. For competition, another Gaussian interaction function with *maxDistance* = 0.3, *maxStrength* = 3.0, *sigma* = 0.1 is used to calculate competition using the *totalNeighbor-Strength* vector of that interaction function to scale an individual's relative fitness as 1.1 - *competition* / *N*. Finally, we establish local dispersal through a *modifyChild* callback, where a newly generated offspring's position is drawn from a Gaussian centered at the location of the maternal individual with standard deviation *sigmaDisp*.

We let this population evolve forward in time for 2000+120 generations during which the processes outlined above lead to spatial population structure, where individuals sampled closely together in 2D space are also more closely related genetically. We do not simulate mutations in SLiM, as this poses a major computational burden. Instead, we use tree sequence recording (Haller et al., 2019; Kelleher et al., 2018) to track the full genealogy of all individuals in the simulation which are either alive at the end of the simulation, or explicitly sampled through time using the *treeSeqRememberIndividuals* function of SLiM. While 2000 generations are enough to establish spatial structure under the parameters we consider, it is by far insufficient for all sampled individuals to fully coalesce. To accurately assess neutral variation however, we need all sampled individuals to have a common ancestor at some point in the past, as mutations may have arisen at any point leading back to this ultimate coalescence event. Therefore, we approximate the deep history of our population with a panmictic population simulated backwards in time using the coalescent with recombination as implemented in msprime (Kelleher et al., 2016). This process has been referred to as 'recapitation' (Haller et al., 2019), where an incomplete genealogy with multiple roots (from SLiM) is 'recapitated' using coalescent simulation backwards in time. This is made possible by using the tree sequence data structure to record and simulate genealogies in both SLiM and msprime.

Since our simulation is only concerned with how processes such as dispersal affect neutral variation across space and through time, we can use the 'recapitated' tree sequence to overlay mutations onto the full genealogy of all sampled individuals, also using msprime. The rationale here is that under neutrality, mutations will not affect the structure of the genealogy, so we can simulate the genealogy without mutations first, and overlay neutral mutations second, thereby greatly reducing computational burden.

We then extracted the resulting genotypes of all individuals from the tree sequence for downstream analysis. We only kept sites segregating in individuals at the end of the simulation ('present day'), and filtered for minor allele frequency of at least 0.01 across the entire dataset, to make downstream analysis of simulated genomes comparable to how the empirical data was ascertained and analyzed.

To assess the relationship between $F_{ST}$ and spatial distance, we split geographic space into a 10-by-10 grid and calculated all pairwise $F_{ST}$ between inhabited grid cells using Hudson's estimator with unequal sample size correction (Bhatia et al., 2013), as well as the average geographic (euclidean) distance between individuals across grid cells. We used this $F_{ST}$ analysis to calibrate the base dispersal *sigmaDisp* as well as the population size *N*, so that $F_{ST}$ at maximum distance ($F_{STmax}$) would approximately match the $F_{STmax}$ we observed at the start of the historical period (~0.03). We used grid search with a range of *sigmaDisp* and *N* values, and found the parameter pair $N=50,000$ & *sigmaDisp* = 0.02 to qualitatively produce the closest match (Figure 7—figure supplement 1). We use this parameter set as our base model of population structure without long-range dispersal, where we allow spatial structure to establish over 2000 generations, and then observe the $F_{ST}$ - Distance relationship over the following 120 generations for a total of 2120 generations simulated in SLiM (Figure 7A).
Given this base model of spatial population structure, we can now start to introduce long-range dispersing individuals. We do this by allowing a specified fraction of individuals to use a higher sigmaD-dispLR than the sigmaDisp used by the rest of the population for the final 120 generations of the simulation, approximately matching the 3000 years since the beginning of the historical period assuming a generation time of 25 years. Since the long-range dispersal is also drawn from a Gaussian, the distribution of dispersal distances will have substantial overlap with the distances produced by base dispersal. To make the fraction of long-range dispersal accurately represent the fraction of individuals that actually disperse longer distances, we thus require a long-range dispersing individual to disperse to a location outside of the 99th percentile of density covered by the short-range base dispersal.

We aimed to choose a sigmaDispLR that approximately matches the empirical distribution of long-range dispersing individuals we observe in the analysis displayed in **Figure 4**. Since the euclidean distances in the simulation are on a different scale than the geodesic distances observed in the data, we aimed to match qualitatively the relationship of long-range dispersal distances to random distances that could be observed if two populated locations were drawn at random. We visually analyzed this relationship from the data, and then performed a search across a range of possible sigmaDistLR to find a qualitative match. This led us to choose a value of sigmaDistLR = 0.20 (Figure 7—figure supplement 2).

Finally, we analyzed simulated 'present-day' genomes (i.e. after 2120 generations of SLiM) using PCA. We used the *sklearn.decomposition.PCA* module (scikit-learn v. 0.24.2) with the `svdreiplier` option to run non-probabilistic PCA to calculate the first 10 principal components. Similarly to how the empirical data was analyzed with *smartpca*, we also did 5 rounds of iterative outlier removal, removing individuals from the PCA that deviated by more than six standard deviations along any of the 10 principal components. The number of variants contributing to these PCA were 624,617, 625,669, and 626,052 for **Figure 7A**, **B** and **C** respectively, and thus comparable to the number of variants contributing to our data analysis.

Appendix 1. Visualization of per-region qqAdm clustering.

## Data availability

All sequence data newly generated for this study are available at the European Nucleotide Archive (ENA) database. Raw sequencing data is available under the accession number PRJEB53565. Sequences mapped to the human reference genome are available under accession number PRJEB53564. Sequences previously reported by Moots et al., 2022 are available under accession number PRJEB49419. All published data used in this study is listed in *Supplementary file 4*, and can be retrieved from primary sources (see Materials and methods, 'Combining new genotypes with ancient and present-day published data') or from the Allen Ancient Data Resource (Mallick et al., 2023).

## Ethics

This study follows ethics guidelines adopted by the ancient DNA field (Alpaslan-Roodenberg et al., 2021). A clear plan of research was laid out before the collection of samples, leading us to focus on sampling from under-sampled historical regions in Eurasia and therefore minimize unnecessary destruction of human remains. The research intent for these samples was clearly communicated to caretakers of the samples prior to collection. Local anthropologists, archaeologists, and museum directors from each geographic region were involved in the sample acquisition, extraction from skeletal material, and interpretation of genetic results. The genetic findings regarding individual samples from each region were communicated to local collaborators, all of whom were included as co-authors on the paper and were supportive of the final results. The involvement of our local collaborators was essential for the interpretation of the genetic results through their input on the historical and archaeological characterization of the specimens. We have supported our local collaborators with immediate access to the raw genetic data, and by communicating results in written and oral forums. Authorities responsible for all archaeological sites provided written documentation for their specimens to be included in this study through collaboration with the Pinasi Lab (Vienna, Austria).

## Acknowledgements
Given this base model of spatial population structure, we can now start to introduce long-range dispersing individuals. We do this by allowing a specified fraction of individuals to use a higher sigmaD-dispLR than the sigmaDisp used by the rest of the population for the final 120 generations of the simulation, approximately matching the 3000 years since the beginning of the historical period assuming a generation time of 25 years. Since the long-range dispersal is also drawn from a Gaussian, the distribution of dispersal distances will have substantial overlap with the distances produced by base dispersal. To make the fraction of long-range dispersal accurately represent the fraction of individuals that actually disperse longer distances, we thus require a long-range dispersing individual to disperse to a location outside of the 99th percentile of density covered by the short-range base dispersal.

We aimed to choose a sigmaDispLR that approximately matches the empirical distribution of long-range dispersing individuals we observe in the analysis displayed in **Figure 4**. Since the euclidean distances in the simulation are on a different scale than the geodesic distances observed in the data, we aimed to match qualitatively the relationship of long-range dispersal distances to random distances that could be observed if two populated locations were drawn at random. We visually analyzed this relationship from the data, and then performed a search across a range of possible sigmaDistLR to find a qualitative match. This led us to choose a value of sigmaDistLR = 0.20 (Figure 7—figure supplement 2).

Finally, we analyzed simulated 'present-day' genomes (i.e. after 2120 generations of SLiM) using PCA. We used the *sklearn.decomposition.PCA* module (scikit-learn v. 0.24.2) with the `svdreiplier` option to run non-probabilistic PCA to calculate the first 10 principal components. Similarly to how the empirical data was analyzed with *smartpca*, we also did 5 rounds of iterative outlier removal, removing individuals from the PCA that deviated by more than six standard deviations along any of the 10 principal components. The number of variants contributing to these PCA were 624,617, 625,669, and 626,052 for **Figure 7A**, **B** and **C** respectively, and thus comparable to the number of variants contributing to our data analysis.

Appendix 1. Visualization of per-region qqAdm clustering.

## Data availability

All sequence data newly generated for this study are available at the European Nucleotide Archive (ENA) database. Raw sequencing data is available under the accession number PRJEB53565. Sequences mapped to the human reference genome are available under accession number PRJEB53564. Sequences previously reported by Moots et al., 2022 are available under accession number PRJEB49419. All published data used in this study is listed in *Supplementary file 4*, and can be retrieved from primary sources (see Materials and methods, 'Combining new genotypes with ancient and present-day published data') or from the Allen Ancient Data Resource (Mallick et al., 2023).

## Ethics

This study follows ethics guidelines adopted by the ancient DNA field (Alpaslan-Roodenberg et al., 2021). A clear plan of research was laid out before the collection of samples, leading us to focus on sampling from under-sampled historical regions in Eurasia and therefore minimize unnecessary destruction of human remains. The research intent for these samples was clearly communicated to caretakers of the samples prior to collection. Local anthropologists, archaeologists, and museum directors from each geographic region were involved in the sample acquisition, extraction from skeletal material, and interpretation of genetic results. The genetic findings regarding individual samples from each region were communicated to local collaborators, all of whom were included as co-authors on the paper and were supportive of the final results. The involvement of our local collaborators was essential for the interpretation of the genetic results through their input on the historical and archaeological characterization of the specimens. We have supported our local collaborators with immediate access to the raw genetic data, and by communicating results in written and oral forums. Authorities responsible for all archaeological sites provided written documentation for their specimens to be included in this study through collaboration with the Pinasi Lab (Vienna, Austria).

## Acknowledgements
We thank Professor Walter Scheidel for helpful discussions and feedback on the historical context, and all members of the Pritchard and Pinhasi labs for their valuable input. We thank Pieter W Faber and the University of Chicago Genomics Facility for sequencing the samples reported here. We thank Benjamin Peter and an anonymous reviewer for their insightful and constructive reviews. This project was partially supported by a National Science Foundation Graduate Research Fellowship (MLA.), a grant from the National Institutes of Health RO1 HG011432 (CLW), the Austrian Science Fund (FWF) M3108-G (SS), and the Howard Hughes Medical Institute (JKP).

## Additional information

### Competing interests

Ziyue Gao: Reviewing editor, eLife. Maja Bilić: Maja Bilić is affiliated with Palisada d.o.o. The author has no financial interests to declare. Josip Burmaz: Josip Burmaz is affiliated with Kaducej Ltd. The author has no financial interests to declare. Dženi Los: Dženi Los is affiliated with Kaducej Ltd. The author has no financial interests to declare. Rene Masaryk: Rene Masaryk is affiliated with Skupina STIK. The author has no financial interests to declare. The other authors declare that no competing interests exist.

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>National Science Foundation</td><td>Graduate Research Fellowship</td><td>Margaret L Antonio</td></tr><tr><td>National Institutes of Health</td><td>HG011432</td><td>Clemens L Weiss</td></tr><tr><td>Howard Hughes Medical Institute</td><td></td><td>Jonathan K Pritchard</td></tr><tr><td>Austrian Science Fund</td><td>M3108-G</td><td>Susanna Sawyer</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

Margaret L Antonio, Clemens L Weiss, Data curation, Formal analysis, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing; Ziyue Gao, Conceptualization, Formal analysis, Writing – original draft, Project administration, Writing – review and editing; Susanna Sawyer, Victoria Oberreiter, Olivia Cheronet, Brina Zagorc, Elisa Praxmarer, Kadir Toykan Özdoğan, Lea Demetz, Resources, Data curation; Hannah M Moots, Resources, Data curation, Writing – review and editing; Jeffrey P Spence, Data curation, Methodology; Pere Gelabert, Daniel Fernandes, Michaela Lucci, Timka Alihodžić, Selma Amrani, Pavel Avetisyan, Christèle Baillif-Ducros, Željka Bedić, Audrey Bertrand, Maja Bilić, Luca Bondioli, Paulina Borówka, Emmanuel Botte, Josip Burmaz, Domagoj Bužanić, Francesca Candilio, Mirna Cvetko, Daniela De Angelis, Ivan Drnić, Kristián Elschek, Mounir Fantar, Andrej Gaspari, Gabriella Gasperetti, Francesco Genchi, Snežana Golubović, Zuzana Huke ľová, Rimantas Jankauskas, Kristina Jelinčić Vučković, Gordana Jeremić, Iva Kaić, Kevin Kazek, Hamazasp Khachatryan, Anahit Khudaverdyan, Sylvia Kirchengast, Miomir Korać, Valérie Kozlowski, Mária Krošláková, Dora Kušan Špalj, Francesco La Pastina, Marie Laguardia, Sandra Legrand, Tino Leleković, Tamara Leskovar, Wiesław Lorkiewicz, Dženi Los, Ana Maria Silva, Rene Masaryk, Vinka Matijević, Yahia Mehdi Seddik Cherifi, Nicholas Meyer, Ilija Mikić, Nataša Miladinović-Radmilović, Branka Milošević Zakić, Lina Nacouzi, Magdalena Natuniewicz-Sekuła, Alessia Nava, Christine Neugebauer-Maresch, Jan Nováček, Anna Osterholtz, Julianne Paige, Lujana Paraman, Dominique Pieri, Karol Pieta, Stefan Pop-Lazić, Matej Ruttkay, Mirjana Sanader, Arkadiusz Sołtysiak, Alessandra Sperduti, Tijana Stan-kovic Pesterac, Maria Teschler-Nicola, Iwona Teul, Domagoj Tončinić, Julien Trapp, Dragana Vulović, Tomasz Waliszewski, Diethard Walter, Milos Zivanovic, Mohamed el Mostefa Filah, Morana Čaušević-Bully, Mario Šlaus, Dušan Borić, Mario Novak, Resources; Alfredo Coppa, Resources, Supervision;
We thank Professor Walter Scheidel for helpful discussions and feedback on the historical context, and all members of the Pritchard and Pinhasi labs for their valuable input. We thank Pieter W Faber and the University of Chicago Genomics Facility for sequencing the samples reported here. We thank Benjamin Peter and an anonymous reviewer for their insightful and constructive reviews. This project was partially supported by a National Science Foundation Graduate Research Fellowship (MLA.), a grant from the National Institutes of Health RO1 HG011432 (CLW), the Austrian Science Fund (FWF) M3108-G (SS), and the Howard Hughes Medical Institute (JKP).

## Additional information

### Competing interests

Ziyue Gao: Reviewing editor, eLife. Maja Bilić: Maja Bilić is affiliated with Palisada d.o.o. The author has no financial interests to declare. Josip Burmaz: Josip Burmaz is affiliated with Kaducej Ltd. The author has no financial interests to declare. Dženi Los: Dženi Los is affiliated with Kaducej Ltd. The author has no financial interests to declare. Rene Masaryk: Rene Masaryk is affiliated with Skupina STIK. The author has no financial interests to declare. The other authors declare that no competing interests exist.

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>National Science Foundation</td><td>Graduate Research Fellowship</td><td>Margaret L Antonio</td></tr><tr><td>National Institutes of Health</td><td>HG011432</td><td>Clemens L Weiss</td></tr><tr><td>Howard Hughes Medical Institute</td><td></td><td>Jonathan K Pritchard</td></tr><tr><td>Austrian Science Fund</td><td>M3108-G</td><td>Susanna Sawyer</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

Margaret L Antonio, Clemens L Weiss, Data curation, Formal analysis, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing; Ziyue Gao, Conceptualization, Formal analysis, Writing – original draft, Project administration, Writing – review and editing; Susanna Sawyer, Victoria Oberreiter, Olivia Cheronet, Brina Zagorc, Elisa Praxmarer, Kadir Toykan Özdoğan, Lea Demetz, Resources, Data curation; Hannah M Moots, Resources, Data curation, Writing – review and editing; Jeffrey P Spence, Data curation, Methodology; Pere Gelabert, Daniel Fernandes, Michaela Lucci, Timka Alihodžić, Selma Amrani, Pavel Avetisyan, Christèle Baillif-Ducros, Željka Bedić, Audrey Bertrand, Maja Bilić, Luca Bondioli, Paulina Borówka, Emmanuel Botte, Josip Burmaz, Domagoj Bužanić, Francesca Candilio, Mirna Cvetko, Daniela De Angelis, Ivan Drnić, Kristián Elschek, Mounir Fantar, Andrej Gaspari, Gabriella Gasperetti, Francesco Genchi, Snežana Golubović, Zuzana Huke ľová, Rimantas Jankauskas, Kristina Jelinčić Vučković, Gordana Jeremić, Iva Kaić, Kevin Kazek, Hamazasp Khachatryan, Anahit Khudaverdyan, Sylvia Kirchengast, Miomir Korać, Valérie Kozlowski, Mária Krošláková, Dora Kušan Špalj, Francesco La Pastina, Marie Laguardia, Sandra Legrand, Tino Leleković, Tamara Leskovar, Wiesław Lorkiewicz, Dženi Los, Ana Maria Silva, Rene Masaryk, Vinka Matijević, Yahia Mehdi Seddik Cherifi, Nicholas Meyer, Ilija Mikić, Nataša Miladinović-Radmilović, Branka Milošević Zakić, Lina Nacouzi, Magdalena Natuniewicz-Sekuła, Alessia Nava, Christine Neugebauer-Maresch, Jan Nováček, Anna Osterholtz, Julianne Paige, Lujana Paraman, Dominique Pieri, Karol Pieta, Stefan Pop-Lazić, Matej Ruttkay, Mirjana Sanader, Arkadiusz Sołtysiak, Alessandra Sperduti, Tijana Stan-kovic Pesterac, Maria Teschler-Nicola, Iwona Teul, Domagoj Tončinić, Julien Trapp, Dragana Vulović, Tomasz Waliszewski, Diethard Walter, Milos Zivanovic, Mohamed el Mostefa Filah, Morana Čaušević-Bully, Mario Šlaus, Dušan Borić, Mario Novak, Resources; Alfredo Coppa, Resources, Supervision;
Ron Pinhasi, Supervision, Funding acquisition; Jonathan K Pritchard, Conceptualization, Supervision, Funding acquisition, Methodology, Writing – review and editing

## Author ORCIDs

Margaret L Antonio https://orcid.org/0000-0003-1049-211X

Clemens L Weiss https://orcid.org/0000-0003-3321-3902

Ziyue Gao https://orcid.org/0000-0001-9244-0238

Victoria Oberreiter https://orcid.org/0000-0003-0766-3782

Jeffrey P Spence https://orcid.org/0000-0002-3199-1447

Brina Zagorc https://orcid.org/0000-0002-7685-1958

Christèle Baillif-Ducros https://orcid.org/0000-0003-3050-3082

Francesca Candilio https://orcid.org/0000-0002-4668-1361

Francesco Genchi https://orcid.org/0000-0002-7696-4207

Kristina Jelinčić Vučković https://orcid.org/0000-0002-1236-734X

Ana Maria Silva https://orcid.org/0000-0002-1912-6581

Matej Ruttkay https://orcid.org/0000-0002-6441-9914

Arkadiusz Sołtysiak https://orcid.org/0000-0002-9040-5022

Alessandra Sperduti https://orcid.org/0000-0001-9338-5891

Mario Šlaus https://orcid.org/0000-0002-4941-2212

Mario Novak https://orcid.org/0000-0002-4567-8742

Alfredo Coppa https://orcid.org/0000-0002-7708-2484

Jonathan K Pritchard https://orcid.org/0000-0002-8828-5236

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.79714.sa1

Author response https://doi.org/10.7554/eLife.79714.sa2

## Additional files

### Supplementary files

* Supplementary file 1. Archaeological context for sampling locations.

* Supplementary file 2. Metadata for all newly reported individuals **Supplementary file 3**.

* Supplementary file 3. AMS and calibration results.

* Supplementary file 4. Published samples that contributed to this study.

* MDAR checklist

### Data availability

All sequence data newly generated for this study are available at the European Nucleotide Archive (ENA) database. Raw sequencing data is available under the accession number PRJEB53565. Sequences mapped to the human reference genome are available under accession number PRJEB53564. Sequences previously reported by Moots et al. are available under accession number PRJEB49419.

The following datasets were generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>This Study</td><td>2022</td><td>Historical Period Genomes from West Eurasia and the Mediterranean - Raw reads</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB53565</td><td>EBI, PRJEB53565</td></tr><tr><td>This Study</td><td>2022</td><td>Historical Period Genomes from West Eurasia and the Mediterranean - Mapped reads</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB53564</td><td>EBI, PRJEB53564</td></tr></tbody></table>
Ron Pinhasi, Supervision, Funding acquisition; Jonathan K Pritchard, Conceptualization, Supervision, Funding acquisition, Methodology, Writing – review and editing

## Author ORCIDs

Margaret L Antonio https://orcid.org/0000-0003-1049-211X

Clemens L Weiss https://orcid.org/0000-0003-3321-3902

Ziyue Gao https://orcid.org/0000-0001-9244-0238

Victoria Oberreiter https://orcid.org/0000-0003-0766-3782

Jeffrey P Spence https://orcid.org/0000-0002-3199-1447

Brina Zagorc https://orcid.org/0000-0002-7685-1958

Christèle Baillif-Ducros https://orcid.org/0000-0003-3050-3082

Francesca Candilio https://orcid.org/0000-0002-4668-1361

Francesco Genchi https://orcid.org/0000-0002-7696-4207

Kristina Jelinčić Vučković https://orcid.org/0000-0002-1236-734X

Ana Maria Silva https://orcid.org/0000-0002-1912-6581

Matej Ruttkay https://orcid.org/0000-0002-6441-9914

Arkadiusz Sołtysiak https://orcid.org/0000-0002-9040-5022

Alessandra Sperduti https://orcid.org/0000-0001-9338-5891

Mario Šlaus https://orcid.org/0000-0002-4941-2212

Mario Novak https://orcid.org/0000-0002-4567-8742

Alfredo Coppa https://orcid.org/0000-0002-7708-2484

Jonathan K Pritchard https://orcid.org/0000-0002-8828-5236

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.79714.sa1

Author response https://doi.org/10.7554/eLife.79714.sa2

## Additional files

### Supplementary files

* Supplementary file 1. Archaeological context for sampling locations.

* Supplementary file 2. Metadata for all newly reported individuals **Supplementary file 3**.

* Supplementary file 3. AMS and calibration results.

* Supplementary file 4. Published samples that contributed to this study.

* MDAR checklist

### Data availability

All sequence data newly generated for this study are available at the European Nucleotide Archive (ENA) database. Raw sequencing data is available under the accession number PRJEB53565. Sequences mapped to the human reference genome are available under accession number PRJEB53564. Sequences previously reported by Moots et al. are available under accession number PRJEB49419.

The following datasets were generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>This Study</td><td>2022</td><td>Historical Period Genomes from West Eurasia and the Mediterranean - Raw reads</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB53565</td><td>EBI, PRJEB53565</td></tr><tr><td>This Study</td><td>2022</td><td>Historical Period Genomes from West Eurasia and the Mediterranean - Mapped reads</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB53564</td><td>EBI, PRJEB53564</td></tr></tbody></table>
The following previously published datasets were used:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>Moots et al.</td><td>2022</td><td>Mobility in the Iron Age Central Mediterranean</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB49419</td><td>EBI, PRJEB49419</td></tr><tr><td>Mallick et al.</td><td>2021</td><td>Allen Ancient DNA Resource (AADR) v44.3</td><td>https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data</td><td>Harvard Dataverse, 10.7910/DVN/FFIDCW</td></tr></tbody></table>

## References

* Abulafia D. 2011. *The Great Sea: A Human History of the Mediterranean* Oxford University Press. DOI: https://doi.org/10.1093/oso/9780195323344.001.0001

* Agranat-Tamir L, Waldman S, Martin MAS, Gokhman D, Mishol N, Eshel T, Cheronet O, Rohland N, Mallick S, Adamski N, Lawson AM, Mah M, Michel M, Oppenheimer J, Stewardson K, Candilio F, Keating D, Gamarra B, Tzur S, Novak M, et al. 2020. The genomic history of the bronze age southern levant. *Cell* **181**:1146–1157. DOI: https://doi.org/10.1016/j.cell.2020.04.024, PMID: 32470400

* Allen Ancient DNA Resource. 2021. Allen ancient DNA resource (AADR): downloadable genotypes of present-day and ancient DNA data. https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data

* Allentoft ME, Sikora M, Sjögren KG, Rasmussen S, Rasmussen M, Stenderup J, Damgaard PB, Schroeder H, Ahlström T, Vinner L, Malaspinas AS, Margaryan A, Higham T, Chivall D, Lynnerup N, Harvig L, Baron J, Della Casa P, Dabrowski P, Duffy PR, et al. 2015. Population genomics of bronze age Eurasia. *Nature* **522**:167–172. DOI: https://doi.org/10.1038/nature14507, PMID: 26062507

* Alpaslan-Roodenberg S, Anthony D, Babiker H, Bánffy E, Booth T, Capone P, Deshpande-Mukherjee A, Eisenmann S, Fehren-Schmitz L, Frachetti M, Fujita R, Frieman CJ, Fu Q, Gibbon V, Haak W, Hajdinjak M, Hofmann KP, Holguin B, Inomata T, Kanzawa-Kiriyama H, et al. 2021. Ethics of DNA research on human remains: five globally applicable guidelines. *Nature* **599**:41–46. DOI: https://doi.org/10.1038/s41586-021-04008-x, PMID: 34671160

* Amorim CEG, Vai S, Posth C, Modi A, Koncz I, Hakenbeck S, La Rocca MC, Mende B, Bobo D, Pohl W, Baricco LP, Bedini E, Francalacci P, Giostra C, Vida T, Winger D, von Freeden U, Ghirotto S, Lari M, Barbujani G, et al. 2018. Understanding 6th-century barbarian social organization and migration through paleogenomics. *Nature Communications* **9**:3547. DOI: https://doi.org/10.1038/s41467-018-06024-4, PMID: 30206220

* Antonio ML, Gao Z, Moots HM, Lucci M, Candilio F, Sawyer S, Oberreiter V, Calderon D, Devitofranceschi K, Aikens RC, Aneli S, Bartoli F, Bedini A, Cheronet O, Cotter DJ, Fernandes DM, Gasperetti G, Grifoni R, Guidi A, La Pastina F, et al. 2019. Ancient Rome: a genetic crossroads of Europe and the Mediterranean. *Science* **366**:708–714. DOI: https://doi.org/10.1126/science.aay6826, PMID: 31699931

* Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR, 1000 Genomes Project Consortium. 2015. A global reference for human genetic variation. *Nature* **526**:68–74. DOI: https://doi.org/10.1038/nature15393, PMID: 26432245

* Beard M. 2015. SPQR: A History of Ancient Rome New York: Liveright Publishing Corporation, a Division of W.W. Norton & Company.

* Bergström A, McCarthy SA, Hui R, Almarri MA, Ayub Q, Danecek P, Chen Y, Felkel S, Hallast P, Kamm J, Blanché H, Deleuze JF, Cann H, Mallick S, Reich D, Sandhu MS, Skoglund P, Scally A, Xue Y, Durbin R, et al. 2020. Insights into human genetic variation and population history from 929 diverse genomes. *Science* **367**:eaay5012. DOI: https://doi.org/10.1126/science.aay5012, PMID: 32193295

* Bhatia G, Patterson N, Sankararaman S, Price AL. 2013. Estimating and interpreting FST: the impact of rare variants. *Genome Research* **23**:1514–1521. DOI: https://doi.org/10.1101/gr.154831.113, PMID: 23861382

* Biagini SA, Solé-Morata N, Matisoo-Smith E, Zalloua P, Comas D, Calafell F. 2019. People from Ibiza: an unexpected isolate in the Western Mediterranean. *European Journal of Human Genetics* **27**:941–951. DOI: https://doi.org/10.1038/s41431-019-0361-1, PMID: 30765884

* Bowes K. 2011. Rural Poverty in the Roman Empire. https://bpb-us-w2.wpmucdn.com/u.osu.edu/dist/a/49661/files/2017/08/bowes-rural-poor-243zk6w.pdf [Accessed February 22, 2021].

* Brace S, Diekmann Y, Booth TJ, van Dorp L, Faltyskova Z, Rohland N, Mallick S, Olalde I, Ferry M, Michel M, Oppenheimer J, Broomandkhoshbacht N, Stewardson K, Martiniano R, Walsh S, Kayser M, Charlton S, Hellenthal G, Armit I, Schulting R, et al. 2019. Ancient genomes indicate population replacement in Early Neolithic Britain. *Nature Ecology & Evolution* **3**:765–771. DOI: https://doi.org/10.1038/s41559-019-0871-9, PMID: 30988490

* Broodbank C. 2013. *The making of the Middle Sea: a history of the Mediterranean from the beginning to the emergence of the classical world* Oxford; New York: Oxford University Press.
The following previously published datasets were used:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>Moots et al.</td><td>2022</td><td>Mobility in the Iron Age Central Mediterranean</td><td>https://www.ebi.ac.uk/ena/browser/view/PRJEB49419</td><td>EBI, PRJEB49419</td></tr><tr><td>Mallick et al.</td><td>2021</td><td>Allen Ancient DNA Resource (AADR) v44.3</td><td>https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data</td><td>Harvard Dataverse, 10.7910/DVN/FFIDCW</td></tr></tbody></table>

## References

* Abulafia D. 2011. *The Great Sea: A Human History of the Mediterranean* Oxford University Press. DOI: https://doi.org/10.1093/oso/9780195323344.001.0001

* Agranat-Tamir L, Waldman S, Martin MAS, Gokhman D, Mishol N, Eshel T, Cheronet O, Rohland N, Mallick S, Adamski N, Lawson AM, Mah M, Michel M, Oppenheimer J, Stewardson K, Candilio F, Keating D, Gamarra B, Tzur S, Novak M, et al. 2020. The genomic history of the bronze age southern levant. *Cell* **181**:1146–1157. DOI: https://doi.org/10.1016/j.cell.2020.04.024, PMID: 32470400

* Allen Ancient DNA Resource. 2021. Allen ancient DNA resource (AADR): downloadable genotypes of present-day and ancient DNA data. https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data

* Allentoft ME, Sikora M, Sjögren KG, Rasmussen S, Rasmussen M, Stenderup J, Damgaard PB, Schroeder H, Ahlström T, Vinner L, Malaspinas AS, Margaryan A, Higham T, Chivall D, Lynnerup N, Harvig L, Baron J, Della Casa P, Dabrowski P, Duffy PR, et al. 2015. Population genomics of bronze age Eurasia. *Nature* **522**:167–172. DOI: https://doi.org/10.1038/nature14507, PMID: 26062507

* Alpaslan-Roodenberg S, Anthony D, Babiker H, Bánffy E, Booth T, Capone P, Deshpande-Mukherjee A, Eisenmann S, Fehren-Schmitz L, Frachetti M, Fujita R, Frieman CJ, Fu Q, Gibbon V, Haak W, Hajdinjak M, Hofmann KP, Holguin B, Inomata T, Kanzawa-Kiriyama H, et al. 2021. Ethics of DNA research on human remains: five globally applicable guidelines. *Nature* **599**:41–46. DOI: https://doi.org/10.1038/s41586-021-04008-x, PMID: 34671160

* Amorim CEG, Vai S, Posth C, Modi A, Koncz I, Hakenbeck S, La Rocca MC, Mende B, Bobo D, Pohl W, Baricco LP, Bedini E, Francalacci P, Giostra C, Vida T, Winger D, von Freeden U, Ghirotto S, Lari M, Barbujani G, et al. 2018. Understanding 6th-century barbarian social organization and migration through paleogenomics. *Nature Communications* **9**:3547. DOI: https://doi.org/10.1038/s41467-018-06024-4, PMID: 30206220

* Antonio ML, Gao Z, Moots HM, Lucci M, Candilio F, Sawyer S, Oberreiter V, Calderon D, Devitofranceschi K, Aikens RC, Aneli S, Bartoli F, Bedini A, Cheronet O, Cotter DJ, Fernandes DM, Gasperetti G, Grifoni R, Guidi A, La Pastina F, et al. 2019. Ancient Rome: a genetic crossroads of Europe and the Mediterranean. *Science* **366**:708–714. DOI: https://doi.org/10.1126/science.aay6826, PMID: 31699931

* Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR, 1000 Genomes Project Consortium. 2015. A global reference for human genetic variation. *Nature* **526**:68–74. DOI: https://doi.org/10.1038/nature15393, PMID: 26432245

* Beard M. 2015. SPQR: A History of Ancient Rome New York: Liveright Publishing Corporation, a Division of W.W. Norton & Company.

* Bergström A, McCarthy SA, Hui R, Almarri MA, Ayub Q, Danecek P, Chen Y, Felkel S, Hallast P, Kamm J, Blanché H, Deleuze JF, Cann H, Mallick S, Reich D, Sandhu MS, Skoglund P, Scally A, Xue Y, Durbin R, et al. 2020. Insights into human genetic variation and population history from 929 diverse genomes. *Science* **367**:eaay5012. DOI: https://doi.org/10.1126/science.aay5012, PMID: 32193295

* Bhatia G, Patterson N, Sankararaman S, Price AL. 2013. Estimating and interpreting FST: the impact of rare variants. *Genome Research* **23**:1514–1521. DOI: https://doi.org/10.1101/gr.154831.113, PMID: 23861382

* Biagini SA, Solé-Morata N, Matisoo-Smith E, Zalloua P, Comas D, Calafell F. 2019. People from Ibiza: an unexpected isolate in the Western Mediterranean. *European Journal of Human Genetics* **27**:941–951. DOI: https://doi.org/10.1038/s41431-019-0361-1, PMID: 30765884

* Bowes K. 2011. Rural Poverty in the Roman Empire. https://bpb-us-w2.wpmucdn.com/u.osu.edu/dist/a/49661/files/2017/08/bowes-rural-poor-243zk6w.pdf [Accessed February 22, 2021].

* Brace S, Diekmann Y, Booth TJ, van Dorp L, Faltyskova Z, Rohland N, Mallick S, Olalde I, Ferry M, Michel M, Oppenheimer J, Broomandkhoshbacht N, Stewardson K, Martiniano R, Walsh S, Kayser M, Charlton S, Hellenthal G, Armit I, Schulting R, et al. 2019. Ancient genomes indicate population replacement in Early Neolithic Britain. *Nature Ecology & Evolution* **3**:765–771. DOI: https://doi.org/10.1038/s41559-019-0871-9, PMID: 30988490

* Broodbank C. 2013. *The making of the Middle Sea: a history of the Mediterranean from the beginning to the emergence of the classical world* Oxford ; New York: Oxford University Press.
Broushaki F, Thomas MG, Link V, López S, van Dorp L, Kirsanow K, Hofmanová Z, Diekmann Y, Cassidy LM, Díez-Del-Molino D, Kousathanas A, Sell C, Robson HK, Martiniano R, Blöcher J, Scheu A, Kreutzer S, Bollongino R, Bobo D, Davudi H, et al. 2016. Early neolithic genomes from the eastern fertile crescent. *Science* **353**:499-503. DOI: https://doi.org/10.1126/science.aaf7943, PMID: 27417496

Brunel S, Bennett EA, Cardin L, Garraud D, Barrand Emam H, Beyler A, Boulestin B, Chenal F, Ciesielski E, Convertini F, Dedet B, Desbrosse-Degobertiere S, Desenne S, Dubouloz J, Duday H, Escalon G, Fabre V, Gailledrat E, Gandelin M, Gleize Y, et al. 2020. Ancient genomes from present-day France unveil 7,000 years of its demographic history. *PNAS* **117**:12791-12798. DOI: https://doi.org/10.1073/pnas.1918034117, PMID: 32457149

Burgess RW. 2007. *The fall of Rome and the end of civilization*, by Bryan Ward-Perkins. *Canadian Journal of History* **42**:83-85. DOI: https://doi.org/10.3138/cjh.42.1.83

Cassidy LM, Martiniano R, Murphy EM, Teasdale MD, Mallory J, Hartwell B, Bradley DG. 2016. Neolithic and Bronze Age migration to Ireland and establishment of the insular Atlantic genome. *PNAS* **113**:368-373. DOI: https://doi.org/10.1073/pnas.1518445113, PMID: 26712024

Cassidy LM, Maoldúin RÓ, Kador T, Lynch A, Jones C, Woodman PC, Murphy E, Ramsey G, Dowd M, Noonan A, Campbell C, Jones ER, Mattiangeli V, Bradley DG. 2020. A dynastic elite in monumental Neolithic society. *Nature* **582**:384-388. DOI: https://doi.org/10.1038/s41586-020-2378-6, PMID: 32555485

Clemente F, Unterländer M, Dolgova O, Amorim CEG, Coroado-Santos F, Neuenschwander S, Ganiatsou E, Cruz Dávalos DI, Anchieri L, Michaud F, Winkelbach L, Blöcher J, Arizmendi Cárdenas YO, Sousa da Mota B, Kalliga E, Souleles A, Kontopoulos I, Karamitrou-Mentessidi G, Philaniotou O, Sampson A, et al. 2021. The genomic history of the Aegean palatial civilizations. *Cell* **184**:2565-2586. DOI: https://doi.org/10.1016/j.cell.2021.03.039, PMID: 33930288

Dabney J, Knapp M, Glocke I, Gansauge MT, Weihmann A, Nickel B, Valdiosera C, García N, Pääbo S, Arsuaga JL, Meyer M. 2013. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. *PNAS* **110**:15758-15763. DOI: https://doi.org/10.1073/pnas.1314445110, PMID: 24019490

Damgaard PB, Marchi N, Rasmussen S, Peyrot M, Renaud G, Korneliussen T, Moreno-Mayar JV, Pedersen MW, Goldberg A, Usmanova E, Baimukhanov N, Loman V, Hedeager L, Pedersen AG, Nielsen K, Afanasiev G, Akmatov K, Aldashev A, Alpaslan A, Baimbetov G, et al. 2018. 137 ancient human genomes from across the Eurasian steppes. *Nature* **557**:369-374. DOI: https://doi.org/10.1038/s41586-018-0094-2, PMID: 29743675

de Barros Damgaard P, Martiniano R, Kamm J, Moreno-Mayar JV, Kroonen G, Peyrot M, Barjamovic G, Rasmussen S, Zacho C, Baimukhanov N, Zaibert V, Merz V, Biddanda A, Merz I, Loman V, Evdokimov V, Usmanova E, Hemphill B, Seguin-Orlando A, Yediyah FE, et al. 2018. The first horse herders and the impact of early bronze age steppe expansions into Asia. *Science* **360**:eaar7711. DOI: https://doi.org/10.1126/science.aar7711, PMID: 29743352

Dey HW. 2015. *Introduction: Urban Living and the "Fall" of the Roman Empire. The Afterlife of the Roman CIT* Cambridge: Cambridge University Press. DOI: https://doi.org/10.1017/CBO9781107706538

Ebenesersdóttir SS, Sandoval-Velasco M, Gunnarsdóttir ED, Jagadeesan A, Guðmundsdóttir VB, Thordardóttir EL, Einarsdóttir MS, Moore KHS, Sigurðsson Á, Magnúsdóttir DN, Jónsson H, Snorradóttir S, Hovig E, Møller P, Kockum I, Olsson T, Alfredsson L, Hansen TF, Werge T, Cavalleri GL, et al. 2018. Ancient genomes from Iceland reveal the making of a human population. *Science* **360**:1028-1032. DOI: https://doi.org/10.1126/science.aar2625, PMID: 29853688

Feldman M, Fernández-Domínguez E, Reynolds L, Baird D, Pearson J, Hershkovitz I, May H, Goring-Morris N, Benz M, Gresky J, Bianco RA, Fairbairn A, Mustafaoğlu G, Stockhammer PW, Posth C, Haak W, Jeong C, Krause J. 2019a. Late Pleistocene human genome suggests a local origin for the first farmers of Central Anatolia. *Nature Communications* **10**:1218. DOI: https://doi.org/10.1038/s41467-019-09209-7, PMID: 30890703

Feldman M, Master DM, Bianco RA, Burri M, Stockhammer PW, Mittnik A, Aja AJ, Jeong C, Krause J. 2019b. Ancient DNA sheds light on the genetic origins of early Iron Age Philistines. *Science Advances* **5**:eaax0061. DOI: https://doi.org/10.1126/sciadv.aax0061, PMID: 31281897

Fernandes DM, Strapagiel D, Borówka P, Marciniak B, Żądzińska E, Sirak K, Siska V, Grygiel R, Carlsson J, Manica A, Lorkiewicz W, Pinhasi R. 2018. A genomic Neolithic time transect of hunter-farmer admixture in Central Poland. *Scientific Reports* **8**:14879. DOI: https://doi.org/10.1038/s41598-018-33067-w, PMID: 30291256

Fernandes DM, Mittnik A, Olalde I, Lazaridis I, Cheronet O, Rohland N, Mallick S, Bernardos R, Broomandkhoshbacht N, Carlsson M, Culleton BJ, Ferry M, Gamarra B, Lari M, Mah M, Michel M, Modi A, Novak M, Oppenheimer J, Sirak KA, et al. 2020. Author correction: the spread of steppe and Iranian-related ancestry in the islands of the Western Mediterranean. *Nature Ecology & Evolution* **4**:764. DOI: https://doi.org/10.1038/s41559-020-1197-3, PMID: 32296139

Fregel R, Méndez FL, Bokbot Y, Martín-Socas D, Camalich-Massieu MD, Santana J, Morales J, Ávila-Arcos MC, Underhill PA, Shapiro B, Wojcik G, Rasmussen M, Soares AER, Kapp J, Sockell A, Rodríguez-Santos FJ, Mikdad A, Trujillo-Mederos A, Bustamante CD. 2018. Ancient genomes from North Africa evidence prehistoric migrations to the Maghreb from both the Levant and Europe. *PNAS* **115**:6774-6779. DOI: https://doi.org/10.1073/pnas.1800851115, PMID: 29895688

Fu Q, Posth C, Hajdinjak M, Petr M, Mallick S, Fernandes D, Furtwängler A, Haak W, Meyer M, Mittnik A, Nickel B, Peltzer A, Rohland N, Slon V, Talamo S, Lazaridis I, Lipson M, Mathieson I, Schiffels S, Skoglund P,
Broushaki F, Thomas MG, Link V, López S, van Dorp L, Kirsanow K, Hofmanová Z, Diekmann Y, Cassidy LM, Díez-Del-Molino D, Kousathanas A, Sell C, Robson HK, Martiniano R, Blöcher J, Scheu A, Kreutzer S, Bollongino R, Bobo D, Davudi H, et al. 2016. Early neolithic genomes from the eastern fertile crescent. *Science* **353**:499-503. DOI: https://doi.org/10.1126/science.aaf7943, PMID: 27417496

Brunel S, Bennett EA, Cardin L, Garraud D, Barrand Emam H, Beyler A, Boulestin B, Chenal F, Ciesielski E, Convertini F, Dedet B, Desbrosse-Degobertiere S, Desenne S, Dubouloz J, Duday H, Escalon G, Fabre V, Gailledrat E, Gandelin M, Gleize Y, et al. 2020. Ancient genomes from present-day France unveil 7,000 years of its demographic history. *PNAS* **117**:12791-12798. DOI: https://doi.org/10.1073/pnas.1918034117, PMID: 32457149

Burgess RW. 2007. *The fall of Rome and the end of civilization*, by Bryan Ward-Perkins. *Canadian Journal of History* **42**:83-85. DOI: https://doi.org/10.3138/cjh.42.1.83

Cassidy LM, Martiniano R, Murphy EM, Teasdale MD, Mallory J, Hartwell B, Bradley DG. 2016. Neolithic and Bronze Age migration to Ireland and establishment of the insular Atlantic genome. *PNAS* **113**:368-373. DOI: https://doi.org/10.1073/pnas.1518445113, PMID: 26712024

Cassidy LM, Maoldúin RÓ, Kador T, Lynch A, Jones C, Woodman PC, Murphy E, Ramsey G, Dowd M, Noonan A, Campbell C, Jones ER, Mattiangeli V, Bradley DG. 2020. A dynastic elite in monumental Neolithic society. *Nature* **582**:384-388. DOI: https://doi.org/10.1038/s41586-020-2378-6, PMID: 32555485

Clemente F, Unterländer M, Dolgova O, Amorim CEG, Coroado-Santos F, Neuenschwander S, Ganiatsou E, Cruz Dávalos DI, Anchieri L, Michaud F, Winkelbach L, Blöcher J, Arizmendi Cárdenas YO, Sousa da Mota B, Kalliga E, Souleles A, Kontopoulos I, Karamitrou-Mentessidi G, Philaniotou O, Sampson A, et al. 2021. The genomic history of the Aegean palatial civilizations. *Cell* **184**:2565-2586. DOI: https://doi.org/10.1016/j.cell.2021.03.039, PMID: 33930288

Dabney J, Knapp M, Glocke I, Gansauge MT, Weihmann A, Nickel B, Valdiosera C, García N, Pääbo S, Arsuaga JL, Meyer M. 2013. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. *PNAS* **110**:15758-15763. DOI: https://doi.org/10.1073/pnas.1314445110, PMID: 24019490

Damgaard PB, Marchi N, Rasmussen S, Peyrot M, Renaud G, Korneliussen T, Moreno-Mayar JV, Pedersen MW, Goldberg A, Usmanova E, Baimukhanov N, Loman V, Hedeager L, Pedersen AG, Nielsen K, Afanasiev G, Akmatov K, Aldashev A, Alpaslan A, Baimbetov G, et al. 2018. 137 ancient human genomes from across the Eurasian steppes. *Nature* **557**:369-374. DOI: https://doi.org/10.1038/s41586-018-0094-2, PMID: 29743675

de Barros Damgaard P, Martiniano R, Kamm J, Moreno-Mayar JV, Kroonen G, Peyrot M, Barjamovic G, Rasmussen S, Zacho C, Baimukhanov N, Zaibert V, Merz V, Biddanda A, Merz I, Loman V, Evdokimov V, Usmanova E, Hemphill B, Seguin-Orlando A, Yediyah FE, et al. 2018. The first horse herders and the impact of early bronze age steppe expansions into Asia. *Science* **360**:eaar7711. DOI: https://doi.org/10.1126/science.aar7711, PMID: 29743352

Dey HW. 2015. *Introduction: Urban Living and the "Fall" of the Roman Empire. The Afterlife of the Roman CIT* Cambridge: Cambridge University Press. DOI: https://doi.org/10.1017/CBO9781107706538

Ebenesersdóttir SS, Sandoval-Velasco M, Gunnarsdóttir ED, Jagadeesan A, Guðmundsdóttir VB, Thordardóttir EL, Einarsdóttir MS, Moore KHS, Sigurðsson Á, Magnúsdóttir DN, Jónsson H, Snorradóttir S, Hovig E, Møller P, Kockum I, Olsson T, Alfredsson L, Hansen TF, Werge T, Cavalleri GL, et al. 2018. Ancient genomes from Iceland reveal the making of a human population. *Science* **360**:1028-1032. DOI: https://doi.org/10.1126/science.aar2625, PMID: 29853688

Feldman M, Fernández-Domínguez E, Reynolds L, Baird D, Pearson J, Hershkovitz I, May H, Goring-Morris N, Benz M, Gresky J, Bianco RA, Fairbairn A, Mustafaoğlu G, Stockhammer PW, Posth C, Haak W, Jeong C, Krause J. 2019a. Late Pleistocene human genome suggests a local origin for the first farmers of Central Anatolia. *Nature Communications* **10**:1218. DOI: https://doi.org/10.1038/s41467-019-09209-7, PMID: 30890703

Feldman M, Master DM, Bianco RA, Burri M, Stockhammer PW, Mittnik A, Aja AJ, Jeong C, Krause J. 2019b. Ancient DNA sheds light on the genetic origins of early Iron Age Philistines. *Science Advances* **5**:eaax0061. DOI: https://doi.org/10.1126/sciadv.aax0061, PMID: 31281897

Fernandes DM, Strapagiel D, Borówka P, Marciniak B, Żądzińska E, Sirak K, Siska V, Grygiel R, Carlsson J, Manica A, Lorkiewicz W, Pinhasi R. 2018. A genomic Neolithic time transect of hunter-farmer admixture in Central Poland. *Scientific Reports* **8**:14879. DOI: https://doi.org/10.1038/s41598-018-33067-w, PMID: 30291256

Fernandes DM, Mittnik A, Olalde I, Lazaridis I, Cheronet O, Rohland N, Mallick S, Bernardos R, Broomandkhoshbacht N, Carlsson M, Culleton BJ, Ferry M, Gamarra B, Lari M, Mah M, Michel M, Modi A, Novak M, Oppenheimer J, Sirak KA, et al. 2020. Author correction: the spread of steppe and Iranian-related ancestry in the islands of the Western Mediterranean. *Nature Ecology & Evolution* **4**:764. DOI: https://doi.org/10.1038/s41559-020-1197-3, PMID: 32296139

Fregel R, Méndez FL, Bokbot Y, Martín-Socas D, Camalich-Massieu MD, Santana J, Morales J, Ávila-Arcos MC, Underhill PA, Shapiro B, Wojcik G, Rasmussen M, Soares AER, Kapp J, Sockell A, Rodríguez-Santos FJ, Mikdad A, Trujillo-Mederos A, Bustamante CD. 2018. Ancient genomes from North Africa evidence prehistoric migrations to the Maghreb from both the Levant and Europe. *PNAS* **115**:6774-6779. DOI: https://doi.org/10.1073/pnas.1800851115, PMID: 29895688

Fu Q, Posth C, Hajdinjak M, Petr M, Mallick S, Fernandes D, Furtwängler A, Haak W, Meyer M, Mittnik A, Nickel B, Peltzer A, Rohland N, Slon V, Talamo S, Lazaridis I, Lipson M, Mathieson I, Schiffels S, Skoglund P,
et al. 2016. The genetic history of Ice age Europe. *Nature* **534**:200–205. DOI: https://doi.org/10.1038/nature17993, PMID: 27135931

**Furtwängler A**, Neukamm J, Böhme L, Reiter E, Vollstedt M, Arora N, Singh P, Cole ST, Knauf S, Calvignac-Spencer S, Krause-Kyora B, Krause J, Schuenemann VJ, Herbig A. 2020. Comparison of target enrichment strategies for ancient pathogen DNA. *BioTechniques* **69**:455–459. DOI: https://doi.org/10.2144/btn-2020-0100, PMID: 33135465

**Gamba C**, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, Domboróczki L, Kövári I. Pap I, Anders A, Whittle A, Dani J, Raczky P, Higham TFG, Hofreiter M, Bradley DG, Pinhasi R. 2014. Genome flux and stasis in a five millennium transect of European prehistory. *Nature Communications* **5**:5257. DOI: https://doi.org/10.1038/ncomms6257, PMID: 25334030

**Gokhman D**, Nissim-Rafinia M, Agranat-Tamir L, Housman G, García-Pérez R, Lizano E, Cheronet O, Mallick S, Nieves-Colón MA, Li H, Alpaslan-Roodenberg S, Novak M, Gu H, Osinski JM, Ferrando-Bernal M, Gelabert P, Lipende I, Mjungu D, Kondova I, Bontrop R, et al. 2020. Differential DNA methylation of vocal and facial anatomy genes in modern humans. *Nature Communications* **11**:1189. DOI: https://doi.org/10.1038/s41467-020-15020-6, PMID: 32132541

**González-Fortes G**, Jones ER, Lightfoot E, Bonsall C, Lazar C, Grandal-d'Anglade A, Garralda MD, Drak L, Siska V, Simalcsik A, Boronean A, Vidal Romaní JR, Vaqueiro Rodríguez M, Arias P, Pinhasi R, Manica A, Hofreiter M. 2017. Paleogenomic evidence for multi-generational mixing between Neolithic farmers and Mesolithic hunter-gatherers in the Lower Danube Basin. *Current Biology* **27**:1801–1810. DOI: https://doi.org/10.1016/j.cub.2017.05.023, PMID: 28552360

**González-Fortes G**, Tassi F, Trucchi E, Henneberger K, Paijmans JLA, Díez-Del-Molino D, Schroeder H, Susca RR, Barroso-Ruíz C, Bermudez FJ, Barroso-Medina C, Bettencourt AMS, Sampaio HA, Grandal-d'Anglade A, Salas A, de Lombera-Hermida A, Fabregas Valcarce R, Vaquero M, Alonso S, Lozano M, et al. 2019. A western route of prehistoric human migration from Africa into the Iberian Peninsula. *Proceedings. Biological Sciences* **286**:20182288. DOI: https://doi.org/10.1098/rspb.2018.2288, PMID: 30963949

**Günther T**, Valdiosera C, Malmström H, Ureña I, Rodriguez-Varela R, Sverrisdóttir ÓÓ, Daskalaki EA, Skoglund P, Naidoo T, Svensson EM, Bermúdez de Castro JM, Carbonell E, Dunn M, Storå J, Iriarte E, Arsuaga JL, Carretero JM, Götherström A, Jakobsson M. 2015. Ancient genomes link early farmers from Atapuerca in Spain to modern-day Basques. *PNAS* **112**:11917–11922. DOI: https://doi.org/10.1073/pnas.1509851112, PMID: 26351665

**Günther T**, Malmström H, Svensson EM, Omrak A, Sánchez-Quinto F, Kilinç GM, Krzewińska M, Eriksson G, Fraser M, Edlund H, Munters AR, Coutinho A, Simões LG, Vicente M, Sjölander A, Jansen Sellevold B, Jørgensen R, Claes P, Shriver MD, Valdiosera C, et al. 2018. Population genomics of Mesolithic Scandinavia: investigating early postglacial migration routes and high-latitude adaptation. *PLOS Biology* **16**:e2003703. DOI: https://doi.org/10.1371/journal.pbio.2003703, PMID: 29315301

**Haak W**, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, Brandt G, Nordenfelt S, Harney E, Stewardson K, Fu Q, Mittnik A, Bánffy E, Economou C, Francken M, Friederich S, Pena RG, Hallgren F, Khartanovich V, Khokhlov A, et al. 2015. Massive migration from the steppe was a source for Indo-European languages in Europe. *Nature* **522**:207–211. DOI: https://doi.org/10.1038/nature14317, PMID: 25731166

**Haber M**, Doumet-Serhal C, Scheib C, Xue Y, Danecek P, Mezzavilla M, Youhanna S, Martiniano R, Prado-Martinez J, Szpak M, Matisoo-Smith E, Schutkowski H, Mikulski R, Zalloua P, Kivisild T, Tyler-Smith C. 2017. Continuity and admixture in the last five millennia of Levantine history from Ancient Canaanite and Present-Day Lebanese genome sequences. *American Journal of Human Genetics* **101**:274–282. DOI: https://doi.org/10.1016/j.ajhg.2017.06.013, PMID: 28757201

**Haber M**, Doumet-Serhal C, Scheib CL, Xue Y, Mikulski R, Martiniano R, Fischer-Genz B, Schutkowski H, Kivisild T, Tyler-Smith C. 2019. A transient pulse of genetic admixture from the Crusaders in the near east identified from ancient genome sequences. *American Journal of Human Genetics* **104**:977–984. DOI: https://doi.org/10.1016/j.ajhg.2019.03.015, PMID: 31006515

**Haber M**, Nassar J, Almarri MA, Saupe T, Saag L, Griffith SJ, Doumet-Serhal C, Chanteau J, Saghieh-Beydoun M, Xue Y, Scheib CL, Tyler-Smith C. 2020. A genetic history of the near east from an aDNA time course sampling eight points in the past 4,000 years. *American Journal of Human Genetics* **107**:149–157. DOI: https://doi.org/10.1016/j.ajhg.2020.05.008, PMID: 32470374

**Haller BC**, Galloway J, Kelleher J, Messer PW, Ralph PL. 2019. Tree-sequence recording in SLiM opens new horizons for forward-time simulation of whole genomes. *Molecular Ecology Resources* **19**:552–566. DOI: https://doi.org/10.1111/1755-0998.12968, PMID: 30565882

**Haller BC**, Messer PW. 2019. SLiM 3: forward genetic simulations beyond the wright-fisher model. *Molecular Biology and Evolution* **36**:632–637. DOI: https://doi.org/10.1093/molbev/msy228, PMID: 30517680

**Harney É**, May H, Shalem D, Rohland N, Mallick S, Lazaridis I, Sarig R, Stewardson K, Nordenfelt S, Patterson N, Hershkovitz I, Reich D. 2018. Publisher correction: ancient DNA from Chalcolithic Israel reveals the role of population mixture in cultural transformation. *Nature Communications* **9**:3913. DOI: https://doi.org/10.1038/s41467-018-06484-8, PMID: 30237547

**Harney É**, Patterson N, Reich D, Wakeley J. 2021. Assessing the performance of qpAdm: a statistical tool for studying population admixture. *Genetics* **217**:iyaa045. DOI: https://doi.org/10.1093/ genetics/iyaa045, PMID: 33772284

**Harper K.** 2017. *The Fate of Rome: Climate, Disease, and the End of an Empire* Princeton: Princeton University Press. DOI: https://doi.org/10.2307/j.ctv9b2txr
et al. 2016. The genetic history of Ice age Europe. *Nature* **534**:200–205. DOI: https://doi.org/10.1038/nature17993, PMID: 27135931

**Furtwängler A**, Neukamm J, Böhme L, Reiter E, Vollstedt M, Arora N, Singh P, Cole ST, Knauf S, Calvignac-Spencer S, Krause-Kyora B, Krause J, Schuenemann VJ, Herbig A. 2020. Comparison of target enrichment strategies for ancient pathogen DNA. *BioTechniques* **69**:455–459. DOI: https://doi.org/10.2144/btn-2020-0100, PMID: 33135465

**Gamba C**, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, Domboróczki L, Kövári I. Pap I, Anders A, Whittle A, Dani J, Raczky P, Higham TFG, Hofreiter M, Bradley DG, Pinhasi R. 2014. Genome flux and stasis in a five millennium transect of European prehistory. *Nature Communications* **5**:5257. DOI: https://doi.org/10.1038/ncomms6257, PMID: 25334030

**Gokhman D**, Nissim-Rafinia M, Agranat-Tamir L, Housman G, García-Pérez R, Lizano E, Cheronet O, Mallick S, Nieves-Colón MA, Li H, Alpaslan-Roodenberg S, Novak M, Gu H, Osinski JM, Ferrando-Bernal M, Gelabert P, Lipende I, Mjungu D, Kondova I, Bontrop R, et al. 2020. Differential DNA methylation of vocal and facial anatomy genes in modern humans. *Nature Communications* **11**:1189. DOI: https://doi.org/10.1038/s41467-020-15020-6, PMID: 32132541

**González-Fortes G**, Jones ER, Lightfoot E, Bonsall C, Lazar C, Grandal-d'Anglade A, Garralda MD, Drak L, Siska V, Simalcsik A, Boronean A, Vidal Romaní JR, Vaqueiro Rodríguez M, Arias P, Pinhasi R, Manica A, Hofreiter M. 2017. Paleogenomic evidence for multi-generational mixing between Neolithic farmers and Mesolithic hunter-gatherers in the Lower Danube Basin. *Current Biology* **27**:1801–1810. DOI: https://doi.org/10.1016/j.cub.2017.05.023, PMID: 28552360

**González-Fortes G**, Tassi F, Trucchi E, Henneberger K, Paijmans JLA, Díez-Del-Molino D, Schroeder H, Susca RR, Barroso-Ruíz C, Bermudez FJ, Barroso-Medina C, Bettencourt AMS, Sampaio HA, Grandal-d'Anglade A, Salas A, de Lombera-Hermida A, Fabregas Valcarce R, Vaquero M, Alonso S, Lozano M, et al. 2019. A western route of prehistoric human migration from Africa into the Iberian Peninsula. *Proceedings. Biological Sciences* **286**:20182288. DOI: https://doi.org/10.1098/rspb.2018.2288, PMID: 30963949

**Günther T**, Valdiosera C, Malmström H, Ureña I, Rodriguez-Varela R, Sverrisdóttir ÓÓ, Daskalaki EA, Skoglund P, Naidoo T, Svensson EM, Bermúdez de Castro JM, Carbonell E, Dunn M, Storå J, Iriarte E, Arsuaga JL, Carretero JM, Götherström A, Jakobsson M. 2015. Ancient genomes link early farmers from Atapuerca in Spain to modern-day Basques. *PNAS* **112**:11917–11922. DOI: https://doi.org/10.1073/pnas.1509851112, PMID: 26351665

**Günther T**, Malmström H, Svensson EM, Omrak A, Sánchez-Quinto F, Kilinç GM, Krzewińska M, Eriksson G, Fraser M, Edlund H, Munters AR, Coutinho A, Simões LG, Vicente M, Sjölander A, Jansen Sellevold B, Jørgensen R, Claes P, Shriver MD, Valdiosera C, et al. 2018. Population genomics of Mesolithic Scandinavia: investigating early postglacial migration routes and high-latitude adaptation. *PLOS Biology* **16**:e2003703. DOI: https://doi.org/10.1371/journal.pbio.2003703, PMID: 29315301

**Haak W**, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, Brandt G, Nordenfelt S, Harney E, Stewardson K, Fu Q, Mittnik A, Bánffy E, Economou C, Francken M, Friederich S, Pena RG, Hallgren F, Khartanovich V, Khokhlov A, et al. 2015. Massive migration from the steppe was a source for Indo-European languages in Europe. *Nature* **522**:207–211. DOI: https://doi.org/10.1038/nature14317, PMID: 25731166

**Haber M**, Doumet-Serhal C, Scheib C, Xue Y, Danecek P, Mezzavilla M, Youhanna S, Martiniano R, Prado-Martinez J, Szpak M, Matisoo-Smith E, Schutkowski H, Mikulski R, Zalloua P, Kivisild T, Tyler-Smith C. 2017. Continuity and admixture in the last five millennia of Levantine history from Ancient Canaanite and Present-Day Lebanese genome sequences. *American Journal of Human Genetics* **101**:274–282. DOI: https://doi.org/10.1016/j.ajhg.2017.06.013, PMID: 28757201

**Haber M**, Doumet-Serhal C, Scheib CL, Xue Y, Mikulski R, Martiniano R, Fischer-Genz B, Schutkowski H, Kivisild T, Tyler-Smith C. 2019. A transient pulse of genetic admixture from the Crusaders in the near east identified from ancient genome sequences. *American Journal of Human Genetics* **104**:977–984. DOI: https://doi.org/10.1016/j.ajhg.2019.03.015, PMID: 31006515

**Haber M**, Nassar J, Almarri MA, Saupe T, Saag L, Griffith SJ, Doumet-Serhal C, Chanteau J, Saghieh-Beydoun M, Xue Y, Scheib CL, Tyler-Smith C. 2020. A genetic history of the near east from an aDNA time course sampling eight points in the past 4,000 years. *American Journal of Human Genetics* **107**:149–157. DOI: https://doi.org/10.1016/j.ajhg.2020.05.008, PMID: 32470374

**Haller BC**, Galloway J, Kelleher J, Messer PW, Ralph PL. 2019. Tree-sequence recording in SLiM opens new horizons for forward-time simulation of whole genomes. *Molecular Ecology Resources* **19**:552–566. DOI: https://doi.org/10.1111/1755-0998.12968, PMID: 30565882

**Haller BC**, Messer PW. 2019. SLiM 3: forward genetic simulations beyond the wright-fisher model. *Molecular Biology and Evolution* **36**:632–637. DOI: https://doi.org/10.1093/molbev/msy228, PMID: 30517680

**Harney É**, May H, Shalem D, Rohland N, Mallick S, Lazaridis I, Sarig R, Stewardson K, Nordenfelt S, Patterson N, Hershkovitz I, Reich D. 2018. Publisher correction: ancient DNA from Chalcolithic Israel reveals the role of population mixture in cultural transformation. *Nature Communications* **9**:3913. DOI: https://doi.org/10.1038/s41467-018-06484-8, PMID: 30237547

**Harney É**, Patterson N, Reich D, Wakeley J. 2021. Assessing the performance of qpAdm: a statistical tool for studying population admixture. *Genetics* **217**:iyaa045. DOI: https://doi.org/10.1093/genetics/iyaa045, PMID: 33772284

**Harper K.** 2017. *The Fate of Rome: Climate, Disease, and the End of an Empire* Princeton: Princeton University Press. DOI: https://doi.org/10.2307/j.ctv9b2txr
Järve M, Saag L, Scheib CL, Pathak AK, Montinaro F, Pagani L, Flores R, Guellil M, Saag L, Tambets K, Kushniarevich A, Solnik A, Varul L, Zadnikov S, Petrauskas O, Avramenko M, Magomedov B, Didenko S, Toshev G, Bruyako I, et al. 2019. Shifts in the genetic landscape of the Western Eurasian Steppe associated with the beginning and end of the Scythian dominance. *Current Biology* **29**:2430-2441. DOI: https://doi.org/10.1016/j.cub.2019.06.019, PMID: 31303491

Jeong C, Balanovsky O, Lukianova E, Kahbatkyzy N, Flegontov P, Zaporozhchenko V, Immel A, Wang CC, Ixan O, Khussainova E, Bekmanov B, Zaibert V, Lavryashina M, Pocheshkhova E, Yusupov Y, Agdzhoyan A, Koshel S, Bukin A, Nymadawa P, Turdikulova S, et al. 2019. The genetic history of admixture across inner Eurasia. *Nature Ecology & Evolution* **3**:966-976. DOI: https://doi.org/10.1038/s41559-019-0878-2, PMID: 31036896

Jones ER, Gonzalez-Fortes G, Connell S, Siska V, Eriksson A, Martiniano R, McLaughlin RL, Gallego Llorente M, Cassidy LM, Gamba C, Meshveliani T, Bar-Yosef O, Müller W, Belfer-Cohen A, Matskevich Z, Jakeli N, Higham TFG, Currat M, Lordkipanidze D, Hofreiter M, et al. 2015. Upper Palaeolithic genomes reveal deep roots of modern Eurasians. *Nature Communications* **6**:8912. DOI: https://doi.org/10.1038/ncomms9912, PMID: 26567969

Jones ER, Zarina G, Moiseyev V, Lightfoot E, Nigst PR, Manica A, Pinhasi R, Bradley DG. 2017. The Neolithic transition in the baltic was not driven by admixture with early European farmers. *Current Biology* **27**:576-582. DOI: https://doi.org/10.1016/j.cub.2016.12.060, PMID: 28162894

Jónsson H, Ginolhac A, Schubert M, Johnson PLF, Orlando L. 2013. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. *Bioinformatics* **29**:1682-1684. DOI: https://doi.org/10.1093/bioinformatics/btt193, PMID: 23613487

Kelleher J, Etheridge AM, McVean G. 2016. Efficient coalescent simulation and genealogical analysis for large sample sizes. *PLOS Computational Biology* **12**:e1004842. DOI: https://doi.org/10.1371/journal.pcbi.1004842, PMID: 27145223

Kelleher J, Thornton KR, Ashander J, Ralph PL. 2018. Efficient pedigree recording for fast population genetics simulation. *PLOS Computational Biology* **14**:e1006581. DOI: https://doi.org/10.1371/journal.pcbi.1006581, PMID: 30383757

Keller A, Graefen A, Ball M, Matzas M, Boisguerin V, Maixner F, Leidinger P, Backes C, Khairat R, Forster M, Stade B, Franke A, Mayer J, Spangler J, McLaughlin S, Shah M, Lee C, Harkins TT, Sartori A, Moreno-Estrada A, et al. 2012. New insights into the Tyrolean Iceman's origin and phenotype as inferred by whole-genome sequencing. *Nature Communications* **3**:698. DOI: https://doi.org/10.1038/ncomms1701, PMID: 22426219

Kılınç GM, Omrak A, Özer F, Günther T, Büyükkarakaya AM, Bıçakçı E, Baird D, Dönerttaş HM, Ghalichi A, Yaka R, Koptekin D, Açan SC, Parvizi P, Krzewińska M, Daskalaki EA, Yüncü E, Dağtaş ND, Fairbairn A, Pearson J, Mustafaoğlu G, et al. 2016. The demographic development of the first farmers in Anatolia. *Current Biology* **26**:2659-2666. DOI: https://doi.org/10.1016/j.cub.2016.07.057, PMID: 27498567

Korneliussen TS, Albrechtsen A, Nielsen R. 2014. ANGSD: analysis of next generation sequencing data. BMC Bioinformatics **15**:356. DOI: https://doi.org/10.1186/s12859-014-0356-4, PMID: 25420514

Kovacevic L, Tambets K, Ilumäe AM, Kushniarevich A, Yunusbayev B, Solnik A, Bego T, Primorac D, Skaro V, Leskovac A, Jakovski Z, Drobnic K, Tolk HV, Kovacevic S, Rudan P, Metspalu E, Marjanovic D. 2014. Standing at the gateway to Europe--the genetic structure of Western balkan populations based on autosomal and haploid markers. *PLOS ONE* **9**:e105090. DOI: https://doi.org/10.1371/journal.pone.0105090, PMID: 25148043

Krzewińska M, Kjellström A, Günther T, Hedenstierna-Jonson C, Zachrisson T, Omrak A, Yaka R, Kılınç GM, Somel M, Sobrado V, Evans J, Knipper C, Jakobsson M, Storå J, Götherström A. 2018a. Genomic and strontium isotope variation reveal immigration patterns in a Viking Age town. *Current Biology* **28**:2730-2738. DOI: https://doi.org/10.1016/j.cub.2018.06.053, PMID: 30146150

Krzewińska M, Kılınç GM, Juras A, Koptekin D, Chyleński M, Nikitin AG, Shcherbakov N, Shuteleva I, Leonova T, Kraeva L, Sungatov FA, Sultanova AN, Potekhina I, Łukasik S, Krenz-Niedbała M, Dalén L, Sinika V, Jakobsson M, Storå J, Götherström A. 2018b. Ancient genomes suggest the eastern Pontic-Caspian steppe as the source of western Iron Age nomads. *Science Advances* **4**:eaat4457. DOI: https://doi.org/10.1126/sciadv.aat4457, PMID: 30417088

Lamnidis TC, Majander K, Jeong C, Salmela E, Wessman A, Moiseyev V, Khartanovich V, Balanovsky O, Ongyerth M, Weihmann A, Sajantila A, Kelso J, Pääbo S, Onkamo P, Haak W, Krause J, Schiffels S. 2018. Ancient Fennoscandian genomes reveal origin and spread of Siberian ancestry in Europe. *Nature Communications* **9**:5018. DOI: https://doi.org/10.1038/s41467-018-07483-5, PMID: 30479341

Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, Sudmant PH, Schraiber JG, Castellano S, Lipson M, Berger B, Economou C, Bollongino R, Fu Q, Bos KI, Nordenfelt S, Li H, de Filippo C, Prüfer K, Sawyer S, et al. 2014. Ancient human genomes suggest three ancestral populations for present-day Europeans. *Nature* **513**:409-413. DOI: https://doi.org/10.1038/nature13673, PMID: 25230663

Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, Fernandes D, Novak M, Gamarra B, Sirak K, Connell S, Stewardson K, Harney E, Fu Q, Gonzalez-Fortes G, Jones ER, Roodenberg SA, Lengyel G, Bocquentin F, Gasparian B, et al. 2016. Genomic insights into the origin of farming in the ancient Near East. *Nature* **536**:419-424. DOI: https://doi.org/10.1038/nature19310, PMID: 27459054

Lazaridis I, Mittnik A, Patterson N, Mallick S, Rohland N, Pfremgle S, Furtwängler A, Peltzer A, Posth C, Vasilakis A, McGeorge PJP, Konsolaki-Yannopoulou E, Korres G, Martlew H, Michalodimitrakis M, Özsait M, Özsait N, Papathanasiou A, Richards M, Roodenberg SA, et al. 2017. Genetic origins of the Minoans and Mycenaean. *Nature* **548**:214-218. DOI: https://doi.org/10.1038/nature23310

Li H, Durbin R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **25**:1754-1760. DOI: https://doi.org/10.1093/bioinformatics/btp324, PMID: 19451168