Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

Strengthening global-change science by integrating aeDNA with paleoecoinformatics

Ancient environmental DNA (aeDNA) data are close to enabling insights into past global-scale biodiversity dynamics at unprecedented taxonomic extent and resolution. However, achieving this potential requires solutions that bridge bioinformatics and paleoecoinformatics. Essential needs include support for dynamic taxonomic inferences, dynamic age inferences, and precise stratigraphic depth. Moreover, aeDNA data are complex and heterogeneous, generated by dispersed researcher networks, with methods advancing rapidly. Hence, expert community governance and curation are essential to building high-value data resources. Immediate recommendations include uploading metabarcoding-based taxonomic inventories into paleoecoinformatic resources, building linkages among open bioinformatic and paleoecoinformatic data resources, harmonizing aeDNA processing workflows, and expanding community data governance. These advances will enable transformative insights into global-scale biodiversity dynamics during large environmental and anthropogenic changes.

Genomics of adaptive evolution in the woolly mammoth

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

High resolution ancient sedimentary DNA shows that alpine plant diversity is associated with human land use and climate change

The European Alps are highly rich in species, but their future may be threatened by ongoing changes in human land use and climate. Here, we reconstructed vegetation, temperature, human impact and livestock over the past ~12,000 years from Lake Sulsseewli, based on sedimentary ancient plant and mammal DNA, pollen, spores, chironomids, and microcharcoal. We assembled a highly-complete local DNA reference library (PhyloAlps, 3923 plant taxa), and used this to obtain an exceptionally rich sedaDNA record of 366 plant taxa. Vegetation mainly responded to climate during the early Holocene, while human activity had an additional influence on vegetation from 6 ka onwards. Land-use shifted from episodic grazing during the Neolithic and Bronze Age to agropastoralism in the Middle Ages. Associated human deforestation allowed the coexistence of plant species typically found at different elevational belts, leading to levels of plant richness that characterise the current high diversity of this region. Our findings indicate a positive association between low intensity agropastoral activities and precipitation with the maintenance of the unique subalpine and alpine plant diversity of the European Alps.

Paleoeconomy more than demography determined prehistoric human impact in Arctic Norway

Population size has increasingly been taken as the driver of past human environmental impact worldwide, and particularly in the Arctic. However, sedimentary ancient DNA (sedaDNA), pollen and archaeological data show that over the last 12,000 years, paleoeconomy and culture determined human impacts on the terrestrial ecology of Arctic Norway. The large Mortensnes site complex (Ceavccageađgi, 70°N) has yielded the most comprehensive multiproxy record in the Arctic to date. The site saw occupation from the Pioneer period (c. 10,000 cal. years BP) with more intensive use from c. 4,200 to 2,000 cal. years BP and after 1,600 cal. years BP. Here, we combine on-site environmental archaeology with a near-site lake record of plant and animal sedaDNA. The rich animal sedaDNA data (42 taxa) and on-site faunal analyses reveal switches in human dietary composition from early-Holocene fish + marine mammals, to mixed marine + reindeer, then finally to marine + reindeer + domesticates (sheep, cattle, pigs), with highest reindeer concentrations in the last millennium. Archaeological evidence suggests these changes are not directly driven by climate or variation in population densities at the site or in the region, but rather are the result of changing socio-economic activities and culture, probably reflecting settlers' origins. This large settlement only had discernable effects on its hinterland in the last 3,600 years (grazing) and more markedly in the last 1,000 years through reindeer keeping/herding and, possibly domestic stock. Near-site sedaDNA can be linked to and validate the faunal record from archaeological excavations, demonstrating that environmental impacts can be assessed at a landscape scale.

Postglacial species arrival and diversity buildup of northern ecosystems took millennia

What drives ecosystem buildup, diversity, and stability? We assess species arrival and ecosystem changes across 16 millennia by combining regional-scale plant sedimentary ancient DNA from Fennoscandia with near-complete DNA and trait databases. We show that postglacial arrival time varies within and between plant growth forms. Further, arrival times were mainly predicted by adaptation to temperature, disturbance, and light. Major break points in ecological trait diversity were seen between 13.9 and 10.8 calibrated thousand years before the present (cal ka BP), as well as break point in functional diversity at 12.0 cal ka BP, shifting from a state of ecosystem buildup to a state where most habitat types and biotic ecosystem components were in place. Trait and functional diversity stabilized around 8 cal ka BP, after which both remained stable, although changes in climate took place and species inflow continued. Our ecosystem reconstruction indicates a millennial-scale time phase of formation to reach stable and resilient levels of diversity and functioning.

