The spread of steppe and Iranian-related ancestry in the islands of the western Mediterranean

Urbanization and genetic homogenization in the medieval Low Countries revealed through a ten-century paleogenomic study of the city of Sint-Truiden

BACKGROUND

Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes.

RESULTS

We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes.

CONCLUSIONS

This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

Performance of qpAdm-based screens for genetic admixture on graph-shaped histories and stepping stone landscapes

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on 2D stepping stone landscapes and in situations with low prestudy odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture graph-shaped and stepping stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: (1) prestudy odds are low and fall rapidly with increasing model complexity; (2) complex migration networks violate the assumptions of the method; hence, there is poor correlation between qpAdm P-values and model optimality, contributing to low but nonzero false-positive rate and low power; and (3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly nonoptimal models have estimates in the same interval, contributing to the false-positive rate. We also reinterpret large sets of qpAdm models from 2 studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: (1) temporal stratification of targets and proxy sources in the case of admixture graph-shaped histories, (2) focused exploration of few models for increasing prestudy odds; and (3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false-positive rate.

High continuity of forager ancestry in the Neolithic period of the eastern Maghreb

Ancient DNA from the Mediterranean region has revealed long-range connections and population transformations associated with the spread of food-producing economies1-6. However, in contrast to Europe, genetic data from this key transition in northern Africa are limited, and have only been available from the far western Maghreb (Morocco)1-3. Here we present genome-wide data for nine individuals from the Later Stone Age through the Neolithic period from Algeria and Tunisia. The earliest individuals cluster with pre-Neolithic people of the western Maghreb (around 15,000-7,600 years before present (BP)), showing that this 'Maghrebi' ancestry profile had a substantial geographic and temporal extent. At least one individual from Djebba (Tunisia), dating to around 8,000 years BP, harboured ancestry from European hunter-gatherers, probably reflecting movement in the Early Holocene across the Strait of Sicily. Later Neolithic people from the eastern Maghreb retained largely local forager ancestry, together with smaller contributions from European farmers (by around 7,000 years BP) and Levantine groups (by around 6,800 years BP), and were thus far less impacted by external gene flow than were populations in other parts of the Neolithic Mediterranean.

Punic people were genetically diverse with almost no Levantine ancestors

The maritime Phoenician civilization from the Levant transformed the entire Mediterranean during the first millennium BCE1-3. However, the extent of human movement between the Levantine Phoenician homeland and Phoenician-Punic settlements in the central and western Mediterranean has been unclear in the absence of comprehensive ancient DNA studies. Here, we generated genome-wide data for 210 individuals, including 196 from 14 sites traditionally identified as Phoenician and Punic in the Levant, North Africa, Iberia, Sicily, Sardinia and Ibiza, and an early Iron Age individual from Algeria. Levantine Phoenicians made little genetic contribution to Punic settlements in the central and western Mediterranean between the sixth and second centuries BCE, despite abundant archaeological evidence of cultural, historical, linguistic and religious links4. Instead, these inheritors of Levantine Phoenician culture derived most of their ancestry from a genetic profile similar to that of Sicily and the Aegean. Much of the remaining ancestry originated from North Africa, reflecting the growing influence of Carthage5. However, this was a minority contributor of ancestry in all of the sampled sites, including in Carthage itself. Different Punic sites across the central and western Mediterranean show similar patterns of high genetic diversity. We also detect genetic relationships across the Mediterranean, reflecting shared demographic processes that shaped the Punic world.

PANE: fast and reliable ancestral reconstruction on ancient genotype data with non-negative least square and principal component analysis

The history of human populations has been strongly shaped by admixture events, contributing to patterns of observed genetic diversity across populations. In this study, we introduce the Principal component Ancestry proportions using NNLS Estimation (PANE) method that leverages principal component analysis and non-negative least squares to assess the ancestral compositions of admixed individuals given a large set of populations. Our results show its ability to reliably estimate ancestry across several scenarios, even those with a significant proportion of missing genotypes, in a fraction of the time required when using other tools.

Performance of qpAdm-based screens for genetic admixture on admixture-graph-shaped histories and stepping-stone landscapes

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on two-dimensional stepping-stone landscapes and in situations with low pre-study odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture-graph-shaped and stepping-stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: 1) pre-study odds are low and fall rapidly with increasing model complexity; 2) complex migration networks violate the assumptions of the method, hence there is poor correlation between qpAdm p-values and model optimality, contributing to low but non-zero false positive rate and low power; 3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly non-optimal models have estimates in the same interval, contributing to the false positive rate. We also re-interpret large sets of qpAdm models from two studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: 1) temporal stratification of targets and proxy sources in the case of admixture-graph-shaped histories; 2) focused exploration of few models for increasing pre-study odds; 3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false positive rate.

Patterns of population structure and genetic variation within the Saudi Arabian population

