We identified a 3000-year-old specimen from the Traditional Territory of the Tr’ondëk Hwëch’in in central Yukon Territory, Canada as the first known mummified remains of an ancient North American porcupine (Erethizon dorsatum), known as “Ts’ey” in the Hän language, using genetic analysis and metagenomic validation. Our analysis of the sample yielded the first-ever complete ancient mitochondrial genome for (E. dorsatum) and only the second full mitogenome for the species. Its Holocene age is considerably younger than the Pleistocene megafauna typically recovered in the Yukon permafrost, demonstrating the potential for these deposits to preserve specimens from interglacial periods. Crucially, this finding confirms the presence of porcupines in the region 3000 years ago, in line with the hypothesis that this species only dispersed into Yukon and Alaska following the establishment of boreal forests after the Last Glacial Period.

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.

Ancient genomic studies have extensively explored human-microbial interactions, yet research on non-human animals remains limited. In this study, we analyzed ancient microbial DNA from 483 mammoth remains spanning over 1 million years, including 440 newly sequenced and unpublished samples from a 1.1-million-year-old steppe mammoth. Using metagenomic screening, contaminant filtering, damage pattern analysis, and phylogenetic inference, we identified 310 microbes associated with different mammoth tissues. While most microbes were environmental or post-mortem colonizers, we recovered genomic evidence of six host-associated microbial clades spanning Actinobacillus, Pasteurella, Streptococcus, and Erysipelothrix. Some of these clades contained putative virulence factors, including a Pasteurella-related bacterium that had previously been linked to the deaths of African elephants. Notably, we reconstructed partial genomes of Erysipelothrix from the oldest mammoth sample, representing the oldest authenticated host-associated microbial DNA to date. This work demonstrates the potential of obtaining ancient animal microbiomes, which can inform further paleoecological and evolutionary research.

The Norwegian lemming (Lemmus lemmus) is a small rodent distributed across the Fennoscandian mountain tundra and the Kola Peninsula. The Norwegian lemming likely evolved during the Late Pleistocene and inhabited Fennoscandia shortly prior to the Last Glacial Maximum. However, the exact timing and origins of the species, and its phylogenetic position relative to the closely related Siberian lemming (Lemmus sibiricus) remain disputed. Moreover, the presence of ancient or contemporary gene flow between both species is largely untested. The Norwegian lemming displays characteristic phenotypic and behavioral adaptations (e.g., coat color, aggression) that are not present in other Lemmus species. We generated a de novo genome assembly for the Norwegian lemming and resequenced nine modern and two ancient Lemmus spp. genomes. We show that all Lemmus species form distinct monophyletic clades, with concordant topology between the mitochondrial and nuclear genome phylogenies. The Siberian lemming is divided into two distinct but paraphyletic clades, one in the east and one in the west, where the western clade represents a sister taxon to the Norwegian lemming. We estimate that the Norwegian and western Siberian lemming diverged shortly before the Last Glacial Maximum, making the Norwegian lemming one of the youngest known mammalian species. We did not find any indication of gene flow between L. lemmus and L. sibiricus, suggesting postglacial isolation of L. lemmus. Furthermore, we identify species-specific genomic differences in genes related to coat color and fat transport, which are likely associated with the distinctive coloration and overwintering behavior observed in the Norwegian lemming.

Climate changes profoundly impact species distributions and can drastically alter dynamics between formerly isolated taxa. The evolution of mammoths within North America was characterized by repeated cycles of dispersal and putative gene flow between woolly and Columbian mammoths. However, as genome-wide studies on mammoths have predominantly focused on Siberia, the consequences of these North American range shifts remain unclear. Here, we generated genome-wide and morphological data for two Late Pleistocene mammoth molars from British Columbia, Canada (BC), and jointly analysed these with previously published data. Our genome-wide analysis (n = 16) revealed gene flow between woolly and Columbian mammoths that would have gone undiscovered based on morphological (n = 48) and mitochondrial analysis (n = 124) alone. Consistent with their hybrid nature, our analyses suggest that these two BC mammoths had elevated genomic diversity. Our results highlight the importance of combining data types to reconstruct past evolutionary events. These findings demonstrate how the geographical range expansion of woolly mammoths resulted in long-term hybridization with local Columbian mammoths and enhance our understanding of the genomic and morphological consequences of climate-mediated dispersal.

Genetic diversity is fundamental to the adaptive potential and survival of species. Although its importance has long been recognized in science, it has a history of neglect within policy, until now. The new Global Biodiversity Framework recently adopted by the Convention on Biological Diversity, states that genetic diversity must be maintained at levels assuring adaptive potential of populations, and includes metrics for systematic monitoring of genetic diversity in so called indicators. Similarly, indicators for genetic diversity are being developed at national levels. Here, we apply new indicators for Swedish national use to one of the northernmost salmonid fishes, the Arctic charr (Salvelinus alpinus). We sequence whole genomes to monitor genetic diversity over four decades in three landlocked populations inhabiting protected alpine lakes in central Sweden. We find levels of genetic diversity, inbreeding and load to differ among lakes but remain stable over time. Effective population sizes are generally small (< 500), suggesting a limited ability to maintain adaptive variability if genetic exchange with nearby populations became eliminated. We identify genomic regions potentially shaped by selection; SNPs exhibiting population divergence exceeding expectations under drift and a putative selective sweep acting within one lake to which the competitive brown trout (Salmo trutta) was introduced during the sampling period. Identified genes appear involved in immunity and salinity tolerance. Present results suggest that genetically vulnerable populations of Arctic charr have maintained neutral and putatively adaptive genetic diversity despite small effective sizes, attesting the importance of continued protection and assurance of gene flow among populations.

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity – so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's ϴ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (F_ROH) spans 10%–30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths’ extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Ancient genomes provide a tool to investigate the genetic basis of adaptations in extinct organisms. However, the identification of species-specific fixed genetic variants requires the analysis of genomes from multiple individuals. Moreover, the long-term scale of adaptive evolution coupled with the short-term nature of tradi-tional time series data has made it difficult to assess when different adaptations evolved. Here, we analyze 23 woolly mammoth genomes, including one of the oldest known specimens at 700,000 years old, to identify fixed derived non-synonymous mutations unique to the species and to obtain estimates of when these mutations evolved. We find that at the time of its origin, the woolly mammoth had already acquired a broad spectrum of positively selected genes, including ones associated with hair and skin development, fat storage and metabolism, and immune system function. Our results also suggest that these phenotypes continued to evolve during the last 700,000 years, but through positive selection on different sets of genes. Finally, we also identify additional genes that underwent comparatively recent positive selection, including multiple genes related to skeletal morphology and body size, as well as one gene that may have contributed to the small ear size in Late Quaternary woolly mammoths.

Woolly mammoths had a set of adaptations that enabled them to thrive in the Arctic environment. Many mammoth-specific single nucleotide polymorphisms (SNPs) responsible for unique mammoth traits have been previously identified from ancient genomes. However, a multitude of other genetic variants likely contributed to woolly mammoth evolution. In this study, we sequenced two woolly mammoth genomes and combined these with previously sequenced mammoth and elephant genomes to conduct a survey of mammoth-specific deletions and indels. We find that deletions are highly enriched in non-coding regions, suggesting selection against structural variants that affect protein sequences. Nonetheless, at least 87 woolly mammoth genes contain deletions or indels that modify the coding sequence, including genes involved in skeletal morphology and hair growth. These results suggest that deletions and indels contributed to the unique phenotypic adaptations of the woolly mammoth, and were potentially critical to surviving in its natural environment.