AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Forensic Microbiology: When, Where and How

Forensic microbiology is a relatively new discipline, born in part thanks to the development of advanced methodologies for the detection, identification and characterization of microorganisms, and also in relation to the growing impact of infectious diseases of iatrogenic origin. Indeed, the increased application of medical practices, such as transplants, which require immunosuppressive treatments, and the growing demand for prosthetic installations, associated with an increasing threat of antimicrobial resistance, have led to a rise in the number of infections of iatrogenic origin, which entails important medico-legal issues. On the other hand, the possibility of detecting minimal amounts of microorganisms, even in the form of residual traces (e.g., their nucleic acids), and of obtaining gene and genomic sequences at contained costs, has made it possible to ask new questions of whether cases of death or illness might have a microbiological origin, with the possibility of also tracing the origin of the microorganisms involved and reconstructing the chain of contagion. In addition to the more obvious applications, such as those mentioned above related to the origin of iatrogenic infections, or to possible cases of infections not properly diagnosed and treated, a less obvious application of forensic microbiology concerns its use in cases of violence or violent death, where the characterization of the microorganisms can contribute to the reconstruction of the case. Finally, paleomicrobiology, e.g., the reconstruction and characterization of microorganisms in historical or even archaeological remnants, can be considered as a sister discipline of forensic microbiology. In this article, we will review these different aspects and applications of forensic microbiology.

Homo sapiens reached the higher latitudes of Europe by 45,000 years ago

The Middle to Upper Palaeolithic transition in Europe is associated with the regional disappearance of Neanderthals and the spread of Homo sapiens. Late Neanderthals persisted in western Europe several millennia after the occurrence of H. sapiens in eastern Europe1. Local hybridization between the two groups occurred2, but not on all occasions3. Archaeological evidence also indicates the presence of several technocomplexes during this transition, complicating our understanding and the association of behavioural adaptations with specific hominin groups4. One such technocomplex for which the makers are unknown is the Lincombian-Ranisian-Jerzmanowician (LRJ), which has been described in northwestern and central Europe5-8. Here we present the morphological and proteomic taxonomic identification, mitochondrial DNA analysis and direct radiocarbon dating of human remains directly associated with an LRJ assemblage at the site Ilsenhöhle in Ranis (Germany). These human remains are among the earliest directly dated Upper Palaeolithic H. sapiens remains in Eurasia. We show that early H. sapiens associated with the LRJ were present in central and northwestern Europe long before the extinction of late Neanderthals in southwestern Europe. Our results strengthen the notion of a patchwork of distinct human populations and technocomplexes present in Europe during this transitional period.

Chronological and genetic analysis of an Upper Palaeolithic female infant burial from Borsuka Cave, Poland

Six infant human teeth and 112 animal tooth pendants from Borsuka Cave were identified as the oldest burial in Poland. However, uncertainties around the dating and the association of the teeth to the pendants have precluded their association with an Upper Palaeolithic archaeological industry. Using <67 mg per tooth, we combined dating and genetic analyses of two human teeth and six herbivore tooth pendants to address these questions. Our interdisciplinary approach yielded informative results despite limited sampling material, and high levels of degradation and contamination. Our results confirm the Palaeolithic origin of the human remains and herbivore pendants, and permit us to identify the infant as female and discuss the association of the assemblage with different Palaeolithic industries. This study exemplifies the progress that has been made toward minimally destructive methods and the benefits of integrating methods to maximize data retrieval from precious but highly degraded and contaminated prehistoric material.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

decOM: similarity-based microbial source tracking of ancient oral samples using k-mer-based methods

BACKGROUND

The analysis of ancient oral metagenomes from archaeological human and animal samples is largely confounded by contaminant DNA sequences from modern and environmental sources. Existing methods for Microbial Source Tracking (MST) estimate the proportions of environmental sources, but do not perform well on ancient metagenomes. We developed a novel method called decOM for Microbial Source Tracking and classification of ancient and modern metagenomic samples using k-mer matrices.

