Bronze age Northern Eurasian genetics in the context of development of metallurgy and Siberian ancestry.

The Eurasian Bronze Age (BA) has been described as a period of substantial human migrations, the emergence of pastoralism, horse domestication, and development of metallurgy. This study focuses on two north Eurasian sites sharing Siberian genetic ancestry. One of the sites, Rostovka, is associated with the Seima-Turbino (ST) phenomenon (~2200-1900 BCE) that is characterized by elaborate metallurgical objects found throughout Northern Eurasia. The genetic profiles of Rostovka individuals vary widely along the forest-tundra Siberian genetic cline represented by many modern Uralic-speaking populations, and the genetic heterogeneity observed is consistent with the current understanding of the ST being a transcultural phenomenon. Individuals from the second site, Bolshoy Oleni Ostrov in Kola, in comparison form a tighter cluster on the Siberian ancestry cline. We further explore this Siberian ancestry profile and assess the role of the ST phenomenon and other contemporaneous BA cultures in the spread of Uralic languages and Siberian ancestry.

Sacrificed Maya boys tied to myth of 'Hero Twins'.

Ancient DNA shows continuity between living and ancient Maya communities.

52,000 years of woolly rhinoceros population dynamics reveal extinction mechanisms.

The extinction of the woolly rhinoceros (Coelodonta antiquitatis) at the onset of the Holocene remains an enigma, with conflicting evidence regarding its cause and spatiotemporal dynamics. This partly reflects challenges in determining demographic responses of late Quaternary megafauna to climatic and anthropogenic causal drivers with available genetic and paleontological techniques. Here, we show that elucidating mechanisms of ancient extinctions can benefit from a detailed understanding of fine-scale metapopulation dynamics, operating over many millennia. Using an abundant fossil record, ancient DNA, and high-resolution simulation models, we untangle the ecological mechanisms and causal drivers that are likely to have been integral in the decline and later extinction of the woolly rhinoceros. Our 52,000-y reconstruction of distribution-wide metapopulation dynamics supports a pathway to extinction that began long before the Holocene, when the combination of cooling temperatures and low but sustained hunting by humans trapped woolly rhinoceroses in suboptimal habitats along the southern edge of their range. Modeling indicates that this ecological trap intensified after the end of the last ice age, preventing colonization of newly formed suitable habitats, weakening stabilizing metapopulation processes, triggering the extinction of the woolly rhinoceros in the early Holocene. Our findings suggest that fragmentation and resultant metapopulation dynamics should be explicitly considered in explanations of late Quaternary megafauna extinctions, sending a clarion call to the fragility of the remaining large-bodied grazers restricted to disjunct fragments of poor-quality habitat due to anthropogenic environmental change.

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.

Bioarchaeology aids the cultural understanding of six characters in search of their agency (Tarquinia, ninth-seventh century BC, central Italy).

Etruria contained one of the great early urban civilisations in the Italian peninsula during the first millennium BC, much studied from a cultural, humanities-based, perspective, but relatively little with scientific data, and rarely in combination. We have addressed the unusual location of twenty inhumations found in the sacred heart of the Etruscan city of Tarquinia, focusing on six of these as illustrative, contrasting with the typical contemporary cremations found in cemeteries on the edge of the city. The cultural evidence suggests that the six skeletons were also distinctive in their ritualization and memorialisation. Focusing on the six, as a representative sample, the scientific evidence of osteoarchaeology, isotopic compositions, and ancient DNA has established that these appear to show mobility, diversity and violence through an integrated bioarchaeological approach. The combination of multiple lines of evidence makes major strides towards a deeper understanding of the role of these extraordinary individuals in the life of the early city of Etruria.

Identification of the 18 World War II executed citizens of Adele, Rethymnon, Crete using an ancient DNA approach and low coverage genomes.

