The genomic history of the Aegean palatial civilizations

Ancient genomics support deep divergence between Eastern and Western Mediterranean Indo-European languages

The Indo-European languages are among the most widely spoken in the world, yet their early diversification remains contentious1-5. It is widely accepted that the spread of this language family across Europe from the 5th millennium BP correlates with the expansion and diversification of steppe-related genetic ancestry from the onset of the Bronze Age6,7. However, multiple steppe-derived populations co-existed in Europe during this period, and it remains unclear how these populations diverged and which provided the demographic channels for the ancestral forms of the Italic, Celtic, Greek, and Armenian languages8,9. To investigate the ancestral histories of Indo-European-speaking groups in Southern Europe, we sequenced genomes from 314 ancient individuals from the Mediterranean and surrounding regions, spanning from 5,200 BP to 2,100 BP, and co-analysed these with published genome data. We additionally conducted strontium isotope analyses on 224 of these individuals. We find a deep east-west divide of steppe ancestry in Southern Europe during the Bronze Age. Specifically, we show that the arrival of steppe ancestry in Spain, France, and Italy was mediated by Bell Beaker (BB) populations of Western Europe, likely contributing to the emergence of the Italic and Celtic languages. In contrast, Armenian and Greek populations acquired steppe ancestry directly from Yamnaya groups of Eastern Europe. These results are consistent with the linguistic Italo-Celtic10,11 and Graeco-Armenian1,12,13 hypotheses accounting for the origins of most Mediterranean Indo-European languages of Classical Antiquity. Our findings thus align with specific linguistic divergence models for the Indo-European language family while contradicting others. This underlines the power of ancient DNA in uncovering prehistoric diversifications of human populations and language communities.

The genomic portrait of the Picene culture provides new insights into the Italic Iron Age and the legacy of the Roman Empire in Central Italy

BACKGROUND

The Italic Iron Age is characterized by the presence of various ethnic groups partially examined from a genomic perspective. To explore the evolution of Iron Age Italic populations and the genetic impact of Romanization, we focus on the Picenes, one of the most fascinating pre-Roman civilizations, who flourished on the Middle Adriatic side of Central Italy between the 9th and the 3rd century BCE, until the Roman colonization.

RESULTS

More than 50 samples are reported, spanning more than 1000 years of history from the Iron Age to Late Antiquity. Despite cultural diversity, our analysis reveals no major differences between the Picenes and other coeval populations, suggesting a shared genetic history of the Central Italian Iron Age ethnic groups. Nevertheless, a slight genetic differentiation between populations along the Adriatic and Tyrrhenian coasts can be observed, possibly due to different population dynamics in the two sides of Italy and/or genetic contacts across the Adriatic Sea. Additionally, we identify several individuals with ancestries deviating from their general population. Lastly, in our Late Antiquity site, we observe a drastic change in the genetic landscape of the Middle Adriatic region, indicating a relevant influx from the Near East, possibly as a consequence of Romanization.

CONCLUSIONS

Our findings, consistently with archeological hypotheses, suggest genetic interactions across the Adriatic Sea during the Bronze/Iron Age and a high level of individual mobility typical of cosmopolitan societies. Finally, we highlight the role of the Roman Empire in shaping genetic and phenotypic changes that greatly impact the Italian peninsula.

An explanation for the sister repulsion phenomenon in Patterson's f-statistics