The Arabian Peninsula is considered the initial site of historic human migration out of Africa. The modern-day indigenous Arabians are believed to be the descendants who remained from the ancient split of the migrants into Eurasia. Here, we investigated how the population history and cultural practices such as endogamy have shaped the genetic variation of the Saudi Arabians. We genotyped 3,352 individuals and identified twelve genetic sub-clusters that corresponded to the geographical distribution of different tribal regions, differentiated by distinct components of ancestry based on comparisons to modern and ancient DNA references. These sub-clusters also showed variation across ranges of the genome covered in runs of homozygosity, as well as differences in population size changes over time. Using 25,488,981 variants found in whole genome sequencing data (WGS) from 302 individuals, we found that the Saudi tend to show proportionally more deleterious alleles than neutral alleles when compared to Africans/African Americans from gnomAD (e.g. a 13% increase of deleterious alleles annotated by AlphaMissense between 0.5 - 5% frequency in Saudi, compared to 7% decrease of the benign alleles; P < 0.001). Saudi sub-clusters with greater inbreeding and lower effective population sizes showed greater enrichment of deleterious alleles as well. Additionally, we found that approximately 10% of the variants discovered in our WGS data are not observed in gnomAD; these variants are also enriched with deleterious annotations. To accelerate studying the population-enriched deleterious alleles and their health consequences in this population, we made available the allele frequency estimates of 25,488,981 variants discovered in our samples. Taken together, our results suggest that Saudi's population history impacts its pattern of genetic variation with potential consequences to the population health. It further highlights the need to sequence diverse and unique populations so to provide a foundation on which to interpret medical- and pharmaco- genomic findings from these populations.

The cranium from the Octagon in Ephesos

During excavations in 1929, a well-preserved skeleton was discovered in a sarcophagus in the Octagon at Ephesos (Turkey). For the following century, archaeologists have speculated about the identity of this obviously notable person. Repeated claim is that the remains could represent Arsinoë IV, daughter of Ptolemy XII, and younger (half-)sister of Cleopatra VII. To address these questions we undertook state-of-the-art morphological, genetic and dating analyses of the cranium and further analyses of bone samples from a femur and a rib of the skeleton found in the same tomb. We confirm based on genetic analyses from the cranium and the femur that they derive from the same person. 14C-dating of the cranium provides a most likely time range between 205-36 BC. The connection with Arsinoë IV can be excluded because we confirmed that the individual is a male. The cranium represents an 11-14-year-old boy who suffered from significant developmental disturbances. Genetics suggest an ancestry from the Italian peninsula or Sardinia. The fate of the body of Arsinoë IV, who reportedly was killed in 41 BC in Ephesos, remains open. In contrast, investigations regarding the fate and social background of the boy from the Octagon can now proceed free of speculation.

A reduction in effective population size has not relaxed purifying selection in the human population of Eivissa (Balearic Islands)

Ibiza (Eivissa) is one of the main Balearic Islands in the western Mediterranean. Recent studies have highlighted the genetic distinctiveness of present-day Eivissans within the region and suggested it could be attributed to the genetic drift caused by recent demographic events. Whether this distinctiveness emerged from a differential demographic history, or rather from a bias for sampling in a small geographic region such as Eivissa, remains an open question, together with the understanding of the functional consequences of demography in the island. In order to clarify these questions and further characterize the distinctiveness of Eivissa within the Balearic and Mediterranean context, we generated whole exome sequences for 31 and 20 individuals from Eivissa and Menorca respectively, a subset of which were also genotyped with the Human Origins array. Our results show that Eivissans present signs of putatively recent genetic isolation that are shared to a lesser extent with Menorca such as more and longer runs of homozygosity and high numbers of intra-population shared IBD segments. Regarding the functional consequences of recent demography, although Eivissans do not present an excess of deleterious alleles or homozygotes comparing to other populations, genetic drift seems to have increased the allele frequencies of neutral and deleterious variants, which can have various medical implications.

Demographic history and genetic variation of the Armenian population

We introduce a sizable (n = 34) whole-genome dataset on Armenians, a population inhabiting the region in West Asia known as the Armenian highlands. Equipped with this genetic data, we conducted a whole-genome study of Armenians and deciphered their fine-scale population structure and complex demographic history. We demonstrated that the Armenian populations from western, central, and eastern parts of the highlands are relatively homogeneous. The Sasun, a population in the south that had been argued to have received a major genetic contribution from Assyrians, was instead shown to have derived its slightly divergent genetic profile from a bottleneck that occurred in the recent past. We also investigated the debated question on the genetic origin of Armenians and failed to find any significant support for historical suggestions by Herodotus of their Balkan-related ancestry. We checked the degree of continuity of modern Armenians with ancient inhabitants of the eastern Armenian highlands and detected a genetic input into the region from a source linked to Neolithic Levantine Farmers at some point after the Early Bronze Age. Additionally, we cataloged an abundance of new mutations unique to the population, including a missense mutation predicted to cause familial Mediterranean fever, an autoinflammatory disorder highly prevalent in Armenians. Thus, we highlight the importance of further genetic and medical studies of this population.

Ancient DNA challenges prevailing interpretations of the Pompeii plaster casts

