The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

Repeated plague infections across six generations of Neolithic Farmers.

In the period between 5,300 and 4,900 calibrated years before present (cal. BP), populations across large parts of Europe underwent a period of demographic decline1,2. However, the cause of this so-called Neolithic decline is still debated. Some argue for an agricultural crisis resulting in the decline3, others for the spread of an early form of plague4. Here we use population-scale ancient genomics to infer ancestry, social structure and pathogen infection in 108 Scandinavian Neolithic individuals from eight megalithic graves and a stone cist. We find that the Neolithic plague was widespread, detected in at least 17% of the sampled population and across large geographical distances. We demonstrate that the disease spread within the Neolithic community in three distinct infection events within a period of around 120 years. Variant graph-based pan-genomics shows that the Neolithic plague genomes retained ancestral genomic variation present in Yersinia pseudotuberculosis, including virulence factors associated with disease outcomes. In addition, we reconstruct four multigeneration pedigrees, the largest of which consists of 38 individuals spanning six generations, showing a patrilineal social organization. Lastly, we document direct genomic evidence for Neolithic female exogamy in a woman buried in a different megalithic tomb than her brothers. Taken together, our findings provide a detailed reconstruction of plague spread within a large patrilineal kinship group and identify multiple plague infections in a population dated to the beginning of the Neolithic decline.

Temporal dynamics of woolly mammoth genome erosion prior to extinction.

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Population Genomic History of the Endangered Anatolian and Cyprian Mouflons in Relation to Worldwide Wild, Feral, and Domestic Sheep Lineages.

