Ancient DNA extracted from the sediments of archaeological sites (sedaDNA) can provide fine-grained information about the composition of past ecosystems and human site use, even in the absence of visible remains. However, the growing amount of available sequencing data and the nature of the data obtained from archaeological sediments pose several computational challenges; among these, the rapid and accurate taxonomic classification of sequences. While alignment-based taxonomic classifiers remain the standard in sedaDNA analysis pipelines, they are too computationally expensive for the processing of large numbers of sedaDNA sequences. In contrast, alignment-free methods offer fast classification but suffer from higher false-positive rates. To address these limits, we developed quicksand, an open-source Nextflow pipeline designed for rapid and accurate taxonomic classification of mammalian mitochondrial DNA in sedaDNA samples. quicksand combines fast alignment-free classification using KrakenUniq with post-classification mapping, filtering, and ancient DNA authentication. Based on simulations and reanalyses of published datasets, we demonstrate that quicksand achieves accuracy and sensitivity comparable to or better than existing methods, while significantly reducing runtime. quicksand offers an easy workflow for large-scale screening of sedaDNA samples for archaeological research and is freely available at https://github.com/mpieva/quicksand.

Denisova Cave in southern Siberia is the type locality of the Denisovans, an archaic hominin group who were related to Neanderthals 1 – 4 . The dozen hominin remains recovered from the deposits also include Neanderthals 5 , 6 and the child of a Neanderthal and a Denisovan 7 , which suggests that Denisova Cave was a contact zone between these archaic hominins. However, uncertainties persist about the order in which these groups appeared at the site, the timing and environmental context of hominin occupation, and the association of particular hominin groups with archaeological assemblages 5 , 8 – 11 . Here we report the analysis of DNA from 728 sediment samples that were collected in a grid-like manner from layers dating to the Pleistocene epoch. We retrieved ancient faunal and hominin mitochondrial (mt)DNA from 685 and 175 samples, respectively. The earliest evidence for hominin mtDNA is of Denisovans, and is associated with early Middle Palaeolithic stone tools that were deposited approximately 250,000 to 170,000 years ago; Neanderthal mtDNA first appears towards the end of this period. We detect a turnover in the mtDNA of Denisovans that coincides with changes in the composition of faunal mtDNA, and evidence that Denisovans and Neanderthals occupied the site repeatedly—possibly until, or after, the onset of the Initial Upper Palaeolithic at least 45,000 years ago, when modern human mtDNA is first recorded in the sediments.

Bones and teeth are important sources of Pleistocene hominin DNA, but are rarely recovered at archaeological sites. Mitochondrial DNA (mtDNA) has been retrieved from cave sediments but provides limited value for studying population relationships. We therefore developed methods for the enrichment and analysis of nuclear DNA from sediments and applied them to cave deposits in western Europe and southern Siberia dated to between 200,000 and 50,000 years ago. We detected a population replacement in northern Spain about 100,000 years ago, which was accompanied by a turnover of mtDNA. We also identified two radiation events in Neanderthal history during the early part of the Late Pleistocene. Our work lays the ground for studying the population history of ancient hominins from trace amounts of nuclear DNA in sediments.

Selenium (Se), a sparse element on earth, is an essential micronutrient in the vertebrate diet and its intake depends on its content in soils and waters worldwide. Selenium is required due to its function in selenoproteins, which contain selenocysteine (Sec), the 21st amino acid in the genetic code, as one of their constituent residues. Selenocysteine is analogous to the amino acid cysteine (Cys), which uses the abounding element sulfur instead. Despite the irregular distribution of Se worldwide, its distinct biochemical properties have made the substitution of Sec for Cys rare in vertebrate proteins. Still, vertebrates inhabited environments with different amounts of Se and may have distinctly adapted to it. To address this question, we compared the evolutionary forces acting on the coding sequences of selenoprotein genes and genes that regulate Se between vertebrate clades and between the Se-dependent genes and their paralogs with Cys. We find that the strength of natural selection in genes that use or regulate Se is distinct between land vertebrates and teleost fishes and more variable than in the Cys paralogs, particularly in genes involved in the preferential supply of Se to some organs and the tissue-specific expression of selenoproteins. This is compatible with vertebrates adapting to Se scarcity in land and its abundance in waters. In agreement, teleost fishes duplicated and subfunctionalized or neofunctionalized selenoprotein genes and maintained their capacity for Se transport in the body, which declined (under neutrality) for millions of years in terrestrial vertebrates. Dietary Se has thus distinctly shaped vertebrate evolution.