The eruption of Somma-Vesuvius in 79 CE buried several nearby Roman towns, killing the inhabitants and burying under pumice lapilli and ash deposits a unique set of civil and private buildings, monuments, sculptures, paintings, and mosaics that provide a rich picture of life in the empire. The eruption also preserved the forms of many of the dying as the ash compacted around their bodies. Although the soft tissue decayed, the outlines of the bodies remained and were recovered by excavators centuries later by filling the cavities with plaster. From skeletal material embedded in the casts, we generated genome-wide ancient DNA and strontium isotopic data to characterize the genetic relationships, sex, ancestry, and mobility of five individuals. We show that the individuals' sexes and family relationships do not match traditional interpretations, exemplifying how modern assumptions about gendered behaviors may not be reliable lenses through which to view data from the past. For example, an adult wearing a golden bracelet with a child on their lap-often interpreted as mother and child-is genetically an adult male biologically unrelated to the child. Similarly, a pair of individuals who were thought to have died in an embrace-often interpreted as sisters-included at least one genetic male. All Pompeiians with genome-wide data consistently derive their ancestry largely from recent immigrants from the eastern Mediterranean, as has also been seen in contemporaneous ancient genomes from the city of Rome, underscoring the cosmopolitanism of the Roman Empire in this period.

Leveraging ancient DNA to uncover signals of natural selection in Europe lost due to admixture or drift

Large ancient DNA (aDNA) studies offer the chance to examine genomic changes over time, providing direct insights into human evolution. While recent studies have used time-stratified aDNA for selection scans, most focus on single-locus methods. We conducted a multi-locus genotype scan on 708 samples spanning 7000 years of European history. We show that the G12 statistic, originally designed for unphased diploid data, can effectively detect selection in aDNA processed to create 'pseudo-haplotypes'. In simulations and at known positive control loci (e.g., lactase persistence), G12 outperforms the allele frequency-based selection statistic, SweepFinder2, previously used on aDNA. Applying our approach, we identified 14 candidate regions of selection across four time periods, with half the signals detectable only in the earliest period. Our findings suggest that selective events in European prehistory, including from the onset of animal domestication, have been obscured by neutral processes like genetic drift and demographic shifts such as admixture.

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.

Low Genetic Impact of the Roman Occupation of Britain in Rural Communities

The Roman period saw the empire expand across Europe and the Mediterranean, including much of what is today Great Britain. While there is written evidence of high mobility into and out of Britain for administrators, traders, and the military, the impact of imperialism on local, rural population structure, kinship, and mobility is invisible in the textual record. The extent of genetic change that occurred in Britain during the Roman military occupation remains underexplored. Here, using genome-wide data from 52 ancient individuals from eight sites in Cambridgeshire covering the period of Roman occupation, we show low levels of genetic ancestry differentiation between Romano-British sites and indications of larger populations than in the Bronze Age and Neolithic. We find no evidence of long-distance migration from elsewhere in the Empire, though we do find one case of possible temporary mobility within a family unit during the Late Romano-British period. We also show that the present-day patterns of genetic ancestry composition in Britain emerged after the Roman period.

Ancient DNA sheds light on the funerary practices of late Neolithic collective burial in southern France

The Aven de la Boucle (Corconne, Gard, southern France) is a karst shaft used as a collective burial between 3600 and 2800 cal BCE. The site encompasses the skeletal remains of approximately 75 individuals comprising a large majority of adult individuals, represented by scattered and commingled remains. To date, few studies have explored the potential of ancient DNA to tackle the documentation of Neolithic collective burials, and the funerary selection rules within such structures remain largely debated. In this study, we combine genomic analysis of 37 individuals with archaeo-anthropological data and Bayesian modelling of radiocarbon dates. Through this multidisciplinary approach, we aim to characterize the identity of the deceased and their relationships, as well as untangle the genetic diversity and funerary dynamics of this community. Genomic results identify 76% of male Neolithic individuals, suggesting a marked sex-biased selection. Available data emphasize the importance of biological relatedness and a male-mediated transmission of social status, as the affiliation to a specific male-lineage appears as a preponderant selection factor. The genomic results argue in favour of 'continuous' deposits between 3600 and 2800 BCE, carried out by the same community, despite cultural changes reflected by the ceramic material.

Archaeology meets environmental genomics: implementing sedaDNA in the study of the human past

Sedimentary ancient DNA (sedaDNA) has become one of the standard applications in the field of paleogenomics in recent years. It has been used for paleoenvironmental reconstructions, detecting the presence of prehistoric species in the absence of macro remains and even investigating the evolutionary history of a few species. However, its application in archaeology has been limited and primarily focused on humans. This article argues that sedaDNA holds significant potential in addressing key archaeological questions concerning the origins, lifestyles, and environments of past human populations. Our aim is to facilitate the integration of sedaDNA into the standard workflows in archaeology as a transformative tool, thereby unleashing its full potential for studying the human past. Ultimately, we not only underscore the challenges inherent in the sedaDNA field but also provide a research agenda for essential enhancements needed for implementing sedaDNA into the archaeological workflow.

Late Neolithic collective burial reveals admixture dynamics during the third millennium BCE and the shaping of the European genome

The third millennium BCE was a pivotal period of profound cultural and genomic transformations in Europe associated with migrations from the Pontic-Caspian steppe, which shaped the ancestry patterns in the present-day European genome. We performed a high-resolution whole-genome analysis including haplotype phasing of seven individuals of a collective burial from ~2500 cal BCE and of a Bell Beaker individual from ~2300 cal BCE in the Paris Basin in France. The collective burial revealed the arrival in real time of steppe ancestry in France. We reconstructed the genome of an unsampled individual through its relatives' genomes, enabling us to shed light on the early-stage admixture patterns, dynamics, and propagation of steppe ancestry in Late Neolithic Europe. We identified two major Neolithic/steppe-related ancestry admixture pulses around 3000/2900 BCE and 2600 BCE. These pulses suggest different population expansion dynamics with striking links to the Corded Ware and Bell Beaker cultural complexes.

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

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.

