Mitogenomic analyses from ancient DNA

First mitogenome phylogeny of the sun bear Helarctos malayanus reveals a deep split between Indochinese and Sundaic lineages

The sun bear Helarctos malayanus is one of the most endangered ursids, and to date classification of sun bear populations has been based almost exclusively on geographic distribution and morphology. The very few molecular studies focussing on this species were limited in geographic scope. Using archival and non-invasively collected sample material, we have added a substantial number of complete or near-complete mitochondrial genome sequences from sun bears of several range countries of the species' distribution. We here report 32 new mitogenome sequences representing sun bears from Cambodia, Thailand, Peninsular Malaysia, Sumatra, and Borneo. Reconstruction of phylogenetic relationships revealed two matrilines that diverged ~295 thousand years ago: one restricted to portions of mainland Indochina (China, Cambodia, Thailand; "Mainland clade"), and one comprising bears from Borneo, Sumatra, Peninsular Malaysia but also Thailand ("Sunda clade"). Generally recent coalescence times in the mitochondrial phylogeny suggest that recent or historical demographic processes have resulted in a loss of mtDNA variation. Additionally, analysis of our data in conjunction with shorter mtDNA sequences revealed that the Bornean sun bear, classified as a distinct subspecies (H. m. euryspilus), does not harbor a distinctive matriline. Further molecular studies of H. malayanus are needed, which should ideally include data from nuclear loci.

A Comprehensive Overview of Baboon Phylogenetic History

Baboons (genus Papio) are an intriguing study system to investigate complex evolutionary processes and the evolution of social systems. An increasing number of studies over the last 20 years has shown that considerable incongruences exist between phylogenies based on morphology, mitochondrial, and nuclear sequence data of modern baboons, and hybridization and introgression have been suggested as the main drivers of these patterns. Baboons, therefore, present an excellent opportunity to study these phenomena and their impact on speciation. Advances both in geographic and genomic coverage provide increasing details on the complexity of the phylogeography of baboons. Here, we compile the georeferenced genetic data of baboons and review the current knowledge on baboon phylogeny, discuss the evolutionary processes that may have shaped the patterns that we observe today, and propose future avenues for research.

Advancements and Challenges in Ancient DNA Research: Bridging the Global North-South Divide

Ancient DNA (aDNA) research first began in 1984 and ever since has greatly expanded our understanding of evolution and migration. Today, aDNA analysis is used to solve various puzzles about the origin of mankind, migration patterns, and the spread of infectious diseases. The incredible findings ranging from identifying the new branches within the human family to studying the genomes of extinct flora and fauna have caught the world by surprise in recent times. However, a closer look at these published results points out a clear Global North and Global South divide. Therefore, through this research, we aim to emphasize encouraging better collaborative opportunities and technology transfer to support researchers in the Global South. Further, the present research also focuses on expanding the scope of the ongoing conversation in the field of aDNA by reporting relevant literature published around the world and discussing the advancements and challenges in the field.

Revised phylogeny from complete mitochondrial genomes of phyllostomid bats resolves subfamilial classification

Classically, molecular phylogenetic trees of Phyllostomidae have been inferred using a combination of a few mitochondrial and nuclear markers. However, there is still uncertainty in the relationships, especially among deep clades within the family. In this study, we provide newly sequenced complete mitochondrial genomes from 26 bat species, including genomes of 23 species reported here for the first time. By carefully analysing these genomes using maximum likelihood and Bayesian methods and different ingroup and outgroup samples, partition schemes and data types, we investigated the robustness and sensitivity of our phylogenetic results. The optimal topologies were those inferred from the complete data matrix of nucleotides, with complex and highly parameterized substitution models and partition schemes. Our results show a statistically robust picture of the evolutionary relationships between phyllostomid subfamilies and clarify hitherto uncertain relationships of Lonchorhininae and Macrotinae.

Mitochondrial Haplogroup Classification of Ancient DNA Samples Using Haplotracker