Patterson's f-statistics are among the most heavily utilized tools for analyzing genome-wide allele frequency data for demographic inference. Beyond studying admixture, f3- and f4-statistics are also used for clustering populations to identify groups with similar histories. However, previous studies have noted an unexpected behavior of f-statistics: multiple populations from a certain region systematically show higher genetic affinity to a more distant population than to their neighbors, a pattern that is mismatched with alternative measures of genetic similarity. We call this counter-intuitive pattern "sister repulsion". We first present a novel instance of sister repulsion, where genomes from Bronze Age East Anatolian sites show higher affinity toward Bronze Age Greece rather than each other. This is observed both using f3- and f4-statistics, contrasts with archaeological/historical expectation, and also contradicts genetic affinity patterns captured using principal components analysis or multidimensional scaling on genetic distances. We then propose a simple demographic model to explain this pattern, where sister populations receive gene flow from a genetically distant source. We calculate f3- and f4-statistics using simulated genetic data with varying population genetic parameters, confirming that low-level gene flow from an external source into populations from 1 region can create sister repulsion in f-statistics. Unidirectional gene flow between the studied regions (without an external source) can likewise create repulsion. Meanwhile, similar to our empirical observations, multidimensional scaling analyses of genetic distances still cluster sister populations together. Overall, our results highlight the impact of low-level admixture events when inferring demographic history using f-statistics.

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

We present a method for detecting evidence of natural selection in ancient DNA time-series data that leverages an opportunity not utilized in previous scans: testing for a consistent trend in allele frequency change over time. By applying this to 8433 West Eurasians who lived over the past 14000 years and 6510 contemporary people, we find an order of magnitude more genome-wide significant signals than previous studies: 347 independent loci with >99% probability of selection. Previous work showed that classic hard sweeps driving advantageous mutations to fixation have been rare over the broad span of human evolution, but in the last ten millennia, many hundreds of alleles have been affected by strong directional selection. Discoveries include an increase from ~0% to ~20% in 4000 years for the major risk factor for celiac disease at HLA-DQB1; a rise from ~0% to ~8% in 6000 years of blood type B; and fluctuating selection at the TYK2 tuberculosis risk allele rising from ~2% to ~9% from ~5500 to ~3000 years ago before dropping to ~3%. We identify instances of coordinated selection on alleles affecting the same trait, with the polygenic score today predictive of body fat percentage decreasing by around a standard deviation over ten millennia, consistent with the "Thrifty Gene" hypothesis that a genetic predisposition to store energy during food scarcity became disadvantageous after farming. We also identify selection for combinations of alleles that are today associated with lighter skin color, lower risk for schizophrenia and bipolar disease, slower health decline, and increased measures related to cognitive performance (scores on intelligence tests, household income, and years of schooling). These traits are measured in modern industrialized societies, so what phenotypes were adaptive in the past is unclear. We estimate selection coefficients at 9.9 million variants, enabling study of how Darwinian forces couple to allelic effects and shape the genetic architecture of complex traits.

Testing times: disentangling admixture histories in recent and complex demographies using ancient DNA

Our knowledge of human evolutionary history has been greatly advanced by paleogenomics. Since the 2020s, the study of ancient DNA has increasingly focused on reconstructing the recent past. However, the accuracy of paleogenomic methods in resolving questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation remains an open question. We evaluated the performance and behavior of two commonly used methods, qpAdm and the f3-statistic, on admixture inference under a diversity of demographic models and data conditions. We performed two complementary simulation approaches-firstly exploring a wide demographic parameter space under four simple demographic models of varying complexities and configurations using branch-length data from two chromosomes-and secondly, we analyzed a model of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudohaploidization. We observe that population differentiation is the primary factor driving qpAdm performance. Notably, while complex gene flow histories influence which models are classified as plausible, they do not reduce overall performance. Under conditions reflective of the historical period, qpAdm most frequently identifies the true model as plausible among a small candidate set of closely related populations. To increase the utility for resolving fine-scaled hypotheses, we provide a heuristic for further distinguishing between candidate models that incorporates qpAdm model P-values and f3-statistics. Finally, we demonstrate a significant performance increase for qpAdm using whole-genome branch-length f2-statistics, highlighting the potential for improved demographic inference that could be achieved with future advancements in f-statistic estimations.

Evaluation of genotype imputation using Glimpse tools on low coverage ancient DNA