RESULTS

We analysed a collection of 360 ancient oral, modern oral, sediment/soil and skin metagenomes, using stratified five-fold cross-validation. decOM estimates the contributions of these source environments in ancient oral metagenomic samples with high accuracy, outperforming two state-of-the-art methods for source tracking, FEAST and mSourceTracker.

CONCLUSIONS

decOM is a high-accuracy microbial source tracking method, suitable for ancient oral metagenomic data sets. The decOM method is generic and could also be adapted for MST of other ancient and modern types of metagenomes. We anticipate that decOM will be a valuable tool for MST of ancient metagenomic studies. Video Abstract.

Artefact Profiling: Panomics Approaches for Understanding the Materiality of Written Artefacts

This review explains the strategies behind genomics, proteomics, metabolomics, metallomics and isotopolomics approaches and their applicability to written artefacts. The respective sub-chapters give an insight into the analytical procedure and the conclusions drawn from such analyses. A distinction is made between information that can be obtained from the materials used in the respective manuscript and meta-information that cannot be obtained from the manuscript itself, but from residues of organisms such as bacteria or the authors and readers. In addition, various sampling techniques are discussed in particular, which pose a special challenge in manuscripts. The focus is on high-resolution, non-targeted strategies that can be used to extract the maximum amount of information about ancient objects. The combination of the various omics disciplines (panomics) especially offers potential added value in terms of the best possible interpretations of the data received. The information obtained can be used to understand the production of ancient artefacts, to gain impressions of former living conditions, to prove their authenticity, to assess whether there is a toxic hazard in handling the manuscripts, and to be able to determine appropriate measures for their conservation and restoration.

Ancient human DNA recovered from a Palaeolithic pendant

Artefacts made from stones, bones and teeth are fundamental to our understanding of human subsistence strategies, behaviour and culture in the Pleistocene. Although these resources are plentiful, it is impossible to associate artefacts to specific human individuals¹ who can be morphologically or genetically characterized, unless they are found within burials, which are rare in this time period. Thus, our ability to discern the societal roles of Pleistocene individuals based on their biological sex or genetic ancestry is limited^2-5. Here we report the development of a non-destructive method for the gradual release of DNA trapped in ancient bone and tooth artefacts. Application of the method to an Upper Palaeolithic deer tooth pendant from Denisova Cave, Russia, resulted in the recovery of ancient human and deer mitochondrial genomes, which allowed us to estimate the age of the pendant at approximately 19,000-25,000 years. Nuclear DNA analysis identifies the presumed maker or wearer of the pendant as a female individual with strong genetic affinities to a group of Ancient North Eurasian individuals who lived around the same time but were previously found only further east in Siberia. Our work redefines how cultural and genetic records can be linked in prehistoric archaeology.

Genome-wide analysis of a collective grave from Mentesh Tepe provides insight into the population structure of early neolithic population in the South Caucasus

Despite the localisation of the southern Caucasus at the outskirt of the Fertile Crescent, the Neolithisation process started there only at the beginning of the sixth millennium with the Shomutepe-Shulaveri culture of yet unclear origins. We present here genomic data for three new individuals from Mentesh Tepe in Azerbaijan, dating back to the beginnings of the Shomutepe-Shulaveri culture. We evidence that two juveniles, buried embracing each other, were brothers. We show that the Mentesh Tepe Neolithic population is the product of a recent gene flow between the Anatolian farmer-related population and the Caucasus/Iranian population, demonstrating that population admixture was at the core of the development of agriculture in the South Caucasus. By comparing Bronze Age individuals from the South Caucasus with Neolithic individuals from the same region, including Mentesh Tepe, we evidence that gene flows between Pontic Steppe populations and Mentesh Tepe-related groups contributed to the makeup of the Late Bronze Age and modern Caucasian populations. Our results show that the high cultural diversity during the Neolithic period of the South Caucasus deserves close genetic analysis.

