Ancient genomes reveal origin and rapid trans-Eurasian migration of 7th century Avar elites

Ancient DNA reveals reproductive barrier despite shared Avar-period culture

After a long-distance migration, Avars with Eastern Asian ancestry arrived in Eastern Central Europe in 567 to 568 CE and encountered groups with very different European ancestry1,2. We used ancient genome-wide data of 722 individuals and fine-grained interdisciplinary analysis of large seventh- to eighth-century CE neighbouring cemeteries south of Vienna (Austria) to address the centuries-long impact of this encounter1,2. We found that even 200 years after immigration, the ancestry at one site (Leobersdorf) remained dominantly East Asian-like, whereas the other site (Mödling) shows local, European-like ancestry. These two nearby sites show little biological relatedness, despite sharing a distinctive late-Avar culture3,4. We reconstructed six-generation pedigrees at both sites including up to 450 closely related individuals, allowing per-generation demographic profiling of the communities. Despite different ancestry, these pedigrees together with large networks of distant relatedness show absence of consanguinity, patrilineal pattern with female exogamy, multiple reproductive partnerships (for example, levirate) and direct correlation of biological connectivity with archaeological markers of social status. The generation-long genetic barrier was maintained by systematically choosing partners with similar ancestry from other sites in the Avar realm. Leobersdorf had more biological connections with the Avar heartlands than with Mödling, which is instead linked to another site from the Vienna Basin with European-like ancestry. Mobility between sites was mostly due to female exogamy pointing to different marriage networks as the main driver of the maintenance of the genetic barrier.

High-resolution genomic history of early medieval Europe

Many known and unknown historical events have remained below detection thresholds of genetic studies because subtle ancestry changes are challenging to reconstruct. Methods based on shared haplotypes1,2 and rare variants3,4 improve power but are not explicitly temporal and have not been possible to adopt in unbiased ancestry models. Here we develop Twigstats, an approach of time-stratified ancestry analysis that can improve statistical power by an order of magnitude by focusing on coalescences in recent times, while remaining unbiased by population-specific drift. We apply this framework to 1,556 available ancient whole genomes from Europe in the historical period. We are able to model individual-level ancestry using preceding genomes to provide high resolution. During the first half of the first millennium CE, we observe at least two different streams of Scandinavian-related ancestry expanding across western, central and eastern Europe. By contrast, during the second half of the first millennium CE, ancestry patterns suggest the regional disappearance or substantial admixture of these ancestries. In Scandinavia, we document a major ancestry influx by approximately 800 CE, when a large proportion of Viking Age individuals carried ancestry from groups related to central Europe not seen in individuals from the early Iron Age. Our findings suggest that time-stratified ancestry analysis can provide a higher-resolution lens for genetic history.

Ancient genomes reveal Avar-Hungarian transformations in the 9th-10th centuries CE Carpathian Basin

During the Early Medieval period, the Carpathian Basin witnessed substantial demographic shifts, notably under the Avar dominance for ~250 years, followed by the settlement of early Hungarians in the region during the late 9th century CE. This study presents the genetic analysis of 296 ancient samples, including 103 shotgun-sequenced genomes, from present-day Western Hungary. By using identity-by-descent segment sharing networks, this research offers detailed insights into the population structure and dynamics of the region from the 5th to 11th centuries CE, with specific focus on certain microregions. Our evaluations reveal spatially different histories in Transdanubia even between communities in close geographical proximity, highlighting the importance of dense sampling and analyses. Our findings highlight extensive homogenization and reorganization processes, as well as discontinuities between Hun, Avar, and Hungarian conquest period immigrant groups, alongside the spread and integration of ancestry related to the Hungarian conquerors.

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.

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 in exploratory archaeogenetic studies for finding optimal admixture models of population history. Despite its popularity, qpAdm remains untested on histories in the form of admixture graphs of random topology or stepping-stone landscapes. We analyzed data from such simulations and found that while for admixture-graph-shaped histories there exist simple solutions (temporal stratification) for minimizing false findings of gene flow, in the case of stepping-stone landscapes the method generates results that do not appear suspect but are misleading: feasible qpAdm models are either accurate but simplistic in the context of landscapes, or highly inaccurate in the case of multi-component models. This is largely is due to two reasons: 1) because of complex migration networks that violate the assumptions of the method, there is poor correlation between qpAdm p -values and model optimality in many sections of the parameter space; 2) admixture fraction estimates between 0 and 1 are largely restricted to symmetric source configurations around targets, hence popular [0, 1] model plausibility criteria confound analyses of landscape-type demographies, unless their interpretations are explicitly spatial. For many species/regions/periods archaeogenetic sampling is very sparse and may be random with respect to population density of ancient individuals. In this situation only a specific combination of landscape properties and feasibility criteria allows to efficiently reject highly asymmetric non-optimal models most abundant in random deme sets. This problem may obscure useful signal (rare optimal models) and might be responsible for some claims about rapid long-distance migrations in the literature.