Mitochondrial DNA haplogroup classification is used to study maternal lineage of ancient human populations. The haplogrouping of ancient DNA is not easy because the DNA is usually found in small pieces in limited quantities. We have developed Haplotracker, a straightforward and efficient high-resolution haplogroup classification tool optimized specifically for ancient DNA samples. Haplotracker offers a user-friendly input interface for multiple mitochondrial DNA sequence fragments in a sample. It provides accurate haplogroup classification with full-length mitochondrial genome sequences and provides high-resolution haplogroup predictions for some fragmented control region sequences using a novel algorithm built on Phylotree mtDNA Build 17 (Phylotree) and our haplotype database (n = 118,869). Its performance for accuracy was demonstrated to be high through haplogroup classification using 8,216 Phylotree full-length and control region mitochondrial DNA sequences compared with HaploGrep 2, one of the most accurate current haplogroup classifiers. Haplotracker provides a novel haplogroup tracking solution for fragmented sequences to track subhaplogroups or verify the haplogroups efficiently. Using Haplotracker, we classified mitochondrial haplogroups to the final subhaplogroup level in nine ancient DNA samples extracted from human skeletal remains found in 2,000-year-old elite Xiongnu cemetery in Northeast Mongolia. Haplotracker can be freely accessed at https://haplotracker.cau.ac.kr.

Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication

Sheep were among the first domesticated animals, but their demographic history is little understood. Here we analyzed nuclear polymorphism and mitochondrial data (mtDNA) from ancient central and west Anatolian sheep dating from Epipaleolithic to late Neolithic, comparatively with modern-day breeds and central Asian Neolithic/Bronze Age sheep (OBI). Analyzing ancient nuclear data, we found that Anatolian Neolithic sheep (ANS) are genetically closest to present-day European breeds relative to Asian breeds, a conclusion supported by mtDNA haplogroup frequencies. In contrast, OBI showed higher genetic affinity to present-day Asian breeds. These results suggest that the east-west genetic structure observed in present-day breeds had already emerged by 6000 BCE, hinting at multiple sheep domestication episodes or early wild introgression in southwest Asia. Furthermore, we found that ANS are genetically distinct from all modern breeds. Our results suggest that European and Anatolian domestic sheep gene pools have been strongly remolded since the Neolithic.

Historical Demographic Processes Dominate Genetic Variation in Ancient Atlantic Cod Mitogenomes

Ancient DNA (aDNA) approaches have been successfully used to infer the long-term impacts of climate change, domestication, and human exploitation in a range of terrestrial species. Nonetheless, studies investigating such impacts using aDNA in marine species are rare. Atlantic cod (Gadus morhua), is an economically important species that has experienced dramatic census population declines during the last century. Here, we investigated 48 ancient mitogenomes from historical specimens obtained from a range of archeological excavations in northern Europe dated up to 6,500 BCE. We compare these mitogenomes to those of 496 modern conspecifics sampled across the North Atlantic Ocean and adjacent seas. Our results confirm earlier observations of high levels of mitogenomic variation and a lack of mutation-drift equilibrium—suggestive of population expansion. Furthermore, our temporal comparison yields no evidence of measurable mitogenomic changes through time. Instead, our results indicate that mitogenomic variation in Atlantic cod reflects past demographic processes driven by major historical events (such as oscillations in sea level) and subsequent gene flow rather than contemporary fluctuations in stock abundance. Our results indicate that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on mitogenomic diversity in a wide-spread marine species with high gene flow such as Atlantic cod. These observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and the importance of genetic variation in implementing conservation strategies. Instead, these observations imply that any measures toward the demographic recovery of Atlantic cod in the eastern Atlantic, will not be constrained by recent loss of historical mitogenomic variation.

Ancient DNA Evidence from China Reveals the Expansion of Pacific Dogs

The ancestral homeland of Australian dingoes and Pacific dogs is proposed to be in South China. However, the location and timing of their dispersal and relationship to dog domestication is unclear. Here, we sequenced 7,000- to 2,000-year-old complete mitochondrial DNA (mtDNA) genomes of 27 ancient canids (one gray wolf and 26 domestic dogs) from the Yellow River and Yangtze River basins (YYRB). These are the first complete ancient mtDNA of Chinese dogs from the cradle of early Chinese civilization. We found that most ancient dogs (18/26) belong to the haplogroup A1b lineage that is found in high frequency in present-day Australian dingoes and precolonial Pacific Island dogs but low frequency in present-day China. Particularly, a 7,000-year-old dog from the Tianluoshan site in Zhejiang province possesses a haplotype basal to the entire haplogroup A1b lineage. We propose that A1b lineage dogs were once widely distributed in the YYRB area. Following their dispersal to South China, and then into Southeast Asia, New Guinea and remote Oceania, they were largely replaced by dogs belonging to other lineages in the last 2,000 years in present-day China, especially North China.

