hapCon: Estimating contamination of ancient genomes by copying from reference haplotypes

Human ancient DNA (aDNA) studies have surged in recent years, revolutionizing the study of the human past. Typically, aDNA is preserved poorly, making such data prone to contamination from other human DNA. Therefore, it is important to rule out substantial contamination before proceeding to downstream analysis. As most aDNA samples can only be sequenced to low coverages (<1× average depth), computational methods that can robustly estimate contamination in the low coverage regime are needed. However, the ultra low-coverage regime (0.1× and below) remains a challenging task for existing approaches.We present a new method to estimate contamination in aDNA for male modern humans. It utilizes a Li&Stephens haplotype copying model for haploid X chromosomes, with mismatches modeled as errors or contamination. We assessed this new approach, hapCon, on simulated and down-sampled empirical aDNA data. Our experiments demonstrate that hapCon outperforms a commonly used tool for estimating male X contamination (ANGSD), with substantially lower variance and narrower confidence intervals, especially in the low coverage regime. We found that hapCon provides useful contamination estimates for coverages as low as 0.1× for SNP capture data (1240k) and 0.02× for whole genome sequencing data, substantially extending the coverage limit of previous male X chromosome-based contamination estimation methods. Our experiments demonstrate that hapCon has little bias for contamination up to 25–30\% as long as the contaminating source is specified within continental genetic variation, and that its application range extends to human aDNA as old as ∼45 000 and various global ancestries.We make hapCon available as part of a python package (hapROH), which is available at the Python Package Index (https://pypi.org/project/hapROH) and can be installed via pip. The documentation provides example use cases as blueprints for custom applications (https://haproh.readthedocs.io/en/latest/hapCon.html). The program can analyze either BAM files or pileup files produced with samtools. An implementation of our software (hapCon) using Python and C is deposited at https://github.com/hyl317/hapROH.Supplementary data are available at Bioinformatics online

Parental relatedness through time revealed by runs of homozygosity in ancient DNA

Parental relatedness of present-day humans varies substantially across the globe, but little is known about the past. Here we analyze ancient DNA, leveraging that parental relatedness leaves genomic traces in the form of runs of homozygosity. We present an approach to identify such runs in low-coverage ancient DNA data aided by haplotype information from a modern phased reference panel. Simulation and experiments show that this method robustly detects runs of homozygosity longer than 4 centimorgan for ancient individuals with at least 0.3 × coverage. Analyzing genomic data from 1,785 ancient humans who lived in the last 45,000 years, we detect low rates of first cousin or closer unions across most ancient populations. Moreover, we find a marked decay in background parental relatedness co-occurring with or shortly after the advent of sedentary agriculture. We observe this signal, likely linked to increasing local population sizes, across several geographic transects worldwide

An extended admixture pulse model reveals the limitations to human-Neandertal introgression dating

Neandertal DNA makes up 2–3% of the genomes of all non-African individuals. The patterns of Neandertal ancestry in modern humans have been used to estimate that this is the result of gene flow that occurred during the expansion of modern humans into Eurasia, but the precise dates of this event remain largely unknown. Here, we introduce an extended admixture pulse model that allows joint estimation of the timing and duration of gene flow. This model leads to simple expressions for both the admixture segment distribution and the decay curve of ancestry linkage disequilibrium, and we show that these two statistics are closely related. In simulations, we find that estimates of the mean time of admixture are largely robust to details in gene flow models, but that the duration of the gene flow can only be recovered if gene flow is very recent and the exact recombination map is known. These results imply that gene flow from Neandertals into modern humans could have happened over hundreds of generations. Ancient genomes from the time around the admixture event are thus likely required to resolve the question when, where, and for how long humans and Neandertals interacted

Reconstructing genetic histories and social organisation in Neolithic and Bronze Age Croatia

Ancient DNA studies have revealed how human migrations from the Neolithic to the Bronze Age transformed the social and genetic structure of European societies. Present-day Croatia lies at the heart of ancient migration routes through Europe, yet our knowledge about social and genetic processes here remains sparse. To shed light on these questions, we report new whole-genome data for 28 individuals dated to between ~ 4700 BCE–400 CE from two sites in present-day eastern Croatia. In the Middle Neolithic we evidence first cousin mating practices and strong genetic continuity from the Early Neolithic. In the Middle Bronze Age community that we studied, we find multiple closely related males suggesting a patrilocal social organisation. We also find in that community an unexpected genetic ancestry profile distinct from individuals found at contemporaneous sites in the region, due to the addition of hunter-gatherer-related ancestry. These findings support archaeological evidence for contacts with communities further north in the Carpathian Basin. Finally, an individual dated to Roman times exhibits an ancestry profile that is broadly present in the region today, adding an important data point to the substantial shift in ancestry that occurred in the region between the Bronze Age and today

Population history from the Neolithic to present on the Mediterranean island of Sardinia: an ancient DNA perspective

Recent ancient DNA studies of western Eurasia have revealed a dynamic history of admixture, with evidence for major migrations during the Neolithic and Bronze Age. The population of the Mediterranean island of Sardinia has been notable in these studies –} Neolithic individuals from mainland Europe cluster more closely with Sardinian individuals than with all other present-day Europeans. The current model to explain this result is that Sardinia received an initial influx of Neolithic ancestry and then remained relatively isolated from expansions in the later Neolithic and Bronze Age that took place in continental Europe. To test this model, we generated genome-wide capture data (approximately 1.2 million variants) for 43 ancient Sardinian individuals spanning the Neolithic through the Bronze Age, including individuals from Sardinia{’}s Nuragic culture, which is known for the construction of numerous large stone towers throughout the island. We analyze these new samples in the context of previously generated genome-wide ancient DNA data from 972 ancient individuals across western Eurasia and whole-genome sequence data from approximately 1,500 modern individuals from Sardinia. The ancient Sardinian individuals show a strong affinity to western Mediterranean Neolithic populations and we infer a high degree of genetic continuity on the island from the Neolithic (around fifth millennium BCE) through the Nuragic period (second millennium BCE). In particular, during the Bronze Age in Sardinia, we do not find significant levels of the {“}Steppe{” ancestry that was spreading in many other parts of Europe at that time. We also characterize subsequent genetic influx between the Nuragic period and the present. We detect novel, modest signals of admixture between 1,000 BCE and present-day, from ancestry sources in the eastern and northern Mediterranean. Within Sardinia, we confirm that populations from the more geographically isolated mountainous provinces have experienced elevated levels of genetic drift and that northern and southwestern regions of the island received more gene flow from outside Sardinia. Overall, our genetic analysis sheds new light on the origin of Neolithic settlement on Sardinia, reinforces models of genetic continuity on the island, and provides enhanced power to detect post-Bronze-Age gene flow. Together, these findings offer a refined demographic model for future medical genetic studies in Sardinia

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century

We report genome-wide data for 33 Ashkenazi Jews (AJ), dated to the 14th century, following a salvageexcavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are geneticallysimilar to modern AJ and have substantial Southern European ancestry, but they show more variabilityin Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried the samenearly-AJ-specific mitochondrial haplogroup and eight carried pathogenic variants known to affect AJtoday. These observations, together with high levels of runs of homozygosity, suggest that the Erfurtcommunity had already experienced the major reduction in size that affected modern AJ. However, theErfurt bottleneck was more severe, implying substructure in medieval AJ. Together, our results suggestthat the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14th centuryand highlight late medieval genetic heterogeneity no longer present in modern AJ

Converting multilevel nonclassicality into genuine multipartite entanglement

Characterizing genuine quantum resources and determining operational rules for their manipulation are crucial steps to appraise possibilities and limitations of quantum technologies. Two such key resources are nonclassicality, manifested as quantum superposition between reference states of a single system, and entanglement, capturing quantum correlations among two or more subsystems. Here we present a general formalism for the conversion of nonclassicality into multipartite entanglement, showing that a faithful reversible transformation between the two resources is always possible within a precise resource-theoretic framework. Specializing to quantum coherence between the levels of a quantum system as an instance of nonclassicality, we introduce explicit protocols for such a mapping. We further show that the conversion relates multilevel coherence and multipartite entanglement not only qualitatively, but also quantitatively, restricting the amount of entanglement achievable in the process and in particular yielding an equality between the two resources when quantified by fidelity-based geometric measures

Extensive pedigrees reveal the social organization of a Neolithic community

Social anthropology and ethnographic studies have described kinship systems and networks of contact and exchange in extant populations 1–4. However, for prehistoric societies, these systems can be studied only indirectly from biological and cultural remains. Stable isotope data, sex and age at death can provide insights into the demographic structure of a burial community and identify local versus non-local childhood signatures, archaeogenetic data can reconstruct the biological relationships between individuals, which enables the reconstruction of pedigrees, and combined evidence informs on kinship practices and residence patterns in prehistoric societies. Here we report ancient DNA, strontium isotope and contextual data from more than 100 individuals from the site Gurgy ‘les Noisats’ (France), dated to the western European Neolithic around 4850–4500 bc. We find that this burial community was genetically connected by two main pedigrees, spanning seven generations, that were patrilocal and patrilineal, with evidence for female exogamy and exchange with genetically close neighbouring groups. The microdemographic structure of individuals linked and unlinked to the pedigrees reveals additional information about the social structure, living conditions and site occupation. The absence of half-siblings and the high number of adult full siblings suggest that there were stable health conditions and a supportive social network, facilitating high fertility and low mortality 5. Age-structure differences and strontium isotope results by generation indicate that the site was used for just a few decades, providing new insights into shifting sedentary farming practices during the European Neolithic

Genetic history from the Middle Neolithic to present on the Mediterranean island of Sardinia

The island of Sardinia has been of particular interest to geneticists for decades. The current model for Sardinia's genetic history describes the island as harboring a founder population that was established largely from the Neolithic peoples of southern Europe and remained isolated from later Bronze Age expansions on the mainland. To evaluate this model, we generate genome-wide ancient DNA data for 70 individuals from 21 Sardinian archaeological sites spanning the Middle Neolithic through the Medieval period. The earliest individuals show a strong affinity to western Mediterranean Neolithic populations, followed by an extended period of genetic continuity on the island through the Nuragic period (second millennium BCE). Beginning with individuals from Phoenician/Punic sites (first millennium BCE), we observe spatially-varying signals of admixture with sources principally from the eastern and northern Mediterranean. Overall, our analysis sheds light on the genetic history of Sardinia, revealing how relationships to mainland populations shifted over time.Joseph H. Marcus ... Wolfgang Haak ... et al

Ancient DNA reveals admixture history and endogamy in the prehistoric Aegean (advance online)

The Neolithic and Bronze Ages were highly transformative periods forthe genetic history of Europe but for the Aegean—a region fundamentalto Europe’s prehistory—the biological dimensions of cultural transitionshave been elucidated only to a limited extent so far. We have analysed newlygenerated genome-wide data from 102 ancient individuals from Crete, theGreek mainland and the Aegean Islands, spanning from the Neolithic tothe Iron Age. We found that the early farmers from Crete shared the sameancestry as other contemporaneous Neolithic Aegeans. In contrast, the endof 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 Europeanancestry in the Greek mainland by the Middle Bronze Age, we additionallyshow that such genetic signatures appeared in Crete gradually from theseventeenth to twelfth centuries bc, a period when the influence of themainland over the island intensified. Biological and cultural connectednesswithin the Aegean is also supported by the finding of consanguineousendogamy practiced at high frequencies, unprecedented in the globalancient DNA record. Our results highlight the potential of archaeogenomicapproaches in the Aegean for unravelling the interplay of genetic admixture,marital and other cultural practice