Complete mitochondrial genome and phylogeny of Pleistocene mammoth Mammuthus primigenius

Do elephants really never forget? What we know about elephant memory and a call for further investigation

Despite popular culture's promotion of the elephant's ability to "never forget," there is remarkably limited empirical research on the memory capacities of any living elephant species (Asian, Elephas maximus; African savanna, Loxodonta africana; African forest, Loxodonta cyclotis). A growing body of literature on elephant cognition and behavioral ecology has provided insight into the elephant's ability to behave flexibly in changing physical and social environments, but little direct evidence of how memory might relate to this flexibility exists. In this paper, we review and discuss the potential relationships between what we know about elephant cognition and behavior and the elephants' memory for the world around them as they navigate their physical, social, and spatial environments. We also discuss future directions for investigating elephant memory and implications for such research on elephant conservation and human-elephant conflict mitigation.

Phylogeography of the woolly mammoth (Mammuthus primigenius) in the Minusinsk Depression of southern Siberia in the Late Pleistocene

To date, a number of studies have been published on the phylogenetics of woolly mammoths (Mammuthus primigenius), ranging from analyses of parts of the mitochondrial genome to studies of complete nuclear genomes. However, until recently nothing was known about the genetic diversity of woolly mammoths in southern Siberia, in the Minusinsk Depression in particular. Within the framework of this effort, libraries for high-throughput sequencing of seven bone samples of woolly mammoths were obtained, two-round enrichment using biotinylated probes of modern mtDNA of Elephas maximus immobilised on magnetic microspheres and sequencing with subsequent bioinformatic analysis were carried out. Phylogenetic reconstructions showed the presence of all studied mammoths in clade I, which expanded its range. The assignment of mammoth mitotypes in the Minusinsk Depression to different clusters within clade I may indicate a sufficiently high diversity of their gene pool. Phylogeographic reconstructions revealed a genetic proximity of mitochondrial lineages of Late Pleistocene mammoths of the Minusinsk Depression and other regions of eastern Siberia and estimated their divergence time in the range of 100-150 thousand years ago, which indicates active migrations of woolly mammoths over vast territories of eastern Siberia in the late Middle Pleistoceneearly Late Pleistocene.

Population genetic structure and historical demography of the population of forest elephants in Côte d'Ivoire

The population of forest elephant (Loxodonta cyclotis) has continuously declined in Côte d'Ivoire and, the remaining population largely consists of subpopulations that are fragmented and isolated. No data actually exist on the level of genetic diversity and population genetic structure of current forest elephant populations in Côte d'Ivoire. In this sense, determining genetic diversity and the underlying mechanisms of population differentiation is crucial for the initiation of effective conservation management. A total of 158 dung samples of forest elephants were collected at stage 1 of decompositions (dung pile intact, very fresh) in three Classified Forests (CF) (Bossématié, Dassioko and Port-Gauthier) in Côte d'Ivoire. A total of 101 sequences of the mitochondrial DNA control region measuring 600 base pair and 26 haplotypes were obtained. A haplotypic diversity ranging from 0.655 ± 0.050 at Bossématié and 0.859 ± 0.088 at Port Gauthier was obtained. Fifteen (15) out of 26 haplotypes observed were singletons and only the Dassioko and Port Gauthier CFs shared the same haplotypes. The strong genetic connectivity between forest elephant populations of the Dassioko and Port Gauthier CFs is supported by the grouping of these populations into a single cluster by Bayesian analysis. Although populations of L. cyclotis exhibit relatively high genetic diversity, habitat fragmentation could affect the genetic variability of current populations. Urgent measures including the reinforcement/establishment of genetic corridors and the strengthening of protection measures need to be undertaken to save the remaining populations of forest elephants in Côte d'Ivoire.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

Genomic resources for Asian (Elephas maximus) and African savannah elephant (Loxodonta africana) conservation and health research

We provide novel genomic resources to help understand the genomic traits involved in elephant health and to aid conservation efforts. We sequence 11 elephant genomes (5 African savannah, 6 Asian) from North American zoos, including 9 de novo assemblies. We estimate elephant germline mutation rates and reconstruct demographic histories. Finally, we provide an in-solution capture assay to genotype Asian elephants. This assay is suitable for analyzing degraded museum and noninvasive samples, such as feces and hair. The elephant genomic resources we present here should allow for more detailed and uniform studies in the future to aid elephant conservation efforts and disease research.