Ancient DNA provides a unique frame for directly studying human population genetics in time and space. Still, since most of the ancient genomic data is low coverage, analysis is confronted with a low number of SNPs, genotype uncertainties, and reference-bias. Here, we for the first time benchmark the two distinct versions of Glimpse tools on 120 ancient human genomes from Eurasia including those largely from previously under-evaluated regions and compare the performance of genotype imputation with de facto analysis approaches for low coverage genomic data analysis. We further investigate the impact of two distinct reference panels on imputation accuracy for low coverage genomic data. We compute accuracy statistics and perform PCA and f4-statistics to explore the behaviour of genotype imputation on low coverage data regarding (i)two versions of Glimpse, (ii)two reference panels, (iii)four post-imputation filters and coverages, as well as (iv)data type and geographical origin of the samples on the analyses. Our results reveal that even for 0.1X coverage ancient human genomes, genotype imputation using Glimpse-v2 is suitable. Additionally, using the 1000 Genomes merged with Human Genome Diversity Panel improves the accuracy of imputation for the rare variants with low MAF, which might be important not only for ancient genomics but also for modern human genomic studies based on low coverage data and for haplotype-based analysis. Most importantly, we reveal that genotype imputation of low coverage ancient human genomes reduces the genetic affinity of the samples towards human reference genome. Through solving one of the most challenging biases in data analysis, so-called reference bias, genotype imputation using Glimpse v2 is promising for low coverage ancient human genomic data analysis and for rare-variant-based and haplotype-based analysis.

Modelling the demographic history of human North African genomes points to a recent soft split divergence between populations

BACKGROUND

North African human populations present a complex demographic scenario due to the presence of an autochthonous genetic component and population substructure, plus extensive gene flow from the Middle East, Europe, and sub-Saharan Africa.

RESULTS

We conducted a comprehensive analysis of 364 genomes to construct detailed demographic models for the North African region, encompassing its two primary ethnic groups, the Arab and Amazigh populations. This was achieved through an Approximate Bayesian Computation with Deep Learning (ABC-DL) framework and a novel algorithm called Genetic Programming for Population Genetics (GP4PG). This innovative approach enabled us to effectively model intricate demographic scenarios, utilizing a subset of 16 whole genomes at > 30X coverage. The demographic model suggested by GP4PG exhibited a closer alignment with the observed data compared to the ABC-DL model. Both point to a back-to-Africa origin of North African individuals and a close relationship with Eurasian populations. Results support different origins for Amazigh and Arab populations, with Amazigh populations originating back in Epipaleolithic times, while GP4PG supports Arabization as the main source of Middle Eastern ancestry. The GP4PG model includes population substructure in surrounding populations (sub-Saharan Africa and Middle East) with continuous decaying gene flow after population split. Contrary to ABC-DL, the best GP4PG model does not require pulses of admixture from surrounding populations into North Africa pointing to soft splits as drivers of divergence in North Africa.

CONCLUSIONS

We have built a demographic model on North Africa that points to a back-to-Africa expansion and a differential origin between Arab and Amazigh populations.

Improving data archiving practices in ancient genomics

Ancient DNA is producing a rich record of past genetic diversity in humans and other species. However, unless the primary data is appropriately archived, its long-term value will not be fully realised. I surveyed publicly archived data from 42 recent ancient genomics studies. Half of the studies archived incomplete datasets, preventing accurate replication and representing a loss of data of potential future use. No studies met all criteria that could be considered best practice. Based on these results, I make six recommendations for data producers: (1) archive all sequencing reads, not just those that aligned to a reference genome, (2) archive read alignments too, but as secondary analysis files, (3) provide correct experiment metadata on samples, libraries and sequencing runs, (4) provide informative sample metadata, (5) archive data from low-coverage and negative experiments, and (6) document archiving choices in papers, and peer review these. Given the reliance on destructive sampling of finite material, ancient genomics studies have a particularly strong responsibility to ensure the longevity and reusability of generated data.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

Ancient genomes illuminate Eastern Arabian population history and adaptation against malaria