Evolutionary consequences of genomic deletions and insertions in the woolly mammoth genome

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.

•

Two new high-quality woolly mammoth genomes have been generated

•

A new method was used to identify deletions and insertions in woolly mammoths

•

At least 87 genes have been affected by deletions or indels in the mammoth lineage

•

Genes involved in skeletal morphology and hair growth are affected by deletions

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

Sedimentary ancient DNA shows terrestrial plant richness continuously increased over the Holocene in northern Fennoscandia

The effects of climate change on species richness are debated but can be informed by the past. Here, we generated a sedimentary ancient DNA dataset covering 10 lakes and applied novel methods for data harmonization. We assessed the impact of Holocene climate changes and nutrients on terrestrial plant richness in northern Fennoscandia. We find that richness increased steeply during the rapidly warming Early Holocene. In contrast to findings from most pollen studies, we show that richness continued to increase thereafter, although the climate was stable, with richness and the regional species pool only stabilizing during the past three millennia. Furthermore, overall increases in richness were greater in catchments with higher soil nutrient availability. We suggest that richness will increase with ongoing warming, especially at localities with high nutrient availability and assuming that human activity remains low in the region, although lags of millennia may be expected.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

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.

Community-curated and standardised metadata of published ancient metagenomic samples with AncientMetagenomeDir

Ancient DNA and RNA are valuable data sources for a wide range of disciplines. Within the field of ancient metagenomics, the number of published genetic datasets has risen dramatically in recent years, and tracking this data for reuse is particularly important for large-scale ecological and evolutionary studies of individual taxa and communities of both microbes and eukaryotes. AncientMetagenomeDir (archived at 10.5281/zenodo.3980833) is a collection of annotated metagenomic sample lists derived from published studies that provide basic, standardised metadata and accession numbers to allow rapid data retrieval from online repositories. These tables are community-curated and span multiple sub-disciplines to ensure adequate breadth and consensus in metadata definitions, as well as longevity of the database. Internal guidelines and automated checks facilitate compatibility with established sequence-read archives and term-ontologies, and ensure consistency and interoperability for future meta-analyses. This collection will also assist in standardising metadata reporting for future ancient metagenomic studies.

Measurement(s)	genome • Metagenome • Metadata • Ancient DNA
Technology Type(s)	digital curation
Factor Type(s)	geographic location • sample age

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13241537

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny

Ancient DNA (aDNA) sequencing has enabled reconstruction of speciation, migration, and admixture events for extinct taxa¹. Outside the permafrost, however, irreversible aDNA post-mortem degradation² has so far limited aDNA recovery to the past ~0.5 million years (Ma)³. Contrarily, tandem mass spectrometry (MS) allowed sequencing ~1.5 million year (Ma) old collagen type I (COL1)⁴ and suggested the presence of protein residues in Cretaceous fossil remains⁵, although with limited phylogenetic use⁶. In the absence of molecular evidence, the speciation of several Early and Middle Pleistocene extinct species remain contentious. In this study, we address the phylogenetic relationships of the Eurasian Pleistocene Rhinocerotidae^7–9 using a ~1.77 Ma old dental enamel proteome of a Stephanorhinus specimen from the Dmanisi archaeological site in Georgia (South Caucasus)¹⁰. Molecular phylogenetic analyses place the Dmanisi Stephanorhinus as a sister group to the woolly (Coelodonta antiquitatis) and Merck’s rhinoceros (S. kirchbergensis) clade. We show that Coelodonta evolved from an early Stephanorhinus lineage and that the latter includes at least two distinct evolutionary lines. As such, the genus Stephanorhinus is currently paraphyletic and its systematic revision is therefore needed. We demonstrate that Early Pleistocene dental enamel proteome sequencing overcomes the limits of ancient collagen- and aDNA-based phylogenetic inference. It also provides additional information about the sex and taxonomic assignment of the specimens analysed. Dental enamel, the hardest tissue in vertebrates¹¹, is highly abundant in the fossil record. Our findings reveal that palaeoproteomic investigation of this material can push biomolecular investigation further back into the Early Pleistocene.

Targeted Amplification and Sequencing of Ancient Environmental and Sedimentary DNA