A unique single nucleotide polymorphism in Agouti Signalling Protein (ASIP) gene changes coat colour of Sri Lankan leopard (Panthera pardus kotiya) to dark black

The Sri Lankan leopard (Panthera pardus kotiya) is an endangered subspecies restricted to isolated and fragmented populations in Sri Lanka. Among them, melanistic leopards have been recorded on a few occasions. Literature suggests the evolution of melanism several times in the Felidae family, with three species having distinct mutations. Nevertheless, the mutations or other variations in the remaining species, including Sri Lankan melanistic leopard, are unknown. We used reference-based assembled nuclear genomes of Sri Lankan wild type and melanistic leopards and de novo assembled mitogenomes of the same to investigate the genetic basis, adaptive significance, and evolutionary history of the Sri Lankan melanistic leopard. Interestingly, we identified a single nucleotide polymorphism in exon-4 Sri Lankan melanistic leopard, which may completely ablate Agouti Signalling Protein (ASIP) function. The wild type leopards in Sri Lanka did not carry this mutation, suggesting the cause for the occurrence of melanistic leopords in the population. Comparative analysis of existing genomic data in the literature suggests it as a P. p. kotiya specific mutation and a novel mutation in the ASIP-gene of the Felidae family, contributing to naturally occurring colour polymorphism. Our data suggested the coalescence time of Sri Lankan leopards at ~0.5 million years, sisters to the Panthera pardus lineage. The genetic diversity was low in Sri Lankan leopards. Further, the P. p. kotiya melanistic leopard is a different morphotype of the P. p. kotiya wildtype leopard resulting from the mutation in the ASIP-gene. The ability of black leopards to camouflage, along with the likelihood of recurrence and transfer to future generations, suggests that this rare mutation could be environment-adaptable.

Genetic features of Sri Lankan elephant, Elephas maximus maximus Linnaeus revealed by high throughput sequencing of mitogenome and ddRAD-seq

Elephas maximus maximus Linnaeus, the Sri Lankan subspecies is the largest and the darkest among Asian elephants. Patches of depigmented areas with no skin color on the ears, face, trunk, and belly morphologically differentiate it from the others. The elephant population in Sri Lanka is now limited to smaller areas and protected under Sri Lankan law. Despite its ecological and evolutionary importance, the relationship between Sri Lankan elephants and their phylogenetic position among Asian elephants remains controversial. While identifying genetic diversity is the key to any conservation and management strategies, limited data is currently available. To address such issues, we analyzed 24 elephants with known parental lineages with high throughput ddRAD-seq. The mitogenome suggested the coalescence time of the Sri Lankan elephant at ~0.2 million years, and sister to Myanmar elephants supporting the hypothesis of the movement of elephants in Eurasia. The ddRAD-seq approach identified 50,490 genome-wide SNPs among Sri Lankan elephants. The genetic diversity within Sri Lankan elephants assessed with identified SNPs suggests a geographical differentiation resulting in three main clusters; north-eastern, mid-latitude, and southern regions. Interestingly, though it was believed that elephants from the Sinharaja rainforest are of an isolated population, the ddRAD-based genetic analysis clustered it with the north-eastern elephants. The effect of habitat fragmentation on genetic diversity could be further assessed with more samples with specific SNPs identified in the current study.

MEANGS: an efficient seed-free tool for de novo assembling animal mitochondrial genome using whole genome NGS data

Advances in next-generation sequencing (NGS) technologies have led to an exponential increase in the number of whole genome sequences (WGS) in databases. This wealth of WGS data has greatly facilitated the recovery of full mitochondrial genomes (mitogenomes), which are vital for phylogenetic, evolutionary and ecological studies. Unfortunately, most existing software cannot easily assemble mitogenome reference sequences conveniently or efficiently. Therefore, we developed a seed-free de novo assembly tool, MEANGS, which applies the trie-search method to extend contigs from self-discovery seeds and assemble a mitogenome from animal WGS data. We then used data from 16 species with different qualities to compare the performance of MEANGS with three other available programs. MEANGS exhibited the best overall performance since it was the only one that completed all tests, and it assembled full or partial mitogenomes for all of the tested samples while the others failed. Furthermore, MEANGS selects superior assembly sequences and annotates protein-coding genes. Thus, MEANGS can be one of the most efficient software for generating high-quality mitogenomes so far, the further use of it will benefit the study on mitogenome based on whole genome NGS data. MEANGS is available at https://github.com/YanCCscu/meangs.