The harsh climate of Arabia has posed challenges in generating ancient DNA from the region, hindering the direct examination of ancient genomes for understanding the demographic processes that shaped Arabian populations. In this study, we report whole-genome sequence data obtained from four Tylos-period individuals from Bahrain. Their genetic ancestry can be modeled as a mixture of sources from ancient Anatolia, Levant, and Iran/Caucasus, with variation between individuals suggesting population heterogeneity in Bahrain before the onset of Islam. We identify the G6PD Mediterranean mutation associated with malaria resistance in three out of four ancient Bahraini samples and estimate that it rose in frequency in Eastern Arabia from 5 to 6 kya onward, around the time agriculture appeared in the region. Our study characterizes the genetic composition of ancient Arabians, shedding light on the population history of Bahrain and demonstrating the feasibility of studies of ancient DNA in the region.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

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

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.

The genomic echoes of the last Green Sahara on the Fulani and Sahelian people

The population history of the Sahara/Sahelian belt is understudied, despite previous work highlighting complex dynamics.1,2,3,4,5,6,7 The Sahelian Fulani, i.e., the largest nomadic pastoral population in the world,8 represent an interesting case because they show a non-negligible proportion of an Eurasian genetic component, usually explained by recent admixture with northern Africans.1,2,5,6,7,9,10,11,12 Nevertheless, their origins are largely unknown, although several hypotheses have been proposed, including a possible link to ancient peoples settled in the Sahara during its last humid phase (Green Sahara, 12,000-5,000 years before present [BP]).13,14,15 To shed light about the Fulani ancient genetic roots, we produced 23 high-coverage (30×) whole genomes from Fulani individuals from 8 Sahelian countries, plus 17 samples from other African groups and 3 from Europeans as controls, for a total of 43 new whole genomes. These data have been compared with 814 published modern whole genomes2,16,17,18 and with relevant published ancient sequences (> 1,800 samples).19 These analyses showed some evidence that the non-sub-Saharan genetic ancestry component of the Fulani might have also been shaped by older events,1,5,6 possibly tracing the Fulani origins to unsampled ancient Green Saharan population(s). The joint analysis of modern and ancient samples allowed us to shed light on the genetic ancestry composition of such ancient Saharans, suggesting a similarity with Late Neolithic Moroccans and possibly pointing to a link with the spread of cattle herding. We also identified two different Fulani clusters whose admixture pattern may be informative about the historical Fulani movements and their later involvement in the western African empires.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f3-statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P-values of qpAdm models, and f3-statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f2-statistics, demonstrating the improved demographic inference that could be achieved with advancements in f-statistic estimations.

The genomic history of the indigenous people of the Canary Islands

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Language trees with sampled ancestors support a hybrid model for the origin of Indo-European languages

The origins of the Indo-European language family are hotly disputed. Bayesian phylogenetic analyses of core vocabulary have produced conflicting results, with some supporting a farming expansion out of Anatolia ~9000 years before present (yr B.P.), while others support a spread with horse-based pastoralism out of the Pontic-Caspian Steppe ~6000 yr B.P. Here we present an extensive database of Indo-European core vocabulary that eliminates past inconsistencies in cognate coding. Ancestry-enabled phylogenetic analysis of this dataset indicates that few ancient languages are direct ancestors of modern clades and produces a root age of ~8120 yr B.P. for the family. Although this date is not consistent with the Steppe hypothesis, it does not rule out an initial homeland south of the Caucasus, with a subsequent branch northward onto the steppe and then across Europe. We reconcile this hybrid hypothesis with recently published ancient DNA evidence from the steppe and the northern Fertile Crescent.

Inundation risk assessment of Eastern Mediterranean Coastal archaeological and historical sites of Türkiye and Greece

