Large-scale DNA screening of palaeontological and archaeological collections remains a limiting and costly factor for ancient DNA studies. Several DNA extraction protocols are routinely used in ancient DNA laboratories and have even been automated on robotic platforms. Robots offer a solution for high-throughput screening but the costs, as well as necessity for trained technicians and engineers, can be prohibitive for some laboratories. Here, we present a high-throughput alternative to robot-based ancient DNA extraction using a 96-column plate. When compared to routine single MinElute columns, we retrieved highly similar endogenous DNA contents, an important metric in ancient DNA screening. Mitogenomes with a coverage depth greater than 0.1× could be generated and allowed for taxonomic assignment. However, average fragment lengths, DNA damage and library complexities significantly differed between methods but these differences became nonsignificant after modification of our library purification protocol. Our high-throughput extraction method allows generation of 96 extracts within approximately 4 hours of laboratory work while bringing the cost down by ~39% compared to using single columns. Additionally, we formally demonstrate that the addition of Tween-20 during the elution step results in higher complexity libraries, thereby enabling higher genome coverage for the same sequencing effort.

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.

Inbreeding increases genome homozygosity within populations, which can exacerbate inbreeding depression by exposing homozygous deleterious alleles that are responsible for declines in fitness traits. In small populations, genetic purging that occurs under pressure of natural selection acts as an opposing force, contributing to a reduction of deleterious alleles. Both inbreeding and genetic purging are paramount in the field of conservation genomics. The Amur tiger (Panthera tigris altaica) lives in small populations in the forests of Northeast Asia and is among the most endangered animals on the planet. Using genome-wide assessment and comparison, we reveal substantially higher and more extensive inbreeding in wild Amur tigers (FROH = 0.50) than in captive individuals (FROH = 0.24). However, a relatively reduced number of loss-of-function mutations in wild Amur tigers is observed when compared with captive individuals, indicating a genetic purging of relatively large-effect inbreeding load. The higher ratio of homozygous mutation load and number of fixed damaging alleles in the wild population indicate a less-efficient genetic purging, with purifying selection also contributing to this process. These findings provide valuable insights for future conservation of Amur tigers.

In the last three decades, DNA sequencing of ancient animal osteological assemblages has become an important tool complementing standard archaeozoological approaches to reconstruct the history of animal domestication. However, osteological assemblages of key archaeological contexts are not always available or do not necessarily preserve enough ancient DNA for a cost-effective genetic analysis. Here, we develop an in-solution target-enrichment approach, based on 80-mer species-specific RNA probes (ranging from 306 to 1686 per species) to characterise (in single experiments) the mitochondrial genetic variation from eight domesticated animal species of major economic interest: cattle, chickens, dogs, donkeys, goats, horses, pigs and sheep. We also illustrate how our design can be adapted to enrich DNA library content and map the Y-chromosomal diversity within Equus caballus. By applying our target-enrichment assay to an extensive panel of ancient osteological remains, farm soil, and cave sediments spanning the last 43 kyrs, we demonstrate that minimal sequencing efforts are necessary to exhaust the DNA library complexity and to characterise mitogenomes to an average depth-of-coverage of 19.4 to 2003.7-fold. Our assay further retrieved horse mitogenome and Y-chromosome data from Late Pleistocene coprolites, as well as bona fide mitochondrial sequences from species that were not part of the probe design, such as bison and cave hyena. Our methodology will prove especially useful to minimise costs related to the genetic analyses of maternal and paternal lineages of a wide range of domesticated and wild animal species, and for mapping their diversity changes over space and time, including from environmental samples.

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Interest continues to grow in Arctic megafaunal ecological engineering, but, since the mass extinction of megafauna ~ 12–15 ka, key physiographic variables and available forage continue to change. Here we sought to assess the extent to which contemporary Arctic ecosystems are conducive to the rewilding of megaherbivores, using a woolly mammoth (M. primigenius) proxy as a model species. We first perform a literature review on woolly mammoth dietary habits. We then leverage Oak Ridge National Laboratories Distributive Active Archive Center Global Aboveground and Belowground Biomass Carbon Density Maps to generate aboveground biomass carbon density estimates in plant functional types consumed by the woolly mammoth at 300 m resolution on Alaska’s North Slope. We supplement these analyses with a NASA Arctic Boreal Vulnerability Experiment dataset to downgrade overall biomass estimates to digestible levels. We further downgrade available forage by using a conversion factor representing the relationship between total biomass and net primary productivity (NPP) for arctic vegetation types. Integrating these estimates with the forage needs of woolly mammoths, we conservatively estimate Alaska’s North Slope could support densities of 0.0–0.38 woolly mammoth km⁻² (mean 0.13) across a variety of habitats. These results may inform innovative rewilding strategies.

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability. Moreover, small population size may promote gene flow with congeners and outbreeding depression. Here, we examine the connection between habitat availability, effective population size (N_e), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa, populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat. However, we found long-term low N_e, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low N_e and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species’ unique traits.

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1–2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

Genomic approaches have the potential to play a pivotal role in conservation, both to detect threats to species and populations and to restore biodiversity through actions. We here separate these approaches into two subdisciplines, vulnerability and restoration genomics, and discuss current applications, outstanding questions, and future potential.