Socio-cultural practices may have affected sex differences in stature in Early Neolithic Europe

The rules and structure of human culture impact health as much as genetics or environment. To study these relationships, we combine ancient DNA (n = 230), skeletal metrics (n = 391), palaeopathology (n = 606) and dietary stable isotopes (n = 873) to analyse stature variation in Early Neolithic Europeans from North Central, South Central, Balkan and Mediterranean regions. In North Central Europe, stable isotopes and linear enamel hypoplasias indicate high environmental stress across sexes, but female stature is low, despite polygenic scores identical to males, and suggests that cultural factors preferentially supported male recovery from stress. In Mediterranean populations, sexual dimorphism is reduced, indicating male vulnerability to stress and no strong cultural preference for males. Our analysis indicates that biological effects of sex-specific inequities can be linked to cultural influences at least as early as 7,000 yr ago, and culture, more than environment or genetics, drove height disparities in Early Neolithic Europe.

Folate gene polymorphisms CBS 844ins68 and RFC1 A80G and risk of Down syndrome offspring in young Iranian women: A cross-sectional study

Background

Cytogenetics and association studies showed that folate gene polymorphisms can increase the risk of chromosomal nondisjunction and aneuploidies. The folate-metabolizing gene polymorphisms in Down syndrome mothers (DSM) have been assessed in a variety of populations. Reduced folate carrier 1 (RFC1) and cystathionine beta-synthase (CBS) are key enzymes in folate metabolism.

Objective

2 common polymorphisms, CBS 844ins68 and RFC1 A80G, were analyzed to determine their probable risk for having Down syndrome (DS) babies in young mothers of Khuzestan province, Iran.

Materials and Methods

This study was conducted on 100 mothers who had trisomy 21 DS children. 100 age- and ethnic-matched mothers with at least 2 healthy children and no history of abnormal pregnancies were considered as control. The samples were collected from all the mothers from June 2019 to April 2021. Genomic DNA was extracted from peripheral blood. The CBS-844ins68 and RFC1-A80G were genotyped using polymerase chain reaction-electrophoresis and restriction fragment length polymorphism, respectively.

Results

The frequency of RFC1 AG and GG genotypes in DSM was significantly higher than the control mothers (odds ratio [OR] of 2.38 and 3.07, respectively). The heterozygote genotype of CBS 844ins68 was significantly more prevalent among DSM than the control (OR: 2.419). The OR was significantly increased to 6.667 when the homozygote of both variants was found together.

Conclusion

Studying polymorphisms possibly increases the susceptibility of having a DS child. However, ethnicity, nutrition, and epistatic interactions are considerable factors to be evaluated in future studies.

An individual with Sarmatian-related ancestry in Roman Britain

In the second century CE the Roman Empire had increasing contact with Sarmatians, nomadic Iranian speakers occupying an area stretching from the Pontic-Caspian steppe to the Carpathian mountains, both in the Caucasus and in the Danubian borders of the empire.1,2,3 In 175 CE, following their defeat in the Marcomannic Wars, emperor Marcus Aurelius drafted Sarmatian cavalry into Roman legions and deployed 5,500 Sarmatian soldiers to Britain, as recorded by contemporary historian Cassius Dio.4,5 Little is known about where the Sarmatian cavalry were stationed, and no individuals connected with this historically attested event have been identified to date, leaving its impact on Britain largely unknown. Here we document Caucasus- and Sarmatian-related ancestry in the whole genome of a Roman-period individual (126-228 calibrated [cal.] CE)-an outlier without traceable ancestry related to local populations in Britain-recovered from a farmstead site in present-day Cambridgeshire, UK. Stable isotopes support a life history of mobility during childhood. Although several scenarios are possible, the historical deployment of Sarmatians to Britain provides a parsimonious explanation for this individual's extraordinary life history. Regardless of the factors behind his migrations, these results highlight how long-range mobility facilitated by the Roman Empire impacted provincial locations outside of urban centers.

Multiproxy bioarchaeological data reveals interplay between growth, diet and population dynamics across the transition to farming in the central Mediterranean

The transition to farming brought on a series of important changes in human society, lifestyle, diet and health. The human bioarchaeology of the agricultural transition has received much attention, however, relatively few studies have directly tested the interrelationship between individual lifestyle factors and their implications for understanding life history changes among the first farmers. We investigate the interplay between skeletal growth, diet, physical activity and population size across 30,000 years in the central Mediterranean through a 'big data' cross-analysis of osteological data related to stature (n = 361), body mass (n = 334) and long bone biomechanics (n = 481), carbon (δ13C) and nitrogen (δ15N) stable isotopes (n = 1986 human, n = 475 animal) and radiocarbon dates (n = 5263). We present the observed trends on a continuous timescale in order to avoid grouping our data into assigned 'time periods', thus achieving greater resolution and chronological control over our analysis. The results identify important changes in human life history strategies associated with the first farmers, but also highlight the long-term nature of these trends in the millennia either side of the agricultural transition. The integration of these different data is an important step towards disentangling the complex relationship between demography, diet and health, and reconstruct life history changes within a southern European context. We believe the methodological approach adopted here has broader global implications for bioarchaeological studies of human adaptation more generally.

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.