Hybridisation capture allows DNA damage analysis of ancient marine eukaryotes

Marine sedimentary ancient DNA (sedaDNA) is increasingly used to study past ocean ecosystems, however, studies have been severely limited by the very low amounts of DNA preserved in the subseafloor, and the lack of bioinformatic tools to authenticate sedaDNA in metagenomic data. We applied a hybridisation capture 'baits' technique to target marine eukaryote sedaDNA (specifically, phyto- and zooplankton, 'Planktonbaits1'; and harmful algal bloom taxa, 'HABbaits1'), which resulted in up to 4- and 9-fold increases, respectively, in the relative abundance of eukaryotes compared to shotgun sequencing. We further used the bioinformatic tool 'HOPS' to authenticate the sedaDNA component, establishing a new proxy to assess sedaDNA authenticity, "% eukaryote sedaDNA damage", that is positively correlated with subseafloor depth. We used this proxy to report the first-ever DNA damage profiles from a marine phytoplankton species, the ubiquitous coccolithophore Emiliania huxleyi. Our approach opens new avenues for the detailed investigation of long-term change and evolution of marine eukaryotes over geological timescales.

Enhancing mitogenomic phylogeny and resolving the relationships of extinct megafaunal placental mammals

Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.

Absence of evidence is not evidence of absence: Nanopore sequencing and complete assembly of the European lobster (Homarus gammarus) mitogenome uncovers the missing nad2 and a new major gene cluster duplication

BACKGROUND

The recently published complete mitogenome of the European lobster (Homarus gammarus) that was generated using long-range PCR exhibits unusual gene composition (missing nad2) and gene rearrangements among decapod crustaceans with strong implications in crustacean phylogenetics. Such atypical mitochondrial features will benefit greatly from validation with emerging long read sequencing technologies such as Oxford Nanopore that can more accurately identify structural variation.

RESULTS

We re-sequenced the H. gammarus mitogenome on an Oxford Nanopore Minion flowcell and performed a long-read only assembly, generating a complete mitogenome assembly for H. gammarus. In contrast to previous reporting, we found an intact mitochondrial nad2 gene in the H. gammarus mitogenome and showed that its gene organization is broadly similar to that of the American lobster (H. americanus) except for the presence of a large tandemly duplicated region with evidence of pseudogenization in one of each duplicated protein-coding genes.

CONCLUSIONS

Using the European lobster as an example, we demonstrate the value of Oxford Nanopore long read technology in resolving problematic mitogenome assemblies. The increasing accessibility of Oxford Nanopore technology will make it an attractive and useful tool for evolutionary biologists to verify new and existing unusual mitochondrial gene rearrangements recovered using first and second generation sequencing technologies, particularly those used to make phylogenetic inferences of evolutionary scenarios.

Resolving the phylogenetic position of Darwin's extinct ground sloth (Mylodon darwinii) using mitogenomic and nuclear exon data

Mylodon darwinii is the extinct giant ground sloth named after Charles Darwin, who first collected its remains in South America. We have successfully obtained a high-quality mitochondrial genome at 99-fold coverage using an Illumina shotgun sequencing of a 12 880-year-old bone fragment from Mylodon Cave in Chile. Low level of DNA damage showed that this sample was exceptionally well preserved for an ancient subfossil, probably the result of the dry and cold conditions prevailing within the cave. Accordingly, taxonomic assessment of our shotgun metagenomic data showed a very high percentage of endogenous DNA with 22% of the assembled metagenomic contigs assigned to Xenarthra. Additionally, we enriched over 15 kb of sequence data from seven nuclear exons, using target sequence capture designed against a wide xenarthran dataset. Phylogenetic and dating analyses of the mitogenomic dataset including all extant species of xenarthrans and the assembled nuclear supermatrix unambiguously place Mylodon darwinii as the sister-group of modern two-fingered sloths, from which it diverged around 22 million years ago. These congruent results from both the mitochondrial and nuclear data support the diphyly of the two modern sloth lineages, implying the convergent evolution of their unique suspensory behaviour as an adaption to arboreality. Our results offer promising perspectives for whole-genome sequencing of this emblematic extinct taxon.