Ancient DNA reveals admixture history and endogamy in the prehistoric Aegean

The Neolithic and Bronze Ages were highly transformative periods for the genetic history of Europe but for the Aegean-a region fundamental to Europe's prehistory-the biological dimensions of cultural transitions have been elucidated only to a limited extent so far. We have analysed newly generated genome-wide data from 102 ancient individuals from Crete, the Greek mainland and the Aegean Islands, spanning from the Neolithic to the Iron Age. We found that the early farmers from Crete shared the same ancestry as other contemporaneous Neolithic Aegeans. In contrast, the end of the Neolithic period and the following Early Bronze Age were marked by 'eastern' gene flow, which was predominantly of Anatolian origin in Crete. Confirming previous findings for additional Central/Eastern European ancestry in the Greek mainland by the Middle Bronze Age, we additionally show that such genetic signatures appeared in Crete gradually from the seventeenth to twelfth centuries BC, a period when the influence of the mainland over the island intensified. Biological and cultural connectedness within the Aegean is also supported by the finding of consanguineous endogamy practiced at high frequencies, unprecedented in the global ancient DNA record. Our results highlight the potential of archaeogenomic approaches in the Aegean for unravelling the interplay of genetic admixture, marital and other cultural practices.

4000-year-old hair from the Middle Nile highlights unusual ancient DNA degradation pattern and a potential source of early eastern Africa pastoralists

Petrous bones and teeth are the skeletal elements most often targeted by researchers for ancient DNA (aDNA) extraction, and the sources of the majority of previously published ancient African genomes. However, the high temperature environments that characterise much of Africa often lead to poor preservation of skeletal remains. Here, we successfully reconstruct and analyse genome-wide data from the naturally mummified hair of a 4000-year-old individual from Sudan in northeastern Africa, after failed attempts at DNA extraction from teeth, petrous, and cranium of this and other individuals from the Kadruka cemeteries. We find that hair DNA extracted with an established single-stranded library protocol is unusually enriched in ultra-short DNA molecules and exhibits substantial interior molecular damage. The aDNA was nonetheless amenable to genetic analyses, which revealed that the genome is genetically indistinguishable from that of early Neolithic eastern African pastoralists located 2500 kms away. Our findings are consistent with established models for the southward dispersal of Middle Nile Valley pastoral populations to the Rift Valley of eastern Africa, and provide a possible genetic source population for this dispersal. Our study highlights the value of mummified hair as an alternate source of aDNA from regions with poor bone preservation.

Genetic analysis of a bronze age individual from Ulug-depe (Turkmenistan)

The Oxus Civilisation (or Bactrio-Margian Archaeological Complex, BMAC) was the main archaeological culture of the Bronze Age in southern Central Asia. Paleogenetic analyses were previously conducted mainly on samples from the eastern part of BMAC. The population associated with BMAC descends from local Chalcolithic populations, with some outliers of steppe or South-Asian descent. Here, we present new genome-wide data for one individual from Ulug-depe (Turkmenistan), one of the main BMAC sites, located at the southwestern edge of the BMAC. We demonstrate that this individual genetically belongs to the BMAC cluster. Using this genome, we confirm that modern Indo-Iranian–speaking populations from Central Asia derive their ancestry from BMAC populations, with additional gene flow from the western and the Altai steppes in higher proportions among the Tajiks than the Yagnobi ethnic group.

A Late Pleistocene human genome from Southwest China