Estimating the dwarfing rate of an extinct Sicilian elephant

Evolution on islands, together with the often extreme phenotypic changes associated with it, has attracted much interest from evolutionary biologists. However, measuring the rate of change of phenotypic traits of extinct animals can be challenging, in part due to the incompleteness of the fossil record. Here, we use combined molecular and fossil evidence to define the minimum and maximum rate of dwarfing in an extinct Mediterranean dwarf elephant from Puntali Cave (Sicily).¹ Despite the challenges associated with recovering ancient DNA from warm climates,² we successfully retrieved a mitogenome from a sample with an estimated age between 175,500 and 50,000 years. Our results suggest that this specific Sicilian elephant lineage evolved from one of the largest terrestrial mammals that ever lived³ to an island species weighing less than 20% of its original mass with an estimated mass reduction between 0.74 and 200.95 kg and height reduction between 0.15 and 41.49 mm per generation. We show that combining ancient DNA with paleontological and geochronological evidence can constrain the timing of phenotypic changes with greater accuracy than could be achieved using any source of evidence in isolation.

Competitive mapping allows for the identification and exclusion of human DNA contamination in ancient faunal genomic datasets

BACKGROUND

After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data.

RESULTS

Here we investigate whether human contaminating DNA can be found in ancient faunal sequencing datasets. We identify variable levels of human contamination, which persists even after the sequence reads have been mapped to the faunal reference genomes. This contamination has the potential to affect a range of downstream analyses.

CONCLUSIONS

We propose a fast and simple method, based on competitive mapping, which allows identifying and removing human contamination from ancient faunal DNA datasets with limited losses of true ancient data. This method could represent an important tool for the ancient DNA field.

Distinguishing extant elephants ivory from mammoth ivory using a short sequence of cytochrome b gene

Trade in ivory from extant elephant species namely Asian elephant (Elephas maximus), African savanna elephant (Loxodonta africana) and African forest elephant (Loxodonta cyclotis) is regulated internationally, while the trade in ivory from extinct species of Elephantidae, including woolly mammoth, is unregulated. This distinction creates opportunity for laundering and trading elephant ivory as mammoth ivory. The existing morphological and molecular genetics methods do not reliably distinguish the source of ivory items that lack clear identification characteristics or for which the quality of extracted DNA cannot support amplification of large gene fragments. We present a PCR-sequencing method based on 116 bp target sequence of the cytochrome b gene to specifically amplify elephantid DNA while simultaneously excluding non-elephantid species and ivory substitutes, and while avoiding contamination by human DNA. The partial Cytochrome b gene sequence enabled accurate association of ivory samples with their species of origin for all three extant elephants and from mammoth. The detection limit of the PCR system was as low as 10 copy numbers of target DNA. The amplification and sequencing success reached 96.7% for woolly mammoth ivory and 100% for African savanna elephant and African forest elephant ivory. This is the first validated method for distinguishing elephant from mammoth ivory and it provides forensic support for investigation of ivory laundering cases.

Signals of positive selection in mitochondrial protein-coding genes of woolly mammoth: Adaptation to extreme environments?

The mammoths originated in warm and equatorial Africa and later colonized cold and high-latitude environments. Studies on nuclear genes suggest that woolly mammoth had evolved genetic variations involved in processes relevant to cold tolerance, including lipid metabolism and thermogenesis, and adaptation to extremely varied light and darkness cycles. The mitochondria is a major regulator of cellular energy metabolism, thus the mitogenome of mammoths may also exhibit adaptive evolution. However, little is yet known in this regard. In this study, we analyzed mitochondrial protein-coding genes (MPCGs) sequences of 75 broadly distributed woolly mammoths (Mammuthus primigenius) to test for signatures of positive selection. Results showed that a total of eleven amino acid sites in six genes, namely ND1, ND4, ND5, ND6, CYTB, and ATP6, displayed strong evidence of positive selection. Two sites were located in close proximity to proton-translocation channels in mitochondrial complex I. Biochemical and homology protein structure modeling analyses demonstrated that five amino acid substitutions in ND1, ND5, and ND6 might have influenced the performance of protein-protein interaction among subunits of complex I, and three substitutions in CYTB and ATP6 might have influenced the performance of metabolic regulatory chain. These findings suggest metabolic adaptations in the mitogenome of woolly mammoths in relation to extreme environments and provide a basis for further tests on the significance of the variations on other systems.

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.

