Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

Ancient pathogens provide a window into health and well-being.

This perspective draws on the record of ancient pathogen genomes and microbiomes illuminating patterns of infectious disease over the course of the Holocene in order to address the following question. How did major changes in living circumstances involving the transition to and intensification of farming alter pathogens and their distributions? Answers to this question via ancient DNA research provide a rapidly expanding picture of pathogen evolution and in concert with archaeological and historical data, give a temporal and behavioral context for heath in the past that is relevant for challenges facing the world today, including the rise of novel pathogens.

Uncovering Signals of Positive Selection in Peruvian Populations from Three Ecological Regions.

Peru hosts extremely diverse ecosystems which can be broadly classified into the following three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide single-nucleotide polymorphism-array data from 286 Peruvians living across the three ecoregions and analyzed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44, and ICAM1 in the coast; CBLB and PRDM1 in the rainforest; and BRD2, HLA-DOA, HLA-DPA1 regions in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

Fostering Responsible Research on Ancient DNA.

Anticipating and addressing the social implications of scientific work is a fundamental responsibility of all scientists. However, expectations for ethically sound practices can evolve over time as the implications of science come to be better understood. Contemporary researchers who work with ancient human remains, including those who conduct ancient DNA research, face precisely this challenge as it becomes clear that practices such as community engagement are needed to address the important social implications of this work. To foster and promote ethical engagement between researchers and communities, we offer five practical recommendations for ancient DNA researchers: (1) formally consult with communities; (2) address cultural and ethical considerations; (3) engage communities and support capacity building; (4) develop plans to report results and manage data; and (5) develop plans for long-term responsibility and stewardship. Ultimately, every member of a research team has an important role in fostering ethical research on ancient DNA.

Mitogenomic phylogeny of Callithrix with special focus on human transferred taxa.

BACKGROUND: Callithrix marmosets are a relatively young primate radiation, whose phylogeny is not yet fully resolved. These primates are naturally para- and allopatric, but three species with highly invasive potential have been introduced into the southeastern Brazilian Atlantic Forest by the pet trade. There, these species hybridize with each other and endangered, native congeners. We aimed here to reconstruct a robust Callithrix phylogeny and divergence time estimates, and identify the biogeographic origins of autochthonous and allochthonous Callithrix mitogenome lineages. We sequenced 49 mitogenomes from four species (C. aurita, C. geoffroyi, C. jacchus, C. penicillata) and anthropogenic hybrids (C. aurita x Callithrix sp., C. penicillata x C. jacchus, Callithrix sp. x Callithrix sp., C. penicillata x C. geoffroyi) via Sanger and whole genome sequencing. We combined these data with previously published Callithrix mitogenomes to analyze five Callithrix species in total. RESULTS: We report the complete sequence and organization of the C. aurita mitogenome. Phylogenetic analyses showed that C. aurita was the first to diverge within Callithrix 3.54 million years ago (Ma), while C. jacchus and C. penicillata lineages diverged most recently 0.5 Ma as sister clades. MtDNA clades of C. aurita, C. geoffroyi, and C. penicillata show intraspecific geographic structure, but C. penicillata clades appear polyphyletic. Hybrids, which were identified by phenotype, possessed mainly C. penicillata or C. jacchus mtDNA haplotypes. The biogeographic origins of mtDNA haplotypes from hybrid and allochthonous Callithrix were broadly distributed across natural Callithrix ranges. Our phylogenetic results also evidence introgression of C. jacchus mtDNA into C. aurita. CONCLUSION: Our robust Callithrix mitogenome phylogeny shows C. aurita lineages as basal and C. jacchus lineages among the most recent within Callithrix. We provide the first evidence that parental mtDNA lineages of anthropogenic hybrid and allochthonous marmosets are broadly distributed inside and outside of the Atlantic Forest. We also show evidence of cryptic hybridization between allochthonous Callithrix and autochthonous C. aurita. Our results encouragingly show that further development of genomic resources will allow to more clearly elucidate Callithrix evolutionary relationships and understand the dynamics of Callithrix anthropogenic introductions into the Brazilian Atlantic Forest.

Ancient DNA Reconstructs the Genetic Legacies of Precontact Puerto Rico Communities.