Long shared haplotypes identify the Southern Urals as a primary source for the 10th century Hungarians

During the Hungarian Conquest in the 10th century CE, the early medieval Magyars, a group of mounted warriors from Eastern Europe, settled in the Carpathian Basin. They likely introduced the Hungarian language to this new settlement area, during an event documented by both written sources and archaeological evidence. Previous archaeogenetic research identified the newcomers as migrants from the Eurasian steppe. However, genome-wide ancient DNA from putative source populations has not been available to test alternative theories of their precise source. We generated genome-wide ancient DNA data for 131 individuals from candidate archaeological contexts in the Circum-Uralic region in present-day Russia. Our results tightly link the Magyars to people of the Early Medieval Karayakupovo archaeological horizon on both the European and Asian sides of the southern Urals. Our analyes show that ancestors of the people of the Karayakupovo archaeological horizon were established in the Southern Urals by the Iron Age and that their descendants persisted locally in the Volga-Kama region until at least the 14th century.

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.

Network of large pedigrees reveals social practices of Avar communities

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

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.

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.

More than one millennium (2nd-16th century CE) of the White Plague in the Carpathian Basin - New cases, expanding knowledge

The causative agent of tuberculosis is still a widespread pathogen, which caused the death of ca. 1.6 million people globally in 2021. The paleopathological study of human remains revealed the antiquity of the disease and its continuous presence throughout the history of humankind. The Carpathian Basin has always been a biocultural melting pot, since it has seen several migrations over the centuries, and served as a location of admixture and interaction for numerous populations of different cultures. Thus, this geographical territory is ideal for the examination of the coevolutionary processes of hosts and their pathogens. We aimed to reveal the spatial and temporal distribution of tuberculosis cases excavated inside the borders of Hungary between the 2nd and 16th centuries CE. We established a comprehensive database by collecting 114 already published cases and introducing 39 new cases. The involved cases include those that have been confirmed by different molecular methods, as well as possible infections that were identified based on the presence of macromorphological and radiological alterations. The progress of future molecular and paleopathological studies can be facilitated by our dataset, as it presents spatial and temporal information concerning the spread of the disease in the Carpathian Basin, as well as the biological profile and detailed paleopathological description of lesions illustrated by photo- and radiographs.

The dahliagram: An interdisciplinary tool for investigation, visualization, and communication of past human-environmental interaction

Investigation into the nexus of human-environmental behavior has seen increasing collaboration of archaeologists, historians, and paleo-scientists. However, many studies still lack interdisciplinarity and overlook incompatibilities in spatiotemporal scaling of environmental and societal data and their uncertainties. Here, we argue for a strengthened commitment to collaborative work and introduce the "dahliagram" as a tool to analyze and visualize quantitative and qualitative knowledge from diverse disciplinary sources and epistemological backgrounds. On the basis of regional cases of past human mobility in eastern Africa, Inner Eurasia, and the North Atlantic, we develop three dahliagrams that illustrate pull and push factors underlying key phases of population movement across different geographical scales and over contrasting periods of time since the end of the last Ice Age. Agnostic to analytical units, dahliagrams offer an effective tool for interdisciplinary investigation, visualization, and communication of complex human-environmental interactions at a diversity of spatiotemporal scales.

Mitogenome analysis reveals predominantly ancient Yellow River origin of population inhabiting Datong agro-pastoral ecotone along Great Wall

The Datong Basin was an important arena for population movement and admixture between the Yellow River Valley and Eastern Steppe. In historical materials, the region was often the setting for a tug-of-war between Han farmers and non-Han nomads. The genetic makeup and population history of this Datong population has, however, remained uncertain. In this study, we analysed 289 mitogenomes from Datong individuals. Our primary findings were: (1) population summary statistics analysis revealed a high level of genetic diversity and strong signals of population expansion in the Datong population; (2) inter-population comparisons (PCA and Fst heatmap) exhibited a close clustering between the Datong population and Northern Han, especially northern frontier groups, such as the Inner Mongolia Han, Heilongjiang Han, Liaoning Han and Tianjin Han; (3) phylogeographic analysis of complete mitogenomes revealed the presence of different components in the maternal gene pools of Datong population-the northern East Asian component was dominant (66.44%), whereas the southern East Asians were the second largest component with 31.49%. We also observed a much reduced west Eurasian (2.07%) component; (4) direct comparisons with ancient groups showed closer relationship between Datong and Yellow River farmers than Eastern Steppe nomads. Despite, therefore, centuries of Eastern Steppe nomadic control over the Datong area, Yellow River farmers had a much more significant impact on the Datong population.