Complete mitochondrial genome of a woolly mammoth (Mammuthus primigenius) from Maly Lyakhovsky Island (New Siberian Islands, Russia) and its phylogenetic assessment

We present a complete sequence and an annotation of the mitochondrial genome of the woolly mammoth (Mammuthus primigenius) found in 2012 on Maly Lyakhovsky Island (North-Eastern Siberia, Russia). The genome was 16,851 bp long and contained 13 protein-coding, 22 tRNA, and 2 rRNA genes. It was AT reach (61.3%) with A = 32.9%, T = 28.4%, C = 25.3%, and G = 13.4%.

Body size downgrading of mammals over the late Quaternary

Since the late Pleistocene, large-bodied mammals have been extirpated from much of Earth. Although all habitable continents once harbored giant mammals, the few remaining species are largely confined to Africa. This decline is coincident with the global expansion of hominins over the late Quaternary. Here, we quantify mammalian extinction selectivity, continental body size distributions, and taxonomic diversity over five time periods spanning the past 125,000 years and stretching approximately 200 years into the future. We demonstrate that size-selective extinction was already under way in the oldest interval and occurred on all continents, within all trophic modes, and across all time intervals. Moreover, the degree of selectivity was unprecedented in 65 million years of mammalian evolution. The distinctive selectivity signature implicates hominin activity as a primary driver of taxonomic losses and ecosystem homogenization. Because megafauna have a disproportionate influence on ecosystem structure and function, past and present body size downgrading is reshaping Earth's biosphere.

Central European Woolly Mammoth Population Dynamics: Insights from Late Pleistocene Mitochondrial Genomes

The population dynamics of the Pleistocene woolly mammoth (Mammuthus primigenius) has been the subject of intensive palaeogenetic research. Although a large number of mitochondrial genomes across Eurasia have been reconstructed, the available data remains geographically sparse and mostly focused on eastern Eurasia. Thus, population dynamics in other regions have not been extensively investigated. Here, we use a multi-method approach utilising proteomic, stable isotope and genetic techniques to identify and generate twenty woolly mammoth mitochondrial genomes, and associated dietary stable isotopic data, from highly fragmentary Late Pleistocene material from central Europe. We begin to address region-specific questions regarding central European woolly mammoth populations, highlighting parallels with a previous replacement event in eastern Eurasia ten thousand years earlier. A high number of shared derived mutations between woolly mammoth mitochondrial clades are identified, questioning previous phylogenetic analysis and thus emphasizing the need for nuclear DNA studies to explicate the increasingly complex genetic history of the woolly mammoth.

Mitogenome evolution in the last surviving woolly mammoth population reveals neutral and functional consequences of small population size

The onset of the Holocene was associated with a global temperature increase, which led to a rise in sea levels and isolation of the last surviving population of woolly mammoths on Wrangel Island. Understanding what happened with the population's genetic diversity at the time of the isolation and during the ensuing 6000 years can help clarify the effects of bottlenecks and subsequent limited population sizes in species approaching extinction. Previous genetic studies have highlighted questions about how the Holocene Wrangel population was established and how the isolation event affected genetic diversity. Here, we generated high‐quality mitogenomes from 21 radiocarbon‐dated woolly mammoths to compare the ancestral large and genetically diverse Late Pleistocene Siberian population and the small Holocene Wrangel population. Our results indicate that mitogenome diversity was reduced to one single haplotype at the time of the isolation, and thus that the Holocene Wrangel Island population was established by a single maternal lineage. Moreover, we show that the ensuing small effective population size coincided with fixation of a nonsynonymous mutation, and a comparative analysis of mutation rates suggests that the evolutionary rate was accelerated in the Holocene population. These results suggest that isolation on Wrangel Island led to an increase in the frequency of deleterious genetic variation, and thus are consistent with the hypothesis that strong genetic drift in small populations leads to purifying selection being less effective in removing deleterious mutations.

Mammoth and Elephant Phylogenetic Relationships: Mammut Americanum, the Missing Outgroup

At the morphological level, the woolly mammoth has most often been considered as the sister-species of Asian elephants, but at the DNA level, different studies have found support for proximity with African elephants. Recent reports have increased the available sequence data and apparently solved the discrepancy, finding mammoths to be most closely related to Asian elephants. However, we demonstrate here that the three competing topologies have similar likelihood, bayesian and parsimony supports. The analysis further suggests the inadequacy of using Sirenia or Hyracoidea as outgroups. We therefore argue that orthologous sequences from the extinct American mastodon will be required to definitively solve this long-standing question.