Indigenous peoples have occupied the island of Puerto Rico since at least 3000 BC. Due to the demographic shifts that occurred after European contact, the origin(s) of these ancient populations, and their genetic relationship to present-day islanders, are unclear. We use ancient DNA to characterize the population history and genetic legacies of precontact Indigenous communities from Puerto Rico. Bone, tooth, and dental calculus samples were collected from 124 individuals from three precontact archaeological sites: Tibes, Punta Candelero, and Paso del Indio. Despite poor DNA preservation, we used target enrichment and high-throughput sequencing to obtain complete mitochondrial genomes (mtDNA) from 45 individuals and autosomal genotypes from two individuals. We found a high proportion of Native American mtDNA haplogroups A2 and C1 in the precontact Puerto Rico sample (40% and 44%, respectively). This distribution, as well as the haplotypes represented, supports a primarily Amazonian South American origin for these populations and mirrors the Native American mtDNA diversity patterns found in present-day islanders. Three mtDNA haplotypes from precontact Puerto Rico persist among Puerto Ricans and other Caribbean islanders, indicating that present-day populations are reservoirs of precontact mtDNA diversity. Lastly, we find similarity in autosomal ancestry patterns between precontact individuals from Puerto Rico and the Bahamas, suggesting a shared component of Indigenous Caribbean ancestry with close affinity to South American populations. Our findings contribute to a more complete reconstruction of precontact Caribbean population history and explore the role of Indigenous peoples in shaping the biocultural diversity of present-day Puerto Ricans and other Caribbean islanders.

The population genomics of within-host Mycobacterium tuberculosis.

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis.

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates.

OBJECTIVES: Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups. MATERIALS AND METHODS: First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6). RESULTS: Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits. DISCUSSION: Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Advancing the understanding of treponemal disease in the past and present.

Syphilis was perceived to be a new disease in Europe in the late 15th century, igniting a debate about its origin that continues today in anthropological, historical, and medical circles. We move beyond this age-old debate using an interdisciplinary approach that tackles broader questions to advance the understanding of treponemal infection (syphilis, yaws, bejel, and pinta). How did the causative organism(s) and humans co-evolve? How did the related diseases caused by Treponema pallidum emerge in different parts of the world and affect people across both time and space? How are T. pallidum subspecies related to the treponeme causing pinta? The current state of scholarship in specific areas is reviewed with recommendations made to stimulate future work. Understanding treponemal biology, genetic relationships, epidemiology, and clinical manifestations is crucial for vaccine development today and for investigating the distribution of infection in both modern and past populations. Paleopathologists must improve diagnostic criteria and use a standard approach for recording skeletal lesions on archaeological human remains. Adequate contextualization of cultural and environmental conditions is necessary, including site dating and justification for any corrections made for marine or freshwater reservoir effects. Biogeochemical analyses may assess aquatic contributions to diet, physiological changes arising from treponemal disease and its treatments (e.g., mercury), or residential mobility of those affected. Shifting the focus from point of origin to investigating who is affected (e.g., by age/sex or socioeconomic status) and disease distribution (e.g., coastal/ inland, rural/urban) will advance our understanding of the treponemal disease and its impact on people through time.

The efficacy of whole human genome capture on ancient dental calculus and dentin.

OBJECTIVES: Dental calculus is among the richest known sources of ancient DNA in the archaeological record. Although most DNA within calculus is microbial, it has been shown to contain sufficient human DNA for the targeted retrieval of whole mitochondrial genomes. Here, we explore whether calculus is also a viable substrate for whole human genome recovery using targeted enrichment techniques. MATERIALS AND METHODS: Total DNA extracted from 24 paired archaeological human dentin and calculus samples was subjected to whole human genome enrichment using in-solution hybridization capture and high-throughput sequencing. RESULTS: Total DNA from calculus exceeded that of dentin in all cases, and although the proportion of human DNA was generally lower in calculus, the absolute human DNA content of calculus and dentin was not significantly different. Whole genome enrichment resulted in up to four-fold enrichment of the human endogenous DNA content for both dentin and dental calculus libraries, albeit with some loss in complexity. Recovering more on-target reads for the same sequencing effort generally improved the quality of downstream analyses, such as sex and ancestry estimation. For nonhuman DNA, comparison of phylum-level microbial community structure revealed few differences between precapture and postcapture libraries, indicating that off-target sequences in human genome-enriched calculus libraries may still be useful for oral microbiome reconstruction. DISCUSSION: While ancient human dental calculus does contain endogenous human DNA sequences, their relative proportion is low when compared with other skeletal tissues. Whole genome enrichment can help increase the proportion of recovered human reads, but in this instance enrichment efficiency was relatively low when compared with other forms of capture. We conclude that further optimization is necessary before the method can be routinely applied to archaeological samples.