Once widespread in their homelands, the Anatolian mouflon (Ovis gmelini anatolica) and the Cyprian mouflon (Ovis gmelini ophion) were driven to near extinction during the 20th century and are currently listed as endangered populations by the International Union for Conservation of Nature. While the exact origins of these lineages remain unclear, they have been suggested to be close relatives of domestic sheep or remnants of proto-domestic sheep. Here, we study whole genome sequences of n = 5 Anatolian mouflons and n = 10 Cyprian mouflons in terms of population history and diversity, comparing them with eight other extant sheep lineages. We find reciprocal genetic affinity between Anatolian and Cyprian mouflons and domestic sheep, higher than all other studied wild sheep genomes, including the Iranian mouflon (O. gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne and mutation load estimates are lower in Cyprian compared with Anatolian mouflons, suggesting the purging of recessive deleterious variants in Cyprian sheep under a small long-term Ne, possibly attributable to founder effects, island isolation, introgression from domestic lineages, or differences in their bottleneck dynamics. Expanding our analyses to worldwide wild and feral Ovis genomes, we observe varying viability metrics among different lineages and a limited consistency between viability metrics and International Union for Conservation of Nature conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Comparison and optimization of protocols and whole-genome capture conditions for ancient DNA samples.

Ancient DNA (aDNA) obtained from human remains is typically fragmented and present in relatively low amounts. Here we investigate a set of optimal methods for producing aDNA data by comparing silica-based DNA extraction and aDNA library preparation protocols. We also test the efficiency of whole-genome enrichment (WGC) on ancient human samples by modifying a number of parameter combinations. We find that the Dabney extraction protocol performs significantly better than alternatives. We further observed a positive trend with the BEST library protocol indicating lower clonality. Notably, our results suggest that WGC is effective at retrieving endogenous DNA, particularly from poorly-preserved human samples, by increasing human endogenous proportions by 5x. Thus, aDNA studies will be most likely to benefit from our results.

Identification of microbial pathogens in Neolithic Scandinavian humans.

With the Neolithic transition, human lifestyle shifted from hunting and gathering to farming. This change altered subsistence patterns, cultural expression, and population structures as shown by the archaeological/zooarchaeological record, as well as by stable isotope and ancient DNA data. Here, we used metagenomic data to analyse if the transitions also impacted the microbiome composition in 25 Mesolithic and Neolithic hunter-gatherers and 13 Neolithic farmers from several Scandinavian Stone Age cultural contexts. Salmonella enterica, a bacterium that may have been the cause of death for the infected individuals, was found in two Neolithic samples from Battle Axe culture contexts. Several species of the bacterial genus Yersinia were found in Neolithic individuals from Funnel Beaker culture contexts as well as from later Neolithic context. Transmission of e.g. Y. enterocolitica may have been facilitated by the denser populations in agricultural contexts.

Metagenomic analysis of Mesolithic chewed pitch reveals poor oral health among stone age individuals.

Prehistoric chewed pitch has proven to be a useful source of ancient DNA, both from humans and their microbiomes. Here we present the metagenomic analysis of three pieces of chewed pitch from Huseby Klev, Sweden, that were dated to 9,890-9,540 before present. The metagenomic profile exposes a Mesolithic oral microbiome that includes opportunistic oral pathogens. We compared the data with healthy and dysbiotic microbiome datasets and we identified increased abundance of periodontitis-associated microbes. In addition, trained machine learning models predicted dysbiosis with 70-80% probability. Moreover, we identified DNA sequences from eukaryotic species such as red fox, hazelnut, red deer and apple. Our results indicate a case of poor oral health during the Scandinavian Mesolithic, and show that pitch pieces have the potential to provide information on material use, diet and oral health.

Population genomics of the muskox' resilience in the near absence of genetic variation.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

aMeta: an accurate and memory-efficient ancient metagenomic profiling workflow.

Analysis of microbial data from archaeological samples is a growing field with great potential for understanding ancient environments, lifestyles, and diseases. However, high error rates have been a challenge in ancient metagenomics, and the availability of computational frameworks that meet the demands of the field is limited. Here, we propose aMeta, an accurate metagenomic profiling workflow for ancient DNA designed to minimize the amount of false discoveries and computer memory requirements. Using simulated data, we benchmark aMeta against a current state-of-the-art workflow and demonstrate its superiority in microbial detection and authentication, as well as substantially lower usage of computer memory.

Patrilocality and hunter-gatherer-related ancestry of populations in East-Central Europe during the Middle Bronze Age.

The demographic history of East-Central Europe after the Neolithic period remains poorly explored, despite this region being on the confluence of various ecological zones and cultural entities. Here, the descendants of societies associated with steppe pastoralists form Early Bronze Age were followed by Middle Bronze Age populations displaying unique characteristics. Particularly, the predominance of collective burials, the scale of which, was previously seen only in the Neolithic. The extent to which this re-emergence of older traditions is a result of genetic shift or social changes in the MBA is a subject of debate. Here by analysing 91 newly generated genomes from Bronze Age individuals from present Poland and Ukraine, we discovered that Middle Bronze Age populations were formed by an additional admixture event involving a population with relatively high proportions of genetic component associated with European hunter-gatherers and that their social structure was based on, primarily patrilocal, multigenerational kin-groups.

Genetic continuity, isolation, and gene flow in Stone Age Central and Eastern Europe.

The genomic landscape of Stone Age Europe was shaped by multiple migratory waves and population replacements, but different regions do not all show similar patterns. To refine our understanding of the population dynamics before and after the dawn of the Neolithic, we generated and analyzed genomic sequence data from human remains of 56 individuals from the Mesolithic, Neolithic, and Eneolithic across Central and Eastern Europe. We found that Mesolithic European populations formed a geographically widespread isolation-by-distance zone ranging from Central Europe to Siberia, which was already established 10,000 years ago. We found contrasting patterns of population continuity during the Neolithic transition: people around the lower Dnipro Valley region, Ukraine, showed continuity over 4000 years, from the Mesolithic to the end of the Neolithic, in contrast to almost all other parts of Europe where population turnover drove this cultural change, including vast areas of Central Europe and around the Danube River.

Iron age genomic data from Althiburos - Tunisia renew the debate on the origins of African taurine cattle.

The Maghreb is a key region for understanding the dynamics of cattle dispersal and admixture with local aurochs following their earliest domestication in the Fertile Crescent more than 10,000 years ago. Here, we present data on autosomal genomes and mitogenomes obtained for four archaeological specimens of Iron Age (∼2,800 cal BP-2,000 cal BP) domestic cattle from the Eastern Maghreb, i.e. Althiburos (El Kef, Tunisia). D-loop sequences were obtained for an additional eight cattle specimens from this site. Maternal lineages were assigned to the elusive R and ubiquitous African-T1 haplogroups found in two and ten Althiburos specimens, respectively. Our results can be explained by post-domestication hybridization of Althiburos cattle with local aurochs. However, we cannot rule out an independent domestication in North Africa considering the shared ancestry of Althiburos cattle with the pre-domestic Moroccan aurochs and present-day African taurine cattle.

Imputed genomes and haplotype-based analyses of the Picts of early medieval Scotland reveal fine-scale relatedness between Iron Age, early medieval and the modern people of the UK.

There are longstanding questions about the origins and ancestry of the Picts of early medieval Scotland (ca. 300-900 CE), prompted in part by exotic medieval origin myths, their enigmatic symbols and inscriptions, and the meagre textual evidence. The Picts, first mentioned in the late 3rd century CE resisted the Romans and went on to form a powerful kingdom that ruled over a large territory in northern Britain. In the 9th and 10th centuries Gaelic language, culture and identity became dominant, transforming the Pictish realm into Alba, the precursor to the medieval kingdom of Scotland. To date, no comprehensive analysis of Pictish genomes has been published, and questions about their biological relationships to other cultural groups living in Britain remain unanswered. Here we present two high-quality Pictish genomes (2.4 and 16.5X coverage) from central and northern Scotland dated from the 5th-7th century which we impute and co-analyse with >8,300 previously published ancient and modern genomes. Using allele frequency and haplotype-based approaches, we can firmly place the genomes within the Iron Age gene pool in Britain and demonstrate regional biological affinity. We also demonstrate the presence of population structure within Pictish groups, with Orcadian Picts being genetically distinct from their mainland contemporaries. When investigating Identity-By-Descent (IBD) with present-day genomes, we observe broad affinities between the mainland Pictish genomes and the present-day people living in western Scotland, Wales, Northern Ireland and Northumbria, but less with the rest of England, the Orkney islands and eastern Scotland-where the political centres of Pictland were located. The pre-Viking Age Orcadian Picts evidence a high degree of IBD sharing across modern Scotland, Wales, Northern Ireland, and the Orkney islands, demonstrating substantial genetic continuity in Orkney for the last ~2,000 years. Analysis of mitochondrial DNA diversity at the Pictish cemetery of Lundin Links (n = 7) reveals absence of direct common female ancestors, with implications for broader social organisation. Overall, our study provides novel insights into the genetic affinities and population structure of the Picts and direct relationships between ancient and present-day groups of the UK.

The genetic history of Scandinavia from the Roman Iron Age to the present.

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

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

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

A 14,000-year-old genome sheds light on the evolution and extinction of a Pleistocene vulture.

The New World Vulture [Coragyps] occidentalis (L. Miller, 1909) is one of many species that were extinct by the end of the Pleistocene. To understand its evolutionary history we sequenced the genome of a 14,000 year old [Coragyps] occidentalis found associated with megaherbivores in the Peruvian Andes. occidentalis has been viewed as the ancestor, or possibly sister, to the extant Black Vulture Coragyps atratus, but genomic data shows occidentalis to be deeply nested within the South American clade of atratus. Coragyps atratus inhabits lowlands, but the fossil record indicates that occidentalis mostly occupied high elevations. Our results suggest that occidentalis evolved from a population of atratus in southwestern South America that colonized the High Andes 300 to 400 kya. The morphological and morphometric differences between occidentalis and atratus may thus be explained by ecological diversification following from the natural selection imposed by this new and extreme, high elevation environment. The sudden evolution of a population with significantly larger body size and different anatomical proportions than atratus thus constitutes an example of punctuated evolution.

Bioarchaeological evidence of one of the earliest Islamic burials in the Levant.

The Middle East plays a central role in human history harbouring a vast diversity of ethnic, cultural and religious groups. However, much remains to be understood about past and present genomic diversity in this region. Here we present a multidisciplinary bioarchaeological analysis of two individuals dated to the late 7th and early 8th centuries, the Umayyad Era, from Tell Qarassa, an open-air site in modern-day Syria. Radiocarbon dates and burial type are consistent with one of the earliest Islamic Arab burials in the Levant. Interestingly, we found genomic similarity to a genotyped group of modern-day Bedouins and Saudi rather than to most neighbouring Levantine groups. This study represents the genomic analysis of a secondary use site with characteristics consistent with an early Islamic burial in the Levant. We discuss our findings and possible historic scenarios in the light of forces such as genetic drift and their possible interaction with religious and cultural processes (including diet and subsistence practices).

The genetic origin of Huns, Avars, and conquering Hungarians.

Huns, Avars, and conquering Hungarians were migration-period nomadic tribal confederations that arrived in three successive waves in the Carpathian Basin between the 5th and 9th centuries. Based on the historical data, each of these groups are thought to have arrived from Asia, although their exact origin and relation to other ancient and modern populations have been debated. Recently, hundreds of ancient genomes were analyzed from Central Asia, Mongolia, and China, from which we aimed to identify putative source populations for the above-mentioned groups. In this study, we have sequenced 9 Hun, 143 Avar, and 113 Hungarian conquest period samples and identified three core populations, representing immigrants from each period with no recent European ancestry. Our results reveal that this "immigrant core" of both Huns and Avars likely originated in present day Mongolia, and their origin can be traced back to Xiongnus (Asian Huns), as suggested by several historians. On the other hand, the "immigrant core" of the conquering Hungarians derived from an earlier admixture of Mansis, early Sarmatians, and descendants of late Xiongnus. We have also shown that a common "proto-Ugric" gene pool appeared in the Bronze Age from the admixture of Mezhovskaya and Nganasan people, supporting genetic and linguistic data. In addition, we detected shared Hun-related ancestry in numerous Avar and Hungarian conquest period genetic outliers, indicating a genetic link between these successive nomadic groups. Aside from the immigrant core groups, we identified that the majority of the individuals from each period were local residents harboring "native European" ancestry.

Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication.

Sheep were among the first domesticated animals, but their demographic history is little understood. Here we analyzed nuclear polymorphism and mitochondrial data (mtDNA) from ancient central and west Anatolian sheep dating from Epipaleolithic to late Neolithic, comparatively with modern-day breeds and central Asian Neolithic/Bronze Age sheep (OBI). Analyzing ancient nuclear data, we found that Anatolian Neolithic sheep (ANS) are genetically closest to present-day European breeds relative to Asian breeds, a conclusion supported by mtDNA haplogroup frequencies. In contrast, OBI showed higher genetic affinity to present-day Asian breeds. These results suggest that the east-west genetic structure observed in present-day breeds had already emerged by 6000 BCE, hinting at multiple sheep domestication episodes or early wild introgression in southwest Asia. Furthermore, we found that ANS are genetically distinct from all modern breeds. Our results suggest that European and Anatolian domestic sheep gene pools have been strongly remolded since the Neolithic.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics.

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.