High-coverage genome of the Tyrolean Iceman reveals unusually high Anatolian farmer ancestry

The Tyrolean Iceman is known as one of the oldest human glacier mummies, directly dated to 3350-3120 calibrated BCE. A previously published low-coverage genome provided novel insights into European prehistory, despite high present-day DNA contamination. Here, we generate a high-coverage genome with low contamination (15.3×) to gain further insights into the genetic history and phenotype of this individual. Contrary to previous studies, we found no detectable Steppe-related ancestry in the Iceman. Instead, he retained the highest Anatolian-farmer-related ancestry among contemporaneous European populations, indicating a rather isolated Alpine population with limited gene flow from hunter-gatherer-ancestry-related populations. Phenotypic analysis revealed that the Iceman likely had darker skin than present-day Europeans and carried risk alleles associated with male-pattern baldness, type 2 diabetes, and obesity-related metabolic syndrome. These results corroborate phenotypic observations of the preserved mummified body, such as high pigmentation of his skin and the absence of hair on his head.

A genetic history of continuity and mobility in the Iron Age central Mediterranean

The Iron Age was a dynamic period in central Mediterranean history, with the expansion of Greek and Phoenician colonies and the growth of Carthage into the dominant maritime power of the Mediterranean. These events were facilitated by the ease of long-distance travel following major advances in seafaring. We know from the archaeological record that trade goods and materials were moving across great distances in unprecedented quantities, but it is unclear how these patterns correlate with human mobility. Here, to investigate population mobility and interactions directly, we sequenced the genomes of 30 ancient individuals from coastal cities around the central Mediterranean, in Tunisia, Sardinia and central Italy. We observe a meaningful contribution of autochthonous populations, as well as highly heterogeneous ancestry including many individuals with non-local ancestries from other parts of the Mediterranean region. These results highlight both the role of local populations and the extreme interconnectedness of populations in the Iron Age Mediterranean. By studying these trans-Mediterranean neighbours together, we explore the complex interplay between local continuity and mobility that shaped the Iron Age societies of the central Mediterranean.

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.

Climbing the Giara: A quantitative reassessment of movement and visibility in the Nuragic landscape of the Gesturi plateau (South-Central Sardinia, Italy)

The built landscape of Nuragic Sardinia is an exceptional case for geostatistical analysis, allowing for a discussion of long-held assumptions and ideas. The function of nuraghi (ca. 1700-1100 BCE), the most prominent settled monument of the Sardinian Bronze Age, has been addressed via a multiplicity of landscape approaches, mainly relying on intuitive assessments of their spatial properties: nuraghi were assumed as means of territorial control. The series of nuraghi crowning the mesa plateau named Giara of Gesturi (South-Central Sardinia) provides a significant case for the study of their relations with visibility and movement. Context-oriented GIS models based on viewsheds and least-cost paths have been devised as targeted tools. The results show a certain correlation between nuraghi and potential movement on the slopes, thanks to the selection of plateau morphologies such as outward crests. Anyway, nuraghi do not stand exactly at the most accessible points of the plateau. Nuraghi offered ample visual control, especially at large distances, but not specifically over the closest accessible ways. This suggests that the function of nuraghi is somehow connected to defense and visibility, but it is not explained directly by local territorial control: a role as landmarks and multifaceted monuments has likely to be envisioned.

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.

Bioarchaeological and paleogenomic profiling of the unusual Neolithic burial from Grotta di Pietra Sant'Angelo (Calabria, Italy)

The Neolithic burial of Grotta di Pietra Sant'Angelo (CS) represents a unique archaeological finding for the prehistory of Southern Italy. The unusual placement of the inhumation at a rather high altitude and far from inhabited areas, the lack of funerary equipment and the prone deposition of the body find limited similarities in coeval Italian sites. These elements have prompted wider questions on mortuary customs during the prehistory of Southern Italy. This atypical case requires an interdisciplinary approach aimed to build an integrated bioarchaeological profile of the individual. The paleopathological investigation of the skeletal remains revealed the presence of numerous markers that could be associated with craft activities, suggesting possible interpretations of the individual's lifestyle. CT analyses, carried out on the maxillary bones, showed the presence of a peculiar type of dental wear, but also a good density of the bone matrix. Biomolecular and micromorphological analyses of dental calculus highlight the presence of a rich Neolithic-like oral microbiome, the composition of which is consistent with the presence pathologies. Finally, paleogenomic data obtained from the individual were compared with ancient and modern Mediterranean populations, including unpublished high-resolution genome-wide data for 20 modern inhabitants of the nearby village of San Lorenzo Bellizzi, which provided interesting insights into the biodemographic landscape of the Neolithic in Southern Italy.

Hunter-gatherer admixture facilitated natural selection in Neolithic European farmers