The influence of positive selection sweeps in human evolution is increasingly debated, although our ability to detect them is hampered by inherent uncertainties in the timing of past events. Ancient genomes provide snapshots of allele frequencies in the past and can help address this question. We combine modern and ancient genomic data in a simple statistic (DAnc) to time allele frequency changes, and investigate the role of drift and adaptation in population differentiation. Only 30% of the most strongly differentiated alleles between Africans and Eurasians changed in frequency during the colonization of Eurasia, but in Europe these alleles are enriched in genic and putatively functional alleles to an extent only compatible with local adaptation. Adaptive alleles--especially those associated with pigmentation--are mostly of hunter-gatherer origin, although lactose persistence arose in a haplotype present in farmers. These results provide evidence for a role of local adaptation in human population differentiation.

As humans migrated around the world, they came to inhabit environments that differ widely in the soil levels of certain micronutrients, including selenium (Se). Coupled with cultural variation in dietary practices, these migrations have led to a wide range of Se intake levels in populations around the world. Both excess and deficiency of Se in the diet can have adverse health consequences in humans, with severe Se deficiency resulting in diseases of the bone and heart. Se is required by humans mainly due to its function in selenoproteins, which contain the amino acid selenocysteine as one of their constituent residues. To understand the evolution of the use of this micronutrient in humans, we surveyed the patterns of polymorphism in all selenoprotein genes and genes involved in their regulation in 50 human populations. We find that single nucleotide polymorphisms from populations in Asia, particularly in populations living in the extreme Se-deficient regions of China, have experienced concerted shifts in their allele frequencies. Such differentiation in allele frequencies across genes is not observed in other regions of the world and is not expected under neutral evolution, being better explained by the action of recent positive selection. Thus, recent changes in the use and regulation of Se may harbor the genetic adaptations that helped humans inhabit environments that do not provide adequate levels of Se in the diet.

SelenoDB (http://www.selenodb.org) aims to provide high-quality annotations of selenoprotein genes, proteins and SECIS elements. Selenoproteins are proteins that contain the amino acid selenocysteine (Sec) and the first release of the database included annotations for eight species. Since the release of SelenoDB 1.0 many new animal genomes have been sequenced. The annotations of selenoproteins in new genomes usually contain many errors in major databases. For this reason, we have now fully annotated selenoprotein genes in 58 animal genomes. We provide manually curated annotations for human selenoproteins, whereas we use an automatic annotation pipeline to annotate selenoprotein genes in other animal genomes. In addition, we annotate the homologous genes containing cysteine (Cys) instead of Sec. Finally, we have surveyed genetic variation in the annotated genes in humans. We use exon capture and resequencing approaches to identify single-nucleotide polymorphisms in more than 50 human populations around the world. We thus present a detailed view of the genetic divergence of Sec- and Cys-containing genes in animals and their diversity in humans. The addition of these datasets into the second release of the database provides a valuable resource for addressing medical and evolutionary questions in selenium biology.

BACKGROUND:The Human Microbiome Project (HMP) is one of the U.S. National Institutes of Health Roadmap for Medical Research. Primary interests of the HMP include the distinctiveness of different gut microbiomes, the factors influencing microbiome diversity, and the functional redundancies of the members of human microbiotas. In this present work, we contribute to these interests by characterizing two extinct human microbiotas.METHODOLOGY/PRINCIPAL FINDINGS:We examine two paleofecal samples originating from cave deposits in Durango Mexico and dating to approximately 1300 years ago. Contamination control is a serious issue in ancient DNA research; we use a novel approach to control contamination. After we determined that each sample originated from a different human, we generated 45 thousand shotgun DNA sequencing reads. The phylotyping and functional analysis of these reads reveals a signature consistent with the modern gut ecology. Interestingly, inter-individual variability for phenotypes but not functional pathways was observed. The two ancient samples have more similar functional profiles to each other than to a recently published profile for modern humans. This similarity could not be explained by a chance sampling of the databases.CONCLUSIONS/SIGNIFICANCE:We conduct a phylotyping and functional analysis of ancient human microbiomes, while providing novel methods to control for DNA contamination and novel hypotheses about past microbiome biogeography. We postulate that natural selection has more of an influence on microbiome functional profiles than it does on the species represented in the microbial ecology. We propose that human microbiomes were more geographically structured during pre-Columbian times than today.