Fetal Hemoglobin Modulators May Be Associated With Symptomology of Football Players with Sickle Cell Trait.

OBJECTIVES: This study investigates whether genetic modifiers previously shown to influence adult fetal hemoglobin (HbF) levels and glucose-6-phosphate dehydrogenase deficiency were associated with variable symptomology in a small sample of collegiate football players with sickle cell trait. METHODS: Survey data on self-assessed symptoms and genotype data from five single nucleotide polymorphisms (SNPs) related to HbF production and two SNPs that cause glucose-6-phosphate dehydrogenase deficiency were collected from current and former college football players. RESULTS: In this sample, SNPs found within the ß-globin gene cluster were found to be associated with a previous diagnosis of exertional sickling and experience of extreme heat during and after training. rs10189857 in the BCL11A gene was associated with body mass index and weight and with experiencing extreme thirst during and after training. No significant correlations were found between the other SNPs and symptoms within this sample. CONCLUSIONS: These findings show that genetic variation known to affect sickle cell disease symptomology may partly explain why some football players with sickle cell trait experience adverse clinical outcomes during periods of extreme physical exertion and others do not.

In-solution Y-chromosome capture-enrichment on ancient DNA libraries.

BACKGROUND: As most ancient biological samples have low levels of endogenous DNA, it is advantageous to enrich for specific genomic regions prior to sequencing. One approach-in-solution capture-enrichment-retrieves sequences of interest and reduces the fraction of microbial DNA. In this work, we implement a capture-enrichment approach targeting informative regions of the Y chromosome in six human archaeological remains excavated in the Caribbean and dated between 200 and 3000 years BP. We compare the recovery rate of Y-chromosome capture (YCC) alone, whole-genome capture followed by YCC (WGC + YCC) versus non-enriched (pre-capture) libraries. RESULTS: The six samples show different levels of initial endogenous content, with very low (< 0.05%, 4 samples) or low (0.1-1.54%, 2 samples) percentages of sequenced reads mapping to the human genome. We recover 12-9549 times more targeted unique Y-chromosome sequences after capture, where 0.0-6.2% (WGC + YCC) and 0.0-23.5% (YCC) of the sequence reads were on-target, compared to 0.0-0.00003% pre-capture. In samples with endogenous DNA content greater than 0.1%, we found that WGC followed by YCC (WGC + YCC) yields lower enrichment due to the loss of complexity in consecutive capture experiments, whereas in samples with lower endogenous content, the libraries' initial low complexity leads to minor proportions of Y-chromosome reads. Finally, increasing recovery of informative sites enabled us to assign Y-chromosome haplogroups to some of the archeological remains and gain insights about their paternal lineages and origins. CONCLUSIONS: We present to our knowledge the first in-solution capture-enrichment method targeting the human Y-chromosome in aDNA sequencing libraries. YCC and WGC + YCC enrichments lead to an increase in the amount of Y-DNA sequences, as compared to libraries not enriched for the Y-chromosome. Our probe design effectively recovers regions of the Y-chromosome bearing phylogenetically informative sites, allowing us to identify paternal lineages with less sequencing than needed for pre-capture libraries. Finally, we recommend considering the endogenous content in the experimental design and avoiding consecutive rounds of capture, as clonality increases considerably with each round.

Novel anelloviruses identified in buccal swabs of Antarctic fur seals.