Ancient DNA has revealed multiple episodes of admixture in human prehistory during geographic expansions associated with cultural innovations. One important example is the expansion of Neolithic agricultural groups out of the Near East into Europe and their consequent admixture with Mesolithic hunter-gatherers.1,2,3,4 Ancient genomes from this period provide an opportunity to study the role of admixture in providing new genetic variation for selection to act upon, and also to identify genomic regions that resisted hunter-gatherer introgression and may thus have contributed to agricultural adaptations. We used genome-wide DNA from 677 individuals spanning Mesolithic and Neolithic Europe to infer ancestry deviations in the genomes of admixed individuals and to test for natural selection after admixture by testing for deviations from a genome-wide null distribution. We find that the region around the pigmentation-associated gene SLC24A5 shows the greatest overrepresentation of Neolithic local ancestry in the genome (|Z| = 3.46). In contrast, we find the greatest overrepresentation of Mesolithic ancestry across the major histocompatibility complex (MHC; |Z| = 4.21), a major immunity locus, which also shows allele frequency deviations indicative of selection following admixture (p = 1 × 10-56). This could reflect negative frequency-dependent selection on MHC alleles common in Neolithic populations or that Mesolithic alleles were positively selected for and facilitated adaptation in Neolithic populations to pathogens or other environmental factors. Our study extends previous results that highlight immune function and pigmentation as targets of adaptation in more recent populations to selection processes in the Stone Age.

Estimating human mobility in Holocene Western Eurasia with large-scale ancient genomic data

The recent increase in openly available ancient human DNA samples allows for large-scale meta-analysis applications. Trans-generational past human mobility is one of the key aspects that ancient genomics can contribute to since changes in genetic ancestry-unlike cultural changes seen in the archaeological record-necessarily reflect movements of people. Here, we present an algorithm for spatiotemporal mapping of genetic profiles, which allow for direct estimates of past human mobility from large ancient genomic datasets. The key idea of the method is to derive a spatial probability surface of genetic similarity for each individual in its respective past. This is achieved by first creating an interpolated ancestry field through space and time based on multivariate statistics and Gaussian process regression and then using this field to map the ancient individuals into space according to their genetic profile. We apply this algorithm to a dataset of 3138 aDNA samples with genome-wide data from Western Eurasia in the last 10,000 y. Finally, we condense this sample-wise record with a simple summary statistic into a diachronic measure of mobility for subregions in Western, Central, and Southern Europe. For regions and periods with sufficient data coverage, our similarity surfaces and mobility estimates show general concordance with previous results and provide a meta-perspective of genetic changes and human mobility.

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.

Phylogeographic and population genetic structure of hound-like native dogs of the Mediterranean Basin

The dog was probably the first domesticated animal. Despite extensive archaeological and genetic investigations, the origin and the evolution of the extant dogs are still being debated. Dog breeds that have over time been selected for hunting share common ancestral traits. This study represents the first comprehensive attempt to survey at the genomic and mitochondrial level eight hound-like dogs breeds indigenous to the Mediterranean Basin to determine if they share common ancient origins. Results from the microsatellite analysis indicate that all the dog populations have a low inbreeding value.The Kelb tal-Fenek has a high divergence from the current Egyptian street population, however there is not enough evidence from this study to exclude completely the potential of an ancient common relationship. Overall, the mitochondrial results indicate high frequencies of haplogroups A and B and a low representation of haplogroup C, while only one Egyptian dog could be assigned to haplogroup D. Results reveal identities and shared clades, suggesting the conservation of ancient European mitotypes in the Mediterranean hound-like breeds, especially in the Egyptian population. Although none of the dog populations/breeds participating in this study indicate to be direct descendants of the Egyptian dogs, they still have a very close morphologically resemblance to those iconic Egyptian dogs often depicted in ancient art forms and share some genetic links with the current Egyptian population. Further research is required with other markers such us complete mitogenomes and SNP panels to confirm the complex history of the Mediterranean dogs involved in this study.

de Pancorbo M, Calafell F. Y-chromosome target enrichment reveals rapid expansion of haplogroup R1b-DF27 in Iberia during the Bronze Age transition

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia.

Neolithic genomic data from southern France showcase intensified interactions with hunter-gatherer communities

Archaeological research shows that the dispersal of the Neolithic took a more complex turn when reaching western Europe, painting a contrasted picture of interactions between autochthonous hunter-gatherers (HGs) and incoming farmers. In order to clarify the mode, the intensity, and the regional variability of biological exchanges implied in these processes, we report new palaeogenomic data from Occitanie, a key region in Southern France. Genomic data from 28 individuals originating from six sites spanning from c. 5,500 to c. 2,500 BCE allow us to characterize regional patterns of ancestries throughout the Neolithic period. Results highlight major differences between the Mediterranean and Continental Neolithic expansion routes regarding both migration and interaction processes. High proportions of HG ancestry in both Early and Late Neolithic groups in Southern France support multiple pulses of inter-group gene flow throughout time and space and confirm the need for regional studies to address the complexity of the processes involved.

•

Genome-wide data from 28 individuals from Southern France (∼5,500–∼2,500 BCE)

•

Small groups associated with the Neolithic expansion along the Mediterranean

•

Early admixture between hunter-gatherers and Neolithic farmers in Southern France

•

Multiple pulses of HG legacy introgression in Western Europe throughout Neolithic

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.

Genetic structure and differentiation from early bronze age in the mediterranean island of sicily: Insights from ancient mitochondrial genomes