The evolutionary and phylogeographic history of woolly mammoths: a comprehensive mitogenomic analysis

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

In-solution hybridization for mammalian mitogenome enrichment: pros, cons and challenges associated with multiplexing degraded DNA

Here, we present a set of RNA-based probes for whole mitochondrial genome in-solution enrichment, targeting a diversity of mammalian mitogenomes. This probes set was designed from seven mammalian orders and tested to determine the utility for enriching degraded DNA. We generated 63 mitogenomes representing five orders and 22 genera of mammals that yielded varying coverage ranging from 0 to >5400X. Based on a threshold of 70% mitogenome recovery and at least 10× average coverage, 32 individuals or 51% of samples were considered successful. The estimated sequence divergence of samples from the probe sequences used to construct the array ranged up to nearly 20%. Sample type was more predictive of mitogenome recovery than sample age. The proportion of reads from each individual in multiplexed enrichments was highly skewed, with each pool having one sample that yielded a majority of the reads. Recovery across each mitochondrial gene varied with most samples exhibiting regions with gaps or ambiguous sites. We estimated the ability of the probes to capture mitogenomes from a diversity of mammalian taxa not included here by performing a clustering analysis of published sequences for 100 taxa representing most mammalian orders. Our study demonstrates that a general array can be cost and time effective when there is a need to screen a modest number of individuals from a variety of taxa. We also address the practical concerns for using such a tool, with regard to pooling samples, generating high quality mitogenomes and detail a pipeline to remove chimeric molecules.

Pros and cons of methylation-based enrichment methods for ancient DNA

The recent discovery that DNA methylation survives in fossil material provides an opportunity for novel molecular approaches in palaeogenomics. Here, we apply to ancient DNA extracts the probe-independent Methylated Binding Domains (MBD)-based enrichment method, which targets DNA molecules containing methylated CpGs. Using remains of a Palaeo-Eskimo Saqqaq individual, woolly mammoths, polar bears and two equine species, we confirm that DNA methylation survives in a variety of tissues, environmental contexts and over a large temporal range (4,000 to over 45,000 years before present). MBD enrichment, however, appears principally biased towards the recovery of CpG-rich and long DNA templates and is limited by the fast post-mortem cytosine deamination rates of methylated epialleles. This method, thus, appears only appropriate for the analysis of ancient methylomes from very well preserved samples, where both DNA fragmentation and deamination have been limited. This work represents an essential step toward the characterization of ancient methylation signatures, which will help understanding the role of epigenetic changes in past environmental and cultural transitions.

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.

Primate phylogenetic relationships and divergence dates inferred from complete mitochondrial genomes

The origins and the divergence times of the most basal lineages within primates have been difficult to resolve mainly due to the incomplete sampling of early fossil taxa. The main source of contention is related to the discordance between molecular and fossil estimates: while there are no crown primate fossils older than 56Ma, most molecule-based estimates extend the origins of crown primates into the Cretaceous. Here we present a comprehensive mitogenomic study of primates. We assembled 87 mammalian mitochondrial genomes, including 62 primate species representing all the families of the order. We newly sequenced eleven mitochondrial genomes, including eight Old World monkeys and three strepsirrhines. Phylogenetic analyses support a strong topology, confirming the monophyly for all the major primate clades. In contrast to previous mitogenomic studies, the positions of tarsiers and colugos relative to strepsirrhines and anthropoids are well resolved. In order to improve our understanding of how fossil calibrations affect age estimates within primates, we explore the effect of seventeen fossil calibrations across primates and other mammalian groups and we select a subset of calibrations to date our mitogenomic tree. The divergence date estimates of the Strepsirrhine/Haplorhine split support an origin of crown primates in the Late Cretaceous, at around 74Ma. This result supports a short-fuse model of primate origins, whereby relatively little time passed between the origin of the order and the diversification of its major clades. It also suggests that the early primate fossil record is likely poorly sampled.

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.

Holarctic genetic structure and range dynamics in the woolly mammoth

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32–34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20–15 kyr BP, when a severe population size decline occurred.

Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads--a baiting and iterative mapping approach