In the Battle of Crete during the World War II occupation of Greece, the German forces faced substantial civilian resistance. To retribute the numerous German losses, a series of mass executions took place in numerous places in Crete; a common practice reported from Greece and elsewhere. In Adele, a village in the regional unit of Rethymnon, 18 male civilians were executed and buried in a burial pit at the Sarakina site. In this study, the first one conducted for a conflict that occurred in Greece, we identified for humanitarian purposes the 18 skulls of the Sarakina victims, following a request from the local community of Adele. The molecular identification of historical human remains via ancient DNA approaches and low coverage whole genome sequencing has only recently been introduced. Here, we performed genome skimming on the living relatives of the victims, as well as high throughput historical DNA analysis on the skulls to infer the kinship degrees among the victims via genetic relatedness analyses. We also conducted targeted anthropological analysis to successfully complete the identification of all Sarakina victims. We demonstrate that our methodological approach constitutes a potentially highly informative forensic tool to identify war victims. It can hence be applied to analogous studies on degraded DNA, thus, paving the path for systematic war victim identification in Greece and beyond.

Ancient genomes revealed the complex human interactions of the ancient western Tibetans.

The western Tibetan Plateau is the crossroad between the Tibetan Plateau, Central Asia, and South Asia, and it is a potential human migration pathway connecting these regions. However, the population history of the western Tibetan Plateau remains largely unexplored due to the lack of ancient genomes covering a long-time interval from this area. Here, we reported genome-wide data of 65 individuals dated to 3,500-300 years before present (BP) in the Ngari prefecture. The ancient western Tibetan Plateau populations share the majority of their genetic components with the southern Tibetan Plateau populations and have maintained genetic continuity since 3,500 BP while maintaining interactions with populations within and outside the Tibetan Plateau. Within the Tibetan Plateau, the ancient western Tibetan Plateau populations were influenced by the additional expansion from the south to the southwest plateau before 1,800 BP. Outside the Tibetan Plateau, the western Tibetan Plateau populations interacted with both South and Central Asian populations at least 2,000 years ago, and the South Asian-related genetic influence, despite being very limited, was from the Indus Valley Civilization (IVC) migrants in Central Asia instead of the IVC populations from the Indus Valley. In light of the new genetic data, our study revealed the complex population interconnections across and within the Tibetan Plateau.

Ancient Sheep Genomes Reveal Four Millennia of North European Short-Tailed Sheep in the Baltic Sea Region.

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

Ancient Mycobacterium leprae genome reveals medieval English red squirrels as animal leprosy host.

Leprosy, one of the oldest recorded diseases in human history, remains prevalent in Asia, Africa, and South America, with over 200,000 cases every year.1,2 Although ancient DNA (aDNA) approaches on the major causative agent, Mycobacterium leprae, have elucidated the disease's evolutionary history,3,4,5 the role of animal hosts and interspecies transmission in the past remains unexplored. Research has uncovered relationships between medieval strains isolated from archaeological human remains and modern animal hosts such as the red squirrel in England.6,7 However, the time frame, distribution, and direction of transmissions remains unknown. Here, we studied 25 human and 12 squirrel samples from two archaeological sites in Winchester, a medieval English city well known for its leprosarium and connections to the fur trade. We reconstructed four medieval M. leprae genomes, including one from a red squirrel, at a 2.2-fold average coverage. Our analysis revealed a phylogenetic placement of all strains on branch 3 as well as a close relationship between the squirrel strain and one newly reconstructed medieval human strain. In particular, the medieval squirrel strain is more closely related to some medieval human strains from Winchester than to modern red squirrel strains from England, indicating a yet-undetected circulation of M. leprae in non-human hosts in the Middle Ages. Our study represents the first One Health approach for M. leprae in archaeology, which is centered around a medieval animal host strain, and highlights the future capability of such approaches to understand the disease's zoonotic past and current potential.

Phylogenetic insights into the genetic legacies of Hungarian-speaking communities in the Carpathian Basin.