Southern East Asia is the dispersal center regarding the prehistoric settlement and migrations of modern humans in Asia-Pacific regions. However, the settlement pattern and population structure of paleolithic humans in this region remain elusive, and ancient DNA can provide direct information. Here, we sequenced the genome of a Late Pleistocene hominin (MZR), dated ∼14.0 thousand years ago from Red Deer Cave located in Southwest China, which was previously reported possessing mosaic features of modern and archaic hominins. MZR is the first Late Pleistocene genome from southern East Asia. Our results indicate that MZR is a modern human who represents an early diversified lineage in East Asia. The mtDNA of MZR belongs to an extinct basal lineage of the M9 haplogroup, reflecting a rich matrilineal diversity in southern East Asia during the Late Pleistocene. Combined with the published data, we detected clear genetic stratification in ancient southern populations of East/Southeast Asia and some degree of south-versus-north divergency during the Late Pleistocene, and MZR was identified as a southern East Asian who exhibits genetic continuity to present day populations. Markedly, MZR is linked deeply to the East Asian ancestry that contributed to First Americans.

hapCon: Estimating Contamination of Ancient Genomes by Copying from Reference Haplotypes

Motivation

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.

Results

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.

Availability and implementation

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 information

Supplementary data are available at Bioinformatics online.

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

Retrieving high-quality endogenous ancient DNA (aDNA) poses several challenges, including low molecular copy number, high rates of fragmentation, damage at read termini, and potential presence of exogenous contaminant DNA. All these factors complicate a reliable reconstruction of consensus aDNA sequences in reads from high-throughput sequencing platforms. Here, we report findings from a thorough evaluation of two alternative tools (ANGSD and schmutzi) aimed at overcoming these issues and constructing high-quality ancient mitogenomes. Raw genomic data (BAM/FASTQ) from a total of 17 previously published whole ancient human genomes ranging from the 14th to the 7th millennium BCE were retrieved and mitochondrial consensus sequences were reconstructed using different quality filters, with their accuracy measured and compared. Moreover, the influence of different sequence parameters (number of reads, sequenced bases, mean coverage, and rate of deamination and contamination) as predictors of derived sequence quality was evaluated. Complete mitogenomes were successfully reconstructed for all ancient samples, and for the majority of them, filtering substantially improved mtDNA consensus calling and haplogroup prediction. Overall, the schmutzi pipeline, which estimates and takes into consideration exogenous contamination, appeared to have the edge over the much faster and user-friendly alternative method (ANGSD) in moderate to high coverage samples (>1,000,000 reads). ANGSD, however, through its read termini trimming filter, showed better capabilities in calling the consensus sequence from low-quality samples. Among all the predictors of overall sample quality examined, the strongest correlation was found for the available number of sequence reads and bases. In the process, we report a previously unassigned haplogroup (U3b) for an Early Chalcolithic individual from Southern Anatolia/Northern Levant.

Environmental paleomicrobiology: using DNA preserved in aquatic sediments to its full potential

In‐depth knowledge about spatial and temporal variation in microbial diversity and function is needed for a better understanding of ecological and evolutionary responses to global change. In particular, the study of microbial ancient DNA preserved in sediment archives from lakes and oceans can help us to evaluate the responses of aquatic microbes in the past and make predictions about future biodiversity change in those ecosystems. Recent advances in molecular genetic methods applied to the analysis of historically deposited DNA in sediments have not only allowed the taxonomic identification of past aquatic microbial communities but also enabled tracing their evolution and adaptation to episodic disturbances and gradual environmental change. Nevertheless, some challenges remain for scientists to take full advantage of the rapidly developing field of paleo‐genetics, including the limited ability to detect rare taxa and reconstruct complete genomes for evolutionary studies. Here, we provide a brief review of some of the recent advances in the field of environmental paleomicrobiology and discuss remaining challenges related to the application of molecular genetic methods to study microbial diversity, ecology, and evolution in sediment archives. We anticipate that, in the near future, environmental paleomicrobiology will shed new light on the processes of microbial genome evolution and microbial ecosystem responses to quaternary environmental changes at an unprecedented level of detail. This information can, for example, aid geological reconstructions of biogeochemical cycles and predict ecosystem responses to environmental perturbations, including in the context of human‐induced global changes.

Ancient DNA Methods Improve Forensic DNA Profiling of Korean War and World War II Unknowns