We present an in silico approach for the reconstruction of complete mitochondrial genomes of non-model organisms directly from next-generation sequencing (NGS) data—mitochondrial baiting and iterative mapping (MITObim). The method is straightforward even if only (i) distantly related mitochondrial genomes or (ii) mitochondrial barcode sequences are available as starting-reference sequences or seeds, respectively. We demonstrate the efficiency of the approach in case studies using real NGS data sets of the two monogenean ectoparasites species Gyrodactylus thymalli and Gyrodactylus derjavinoides including their respective teleost hosts European grayling (Thymallus thymallus) and Rainbow trout (Oncorhynchus mykiss). MITObim appeared superior to existing tools in terms of accuracy, runtime and memory requirements and fully automatically recovered mitochondrial genomes exceeding 99.5% accuracy from total genomic DNA derived NGS data sets in <24 h using a standard desktop computer. The approach overcomes the limitations of traditional strategies for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information at hand and represents a fast and highly efficient in silico alternative to laborious conventional strategies relying on initial long-range PCR. We furthermore demonstrate the applicability of MITObim for metagenomic/pooled data sets using simulated data. MITObim is an easy to use tool even for biologists with modest bioinformatics experience. The software is made available as open source pipeline under the MIT license at https://github.com/chrishah/MITObim.

Mitochondrial genome sequences of Artemia tibetiana and Artemia urmiana: assessing molecular changes for high plateau adaptation

Brine shrimps, Artemia (Crustacea, Anostraca), inhabit hypersaline environments and have a broad geographical distribution from sea level to high plateaus. Artemia therefore possess significant genetic diversity, which gives them their outstanding adaptability. To understand this remarkable plasticity, we sequenced the mitochondrial genomes of two Artemia tibetiana isolates from the Tibetan Plateau in China and one Artemia urmiana isolate from Lake Urmia in Iran and compared them with the genome of a low-altitude Artemia, A. franciscana. We compared the ratio of the rate of nonsynonymous (Ka) and synonymous (Ks) substitutions (Ka/Ks ratio) in the mitochondrial protein-coding gene sequences and found that atp8 had the highest Ka/Ks ratios in comparisons of A. franciscana with either A. tibetiana or A. urmiana and that atp6 had the highest Ka/Ks ratio between A. tibetiana and A. urmiana. Atp6 may have experienced strong selective pressure for high-altitude adaptation because although A. tibetiana and A. urmiana are closely related they live at different altitudes. We identified two extended termination-associated sequences and three conserved sequence blocks in the D-loop region of the mitochondrial genomes. We propose that sequence variations in the D-loop region and in the subunits of the respiratory chain complexes independently or collectively contribute to the adaptation of Artemia to different altitudes.

Phylogenetic analysis of ancient DNA using BEAST

Under exceptional circumstances, it is possible to obtain DNA sequences from samples that are up to hundreds of thousands of years old. These data provide an opportunity to look directly at past genetic diversity, to trace the evolutionary process through time, and to infer demographic and phylogeographic trends. Ancient DNA (aDNA) data sets have some degree of intrinsic temporal structure because the sequences have been obtained from samples of different ages. When analyzing these data sets, it is usually necessary to take the sampling times into account. A number of phylogenetic methods have been designed with this purpose in mind. Here I describe the steps involved in Bayesian phylogenetic analysis of aDNA data. I outline a procedure that can be used to co-estimate the genealogical relationships, mutation rate, evolutionary timescale, and demographic history of the study species in a single analytical framework. A number of modifications to the methodology can be made in order to deal with complicating factors such as postmortem damage, sequences from undated samples, and data sets with low information content.

Forest elephant mitochondrial genomes reveal that elephantid diversification in Africa tracked climate transitions

Among elephants, the phylogeographic patterns of mitochondrial (mt) and nuclear markers are often incongruent. One hypothesis attributes this to sex differences in dispersal and in the variance of reproductive success. We tested this hypothesis by examining the coalescent dates of genetic markers within elephantid lineages, predicting that lower dispersal and lower variance in reproductive success among females would have increased mtDNA relative to nuclear coalescent dates. We sequenced the mitochondrial genomes of two forest elephants, aligning them to mitogenomes of African savanna and Asian elephants, and of woolly mammoths, including the most divergent mitogenomes within each lineage. Using fossil calibrations, the divergence between African elephant F and S clade mitochondrial genomes (originating in forest and savanna elephant lineages, respectively) was estimated as 5.5 Ma. We estimated that the (African) ancestor of the mammoth and Asian elephant lineages diverged 6.0 Ma, indicating that four elephantid lineages had differentiated in Africa by the Miocene-Pliocene transition, concurrent with drier climates. The coalescent date for forest elephant mtDNAs was c. 2.4 Ma, suggesting that the decrease in tropical forest cover during the Pleistocene isolated distinct African forest elephant lineages. For all elephantid lineages, the ratio of mtDNA to nuclear coalescent dates was much greater than 0.25. This is consistent with the expectation that sex differences in dispersal and in variance of reproductive success would have increased the effective population size of mtDNA relative to nuclear markers in elephantids, contributing to the persistence of incongruent mtDNA phylogeographic patterns.