This study focuses on exploring the uniparental genetic lineages of Hungarian-speaking minorities residing in rural villages of Baranja (Croatia) and the Zobor region (Slovakia). We aimed to identify ancestral lineages by examining genetic markers distributed across the entire mitogenome and on the Y-chromosome. This allowed us to discern disparities in regional genetic structures within these communities. By integrating our newly acquired genetic data from a total of 168 participants with pre-existing Eurasian and ancient DNA datasets, our goal was to enrich the understanding of the genetic history trajectories of Carpathian Basin populations. Our findings suggest that while population-based analyses may not be sufficiently robust to detect fine-scale uniparental genetic patterns with the sample sizes at hand, phylogenetic analysis of well-characterized Y-chromosomal Short Tandem Repeat (STR) data and entire mitogenome sequences did uncover multiple lineage ties to far-flung regions and eras. While the predominant portions of both paternal and maternal DNA align with the East-Central European spectrum, rarer subhaplogroups and lineages have unveiled ancient ties to both prehistoric and historic populations spanning Europe and Eastern Eurasia. This research augments the expansive field of phylogenetics, offering critical perspectives on the genetic constitution and heritage of the communities in East-Central Europe.

New Canary Islands Roman mediated settlement hypothesis deduced from coalescence ages of curated maternal indigenous lineages.

Numerous genetic studies have contributed to reconstructing the human history of the Canary Islands population. The recent use of new ancient DNA targeted enrichment and next-generation sequencing techniques on new Canary Islands samples have greatly improved these molecular results. However, the bulk of the available data is still provided by the classic mitochondrial DNA phylogenetic and phylogeographic studies carried out on the indigenous, historical, and extant human populations of the Canary Islands. In the present study, making use of all the accumulated mitochondrial information, the existence of DNA contamination and archaeological sample misidentification in those samples is evidenced. Following a thorough review of these cases, the new phylogeographic analysis revealed the existence of a heterogeneous indigenous Canarian population, asymmetrically distributed across the various islands, which most likely descended from a unique mainland settlement. These new results and new proposed coalescent ages are compatible with a Roman-mediated arrival driven by the exploitation of the purple dye manufacture in the Canary Islands.

A 1 Ma sedimentary ancient DNA (sedaDNA) record of catchment vegetation changes and the developmental history of tropical Lake Towuti (Sulawesi, Indonesia).

Studying past ecosystems from ancient environmental DNA preserved in lake sediments (sedaDNA) is a rapidly expanding field. This research has mainly involved Holocene sediments from lakes in cool climates, with little known about the suitability of sedaDNA to reconstruct substantially older ecosystems in the warm tropics. Here, we report the successful recovery of chloroplast trnL (UAA) sequences (trnL-P6 loop) from the sedimentary record of Lake Towuti (Sulawesi, Indonesia) to elucidate changes in regional tropical vegetation assemblages during the lake's Late Quaternary paleodepositional history. After the stringent removal of contaminants and sequence artifacts, taxonomic assignment of the remaining genuine trnL-P6 reads showed that native nitrogen-fixing legumes, C3 grasses, and shallow wetland vegetation (Alocasia) were most strongly associated with >1-million-year-old (>1 Ma) peats and silts (114-98.8 m composite depth; mcd), which were deposited in a landscape of active river channels, shallow lakes, and peat-swamps. A statistically significant shift toward partly submerged shoreline vegetation that was likely rooted in anoxic muddy soils (i.e., peatland forest trees and wetland C3 grasses (Oryzaceae) and nutrient-demanding aquatic herbs (presumably Oenanthe javanica)) occurred at 76 mcd (~0.8 Ma), ~0.2 Ma after the transition into a permanent lake. This wetland vegetation was most strongly associated with diatom ooze (46-37 mcd), thought to be deposited during maximum nutrient availability and primary productivity. Herbs (Brassicaceae), trees/shrubs (Fabaceae and Theaceae), and C3 grasses correlated with inorganic parameters, indicating increased drainage of ultramafic sediments and laterite soils from the lakes' catchment, particularly at times of inferred drying. Downcore variability in trnL-P6 from tropical forest trees (Toona), shady ground cover herbs (Zingiberaceae), and tree orchids (Luisia) most strongly correlated with sediments of a predominantly felsic signature considered to be originating from the catchment of the Loeha River draining into Lake Towuti during wetter climate conditions. However, the co-correlation with dry climate-adapted trees (i.e., Castanopsis or Lithocarpus) plus C4 grasses suggests that increased precipitation seasonality also contributed to the increased drainage of felsic Loeha River sediments. This multiproxy approach shows that despite elevated in situ temperatures, tropical lake sediments potentially comprise long-term archives of ancient environmental DNA for reconstructing ecosystems, which warrants further exploration.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations.

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