Viral diversity associated with Antarctic wildlife remains poorly studied. Nonetheless, over the past 5 years, there has been a concerted effort using viral metagenomics approaches to identify and characterize viruses associated with Antarctic pinniped and avian species. Here we used a viral metagenomics approach to identify circular DNA viruses in buccal swab samples from Antarctic fur seals (Arctocephalus gazella) breeding on Livingston Island, Antarctica during the 2016/2017 field season. We identified two new lineages of anelloviruses, torque teno Arctocephalus gazella virus (TTAgV)-1 (2083 nts) and -2 (2127-2143 nts), which are most closely related to and cluster with a previously identified anellovirus associated with California sea lions (Zalophus californianus) sharing ~ 60% genome-wide pairwise identity. The ORF1 of TTAgVs share 26-41% amino acid similarity to the ORF1 of other previously identified pinniped-associated anelloviruses.

Comparison of two ancient DNA extraction protocols for skeletal remains from tropical environments.

OBJECTIVES: The tropics harbor a large part of the world's biodiversity and have a long history of human habitation. However, paleogenomics research in these climates has been constrained so far by poor ancient DNA yields. Here we compare the performance of two DNA extraction methods on ancient samples of teeth and petrous portions excavated from tropical and semi-tropical sites in Tanzania, Mexico, and Puerto Rico (N = 12). MATERIALS AND METHODS: All samples were extracted twice, built into double-stranded sequencing libraries, and shotgun sequenced on the Illumina HiSeq 2500. The first extraction protocol, Method D, was previously designed for recovery of ultrashort DNA fragments from skeletal remains. The second, Method H, modifies the first by adding an initial EDTA wash and an extended digestion and decalcification step. RESULTS: No significant difference was found in overall ancient DNA yields or post-mortem damage patterns recovered from samples extracted with either method, irrespective of tissue type. However, Method H samples had higher endogenous content and more mapped reads after quality-filtering, but also higher clonality. In contrast, samples extracted with Method D had shorter average DNA fragments. DISCUSSION: Both methods successfully recovered endogenous ancient DNA. But, since surviving DNA in ancient or historic remains from tropical contexts is extremely fragmented, our results suggest that Method D is the optimal choice for working with samples from warm and humid environments. Additional optimization of extraction conditions and further testing of Method H with different types of samples may allow for improvement of this protocol in the future.

Facing the facts: The Runx2 gene is associated with variation in facial morphology in primates.

The phylogenetic and adaptive factors that cause variation in primate facial form-including differences among the major primate clades and variation related to feeding and/or social behavior-are relatively well understood. However, comparatively little is known about the genetic mechanisms that underlie diversity in facial form in primates. Because it is essential for osteoblastic differentiation and skeletal development, the runt-related transcription factor 2 (Runx2) is one gene that may play a role in these genetic mechanisms. Specifically, polymorphisms in the QA ratio (determined by the ratio of the number of polyglutamines to polyalanines in one functional domain of Runx2) have been shown to be correlated with variation in facial length and orientation in other mammal groups. However, to date, the relationship between variation in this gene and variation in facial form in primates has not been explicitly tested. To test the hypothesis that the QA ratio is correlated with facial form in primates, the current study quantified the QA ratio, facial length, and facial angle in a sample of 33 primate species and tested for correlation using phylogenetic generalized least squares. The results indicate that the QA ratio of the Runx2 gene is positively correlated with variation in relative facial length in anthropoid primates. However, no correlation was found in strepsirrhines, and there was no correlation between facial angle and the QA ratio in any groups. These results suggest that, in primates, the QA ratio of the Runx2 gene may play a role in modulating facial size, but not facial orientation. This study therefore provides important clues about the genetic and developmental mechanisms that may underlie variation in facial form in primates.

Phylogeny of the jacchus group of Callithrix marmosets based on complete mitochondrial genomes.