Sicily is one of the main islands of the Mediterranean Sea, and it is characterized by a variety of archaeological records, material culture and traditions, reflecting the history of migrations and populations' interaction since its first colonization, during the Paleolithic. These deep and complex demographic and cultural dynamics should have affected the genomic landscape of Sicily at different levels; however, the relative impact of these migrations on the genomic structure and differentiation within the island remains largely unknown. The available Sicilian modern genetic data gave a picture of the current genetic structure, but the paucity of ancient data did not allow so far to make predictions about the level of historical variation. In this work, we sequenced and analyzed the complete mitochondrial genomes of 36 individuals from five different locations in Sicily, spanning from Early Bronze Age to Iron Age, and with different cultural backgrounds. The comparison with coeval groups from the Mediterranean Basin highlighted structured genetic variation in Sicily since Early Bronze Age, thus supporting a demic impact of the cultural transitions within the Island. Explicit model testing through Approximate Bayesian Computation allowed us to make predictions about the origin of Sicanians, one of the three indigenous peoples of Sicily, whose foreign origin from Spain, historically attributed, was not confirmed by our analysis of genetic data. Sicilian modern mitochondrial data show a different, more homogeneous, genetic composition, calling for a recent genetic replacement in the Island of pre-Iron Age populations, that should be further investigated.

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.

Temporal population structure, a genetic dating method for ancient Eurasian genomes from the past 10,000 years

Radiocarbon dating is the gold standard in archeology to estimate the age of skeletons, a key to studying their origins. Many published ancient genomes lack reliable and direct dates, which results in obscure and contradictory reports. We developed the temporal population structure (TPS), a DNA-based dating method for genomes ranging from the Late Mesolithic to today, and applied it to 3,591 ancient and 1,307 modern Eurasians. TPS predictions aligned with the known dates and correctly accounted for kin relationships. TPS dating of poorly dated Eurasian samples resolved conflicting reports in the literature, as illustrated by one test case. We also demonstrated how TPS improved the ability to study phenotypic traits over time. TPS can be used when radiocarbon dating is unfeasible or uncertain or to develop alternative hypotheses for samples younger than 10,000 years ago, a limitation that may be resolved over time as ancient data accumulate.

Dairying, diseases and the evolution of lactase persistence in Europe

In European and many African, Middle Eastern and southern Asian populations, lactase persistence (LP) is the most strongly selected monogenic trait to have evolved over the past 10,000 years1. Although the selection of LP and the consumption of prehistoric milk must be linked, considerable uncertainty remains concerning their spatiotemporal configuration and specific interactions2,3. Here we provide detailed distributions of milk exploitation across Europe over the past 9,000 years using around 7,000 pottery fat residues from more than 550 archaeological sites. European milk use was widespread from the Neolithic period onwards but varied spatially and temporally in intensity. Notably, LP selection varying with levels of prehistoric milk exploitation is no better at explaining LP allele frequency trajectories than uniform selection since the Neolithic period. In the UK Biobank4,5 cohort of 500,000 contemporary Europeans, LP genotype was only weakly associated with milk consumption and did not show consistent associations with improved fitness or health indicators. This suggests that other reasons for the beneficial effects of LP should be considered for its rapid frequency increase. We propose that lactase non-persistent individuals consumed milk when it became available but, under conditions of famine and/or increased pathogen exposure, this was disadvantageous, driving LP selection in prehistoric Europe. Comparison of model likelihoods indicates that population fluctuations, settlement density and wild animal exploitation-proxies for these drivers-provide better explanations of LP selection than the extent of milk exploitation. These findings offer new perspectives on prehistoric milk exploitation and LP evolution.

Ancient Maltese genomes and the genetic geography of Neolithic Europe

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

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.

The spatiotemporal patterns of major human admixture events during the European Holocene

Recent studies have shown that admixture has been pervasive throughout human history. While several methods exist for dating admixture in contemporary populations, they are not suitable for sparse, low coverage ancient genomic data. Thus, we developed DATES (Distribution of Ancestry Tracts of Evolutionary Signals) that leverages ancestry covariance patterns across the genome of a single individual to infer the timing of admixture. DATES provides reliable estimates under various demographic scenarios and outperforms available methods for ancient DNA applications. Using DATES on~1100 ancient genomes from sixteen regions in Europe and west Asia, we reconstruct the chronology of the formation of the ancestral populations and the fine-scale details of the spread of Neolithic farming and Steppe pastoralist-related ancestry across Europe. By studying the genetic formation of Anatolian farmers, we infer that gene flow related to Iranian Neolithic farmers occurred before 9600 BCE, predating the advent of agriculture in Anatolia. Contrary to the archaeological evidence, we estimate that early Steppe pastoralist groups (Yamnaya and Afanasievo) were genetically formed more than a millennium before the start of Steppe pastoralism. Our analyses provide new insights on the origins and spread of farming and Indo-European languages, highlighting the power of genomic dating methods to elucidate the legacy of human migrations.

Bioarchaeological and palaeogenomic portrait of two Pompeians that died during the eruption of Vesuvius in 79 AD