The integration of massively parallel sequencing (MPS) technology into forensic casework has been of particular benefit to the identification of unknown military service members. However, highly degraded or chemically treated skeletal remains often fail to provide usable DNA profiles, even with sensitive mitochondrial (mt) DNA capture and MPS methods. In parallel, the ancient DNA field has developed workflows specifically for degraded DNA, resulting in the successful recovery of nuclear DNA and mtDNA from skeletal remains as well as sediment over 100,000 years old. In this study we use a set of disinterred skeletal remains from the Korean War and World War II to test if ancient DNA extraction and library preparation methods improve forensic DNA profiling. We identified an ancient DNA extraction protocol that resulted in the recovery of significantly more human mtDNA fragments than protocols previously used in casework. In addition, utilizing single-stranded rather than double-stranded library preparation resulted in increased attainment of reportable mtDNA profiles. This study emphasizes that the combination of ancient DNA extraction and library preparation methods evaluated here increases the success rate of DNA profiling, and likelihood of identifying historical remains.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

A curated dataset of modern and ancient high-coverage shotgun human genomes

Over the last few years, genome-wide data for a large number of ancient human samples have been collected. Whilst datasets of captured SNPs have been collated, high coverage shotgun genomes (which are relatively few but allow certain types of analyses not possible with ascertained captured SNPs) have to be reprocessed by individual groups from raw reads. This task is computationally intensive. Here, we release a dataset including 35 whole-genome sequenced samples, previously published and distributed worldwide, together with the genetic pipeline used to process them. The dataset contains 72,041,355 sites called across 19 ancient and 16 modern individuals and includes sequence data from four previously published ancient samples which we sequenced to higher coverage (10-18x). Such a resource will allow researchers to analyse their new samples with the same genetic pipeline and directly compare them to the reference dataset without re-processing published samples. Moreover, this dataset can be easily expanded to increase the sample distribution both across time and space.

Human mobility at Tell Atchana (Alalakh), Hatay, Turkey during the 2nd millennium BC: Integration of isotopic and genomic evidence

The Middle and Late Bronze Age, a period roughly spanning the 2^nd millennium BC (ca. 2000–1200 BC) in the Near East, is frequently referred to as the first ‘international age’, characterized by intense and far-reaching contacts between different entities from the eastern Mediterranean to the Near East and beyond. In a large-scale tandem study of stable isotopes and ancient DNA of individuals excavated at Tell Atchana (Alalakh, located in Hatay, Turkey), we explored the role of mobility at the capital of a regional kingdom, named Mukish during the Late Bronze Age, which spanned the Amuq Valley and some areas beyond. We generated strontium and oxygen isotope data from dental enamel for 53 individuals and 77 individuals, respectively, and added ancient DNA data of 10 newly sequenced individuals to a dataset of 27 individuals published in 2020. Additionally, we improved the DNA coverage of one individual from this 2020 dataset. The DNA data revealed a very homogeneous gene pool. This picture of an overwhelmingly local ancestry was consistent with the evidence of local upbringing in most of the individuals indicated by the isotopic data, where only five were found to be non-local. High levels of contact, trade, and exchange of ideas and goods in the Middle and Late Bronze Ages, therefore, seem not to have translated into high levels of individual mobility detectable at Tell Atchana.

A genome sequence from a modern human skull over 45,000 years old from Zlatý kůň in Czechia

Modern humans expanded into Eurasia more than 40,000 years ago following their dispersal out of Africa. These Eurasians carried ~2–3% Neanderthal ancestry in their genomes, originating from admixture with Neanderthals that took place sometime between 50,000 and 60,000 years ago, probably in the Middle East. In Europe, the modern human expansion preceded the disappearance of Neanderthals from the fossil record by 3,000–5,000 years. The genetic makeup of the first Europeans who colonized the continent more than 40,000 years ago remains poorly understood since few specimens have been studied. Here, we analyse a genome generated from the skull of a female individual from Zlatý kůň, Czechia. We found that she belonged to a population that appears to have contributed genetically neither to later Europeans nor to Asians. Her genome carries ~3% Neanderthal ancestry, similar to those of other Upper Palaeolithic hunter-gatherers. However, the lengths of the Neanderthal segments are longer than those observed in the currently oldest modern human genome of the ~45,000-year-old Ust’-Ishim individual from Siberia, suggesting that this individual from Zlatý kůň is one of the earliest Eurasian inhabitants following the expansion out of Africa.