All organisms release their DNA into the environment through processes such as excretion and the senescence of tissues and limbs. This DNA, often referred to as environmental DNA (eDNA) or sedimentary ancient DNA (sedaDNA), can be recovered from both present-day and ancient soils, fecal samples, bodies of water and lake cores, and even air. While eDNA is a potentially useful record of past and present biodiversity, several challenges complicate data generation and interpretation of results. Most importantly, eDNA samples tend to be highly taxonomically mixed, and the target organism or group of organisms may be present at very low abundance within this mixture. To overcome this challenge, enrichment approaches are often used to target specific taxa of interest. Here, we describe a protocol to amplify metabarcodes or short, variable loci that identify lineages within broad taxonomic groups (e.g., plants, mammals), using the polymerase chain reaction (PCR) with established generic "barcode" primers. We also provide a catalog of animal and plant barcode primers that, because they target relatively short fragments of DNA, are potentially suitable for use with degraded DNA.

Minimizing polymerase biases in metabarcoding

DNA metabarcoding is an increasingly popular method to characterize and quantify biodiversity in environmental samples. Metabarcoding approaches simultaneously amplify a short, variable genomic region, or "barcode," from a broad taxonomic group via the polymerase chain reaction (PCR), using universal primers that anneal to flanking conserved regions. Results of these experiments are reported as occurrence data, which provide a list of taxa amplified from the sample, or relative abundance data, which measure the relative contribution of each taxon to the overall composition of amplified product. The accuracy of both occurrence and relative abundance estimates can be affected by a variety of biological and technical biases. For example, taxa with larger biomass may be better represented in environmental samples than those with smaller biomass. Here, we explore how polymerase choice, a potential source of technical bias, might influence results in metabarcoding experiments. We compared potential biases of six commercially available polymerases using a combination of mixtures of amplifiable synthetic sequences and real sedimentary DNA extracts. We find that polymerase choice can affect both occurrence and relative abundance estimates and that the main source of this bias appears to be polymerase preference for sequences with specific GC contents. We further recommend an experimental approach for metabarcoding based on results of our synthetic experiments.

Natural selection shaped the rise and fall of passenger pigeon genomic diversity

The extinct passenger pigeon was once the most abundant bird in North America, and possibly the world. Although theory predicts that large populations will be more genetically diverse, passenger pigeon genetic diversity was surprisingly low. To investigate this disconnect, we analyzed 41 mitochondrial and 4 nuclear genomes from passenger pigeons and 2 genomes from band-tailed pigeons, which are passenger pigeons' closest living relatives. Passenger pigeons' large population size appears to have allowed for faster adaptive evolution and removal of harmful mutations, driving a huge loss in their neutral genetic diversity. These results demonstrate the effect that selection can have on a vertebrate genome and contradict results that suggested that population instability contributed to this species's surprisingly rapid extinction.

A new genus of horse from Pleistocene North America

The extinct ‘New World stilt-legged’, or NWSL, equids constitute a perplexing group of Pleistocene horses endemic to North America. Their slender distal limb bones resemble those of Asiatic asses, such as the Persian onager. Previous palaeogenetic studies, however, have suggested a closer relationship to caballine horses than to Asiatic asses. Here, we report complete mitochondrial and partial nuclear genomes from NWSL equids from across their geographic range. Although multiple NWSL equid species have been named, our palaeogenomic and morphometric analyses support the idea that there was only a single species of middle to late Pleistocene NWSL equid, and demonstrate that it falls outside of crown group Equus. We therefore propose a new genus, Haringtonhippus, for the sole species H. francisci. Our combined genomic and phenomic approach to resolving the systematics of extinct megafauna will allow for an improved understanding of the full extent of the terminal Pleistocene extinction event.

The horse family – which also includes zebras, donkeys and asses – is often featured on the pages of textbooks about evolution. All living horses belong to a group, or genus, called Equus. The fossil record shows how the ancestors of these animals evolved from dog-sized, three-toed browsers to larger, one-toed grazers. This process took around 55 million years, and many members of the horse family tree went extinct along the way.

Nevertheless, the details of the horse family tree over the past 2.5 million years remain poorly understood. In North America, horses from this period – which is referred to as the Pleistocene – have been classed into two major groups: stout-legged horses and stilt-legged horses. Both groups became extinct near the end of the Pleistocene in North America, and it was not clear how they relate to one another. Based on their anatomy, many scientists suggested that stilt-legged horses were most closely related to modern-day asses living in Asia. Yet, other studies using ancient DNA placed the stilt-legged horses closer to the stout-legged horses.