Genetic Analysis of Mingrelians Reveals Long-Term Continuity of Populations in Western Georgia (Caucasus)

To elucidate the population history of the Caucasus, we conducted a survey of genetic diversity in Samegrelo (Mingrelia), western Georgia. We collected DNA samples and genealogical information from 485 individuals residing in 30 different locations, the vast majority of whom being Mingrelian speaking. From these DNA samples, we generated mitochondrial DNA (mtDNA) control region sequences for all 485 participants (female and male), Y-short tandem repeat haplotypes for the 372 male participants, and analyzed all samples at nearly 590,000 autosomal single nucleotide polymorphisms (SNPs) plus around 33,000 on the sex chromosomes, with 27,000 SNP removed for missingness, using the GenoChip 2.0+ microarray. The resulting data were compared with those from populations from Anatolia, the Caucasus, the Near East, and Europe. Overall, Mingrelians exhibited considerable mtDNA haplogroup diversity, having high frequencies of common West Eurasian haplogroups (H, HV, I, J, K, N1, R1, R2, T, U, and W. X2) and low frequencies of East Eurasian haplogroups (A, C, D, F, and G). From a Y-chromosome standpoint, Mingrelians possessed a variety of haplogroups, including E1b1b, G2a, I2, J1, J2, L, Q, R1a, and R1b. Analysis of autosomal SNP data further revealed that Mingrelians are genetically homogeneous and cluster with other modern-day South Caucasus populations. When compared with ancient DNA samples from Bronze Age archaeological contexts in the broader region, these data indicate that the Mingrelian gene pool began taking its current form at least by this period, probably in conjunction with the formation of a distinct linguistic community.

Simulation-based Benchmarking of Ancient Haplotype Inference for Detecting Population Structure

Paleogenomic data has informed us about the movements, growth, and relationships of ancient populations. It has also given us context for medically relevant adaptations that appear in present-day humans due to introgression from other hominids, and it continues to help us characterize the evolutionary history of humans. However, ancient DNA (aDNA) presents several practical challenges as various factors such as deamination, high fragmentation, environmental contamination of aDNA, and low amounts of recoverable endogenous DNA, make aDNA recovery and analysis more difficult than modern DNA. Most studies with aDNA leverage only SNP data, and only a few studies have made inferences on human demographic history based on haplotype data, possibly because haplotype estimation (or phasing) has not yet been systematically evaluated in the context of aDNA. Here, we evaluate how the unique challenges of aDNA can impact phasing quality. We also develop a software tool that simulates aDNA taking into account the features of aDNA as well as the evolutionary history of the population. We measured phasing error as a function of aDNA quality and demographic history, and found that low phasing error is achievable even for very ancient individuals (~ 400 generations in the past) as long as contamination and read depth are adequate. Our results show that population splits or bottleneck events occurring between the reference and phased populations affect phasing quality, with bottlenecks resulting in the highest average error rates. Finally, we found that using estimated haplotypes, even if not completely accurate, is superior to using the simulated genotype data when reconstructing changes in population structure after population splits between present-day and ancient populations.

Fine-scale sampling uncovers the complexity of migrations in 5th-6th century Pannonia

As the collapse of the Western Roman Empire accelerated during the 4th and 5th centuries, arriving "barbarian" groups began to establish new communities in the border provinces of the declining (and eventually former) empire. This was a time of significant cultural and political change throughout not only these border regions but Europe as a whole.1,2 To better understand post-Roman community formation in one of these key frontier zones after the collapse of the Hunnic movement, we generated new paleogenomic data for a set of 38 burials from a time series of three 5th century cemeteries3,4,5 at Lake Balaton, Hungary. We utilized a comprehensive sampling approach to characterize these cemeteries along with data from 38 additional burials from a previously published mid-6th century site6 and analyzed them alongside data from over 550 penecontemporaneous individuals.7,8,9,10,11,12,13,14,15,16,17,18,19 The range of genetic diversity in all four of these local burial communities is extensive and wider ranging than penecontemporaneous Europeans sequenced to date. Despite many commonalities in burial customs and demography, we find that there were substantial differences in genetic ancestry between the sites. We detect evidence of northern European gene flow into the Lake Balaton region. Additionally, we observe a statistically significant association between dress artifacts and genetic ancestry among 5th century genetically female burials. Our analysis shows that the formation of early Medieval communities was a multifarious process even at a local level, consisting of genetically heterogeneous groups.

The Rurikids: The First Experience of Reconstructing the Genetic Portrait of the Ruling Family of Medieval Rus' Based on Paleogenomic Data