The authors present the genome sequence of a >45,000-year-old female Homo sapiens individual from the site of Zlatý kůň, Czechia. Although radiometric dating of the human remains was inconclusive, the authors were able to use molecular methods to demonstrate that she was probably among the earliest Eurasian inhabitants following expansion out of Africa.

Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.

New Insights Into Mitochondrial DNA Reconstruction and Variant Detection in Ancient Samples

Ancient DNA (aDNA) studies are frequently focused on the analysis of the mitochondrial DNA (mtDNA), which is much more abundant than the nuclear genome, hence can be better retrieved from ancient remains. However, postmortem DNA damage and contamination make the data analysis difficult because of DNA fragmentation and nucleotide alterations. In this regard, the assessment of the heteroplasmic fraction in ancient mtDNA has always been considered an unachievable goal due to the complexity in distinguishing true endogenous variants from artifacts. We implemented and applied a computational pipeline for mtDNA analysis to a dataset of 30 ancient human samples from an Iron Age necropolis in Polizzello (Sicily, Italy). The pipeline includes several modules from well-established tools for aDNA analysis and a recently released variant caller, which was specifically conceived for mtDNA, applied for the first time to aDNA data. Through a fine-tuned filtering on variant allele sequencing features, we were able to accurately reconstruct nearly complete (>88%) mtDNA genome for almost all the analyzed samples (27 out of 30), depending on the degree of preservation and the sequencing throughput, and to get a reliable set of variants allowing haplogroup prediction. Additionally, we provide guidelines to deal with possible artifact sources, including nuclear mitochondrial sequence (NumtS) contamination, an often-neglected issue in ancient mtDNA surveys. Potential heteroplasmy levels were also estimated, although most variants were likely homoplasmic, and validated by data simulations, proving that new sequencing technologies and software are sensitive enough to detect partially mutated sites in ancient genomes and discriminate true variants from artifacts. A thorough functional annotation of detected and filtered mtDNA variants was also performed for a comprehensive evaluation of these ancient samples.

Genome of a middle Holocene hunter-gatherer from Wallacea

Much remains unknown about the population history of early modern humans in southeast Asia, where the archaeological record is sparse and the tropical climate is inimical to the preservation of ancient human DNA¹. So far, only two low-coverage pre-Neolithic human genomes have been sequenced from this region. Both are from mainland Hòabìnhian hunter-gatherer sites: Pha Faen in Laos, dated to 7939-7751 calibrated years before present (yr cal BP; present taken as AD 1950), and Gua Cha in Malaysia (4.4-4.2 kyr cal BP)¹. Here we report, to our knowledge, the first ancient human genome from Wallacea, the oceanic island zone between the Sunda Shelf (comprising mainland southeast Asia and the continental islands of western Indonesia) and Pleistocene Sahul (Australia-New Guinea). We extracted DNA from the petrous bone of a young female hunter-gatherer buried 7.3-7.2 kyr cal BP at the limestone cave of Leang Panninge² in South Sulawesi, Indonesia. Genetic analyses show that this pre-Neolithic forager, who is associated with the 'Toalean' technocomplex^3,4, shares most genetic drift and morphological similarities with present-day Papuan and Indigenous Australian groups, yet represents a previously unknown divergent human lineage that branched off around the time of the split between these populations approximately 37,000 years ago⁵. We also describe Denisovan and deep Asian-related ancestries in the Leang Panninge genome, and infer their large-scale displacement from the region today.