Heintzman et al. set out to resolve where the stilt-legged horses sit within the horse family tree by examining more ancient DNA than the previous studies. The analyses showed that the stilt-legged horses were much more distinct than previously thought. In fact, contrary to all previous findings, these animals actually belonged outside of the genus Equus. Heintzman et al. named the new genus for the stilt-legged horses Haringtonhippus, and showed that all stilt-legged horses belonged to a single species within this genus, Haringtonhippus francisci.

Together these new findings provide a benchmark for reclassifying problematic fossil groups across the tree of life. A similar approach could be used to resolve the relationships in other problematic groups of Pleistocene animals, such as mammoths and bison. This would give scientists a more nuanced understanding of evolution and extinction during this period.

Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe’s largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

Genomic Data from Extinct North American Camelops Revise Camel Evolutionary History

Recent advances in paleogenomic technologies have enabled an increasingly detailed understanding of the evolutionary relationships of now-extinct mammalian taxa. However, a number of enigmatic Quaternary species have never been characterized with molecular data, often because available fossils are rare or are found in environments that are not optimal for DNA preservation. Here, we analyze paleogenomic data extracted from bones attributed to the late Pleistocene western camel, Camelops cf. hesternus, a species that was distributed across central and western North America until its extinction approximately 13,000 years ago. Despite a modal sequence length of only around 35 base pairs, we reconstructed high-coverage complete mitochondrial genomes and low-coverage partial nuclear genomes for each specimen. We find that Camelops is sister to African and Asian bactrian and dromedary camels, to the exclusion of South American camelids (llamas, guanacos, alpacas, and vicuñas). These results contradict previous morphology-based phylogenetic models for Camelops, which suggest instead a closer relationship between Camelops and the South American camelids. The molecular data imply a Late Miocene divergence of the Camelops clade from lineages that separately gave rise to the extant camels of Eurasia. Our results demonstrate the increasing capacity of modern paleogenomic methods to resolve evolutionary relationships among distantly related lineages.

Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch

Viruses preserved in ancient materials provide snapshots of past viral diversity and a means to trace viral evolution through time. Here, we use a metagenomics approach to identify filterable and nuclease-resistant nucleic acids preserved in 700-y-old caribou feces frozen in a permanent ice patch. We were able to recover and characterize two viruses in replicated experiments performed in two different laboratories: a small circular DNA viral genome (ancient caribou feces associated virus, or aCFV) and a partial RNA viral genome (Ancient Northwest Territories cripavirus, or aNCV). Phylogenetic analysis identifies aCFV as distantly related to the plant-infecting geminiviruses and the fungi-infecting Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 and aNCV as within the insect-infecting Cripavirus genus. We hypothesize that these viruses originate from plant material ingested by caribou or from flying insects and that their preservation can be attributed to protection within viral capsids maintained at cold temperatures. To investigate the tropism of aCFV, we used the geminiviral reverse genetic system and introduced a multimeric clone into the laboratory model plant Nicotiana benthamiana. Evidence for infectivity came from the detection of viral DNA in newly emerged leaves and the precise excision of the viral genome from the multimeric clones in inoculated leaves. Our findings indicate that viral genomes may in some circumstances be protected from degradation for centuries.

Characterizing DNA preservation in degraded specimens of Amara alpina (Carabidae: Coleoptera)

DNA preserved in degraded beetle (Coleoptera) specimens, including those derived from dry-stored museum and ancient permafrost-preserved environments, could provide a valuable resource for researchers interested in species and population histories over timescales from decades to millenia. However, the potential of these samples as genetic resources is currently unassessed. Here, using Sanger and Illumina shotgun sequence data, we explored DNA preservation in specimens of the ground beetle Amara alpina, from both museum and ancient environments. Nearly all museum specimens had amplifiable DNA, with the maximum amplifiable fragment length decreasing with age. Amplification of DNA was only possible in 45% of ancient specimens. Preserved mitochondrial DNA fragments were significantly longer than those of nuclear DNA in both museum and ancient specimens. Metagenomic characterization of extracted DNA demonstrated that parasite-derived sequences, including Wolbachia and Spiroplasma, are recoverable from museum beetle specimens. Ancient DNA extracts contained beetle DNA in amounts comparable to museum specimens. Overall, our data demonstrate that there is great potential for both museum and ancient specimens of beetles in future genetic studies, and we see no reason why this would not be the case for other orders of insect.