Complete Columbian mammoth mitogenome suggests interbreeding with woolly mammoths

BACKGROUND

Late Pleistocene North America hosted at least two divergent and ecologically distinct species of mammoth: the periglacial woolly mammoth (Mammuthus primigenius) and the subglacial Columbian mammoth (Mammuthus columbi). To date, mammoth genetic research has been entirely restricted to woolly mammoths, rendering their genetic evolution difficult to contextualize within broader Pleistocene paleoecology and biogeography. Here, we take an interspecific approach to clarifying mammoth phylogeny by targeting Columbian mammoth remains for mitogenomic sequencing.

RESULTS

We sequenced the first complete mitochondrial genome of a classic Columbian mammoth, as well as the first complete mitochondrial genome of a North American woolly mammoth. Somewhat contrary to conventional paleontological models, which posit that the two species were highly divergent, the M. columbi mitogenome we obtained falls securely within a subclade of endemic North American M. primigenius.

CONCLUSIONS

Though limited, our data suggest that the two species interbred at some point in their evolutionary histories. One potential explanation is that woolly mammoth haplotypes entered Columbian mammoth populations via introgression at subglacial ecotones, a scenario with compelling parallels in extant elephants and consistent with certain regional paleontological observations. This highlights the need for multi-genomic data to sufficiently characterize mammoth evolutionary history. Our results demonstrate that the use of next-generation sequencing technologies holds promise in obtaining such data, even from non-cave, non-permafrost Pleistocene depositional contexts.

The complete mitochondrial genome of an 11,450-year-old aurochsen (Bos primigenius) from Central Italy

Background

Bos primigenius, the aurochs, is the wild ancestor of modern cattle breeds and was formerly widespread across Eurasia and northern Africa. After a progressive decline, the species became extinct in 1627. The origin of modern taurine breeds in Europe is debated. Archaeological and early genetic evidence point to a single Near Eastern origin and a subsequent spread during the diffusion of herding and farming. More recent genetic data are instead compatible with local domestication events or at least some level of local introgression from the aurochs. Here we present the analysis of the complete mitochondrial genome of a pre-Neolithic Italian aurochs.

Results

In this study, we applied a combined strategy employing both multiplex PCR amplifications and 454 pyrosequencing technology to sequence the complete mitochondrial genome of an 11,450-year-old aurochs specimen from Central Italy. Phylogenetic analysis of the aurochs mtDNA genome supports the conclusions from previous studies of short mtDNA fragments - namely that Italian aurochsen were genetically very similar to modern cattle breeds, but highly divergent from the North-Central European aurochsen.

Conclusions

Complete mitochondrial genome sequences are now available for several modern cattle and two pre-Neolithic mtDNA genomes from very different geographic areas. These data suggest that previously identified sub-groups within the widespread modern cattle mitochondrial T clade are polyphyletic, and they support the hypothesis that modern European breeds have multiple geographic origins.

Genomic DNA sequences from mastodon and woolly mammoth reveal deep speciation of forest and savanna elephants

To elucidate the history of living and extinct elephantids, we generated 39,763 bp of aligned nuclear DNA sequence across 375 loci for African savanna elephant, African forest elephant, Asian elephant, the extinct American mastodon, and the woolly mammoth. Our data establish that the Asian elephant is the closest living relative of the extinct mammoth in the nuclear genome, extending previous findings from mitochondrial DNA analyses. We also find that savanna and forest elephants, which some have argued are the same species, are as or more divergent in the nuclear genome as mammoths and Asian elephants, which are considered to be distinct genera, thus resolving a long-standing debate about the appropriate taxonomic classification of the African elephants. Finally, we document a much larger effective population size in forest elephants compared with the other elephantid taxa, likely reflecting species differences in ancient geographic structure and range and differences in life history traits such as variance in male reproductive success.