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

A Million Years of Mammoth Mitogenome Evolution

The genomic study of specimens dating to the Early and Middle Pleistocene (EP and MP), a period spanning from 2.6 million years ago (Ma) to 126 thousand years ago (ka), has the potential to elucidate the evolutionary processes that shaped present-day biodiversity. Obtaining genomic data from this period is challenging, but mitochondrial DNA, given its higher abundance compared to nuclear DNA, could play an important role to understand evolutionary processes at this time scale. In this study, we report 34 new mitogenomes, including two EP and nine MP mammoth (Mammuthus spp.) specimens from Siberia and North America and analyze them jointly with >200 publicly available mitogenomes to reconstruct a transect of mammoth mitogenome diversity throughout the last million years. We find that our EP mitogenomes fall outside the diversity of all Late Pleistocene (LP) mammoths, while those derived from MP mammoths are basal to LP mammoth Clades 2 and 3, supporting an ancient Siberian origin of these lineages. In contrast, the geographical origin of Clade 1 remains unresolved. With these new deep-time mitogenomes, we observe diversification events across all clades that appear consistent with previously hypothesized MP and LP demographic changes. Furthermore, we improve upon an existing methodology for molecular clock dating of specimens >50 ka, demonstrating that specimens need to be individually dated to avoid biases in their age estimates. Both the molecular and analytical improvements presented here highlight the importance of deep-time genomic data to discover long-lost genetic diversity, enabling better assessments of evolutionary histories.

Bayesian Inference Under the Multispecies Coalescent with Ancient DNA Sequences

Ancient DNA (aDNA) is increasingly being used to investigate questions such as the phylogenetic relationships and divergence times of extant and extinct species. If aDNA samples are sufficiently old, expected branch lengths (in units of nucleotide substitutions) are reduced relative to contemporary samples. This can be accounted for by incorporating sample ages into phylogenetic analyses. Existing methods that use tip (sample) dates infer gene trees rather than species trees, which can lead to incorrect or biased inferences of the species tree. Methods using a multispecies coalescent (MSC) model overcome these issues. We developed an MSC model with tip dates and implemented it in the program BPP. The method performed well for a range of biologically realistic scenarios, estimating calibrated divergence times and mutation rates precisely. Simulations suggest that estimation precision can be best improved by prioritizing sampling of many loci and more ancient samples. Incorrectly treating ancient samples as contemporary in analyzing simulated data, mimicking a common practice of empirical analyses, led to large systematic biases in model parameters, including divergence times. Two genomic datasets of mammoths and elephants were analyzed, demonstrating the method's empirical utility.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.¹ In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.² Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia³-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.^4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,^6-12 in paleodietary reconstructions,^13-19 and facilitated genomic reconstructions from cave contexts.^8,20-22 Murchie et al.^6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.²⁴ Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 30,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Pleistocene allopatric differentiation followed by recent range expansion explains the distribution and molecular diversity of two congeneric crustacean species in the Palaearctic

Pleistocene glaciations had a tremendous impact on the biota across the Palaearctic, resulting in strong phylogeographic signals of range contraction and rapid postglacial recolonization of the deglaciated areas. Here, we explore the diversity patterns and history of two sibling species of passively dispersing taxa typical of temporary ponds, fairy shrimps (Anostraca). We combine mitochondrial (COI) and nuclear (ITS2 and 18S) markers to conduct a range-wide phylogeographic study including 56 populations of Branchinecta ferox and Branchinecta orientalis in the Palaearctic. Specifically, we investigate whether their largely overlapping ranges in Europe resulted from allopatric differentiation in separate glacial refugia followed by a secondary contact and reconstruct their postglacial recolonization from the inhabited refugia. Our results suggest the existence of distinct refugia for the two species, with genetic divergence among intraspecific lineages consistent with late Pleistocene glacial cycles. While B. ferox lineages originated from Mediterranean refugia, the origin of B. orientalis lineages was possibly located on the Pannonian Plain. We showed that most dispersal events predominantly happened within 100 km, coupled with several recent long-distance events (> 1000 km). Hence the regional habitat density of suitable habitats in Central Europe is possibly a key to the co-existence of the two species. Overall, our study illustrates how isolation in combination with stochastic effects linked to glacial periods are important drivers of the allopatric differentiation of Palaearctic taxa.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

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.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge

The Bering Land Bridge (BLB) last connected Eurasia and North America during the Late Pleistocene. Although the BLB would have enabled transfers of terrestrial biota in both directions, it also acted as an ecological filter whose permeability varied considerably over time. Here we explore the possible impacts of this ecological corridor on genetic diversity within, and connectivity among, populations of a once wide-ranging group, the caballine horses (Equus spp.). Using a panel of 187 mitochondrial and eight nuclear genomes recovered from present-day and extinct caballine horses sampled across the Holarctic, we found that Eurasian horse populations initially diverged from those in North America, their ancestral continent, around 1.0-0.8 million years ago. Subsequent to this split our mitochondrial DNA analysis identified two bidirectional long-range dispersals across the BLB ~875-625 and ~200-50 thousand years ago, during the Middle and Late Pleistocene. Whole genome analysis indicated low levels of gene flow between North American and Eurasian horse populations, which probably occurred as a result of these inferred dispersals. Nonetheless, mitochondrial and nuclear diversity of caballine horse populations retained strong phylogeographical structuring. Our results suggest that barriers to gene flow, currently unidentified but possibly related to habitat distribution across Beringia or ongoing evolutionary divergence, played an important role in shaping the early genetic history of caballine horses, including the ancestors of living horses within Equus ferus.

Global phylogeography and invasion history of the spotted lanternfly revealed by mitochondrial phylogenomics

Biological invasion has been a serious global threat due to increasing international trade and population movements. Tracking the source and route of invasive species and evaluating the genetic differences in their native regions have great significance for the effective monitoring and management, and further resolving the invasive mechanism. The spotted lanternfly Lycorma delicatula is native to China and invaded South Korea, Japan, and the United States during the last decade, causing severe damages to the fruits and timber industries. However, its global phylogeographic pattern and invasion history are not clearly understood. We applied high-throughput sequencing to obtain 392 whole mitochondrial genome sequences from four countries to ascertain the origin, dispersal, and invasion history of the spotted lanternfly. Phylogenomic analyses revealed that the spotted lanternfly originated from southwestern China, diverged into six phylogeographic lineages, and experienced northward expansion across the Yangtze River in the late Pleistocene. South Korea populations were derived from multiple invasions from eastern China and Japan with two different genetic sources of northwestern (Loess Plateau) and eastern (East Plain) lineages in China, whereas the each of Japan and the United States had only one. The United States populations originated through single invasive event from South Korea, which served as a bridgehead of invasion. The environmental conditions, especially the distribution of host Ailanthus trees, and adaptability possibly account for the rapid spread of the spotted lanternfly in the native and introduced regions.

Molecular sexing of degraded DNA from elephants and mammoths: a genotyping assay relevant both to conservation biology and to paleogenetics

It is important to determine the sex of elephants from their samples-faeces from the field or seized ivory-for forensic reasons or to understand population demography and genetic structure. Molecular sexing methods developed in the last two decades have often shown limited efficiency, particularly in terms of sensitivity and specificity, due to the degradation of DNA in these samples. These limitations have also prevented their use with ancient DNA samples of elephants or mammoths. Here we propose a novel TaqMan-MGB qPCR assay to address these difficulties. We designed it specifically to allow the characterization of the genetic sex for highly degraded samples of all elephantine taxa (elephants and mammoths). In vitro experiments demonstrated a high level of sensitivity and low contamination risks. We applied this assay in two actual case studies where it consistently recovered the right genotype for specimens of known sex a priori. In the context of a modern conservation survey of African elephants, it allowed determining the sex for over 99% of fecal samples. In a paleogenetic analysis of woolly mammoths, it produced a robust hypothesis of the sex for over 65% of the specimens out of three PCR replicates. This simple, rapid, and cost-effective procedure makes it readily applicable to large sample sizes.

Examining Natural History through the Lens of Palaeogenomics

The many high-resolution tools that are uniquely applicable to specimens from the Quaternary period (the past ~2.5 Ma) provide an opportunity to cross-validate data and test hypotheses based on the morphology and distribution of fossils. Among these tools is palaeogenomics - the genome-scale sequencing of genetic material from ancient specimens - that can provide direct insight into ecology and evolution, potentially improving the accuracy of inferences about past ecological communities over longer timescales. Palaeogenomics has revealed instances of over- and underestimation of extinct diversity, detected cryptic faunal migration and turnover, allowed quantification of widespread sex biases and sexual dimorphism in the fossil record, revealed past hybridisation events and hybrid individuals, and has highlighted previously unrecognised routes of zoonotic disease transfer.

An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole-genome sequencing data