Assembly of Ancient Mitochondrial Genomes Without a Closely Related Reference Sequence

Recent methodological advances have transformed the field of ancient DNA (aDNA). Basic bioinformatics skills are becoming essential requirements to process and analyze the sheer amounts of data generated by current aDNA studies and in biomedical research in general. This chapter is intended as a practical guide to the assembly of ancient mitochondrial genomes, directly from genomic DNA-derived next-generation sequencing (NGS) data, specifically in the absence of closely related reference genomes. In a hands-on tutorial suitable for readers with little to no prior bioinformatics experience, we reconstruct the mitochondrial genome of a woolly mammoth deposited ~45,000 years ago. We introduce key software tools and outline general strategies for mitogenome assembly, including the critical quality assessment of assembly results without a reference genome.

Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations

Background

Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts?

Results

In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these.

Conclusions

The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1268-0) contains supplementary material, which is available to authorized users.

Complex history of dog (Canis familiaris) origins and translocations in the Pacific revealed by ancient mitogenomes

Archaeological evidence suggests that dogs were introduced to the islands of Oceania via Island Southeast Asia around 3,300 years ago, and reached the eastern islands of Polynesia by the fourteenth century AD. This dispersal is intimately tied to human expansion, but the involvement of dogs in Pacific migrations is not well understood. Our analyses of seven new complete ancient mitogenomes and five partial mtDNA sequences from archaeological dog specimens from Mainland and Island Southeast Asia and the Pacific suggests at least three dog dispersal events into the region, in addition to the introduction of dingoes to Australia. We see an early introduction of dogs to Island Southeast Asia, which does not appear to extend into the islands of Oceania. A shared haplogroup identified between Iron Age Taiwanese dogs, terminal-Lapita and post-Lapita dogs suggests that at least one dog lineage was introduced to Near Oceania by or as the result of interactions with Austronesian language speakers associated with the Lapita Cultural Complex. We did not find any evidence that these dogs were successfully transported beyond New Guinea. Finally, we identify a widespread dog clade found across the Pacific, including the islands of Polynesia, which likely suggests a post-Lapita dog introduction from southern Island Southeast Asia.

A comparison of methods for estimating substitution rates from ancient DNA sequence data

Background

Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species.

Results

Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception.

Conclusions

We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1192-3) contains supplementary material, which is available to authorized users.

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.

Effective Characterisation of the Complete Orang-Utan Mitochondrial DNA Control Region, in the Face of Persistent Focus in Many Taxa on Shorter Hypervariable Regions

The hypervariable region I (HVRI) is persistently used to discern haplotypes, to distinguish geographic subpopulations, and to infer taxonomy in a range of organisms. Numerous studies have highlighted greater heterogeneity elsewhere in the mitochondrial DNA control region, however-particularly, in some species, in other understudied hypervariable regions. To assess the abundance and utility of such potential variations in orang-utans, we characterised 36 complete control-region haplotypes, of which 13 were of Sumatran and 23 of Bornean maternal ancestry, and compared polymorphisms within these and within shorter HVRI segments predominantly analysed in prior phylogenetic studies of Sumatran (~385 bp) and Bornean (~323 bp) orang-utans. We amplified the complete control region in a single PCR that proved successful even with highly degraded, non-invasive samples. By using species-specific primers to produce a single large amplicon (~1600 bp) comprising flanking coding regions, our method also serves to better avoid amplification of nuclear mitochondrial insertions (numts). We found the number, length and position of hypervariable regions is inconsistent between orang-utan species, and that prior definitions of the HVRI were haphazard. Polymorphisms occurring outside the predominantly analysed segments were phylogeographically informative in isolation, and could be used to assign haplotypes to comparable clades concordant with geographic subpopulations. The predominantly analysed segments could discern only up to 76% of all haplotypes, highlighting the forensic utility of complete control-region sequences. In the face of declining sequencing costs and our proven application to poor-quality DNA extracts, we see no reason to ever amplify only specific 'hypervariable regions' in any taxa, particularly as their lengths and positions are inconsistent and cannot be reliably defined-yet this strategy predominates widely. Given their greater utility and consistency, we instead advocate analysis of complete control-region sequences in future studies, where any shorter segment might otherwise have proven the region of choice.

