Typing single polymorphic nucleotides in mitochondrial DNA as a way to access Middle Pleistocene DNA

How the Sima de los Huesos was won

Although the first discovery of a human fossil in the Sima de los Huesos took place in 1976, systematic excavations did not begin there until 1984. Since then, this site has been continuously excavated in month-long camps. The site is dated by different radiometric techniques to between 430,000 and 300,000 years ago. Until the 2023 campaign, just over 7000 human fossils have been recovered, constituting the largest collection of fossils prior to Homo sapiens ever discovered. The fossils correspond to a minimum of 29 individuals of both sexes and different ages at death, from preadolescents to a specimen of advanced age. Comparative anatomy and ancient DNA studies both suggest that this is a population closely related to Homo neanderthalensis. The great variety and extraordinary quality of the fossils recovered have allowed us to carry out a series of investigations that have greatly increased our knowledge about the evolution of Homo in the Middle Pleistocene. Among the most important discoveries, it has been possible to establish body size and proportions, the confirmation that the origin of the accumulation of human fossils was of an anthropic nature, that those past humans took care of disabled individuals and who were capable of having an oral language almost as complex and efficient as that of our own species.

Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

Ancient genomics

The past decade has witnessed a revolution in ancient DNA (aDNA) research. Although the field's focus was previously limited to mitochondrial DNA and a few nuclear markers, whole genome sequences from the deep past can now be retrieved. This breakthrough is tightly connected to the massive sequence throughput of next generation sequencing platforms and the ability to target short and degraded DNA molecules. Many ancient specimens previously unsuitable for DNA analyses because of extensive degradation can now successfully be used as source materials. Additionally, the analytical power obtained by increasing the number of sequence reads to billions effectively means that contamination issues that have haunted aDNA research for decades, particularly in human studies, can now be efficiently and confidently quantified. At present, whole genomes have been sequenced from ancient anatomically modern humans, archaic hominins, ancient pathogens and megafaunal species. Those have revealed important functional and phenotypic information, as well as unexpected adaptation, migration and admixture patterns. As such, the field of aDNA has entered the new era of genomics and has provided valuable information when testing specific hypotheses related to the past.

A 400,000-year-old mitochondrial genome questions phylogenetic relationships amongst archaic hominins: using the latest advances in ancient genomics, the mitochondrial genome sequence of a 400,000-year-old hominin has been deciphered

By combining state-of-the-art approaches in ancient genomics, Meyer and co-workers have reconstructed the mitochondrial sequence of an archaic hominin that lived at Sierra de Atapuerca, Spain about 400,000 years ago. This achievement follows recent advances in molecular anthropology that delivered the genome sequence of younger archaic hominins, such as Neanderthals and Denisovans. Molecular phylogenetic reconstructions placed the Atapuercan as a sister group to Denisovans, although its morphology suggested closer affinities with Neanderthals. In addition to possibly challenging our interpretation of the fossil record, this study confirms that genomic information can be recovered from extremely damaged DNA molecules, even in the presence of significant levels of human contamination. Together with the recent characterization of a 700,000-year-old horse genome, this study opens the Middle Pleistocene to genomics, thereby extending the scope of ancient DNA to the last million years.

A mitochondrial genome sequence of a hominin from Sima de los Huesos

Excavations of a complex of caves in the Sierra de Atapuerca in northern Spain have unearthed hominin fossils that range in age from the early Pleistocene to the Holocene. One of these sites, the 'Sima de los Huesos' ('pit of bones'), has yielded the world's largest assemblage of Middle Pleistocene hominin fossils, consisting of at least 28 individuals dated to over 300,000 years ago. The skeletal remains share a number of morphological features with fossils classified as Homo heidelbergensis and also display distinct Neanderthal-derived traits. Here we determine an almost complete mitochondrial genome sequence of a hominin from Sima de los Huesos and show that it is closely related to the lineage leading to mitochondrial genomes of Denisovans, an eastern Eurasian sister group to Neanderthals. Our results pave the way for DNA research on hominins from the Middle Pleistocene.

Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments

Although an inverse relationship is expected in ancient DNA samples between the number of surviving DNA fragments and their length, ancient DNA sequencing libraries are strikingly deficient in molecules shorter than 40 bp. We find that a loss of short molecules can occur during DNA extraction and present an improved silica-based extraction protocol that enables their efficient retrieval. In combination with single-stranded DNA library preparation, this method enabled us to reconstruct the mitochondrial genome sequence from a Middle Pleistocene cave bear (Ursus deningeri) bone excavated at Sima de los Huesos in the Sierra de Atapuerca, Spain. Phylogenetic reconstructions indicate that the U. deningeri sequence forms an early diverging sister lineage to all Western European Late Pleistocene cave bears. Our results prove that authentic ancient DNA can be preserved for hundreds of thousand years outside of permafrost. Moreover, the techniques presented enable the retrieval of phylogenetically informative sequences from samples in which virtually all DNA is diminished to fragments shorter than 50 bp.

Improving the performance of true single molecule sequencing for ancient DNA

BACKGROUND

Second-generation sequencing technologies have revolutionized our ability to recover genetic information from the past, allowing the characterization of the first complete genomes from past individuals and extinct species. Recently, third generation Helicos sequencing platforms, which perform true Single-Molecule DNA Sequencing (tSMS), have shown great potential for sequencing DNA molecules from Pleistocene fossils. Here, we aim at improving even further the performance of tSMS for ancient DNA by testing two novel tSMS template preparation methods for Pleistocene bone fossils, namely oligonucleotide spiking and treatment with DNA phosphatase.

RESULTS

We found that a significantly larger fraction of the horse genome could be covered following oligonucleotide spiking however not reproducibly and at the cost of extra post-sequencing filtering procedures and skewed %GC content. In contrast, we showed that treating ancient DNA extracts with DNA phosphatase improved the amount of endogenous sequence information recovered per sequencing channel by up to 3.3-fold, while still providing molecular signatures of endogenous ancient DNA damage, including cytosine deamination and fragmentation by depurination. Additionally, we confirmed the existence of molecular preservation niches in large bone crystals from which DNA could be preferentially extracted.

CONCLUSIONS

We propose DNA phosphatase treatment as a mechanism to increase sequence coverage of ancient genomes when using Helicos tSMS as a sequencing platform. Together with mild denaturation temperatures that favor access to endogenous ancient templates over modern DNA contaminants, this simple preparation procedure can improve overall Helicos tSMS performance when damaged DNA templates are targeted.

Using classical population genetics tools with heterochroneous data: time matters!

BACKGROUND

New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.

METHODOLOGY/PRINCIPAL FINDINGS

We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)).

CONCLUSIONS/SIGNIFICANCE

Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.

Surprising migration and population size dynamics in ancient Iberian brown bears (Ursus arctos)

The endangered brown bear populations (Ursus arctos) in Iberia have been suggested to be the last fragments of the brown bear population that served as recolonization stock for large parts of Europe during the Pleistocene. Conservation efforts are intense, and results are closely monitored. However, the efforts are based on the assumption that the Iberian bears are a unique unit that has evolved locally for an extended period. We have sequenced mitochondrial DNA (mtDNA) from ancient Iberian bear remains and analyzed them as a serial dataset, monitoring changes in diversity and occurrence of European haplogroups over time. Using these data, we show that the Iberian bear population has experienced a dynamic, recent evolutionary history. Not only has the population undergone mitochondrial gene flow from other European brown bears, but the effective population size also has fluctuated substantially. We conclude that the Iberian bear population has been a fluid evolutionary unit, developed by gene flow from other populations and population bottlenecks, far from being in genetic equilibrium or isolated from other brown bear populations. Thus, the current situation is highly unusual and the population may in fact be isolated for the first time in its history.

Does the 43 bp sequence from an 800,000 year old cretan dwarf elephantid really rewrite the textbook on mammoths?

Pigmy elephants inhabited the islands from the Mediterranean region during the Pleistocene period but became extinct in the course of the Holocene. Despite striking distinctive anatomical characteristics related to insularity, some similarities with the lineage of extant Asian elephants have suggested that pigmy elephants could be most probably seen as members of the genus Elephas. Poulakakis et al (2006) have recently challenged this view by recovering a short mtDNA sequence from an 800 000 year old fossil of the Cretan pigmy elephant (Elephas creticus). According to the authors of this study, a deep taxonomic revision of Cretan dwarf elephants would be needed, as the sequence exhibits clear affinities with woolly mammoth haplotypes. However, we point here many aspects that seriously weaken the strength of the ancient DNA evidence reported.