Ancient DNA research has developed rapidly over the past few decades due to improvements in PCR and next-generation sequencing (NGS) technologies, but challenges still exist. One major challenge in relation to ancient DNA research is to recover genuine endogenous ancient DNA sequences from raw sequencing data. This is often difficult due to degradation of ancient DNA and high levels of contamination, especially homologous contamination that has extremely similar genetic background with that of the real ancient DNA. In this study, we collected whole-genome sequencing (WGS) data from 6 ancient samples to compare different mapping algorithms. To further explore more effective methods to separate endogenous DNA from homologous contaminations, we attempted to recover reads based on ancient DNA specific characteristics of deamination, depurination, and DNA fragmentation with different parameters. We propose a quick and improved pipeline for separating endogenous ancient DNA while simultaneously decreasing homologous contaminations to very low proportions. Our goal in this research was to develop useful recommendations for ancient DNA mapping and for separation of endogenous DNA to facilitate future studies of ancient DNA.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

Pleistocene population dynamics can inform the consequences of current climate change. This phylogeography of 35 complete American mastodon mitochondrial genomes suggests distinct lineages in this species repeatedly expanded northwards and then went locally extinct in response to glacial cycles.

Molecular Clocks without Rocks: New Solutions for Old Problems

Molecular data have been used to date species divergences ever since they were described as documents of evolutionary history in the 1960s. Yet, an inadequate fossil record and discordance between gene trees and species trees are persistently problematic. We examine how, by accommodating gene tree discordance and by scaling branch lengths to absolute time using mutation rate and generation time, multispecies coalescent (MSC) methods can potentially overcome these challenges. We find that time estimates can differ - in some cases, substantially - depending on whether MSC methods or traditional phylogenetic methods that apply concatenation are used, and whether the tree is calibrated with pedigree-based mutation rates or with fossils. We discuss the advantages and shortcomings of both approaches and provide practical guidance for data analysis when using these methods.

Museomics of tree squirrels: a dense taxon sampling of mitogenomes reveals hidden diversity, phenotypic convergence, and the need of a taxonomic overhaul

BACKGROUND

Tree squirrels (Sciuridae, Sciurini), in particular the highly diverse Neotropical lineages, are amongst the most rapidly diversifying branches of the mammal tree of life but also some of the least known. Negligence of this group by systematists is likely a product of the difficulties in assessing morphological informative traits and of the scarcity or unavailability of fresh tissue samples for DNA sequencing. The highly discrepant taxonomic arrangements are a consequence of the lack of phylogenies and the exclusive phenotypic-based classifications, which can be misleading in a group with conservative morphology. Here we used high-throughput sequencing and an unprecedented sampling of museum specimens to provide the first comprehensive phylogeny of tree squirrels, with a special emphasis on Neotropical taxa.

RESULTS

We obtained complete or partial mitochondrial genomes from 232 historical and modern samples, representing 40 of the 43 currently recognized species of Sciurini. Our phylogenetic analyses-performed with datasets differing on levels of missing data and taxa under distinct analytical methods-strongly support the monophyly of Sciurini and consistently recovered 12 major clades within the tribe. We found evidence that the diversity of Neotropical tree squirrels is underestimated, with at least six lineages that represent taxa to be named or revalidated. Ancestral state reconstructions of number of upper premolars and number of mammae indicated that alternative conditions of both characters must have evolved multiple times throughout the evolutionary history of tree squirrels.

CONCLUSIONS

Complete mitogenomes were obtained from museum specimens as old as 120 years, reinforcing the potential of historical samples for phylogenetic inferences of elusive lineages of the tree of life. None of the taxonomic arrangements ever proposed for tree squirrels fully corresponded to our phylogenetic reconstruction, with only a few of the currently recognized genera recovered as monophyletic. By investigating the evolution of two morphological traits widely employed in the taxonomy of the group, we revealed that their homoplastic nature can help explain the incongruence between phylogenetic results and the classification schemes presented so far. Based on our phylogenetic results we suggest a tentative supraspecific taxonomic arrangement for Sciurini, employing 13 generic names used in previous taxonomic classifications.

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Many factors complicate the estimation of time scales for phylogenetic histories, requiring increasingly complex evolutionary models and inference procedures. The widespread application of molecular clock dating has led to the insight that evolutionary rate estimates may vary with the time frame of measurement. This is particularly well established for rapidly evolving viruses that can accumulate sequence divergence over years or even months. However, this rapid evolution stands at odds with a relatively high degree of conservation of viruses or endogenous virus elements over much longer time scales. Building on recent insights into time-dependent evolutionary rates, we develop a formal and flexible Bayesian statistical inference approach that accommodates rate variation through time. We evaluate the novel molecular clock model on a foamy virus cospeciation history and a lentivirus evolutionary history and compare the performance to other molecular clock models. For both virus examples, we estimate a similarly strong time-dependent effect that implies rates varying over four orders of magnitude. The application of an analogous codon substitution model does not implicate long-term purifying selection as the cause of this effect. However, selection does appear to affect divergence time estimates for the less deep evolutionary history of the Ebolavirus genus. Finally, we explore the application of our approach on woolly mammoth ancient DNA data, which shows a much weaker, but still important, time-dependent rate effect that has a noticeable impact on node age estimates. Future developments aimed at incorporating more complex evolutionary processes will further add to the broad applicability of our approach.