Palaeogenomic insights into the origins of early settlers on the island of Cyprus.

Archaeological evidence supports sporadic seafaring visits to the Eastern Mediterranean island of Cyprus by Epipaleolithic hunter-gatherers over 12,000 years ago, followed by permanent settlements during the early Neolithic. The geographical origins of these early seafarers have so far remained elusive. By systematically analysing all available genomes from the late Pleistocene to early Holocene Near East (c. 14,000-7000 cal BCE), we provide a comprehensive overview of the genetic landscape of the early Neolithic Fertile Crescent and Anatolia and infer the likely origins of three recently published genomes from Kissonerga-Mylouthkia (Cypriot Late Pre-Pottery Neolithic B, c. 7600-6800 cal BCE). These appear to derive roughly 80% of their ancestry from Aceramic Neolithic Central Anatolians residing in or near the Konya plain, and the remainder from a genetically basal Levantine population. Based on genome-wide weighted ancestry covariance analysis, we infer that this admixture event took place roughly between 14,000 and 10,000 BCE, coinciding with the transition from the Cypriot late Epipaleolithic to the Pre-Pottery Neolithic A (PPNA). Additionally, we identify strong genetic affinities between the examined Cypro-LPPNB individuals and later northwestern Anatolians and the earliest European Neolithic farmers. Our results inform archaeological evidence on prehistoric demographic processes in the Eastern Mediterranean, providing important insights into early seafaring, maritime connections, and insular settlement.

Assessing autosomal aneuploidy in ancient genomes.

Using genetic methods, aneuploidies can be detected in ancient human remains, which is so far the only way to reliably prove their existence in the past. As highlighted in recent studies by Rohrlach et al. and by Anastasiadou et al., this initial step enables a deeper exploration of the history of rare diseases, encompassing the social and historical contexts of the afflicted individuals.

Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change.

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

Ancient DNA indicates a century of overhunting did not reduce genetic diversity in Pacific Walruses (Odobenus rosmarus divergens).

Pacific Walruses (Odobenus rosmarus divergens [Illiger 1815]) are gregarious marine mammals considered to be sentinels of the Arctic because of their dependence on sea ice for feeding, molting, and parturition. Like many other marine mammal species, their population sizes were decimated by historical overhunting in the nineteenth and twentieth centuries. Although they have since been protected from nearly all commercial hunting pressure, they now face rapidly accelerating habitat loss as global warming reduces the extent of summer sea ice in the Arctic. To investigate how genetic variation was impacted by overhunting, we obtained mitochondrial DNA sequences from historic Pacific Walrus samples in Alaska that predate the period of overhunting, as well as from extant populations. We found that genetic variation was unchanged over this period, suggesting Pacific Walruses are resilient to genetic attrition in response to reduced population size, and that this may be related to their high vagility and lack of population structure. Although Pacific Walruses will almost certainly continue to decline in number as the planet warms and summer sea ice is further reduced, they may be less susceptible to the ratcheting effects of inbreeding that typically accompany shrinking populations.

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

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

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.