The archaeological site of Pompeii is one of the 54 UNESCO World Heritage sites in Italy, thanks to its uniqueness: the town was completely destroyed and buried by a Vesuvius’ eruption in 79 AD. In this work, we present a multidisciplinary approach with bioarchaeological and palaeogenomic analyses of two Pompeian human remains from the Casa del Fabbro. We have been able to characterize the genetic profile of the first Pompeian’ genome, which has strong affinities with the surrounding central Italian population from the Roman Imperial Age. Our findings suggest that, despite the extensive connection between Rome and other Mediterranean populations, a noticeable degree of genetic homogeneity exists in the Italian peninsula at that time. Moreover, palaeopathological analyses identified the presence of spinal tuberculosis and we further investigated the presence of ancient DNA from Mycobacterium tuberculosis. In conclusion, our study demonstrates the power of a combined approach to investigate ancient humans and confirms the possibility to retrieve ancient DNA from Pompeii human remains. Our initial findings provide a foundation to promote an intensive and extensive paleogenetic analysis in order to reconstruct the genetic history of population from Pompeii, a unique archaeological site.

Ancient DNA gives new insights into a Norman Neolithic monumental cemetery dedicated to male elites

By integrating genomic and archaeological data, we provide new insights into the Neolithic French monumental site of Fleury-sur-Orne in Normandy, where a group of selected individuals was buried in impressively long monuments. The earliest individuals buried at Fleury-sur-Orne match the expected western European Neolithic genetic diversity, while three individuals, designated as genetic outliers, were buried after 4,000 calibrated BCE. We hypothesize that different, unrelated families or clans used the site over several centuries. Thirteen of 14 of the analyzed individuals were male, indicating an overarching patrilineal system. However, one exception, a female buried with a symbolically male artifact, suggests that the embodiment of the male gender in death was required to access burial at the monumental structures.

The Middle Neolithic in western Europe is characterized by monumental funerary structures, known as megaliths, along the Atlantic façade. The first manifestations of this phenomenon occurred in modern-day France with the long mounds of the Cerny culture. Here, we present genome-wide data from the fifth-millennium BCE site of Fleury-sur-Orne in Normandy (France), famous for its impressively long monuments built for selected individuals. The site encompasses 32 monuments of variable sizes, containing the burials of 19 individuals from the Neolithic period. To address who was buried at the site, we generated genome-wide data for 14 individuals, of whom 13 are males, completing previously published data [M. Rivollat et al., Sci. Adv. 6, eaaz5344 (2020)]. Population genetic and Y chromosome analyses show that the Fleury-sur-Orne group fits within western European Neolithic genetic diversity and that the arrival of a new group is detected after 4,000 calibrated BCE. The results of analyzing uniparentally inherited markers and an overall low number of long runs of homozygosity suggest a patrilineal group practicing female exogamy. We find two pairs of individuals to be father and son, buried together in the same monument/grave. No other biological relationship can link monuments together, suggesting that each monument was dedicated to a genetically independent lineage. The combined data and documented father–son line of descent suggest a male-mediated transmission of sociopolitical authority. However, a single female buried with an arrowhead, otherwise considered a symbol of power of the male elite of the Cerny culture, questions a strictly biological sex bias in the burial rites of this otherwise “masculine” monumental cemetery.

Genomic and dietary discontinuities during the Mesolithic and Neolithic in Sicily

Sicily is a key region for understanding the agricultural transition in the Mediterranean because of its central position. Here, we present genomic and stable isotopic data for 19 prehistoric Sicilians covering the Mesolithic to Bronze Age periods (10,700–4,100 yBP). We find that Early Mesolithic hunter-gatherers (HGs) from Sicily are a highly drifted lineage of the Early Holocene western European HGs, whereas Late Mesolithic HGs carry ∼20% ancestry related to northern and (south) eastern European HGs, indicating substantial gene flow. Early Neolithic farmers are genetically most similar to farmers from the Balkans and Greece, with only ∼7% of ancestry from local Mesolithic HGs. The genetic discontinuities during the Mesolithic and Early Neolithic match the changes in material culture and diet. Three outlying individuals dated to ∼8,000 yBP; however, suggest that hunter-gatherers interacted with incoming farmers at Grotta dell’Uzzo, resulting in a mixed economy and diet for a brief interlude at the Mesolithic-Neolithic transition.

•

Genetic transition between Early Mesolithic and Late Mesolithic hunter-gatherers

•

A near-complete genetic turnover during the Mesolithic-Neolithic transition

•

Exchange of subsistence practices between hunter-gatherers and early farmers

Origin and mobility of Iron Age Gaulish groups in present-day France revealed through archaeogenomics

The Iron Age period occupies an important place in French history because the Gauls are regularly presented as the direct ancestors of the extant French population. We documented here the genomic diversity of Iron Age communities originating from six French regions. The 49 acquired genomes permitted us to highlight an absence of discontinuity between Bronze Age and Iron Age groups in France, lending support to a cultural transition linked to progressive local economic changes rather than to a massive influx of allochthone groups. Genomic analyses revealed strong genetic homogeneity among the regional groups associated with distinct archaeological cultures. This genomic homogenization appears to be linked to individuals' mobility between regions and gene flow with neighbouring groups from England and Spain. Thus, the results globally support a common genomic legacy for the Iron Age population of modern-day France that could be linked to recurrent gene flow between culturally differentiated communities.