Large-scale mitogenomics enables insights into Schizophora (Diptera) radiation and population diversity

True flies are insects of the order Diptera and encompass one of the most diverse groups of animals on Earth. Within dipterans, Schizophora represents a recent radiation of insects that was used as a model to develop a pipeline for generating complete mitogenomes using various sequencing platforms and strategies. 91 mitogenomes from 32 different species were sequenced and assembled with high fidelity, using amplicon, whole genome shotgun or single molecule sequencing approaches. Based on the novel mitogenomes, we estimate the origin of Schizophora within the Cretaceous-Paleogene (K-Pg) boundary, about 68.3 Ma. Detailed analyses of the blowfly family (Calliphoridae) place its origin at 22 Ma, concomitant with the radiation of grazing mammals. The emergence of ectoparasitism within calliphorids was dated 6.95 Ma for the screwworm fly and 2.3 Ma for the Australian sheep blowfly. Varying population histories were observed for the blowfly Chrysomya megacephala and the housefly Musca domestica samples in our dataset. Whereas blowflies (n = 50) appear to have undergone selective sweeps and/or severe bottlenecks in the New World, houseflies (n = 14) display variation among populations from different zoogeographical zones and low levels of gene flow. The reported high-throughput mitogenomics approach for insects enables new insights into schizophoran diversity and population history of flies.

Natural history collections-based research: progress, promise, and best practices

Specimens and associated data in natural history collections (NHCs) foster substantial scientific progress. In this paper, we explore recent contributions of NHCs to the study of systematics and biogeography, genomics, morphology, stable isotope ecology, and parasites and pathogens of mammals. To begin to assess the magnitude and scope of these contributions, we analyzed publications in the Journal of Mammalogy over the last decade, as well as recent research supported by a single university mammal collection (Museum of Southwestern Biology, Division of Mammals). Using these datasets, we also identify weak links that may be hindering the development of crucial NHC infrastructure. Maintaining the vitality and growth of this foundation of mammalogy depends on broader engagement and support from across the scientific community and is both an ethical and scientific imperative given the rapidly changing environmental conditions on our planet.

Shotgun Mitogenomics Provides a Reference Phylogenetic Framework and Timescale for Living Xenarthrans

Xenarthra (armadillos, sloths, and anteaters) constitutes one of the four major clades of placental mammals. Despite their phylogenetic distinctiveness in mammals, a reference phylogeny is still lacking for the 31 described species. Here we used Illumina shotgun sequencing to assemble 33 new complete mitochondrial genomes, establishing Xenarthra as the first major placental clade to be fully sequenced at the species level for mitogenomes. The resulting data set allowed the reconstruction of a robust phylogenetic framework and timescale that are consistent with previous studies conducted at the genus level using nuclear genes. Incorporating the full species diversity of extant xenarthrans points to a number of inconsistencies in xenarthran systematics and species definition. We propose to split armadillos into two distinct families Dasypodidae (dasypodines) and Chlamyphoridae (euphractines, chlamyphorines, and tolypeutines) to better reflect their ancient divergence, estimated around 42 Ma. Species delimitation within long-nosed armadillos (genus Dasypus) appeared more complex than anticipated, with the discovery of a divergent lineage in French Guiana. Diversification analyses showed Xenarthra to be an ancient clade with a constant diversification rate through time with a species turnover driven by high but constant extinction. We also detected a significant negative correlation between speciation rate and past temperature fluctuations with an increase in speciation rate corresponding to the general cooling observed during the last 15 My. Biogeographic reconstructions identified the tropical rainforest biome of Amazonia and the Guiana Shield as the cradle of xenarthran evolutionary history with subsequent dispersions into more open and dry habitats.

Analysis of Whole Mitogenomes from Ancient Samples

Ancient mitochondrial DNA has been used in a wide variety of paleontological and archeological studies, ranging from population dynamics of extinct species to patterns of domestication. Most of these studies have traditionally been based on the analysis of short fragments from the mitochondrial control region, analyzed using PCR coupled with Sanger sequencing. With the introduction of high-throughput sequencing, as well as new enrichment technologies, the recovery of full mitochondrial genomes (mitogenomes) from ancient specimens has become significantly less complicated. Here we present a protocol to build ancient extracts into Illumina high-throughput sequencing libraries, and subsequent Agilent array-based capture to enrich for the desired mitogenome. Both are based on previously published protocols, with the introduction of several improvements aimed to increase the recovery of short DNA fragments, while keeping the cost and effort requirements low. This protocol was designed for enrichment of mitochondrial DNA in ancient or other degraded samples. However, the protocols can be easily adapted for using for building libraries for shotgun-sequencing of whole genomes, or enrichment of other genomic regions.