The Rurikids were the reigning house of Rus', its principalities and, ultimately the Tsardom of Russia, for seven centuries: from the IX to the end of the XVI century. According to the Primary Chronicle (the Tale of Bygone Years), the main chronicle of Rus', the Rurik dynasty was founded by the Varangian prince Rurik, invited to reign in Novgorod in 862, but still there is no direct genetic evidence of the origin of the early Rurikids. This research, for the first time, provides a genome-wide paleogenetic analysis of bone remains belonging to one of the Rurikids, Prince Dmitry Alexandrovich (?-1294), the son of the Grand Prince of Vladimir Alexander Yaroslavich Nevsky (1221-1263). It has been established that his Y chromosome belongs to the N1a haplogroup. Most of the modern Rurikids, according to their genealogies, belonging to the N1a haplogroup, have the most similar variants of Y chromosomes to each other, as well as to the Y chromosome of Prince Dmitry Alexandrovich. Genome-wide data of the medieval and modern Rurikids unequivocally indicates that they belong to the N1a haplogroup of the Y chromosome, starting at least from the XI century (since the time of Prince Yaroslav the Wise). All the other alleged Rurikids, both ancient and modern, being carriers of other haplogroups (R1a, I2a), possess high heterogeneity of the sequence of Y chromosomes, meaning that we cannot confirm their common ancestry. The most probable ancestors of Prince Dmitry Alexandrovich in the male line were the men who left the burial ground Bolshoy Oleny Island on the coast of the Kola Peninsula about 3,600 years ago. The reconstruction of the genome of Prince Dmitry Alexandrovich indicates the contribution of three ancestral components to his origin: (1) the early medieval population of the east of Scandinavia from the island of Oland, (2) representatives of the steppe nomadic peoples of the Eurasian steppes of the Iron Age or the early medieval population of central Europe (steppe nomads from the territory of Hungary), and (3) the ancient East-Eurasian component. Reliable statistics were also obtained when the Scandinavians were replaced with the Medieval Russian Slavic populations of the XI century. Thus, for the first time, we have shown the complex nature of interethnic interactions in the formation of the nobility of medieval Rus' on the example of the ancient Rurikid.

Integrating Linguistic, Archaeological and Genetic Perspectives Unfold the Origin of Ugrians

In the last year two publications shed new light on the linguistic and genomic history of ancient Uralic speakers. Here I show that these novel genetic and linguistic data are compatible with each-other and with the archaeological inferences, allowing us to formulate a very plausible hypothesis about the prehistory of Ugric speakers. Both genetic and archaeological data indicate the admixture of the Mezhovskaya population with northern forest hunters in the late Bronze Age, which gave rise to a "proto-Ugric" community. This finding is consistent with the linguistic reconstruction of the proto-Ugric language. Genetic data indicate an admixture of proto-Hungarians with early Sarmatians and early Huns, and I show that the first admixture can be reconciled with the formation of the Gorokhovo culture and its integration into the early Sarmatian Prokhorovka culture, while the second admixture corresponds to the transformation of the Sargat and Sarmatian cultures due to Xiongnu invasions.

High Coverage Mitogenomes and Y-Chromosomal Typing Reveal Ancient Lineages in the Modern-Day Székely Population in Romania

Here we present 115 whole mitogenomes and 92 Y-chromosomal Short Tandem Repeat (STR) and Single Nucleotide Polymorphism (SNP) profiles from a Hungarian ethnic group, the Székelys (in Romanian: Secuii, in German: Sekler), living in southeast Transylvania (Romania). The Székelys can be traced back to the 12th century in the region, and numerous scientific theories exist as to their origin. We carefully selected sample providers that had local ancestors inhabiting small villages in the area of Odorheiu Secuiesc/Székelyudvarhely in Romania. The results of our research and the reported data signify a qualitative leap compared to previous studies since it presents the first complete mitochondrial DNA sequences and Y-chromosomal profiles of 23 STRs from the region. We evaluated the results with population genetic and phylogenetic methods in the context of the modern and ancient populations that are either geographically or historically related to the Székelys. Our results demonstrate a predominantly local uniparental make-up of the population that also indicates limited admixture with neighboring populations. Phylogenetic analyses confirmed the presumed eastern origin of certain maternal (A, C, D) and paternal (Q, R1a) lineages, and, in some cases, they could also be linked to ancient DNA data from the Migration Period (5th-9th centuries AD) and Hungarian Conquest Period (10th century AD) populations.

Historical human migrations: From the steppe to the basin

Many migrations during human history have made the Carpathian Basin the melting pot of Europe. New ancient genomes confirm the Asian origin of European Huns, Avars and Magyars and huge within-group variability that is linked with social structure.

The origin of the Carpathian Avar elites revealed

In the sixth century, the Avar elites established themselves in the Carpathian Basin as evidenced today by the astounding archaeological sites with exquisite grave goods made of silver and gold. Recently reported in Cell, Gnecchi-Ruscone et al.¹ obtained paleogenomic evidence from Avar sites in Hungary, placing the Avar elites' origin in Mongolia and confirming their own historical claims.