Demographic reconstruction from ancient DNA supports rapid extinction of the great auk

The great auk was once abundant and distributed across the North Atlantic. It is now extinct, having been heavily exploited for its eggs, meat, and feathers. We investigated the impact of human hunting on its demise by integrating genetic data, GPS-based ocean current data, and analyses of population viability. We sequenced complete mitochondrial genomes of 41 individuals from across the species' geographic range and reconstructed population structure and population dynamics throughout the Holocene. Taken together, our data do not provide any evidence that great auks were at risk of extinction prior to the onset of intensive human hunting in the early 16th century. In addition, our population viability analyses reveal that even if the great auk had not been under threat by environmental change, human hunting alone could have been sufficient to cause its extinction. Our results emphasise the vulnerability of even abundant and widespread species to intense and localised exploitation.

Large-scale mitogenomic analysis of the phylogeography of the Late Pleistocene cave bear

The cave bear (Ursus spelaeus) is one of the Late Pleistocene megafauna species that faced extinction at the end of the last ice age. Although it is represented by one of the largest fossil records in Europe and has been subject to several interdisciplinary studies including palaeogenetic research, its fate remains highly controversial. Here, we used a combination of hybridisation capture and next generation sequencing to reconstruct 59 new complete cave bear mitochondrial genomes (mtDNA) from 14 sites in Western, Central and Eastern Europe. In a Bayesian phylogenetic analysis, we compared them to 64 published cave bear mtDNA sequences to reconstruct the population dynamics and phylogeography during the Late Pleistocene. We found five major mitochondrial DNA lineages resulting in a noticeably more complex biogeography of the European lineages during the last 50,000 years than previously assumed. Furthermore, our calculated effective female population sizes suggest a drastic cave bear population decline starting around 40,000 years ago at the onset of the Aurignacian, coinciding with the spread of anatomically modern humans in Europe. Thus, our study supports a potential significant human role in the general extinction and local extirpation of the European cave bear and illuminates the fate of this megafauna species.

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.

Climate-driven ecological stability as a globally shared cause of Late Quaternary megafaunal extinctions: the Plaids and Stripes Hypothesis

Controversy persists about why so many large-bodied mammal species went extinct around the end of the last ice age. Resolving this is important for understanding extinction processes in general, for assessing the ecological roles of humans, and for conserving remaining megafaunal species, many of which are endangered today. Here we explore an integrative hypothesis that asserts that an underlying cause of Late Quaternary megafaunal extinctions was a fundamental shift in the spatio-temporal fabric of ecosystems worldwide. This shift was triggered by the loss of the millennial-scale climate fluctuations that were characteristic of the ice age but ceased approximately 11700 years ago on most continents. Under ice-age conditions, which prevailed for much of the preceding 2.6 Ma, these radical and rapid climate changes prevented many ecosystems from fully equilibrating with their contemporary climates. Instead of today's 'striped' world in which species' ranges have equilibrated with gradients of temperature, moisture, and seasonality, the ice-age world was a disequilibrial 'plaid' in which species' ranges shifted rapidly and repeatedly over time and space, rarely catching up with contemporary climate. In the transient ecosystems that resulted, certain physiological, anatomical, and ecological attributes shared by megafaunal species pre-adapted them for success. These traits included greater metabolic and locomotory efficiency, increased resistance to starvation, longer life spans, greater sensory ranges, and the ability to be nomadic or migratory. When the plaid world of the ice age ended, many of the advantages of being large were either lost or became disadvantages. For instance in a striped world, the low population densities and slow reproductive rates associated with large body size reduced the resiliency of megafaunal species to population bottlenecks. As the ice age ended, the downsides of being large in striped environments lowered the extinction thresholds of megafauna worldwide, which then increased the vulnerability of individual species to a variety of proximate threats they had previously tolerated, such as human predation, competition with other species, and habitat loss. For many megafaunal species, the plaid-to-stripes transition may have been near the base of a hierarchy of extinction causes whose relative importances varied geographically, temporally, and taxonomically.