Archaeological and historical sites, including World Heritage Sites, archaeological sites, castles, fortresses, towers, bridges, necropolises, and ruins, which are a shared legacy of humanity, are at risk of inundation due to sea level rise resulting from global climate change and increased carbon emissions, much like settlements, transportation systems, agricultural areas, and animal habitats. These locations are significant hot spots for connecting the past to the present regarding human history. Throughout history, settlements have been established around wetlands and shorelines. For this reason, the most significant danger threatening coastal archaeological and historical sites is the risk of inundation due to sea level rise. Creating an inundation risk map could be beneficial for planning and conservation studies due to reasons such as the long duration of archaeological studies, lack of a qualified workforce, and budget constraints. To this end, inundation risk maps were created for five different scenarios with risk classes varying from very high to very low, considering global and local sea level rise predictions. The study evaluated 464 historic areas located on the Eastern Mediterranean coasts of Türkiye and Greece. The findings indicate that a significant number of archaeological and historical sites, such as Knidos and Delos, will be inundated even in the lowest sea level rise scenarios. According to the inundation risk map, 34 sites are at very high risk, 21 at high risk, 25 at medium risk, 26 at low risk, and 41 at very low risk of inundation.

Imputation of ancient human genomes

Due to postmortem DNA degradation and microbial colonization, most ancient genomes have low depth of coverage, hindering genotype calling. Genotype imputation can improve genotyping accuracy for low-coverage genomes. However, it is unknown how accurate ancient DNA imputation is and whether imputation introduces bias to downstream analyses. Here we re-sequence an ancient trio (mother, father, son) and downsample and impute a total of 43 ancient genomes, including 42 high-coverage (above 10x) genomes. We assess imputation accuracy across ancestries, time, depth of coverage, and sequencing technology. We find that ancient and modern DNA imputation accuracies are comparable. When downsampled at 1x, 36 of the 42 genomes are imputed with low error rates (below 5%) while African genomes have higher error rates. We validate imputation and phasing results using the ancient trio data and an orthogonal approach based on Mendel's rules of inheritance. We further compare the downstream analysis results between imputed and high-coverage genomes, notably principal component analysis, genetic clustering, and runs of homozygosity, observing similar results starting from 0.5x coverage, except for the African genomes. These results suggest that, for most populations and depths of coverage as low as 0.5x, imputation is a reliable method that can improve ancient DNA studies.

Mapache: a flexible pipeline to map ancient DNA

SUMMARY

We introduce mapache, a flexible, robust and scalable pipeline to map, quantify and impute ancient and present-day DNA in a reproducible way. Mapache is implemented in the workflow manager Snakemake and is optimized for low-space consumption, allowing to efficiently (re)map large datasets-such as reference panels and multiple extracts and libraries per sample - to one or several genomes. Mapache can easily be customized or combined with other Snakemake tools.

AVAILABILITY AND IMPLEMENTATION