The future of ancient DNA: Technical advances and conceptual shifts

Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7× coverage (mammoth) in 2008 to more than 50× coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.

Morphological and genetic evidence for early Holocene cattle management in northeastern China

The domestication of cattle is generally accepted to have taken place in two independent centres: around 10,500 years ago in the Near East, giving rise to modern taurine cattle, and two millennia later in southern Asia, giving rise to zebu cattle. Here we provide firmly dated morphological and genetic evidence for early Holocene management of taurine cattle in northeastern China. We describe conjoining mandibles from this region that show evidence of oral stereotypy, dated to the early Holocene by two independent ¹⁴C dates. Using Illumina high-throughput sequencing coupled with DNA hybridization capture, we characterize 15,406 bp of the mitogenome with on average 16.7-fold coverage. Phylogenetic analyses reveal a hitherto unknown mitochondrial haplogroup that falls outside the known taurine diversity. Our data suggest that the first attempts to manage cattle in northern China predate the introduction of domestic cattle that gave rise to the current stock by several thousand years.

Attempted DNA extraction from a Rancho La Brea Columbian mammoth (Mammuthus columbi): prospects for ancient DNA from asphalt deposits

Fossil-bearing asphalt deposits are an understudied and potentially significant source of ancient DNA. Previous attempts to extract DNA from skeletons preserved at the Rancho La Brea tar pits in Los Angeles, California, have proven unsuccessful, but it is unclear whether this is due to a lack of endogenous DNA, or if the problem is caused by asphalt-mediated inhibition. In an attempt to test these hypotheses, a recently recovered Columbian mammoth (Mammuthus columbi) skeleton with an unusual pattern of asphalt impregnation was studied. Ultimately, none of the bone samples tested successfully amplified M. columbi DNA. Our work suggests that reagents typically used to remove asphalt from ancient samples also inhibit DNA extraction. Ultimately, we conclude that the probability of recovering ancient DNA from fossils in asphalt deposits is strongly (perhaps fatally) hindered by the organic compounds that permeate the bones and that at the Rancho La Brea tar pits, environmental conditions might not have been ideal for the general preservation of genetic material.

Improving evolutionary models for mitochondrial protein data with site-class specific amino acid exchangeability matrices

Adequate modeling of mitochondrial sequence evolution is an essential component of mitochondrial phylogenomics (comparative mitogenomics). There is wide recognition within the field that lineage-specific aspects of mitochondrial evolution should be accommodated through lineage-specific amino-acid exchangeability matrices (e.g., mtMam for mammalian data). However, such a matrix must be applied to all sites and this implies that all sites are subject to the same, or largely similar, evolutionary constraints. This assumption is unjustified. Indeed, substantial differences are expected to arise from three-dimensional structures that impose different physiochemical environments on individual amino acid residues. The objectives of this paper are (1) to investigate the extent to which amino acid evolution varies among sites of mitochondrial proteins, and (2) to assess the potential benefits of explicitly modeling such variability. To achieve this, we developed a novel method for partitioning sites based on amino acid physiochemical properties. We apply this method to two datasets derived from complete mitochondrial genomes of mammals and fish, and use maximum likelihood to estimate amino acid exchangeabilities for the different groups of sites. Using this approach we identified large groups of sites evolving under unique physiochemical constraints. Estimates of amino acid exchangeabilities differed significantly among such groups. Moreover, we found that joint estimates of amino acid exchangeabilities do not adequately represent the natural variability in evolutionary processes among sites of mitochondrial proteins. Significant improvements in likelihood are obtained when the new matrices are employed. We also find that maximum likelihood estimates of branch lengths can be strongly impacted. We provide sets of matrices suitable for groups of sites subject to similar physiochemical constraints, and discuss how they might be used to analyze real data. We also discuss how the general approach might be employed to improve a variety of mitogenomic-based research activities.