Accessible molecular phylogenomics at no cost: obtaining 14 new mitogenomes for the ant subfamily Pseudomyrmecinae from public data

The advent of Next Generation Sequencing has reduced sequencing costs and increased genomic projects from a huge amount of organismal taxa, generating an unprecedented amount of genomic datasets publicly available. Often, only a tiny fraction of outstanding relevance of the genomic data produced by researchers is used in their works. This fact allows the data generated to be recycled in further projects worldwide. The assembly of complete mitogenomes is frequently overlooked though it is useful to understand evolutionary relationships among taxa, especially those presenting poor mtDNA sampling at the level of genera and families. This is exactly the case for ants (Hymenoptera:Formicidae) and more specifically for the subfamily Pseudomyrmecinae, a group of arboreal ants with several cases of convergent coevolution without any complete mitochondrial sequence available. In this work, we assembled, annotated and performed comparative genomics analyses of 14 new complete mitochondria from Pseudomyrmecinae species relying solely on public datasets available from the Sequence Read Archive (SRA). We used all complete mitogenomes available for ants to study the gene order conservation and also to generate two phylogenetic trees using both (i) concatenated set of 13 mitochondrial genes and (ii) the whole mitochondrial sequences. Even though the tree topologies diverged subtly from each other (and from previous studies), our results confirm several known relationships and generate new evidences for sister clade classification inside Pseudomyrmecinae clade. We also performed a synteny analysis for Formicidae and identified possible sites in which nucleotidic insertions happened in mitogenomes of pseudomyrmecine ants. Using a data mining/bioinformatics approach, the current work increased the number of complete mitochondrial genomes available for ants from 15 to 29, demonstrating the unique potential of public databases for mitogenomics studies. The wide applications of mitogenomes in research and presence of mitochondrial data in different public dataset types makes the "no budget mitogenomics" approach ideal for comprehensive molecular studies, especially for subsampled taxa.

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%.

Ancient DNA Resolves the History of Tetragnatha (Araneae, Tetragnathidae) Spiders on Rapa Nui

Rapa Nui is one of the most remote islands in the world. As a young island, its biota is a consequence of both natural dispersals over the last ~1 million years and recent human introductions. It therefore provides an opportunity to study a unique community assemblage. Here, we extract DNA from museum-preserved and newly field-collected spiders from the genus Tetragnatha to explore their history on Rapa Nui. Using an optimized protocol to recover ancient DNA from museum-preserved spiders, we sequence and assemble partial mitochondrial genomes from nine Tetragnatha species, two of which were found on Rapa Nui, and estimate the evolutionary relationships between these and other Tetragnatha species. Our phylogeny shows that the two Rapa Nui species are not closely related. One, the possibly extinct, T. paschae, is nested within a circumtropical species complex (T. nitens), and the other (Tetragnatha sp. Rapa Nui) appears to be a recent human introduction. Our results highlight the power of ancient DNA approaches in identifying cryptic and rare species, which can contribute to our understanding of the global distribution of biodiversity in all taxonomic lineages.

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.

The Small and the Dead: A Review of Ancient DNA Studies Analysing Micromammal Species

The field of ancient DNA (aDNA) has recently been in a state of exponential growth, largely driven by the uptake of Next Generation Sequencing (NGS) techniques. Much of this work has focused on the mammalian megafauna and ancient humans, with comparatively less studies looking at micromammal fauna, despite the potential of these species in testing evolutionary, environmental and taxonomic theories. Several factors make micromammal fauna ideally suited for aDNA extraction and sequencing. Micromammal subfossil assemblages often include the large number of individuals appropriate for population level analyses, and, furthermore, the assemblages are frequently found in cave sites where the constant temperature and sheltered environment provide favourable conditions for DNA preservation. This review looks at studies that include the use of aDNA in molecular analysis of micromammal fauna, in order to examine the wide array of questions that can be answered in the study of small mammals using new palaeogenetic techniques. This study highlights the bias in current aDNA studies and assesses the future use of aDNA as a tool for the study of micromammal fauna.