OBJECTIVES: Two subgroups make up the marmoset genus Callithrix. The "aurita" group is composed of two species, whereas evolutionary relationships among the four species of the "jacchus" group remain unclear. To uncover these relationships, we first sequenced mitochondrial genomes for C. kuhlii and C. penicillata to complement data available for congeners. We then constructed a phylogenetic tree based on mtDNA heavy chain protein coding genes from several primates to untangle species relationships and estimate divergence times of the jacchus group. MATERIALS AND METHODS: MtDNA genomes of C. kuhlii and C. penicillata were Sanger sequenced. These Callithrix mitogenomes were combined with other publically available primate mtDNA genomes. Phylogenies were produced using maximum likelihood and Bayesian inference. Finally, divergence times within the jacchus group of marmosets were estimated with Bayesian inference. RESULTS: In our phylogenetic tree, C. geoffroyi was the sister to all other jacchus group species, followed by C. kuhlii, while C. jacchus and C. penicillata diverged most recently. Bayesian inference showed that C. jacchus and C. penicillata diverged approximately 0.70 MYA and that the jacchus group radiated approximately 1.30 MYA. DISCUSSION: Callithrix nuclear and mtDNA phylogenies frequently result in polytomies and paraphyly. Here, we present a well-supported phylogenetic tree based on mitochondrial genome sequences, which facilitates the understanding of the divergence of the jacchus marmosets. Our results demonstrate how mitochondrial genomes can enrich Callithrix phylogenetic studies by alleviating some of the difficulties faced by previous mtDNA studies and allow formulation of hypotheses to test further under larger genomic-scale analyses.

Application of PE-RADSeq to the study of genomic diversity and divergence of two Brazilian marmoset species (Callithrix jacchus and C. penicillata).

Callithrix jacchus and C. penicillata are among the smallest anthropoid primates, are highly specialized tree gougers, and largely occupy Brazil's most extreme, semi-arid biomes. However, the underlying genomic factors that underpin the evolution of these species and their unique traits are under-investigated. Additionally, exotic populations of these two species are widely established throughout Brazil and hybridize with threatened native congers. Thus, both genomic and conservation factors call for a better understanding of C. jacchus and C. penicillata evolution. Here, we applied PE-RADseq to characterize genomic variation in these two species, using six C. jacchus and seven C. penicillata individuals. We identified an average of 7,463 and 5,180 SNPs/individual in C. penicillata and C. jacchus, respectively, and also found 1,395 variable sites that were represented in both species. C. penicillata showed overall higher levels of genetic diversity than C. jacchus at the variable sites present in both species. Additionally, among these variable sites, 106 showed relative interspecific divergence levels that were significantly higher than the genome-wide average. We further compared relative and absolute divergence for C. penicillata and C. jacchus between RAD loci associated with the 106 significantly diverged variable sites and all other RAD loci present in both species. The former RAD loci set showed significantly elevated relative and absolute divergence measures in comparison to the latter set. This convergence suggests that C. jacchus and C. penicillata may have diverged under a scenario of gene flow under secondary contact. Here, we demonstrate that RADseq is an efficient method to simultaneously discover and genotype a large number of markers and validate the utility of RADseq for examining Callithrix evolution.

Successful enrichment and recovery of whole mitochondrial genomes from ancient human dental calculus.

OBJECTIVES: Archaeological dental calculus is a rich source of host-associated biomolecules. Importantly, however, dental calculus is more accurately described as a calcified microbial biofilm than a host tissue. As such, concerns regarding destructive analysis of human remains may not apply as strongly to dental calculus, opening the possibility of obtaining human health and ancestry information from dental calculus in cases where destructive analysis of conventional skeletal remains is not permitted. Here we investigate the preservation of human mitochondrial DNA (mtDNA) in archaeological dental calculus and its potential for full mitochondrial genome (mitogenome) reconstruction in maternal lineage ancestry analysis. MATERIALS AND METHODS: Extracted DNA from six individuals at the 700-year-old Norris Farms #36 cemetery in Illinois was enriched for mtDNA using in-solution capture techniques, followed by Illumina high-throughput sequencing. RESULTS: Full mitogenomes (7-34×) were successfully reconstructed from dental calculus for all six individuals, including three individuals who had previously tested negative for DNA preservation in bone using conventional PCR techniques. Mitochondrial haplogroup assignments were consistent with previously published findings, and additional comparative analysis of paired dental calculus and dentine from two individuals yielded equivalent haplotype results. All dental calculus samples exhibited damage patterns consistent with ancient DNA, and mitochondrial sequences were estimated to be 92-100% endogenous. DNA polymerase choice was found to impact error rates in downstream sequence analysis, but these effects can be mitigated by greater sequencing depth. DISCUSSION: Dental calculus is a viable alternative source of human DNA that can be used to reconstruct full mitogenomes from archaeological remains. Am J Phys Anthropol 160:220-228, 2016. © 2016 The Authors American Journal of Physical Anthropology Published by Wiley Periodicals, Inc.