Mapache is freely available on GitHub (https://github.com/sneuensc/mapache). An extensive manual is provided at https://github.com/sneuensc/mapache/wiki.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ancient DNA reveals admixture history and endogamy in the prehistoric Aegean

The Neolithic and Bronze Ages were highly transformative periods for the genetic history of Europe but for the Aegean-a region fundamental to Europe's prehistory-the biological dimensions of cultural transitions have been elucidated only to a limited extent so far. We have analysed newly generated genome-wide data from 102 ancient individuals from Crete, the Greek mainland and the Aegean Islands, spanning from the Neolithic to the Iron Age. We found that the early farmers from Crete shared the same ancestry as other contemporaneous Neolithic Aegeans. In contrast, the end of the Neolithic period and the following Early Bronze Age were marked by 'eastern' gene flow, which was predominantly of Anatolian origin in Crete. Confirming previous findings for additional Central/Eastern European ancestry in the Greek mainland by the Middle Bronze Age, we additionally show that such genetic signatures appeared in Crete gradually from the seventeenth to twelfth centuries BC, a period when the influence of the mainland over the island intensified. Biological and cultural connectedness within the Aegean is also supported by the finding of consanguineous endogamy practiced at high frequencies, unprecedented in the global ancient DNA record. Our results highlight the potential of archaeogenomic approaches in the Aegean for unravelling the interplay of genetic admixture, marital and other cultural practices.

Spatial and temporal heterogeneity in human mobility patterns in Holocene Southwest Asia and the East Mediterranean

We present a spatiotemporal picture of human genetic diversity in Anatolia, Iran, Levant, South Caucasus, and the Aegean, a broad region that experienced the earliest Neolithic transition and the emergence of complex hierarchical societies. Combining 35 new ancient shotgun genomes with 382 ancient and 23 present-day published genomes, we found that genetic diversity within each region steadily increased through the Holocene. We further observed that the inferred sources of gene flow shifted in time. In the first half of the Holocene, Southwest Asian and the East Mediterranean populations homogenized among themselves. Starting with the Bronze Age, however, regional populations diverged from each other, most likely driven by gene flow from external sources, which we term "the expanding mobility model." Interestingly, this increase in inter-regional divergence can be captured by outgroup-f₃-based genetic distances, but not by the commonly used F_ST statistic, due to the sensitivity of F_ST, but not outgroup-f₃, to within-population diversity. Finally, we report a temporal trend of increasing male bias in admixture events through the Holocene.

A genomic snapshot of demographic and cultural dynamism in Upper Mesopotamia during the Neolithic Transition

Upper Mesopotamia played a key role in the Neolithic Transition in Southwest Asia through marked innovations in symbolism, technology, and diet. We present 13 ancient genomes (c. 8500 to 7500 cal BCE) from Pre-Pottery Neolithic Çayönü in the Tigris basin together with bioarchaeological and material culture data. Our findings reveal that Çayönü was a genetically diverse population, carrying mixed ancestry from western and eastern Fertile Crescent, and that the community received immigrants. Our results further suggest that the community was organized along biological family lines. We document bodily interventions such as head shaping and cauterization among the individuals examined, reflecting Çayönü's cultural ingenuity. Last, we identify Upper Mesopotamia as the likely source of eastern gene flow into Neolithic Anatolia, in line with material culture evidence. We hypothesize that Upper Mesopotamia's cultural dynamism during the Neolithic Transition was the product not only of its fertile lands but also of its interregional demographic connections.

The diverse genetic origins of a Classical period Greek army

Trade and colonization caused an unprecedented increase in Mediterranean human mobility in the first millennium BCE. Often seen as a dividing force, warfare is in fact another catalyst of culture contact. We provide insight into the demographic dynamics of ancient warfare by reporting genome-wide data from fifth-century soldiers who fought for the army of the Greek Sicilian colony of Himera, along with representatives of the civilian population, nearby indigenous settlements, and 96 present-day individuals from Italy and Greece. Unlike the rest of the sample, many soldiers had ancestral origins in northern Europe, the Steppe, and the Caucasus. Integrating genetic, archaeological, isotopic, and historical data, these results illustrate the significant role mercenaries played in ancient Greek armies and highlight how participation in war contributed to continental-scale human mobility in the Classical world.

The genetic history of the Southern Arc: A bridge between West Asia and Europe

By sequencing 727 ancient individuals from the Southern Arc (Anatolia and its neighbors in Southeastern Europe and West Asia) over 10,000 years, we contextualize its Chalcolithic period and Bronze Age (about 5000 to 1000 BCE), when extensive gene flow entangled it with the Eurasian steppe. Two streams of migration transmitted Caucasus and Anatolian/Levantine ancestry northward, and the Yamnaya pastoralists, formed on the steppe, then spread southward into the Balkans and across the Caucasus into Armenia, where they left numerous patrilineal descendants. Anatolia was transformed by intra-West Asian gene flow, with negligible impact of the later Yamnaya migrations. This contrasts with all other regions where Indo-European languages were spoken, suggesting that the homeland of the Indo-Anatolian language family was in West Asia, with only secondary dispersals of non-Anatolian Indo-Europeans from the steppe.

A genetic probe into the ancient and medieval history of Southern Europe and West Asia

Literary and archaeological sources have preserved a rich history of Southern Europe and West Asia since the Bronze Age that can be complemented by genetics. Mycenaean period elites in Greece did not differ from the general population and included both people with some steppe ancestry and others, like the Griffin Warrior, without it. Similarly, people in the central area of the Urartian Kingdom around Lake Van lacked the steppe ancestry characteristic of the kingdom's northern provinces. Anatolia exhibited extraordinary continuity down to the Roman and Byzantine periods, with its people serving as the demographic core of much of the Roman Empire, including the city of Rome itself. During medieval times, migrations associated with Slavic and Turkic speakers profoundly affected the region.

Ancient mitochondrial diversity reveals population homogeneity in Neolithic Greece and identifies population dynamics along the Danubian expansion axis

The aim of the study is to investigate mitochondrial diversity in Neolithic Greece and its relation to hunter-gatherers and farmers who populated the Danubian Neolithic expansion axis. We sequenced 42 mitochondrial palaeogenomes from Greece and analysed them together with European set of 328 mtDNA sequences dating from the Early to the Final Neolithic and 319 modern sequences. To test for population continuity through time in Greece, we use an original structured population continuity test that simulates DNA from different periods by explicitly considering the spatial and temporal dynamics of populations. We explore specific scenarios of the mode and tempo of the European Neolithic expansion along the Danubian axis applying spatially explicit simulations coupled with Approximate Bayesian Computation. We observe a striking genetic homogeneity for the maternal line throughout the Neolithic in Greece whereas population continuity is rejected between the Neolithic and present-day Greeks. Along the Danubian expansion axis, our best-fitting scenario supports a substantial decrease in mobility and an increasing local hunter-gatherer contribution to the gene-pool of farmers following the initial rapid Neolithic expansion. Οur original simulation approach models key demographic parameters rather than inferring them from fragmentary data leading to a better understanding of this important process in European prehistory.

High-throughput degraded DNA sequencing of subfossil shells of a critically endangered stenoendemic land snail in the Aegean

High-throughput sequencing has enabled the comprehensive genetic exploration of biological diversity, especially by using natural history collections to study hard-to-find, threatened or even extinct-in-the-wild taxa. Mollusk shells are under-exploited as a source for DNA-based approaches, despite their apparent advantages in the field of conservation genetics. More particularly, degraded DNA techniques combined with high-throughput sequencing have never been used to gain insights about the DNA preservation in land snail subfossil or historical shells. Here, we applied degraded DNA analysis on two historical shells of Levantina rechingeri, a stenoendemic Critically Endangered species that has never been found alive, in order to explore the patterns of DNA preservation on land snail shells originating from the eastern Mediterranean, as well as to infer its molecular phylogenetic placement. Our results showed that centuries to decades-old DNA from an empty shell collected in an Aegean island exhibits characteristic post-mortem damage patterns similar to those observed in ancient DNA from eastern Mediterranean terrestrial animals, setting a precedent for future museomics studies on taxa distributed in areas with similar climate. Finally, genome skimming of the empty shell allowed high coverage of multiple nuclear and mitochondrial loci, enabling the phylogenetic placement of the focal taxon, the re-evaluation of its taxonomic classification, and the revealing of a new Aegean land snail lineage, Aristina genus novum. This approach is a non-invasive way to sample DNA from threatened land snail species and suitable to study the evolutionary history of taxa with cryptic ecology, stenoendemics, or extinct-in-the-wild, as well as old museum specimens.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

Airborne magnetic nanoparticles may contribute to COVID-19 outbreak: Relationships in Greece and Iran

This work attempts to shed light on whether the COVID-19 pandemic rides on airborne pollution. In particular, a two-city study provides evidence that PM_2.5 contributes to the timing and severity of the epidemic, without adjustment for confounders. The publicly available data of deaths between March and October 2020, updated it on May 30, 2021, and the average seasonal concentrations of PM_2.5 pollution over the previous years in Thessaloniki, the second-largest city of Greece, were investigated. It was found that changes in coronavirus-related deaths follow changes in air pollution and that the correlation between the two data sets is maximized at the lag time of one month. Similar data from Tehran were gathered for comparison. The results of this study underscore that it is possible, if not likely, that pollution nanoparticles are related to COVID-19 fatalities (Granger causality, p < 0.05), contributing to the understanding of the environmental impact on pandemics.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.