Wide-spread dispersal in a deep-sea brooding polychaete: the role of natural history collections in assessing the distribution in quill worms (Onuphidae, Annelida).

BACKGROUND: Modern integrative taxonomy-based annelid species descriptions are detailed combining morphological data and, since the last decades, also molecular information. Historic species descriptions are often comparatively brief lacking such detail. Adoptions of species names from western literature in the past led to the assumption of cosmopolitan ranges for many species, which, in many cases, were later found to include cryptic or pseudocryptic lineages with subtle morphological differences. Natural history collections and databases can aid in assessing the geographic ranges of species but depend on correct species identification. Obtaining DNA sequence information from wet-collection museum specimens of marine annelids is often impeded by the use of  formaldehyde and/or long-term storage in ethanol resulting in DNA degradation and cross-linking. RESULTS: The application of ancient DNA extraction methodology in combination with single-stranded DNA library preparation and target gene capture resulted in successful sequencing of a 110-year-old collection specimen of quill worms. Furthermore, a 40-year-old specimen of quill worms was successfully sequenced using a standard extraction protocol for modern samples, PCR and Sanger sequencing. Our study presents the first molecular analysis of Hyalinoecia species including the previously known species Hyalinoecia robusta, H. tubicloa, H. artifex, and H. longibranchiata, and a potentially undescribed species from equatorial western Africa morphologically indistinguishable from H. tubicola. The study also investigates the distribution of these five Hyalinoecia species. Reassessing the distribution of H. robusta reveals a geographical range covering both the Atlantic and the Indian Oceans as indicated by molecular data obtained from recent and historical specimens. CONCLUSION: Our results represent an example of a very wide geographical distribution of a brooding deep-sea annelid with a complex reproduction strategy and seemingly very limited dispersal capacity of its offspring, and highlights the importance of molecular information from museum specimens for integrative annelid taxonomy and biogeography.

Genome Sequencing of up to 6,000-Year-Old Citrullus Seeds Reveals Use of a Bitter-Fleshed Species Prior to Watermelon Domestication.

Iconographic evidence from Egypt suggests that watermelon pulp was consumed there as a dessert by 4,360 BP. Earlier archaeobotanical evidence comes from seeds from Neolithic settlements in Libya, but whether these were watermelons with sweet pulp or other forms is unknown. We generated genome sequences from 6,000- and 3,300-year-old seeds from Libya and Sudan, and from worldwide herbarium collections made between 1824 and 2019, and analyzed these data together with resequenced genomes from important germplasm collections for a total of 131 accessions. Phylogenomic and population-genomic analyses reveal that (1) much of the nuclear genome of both ancient seeds is traceable to West African seed-use "egusi-type" watermelon (Citrullus mucosospermus) rather than domesticated pulp-use watermelon (Citrullus lanatus ssp. vulgaris); (2) the 6,000-year-old watermelon likely had bitter pulp and greenish-white flesh as today found in C. mucosospermus, given alleles in the bitterness regulators ClBT and in the red color marker LYCB; and (3) both ancient genomes showed admixture from C. mucosospermus, C. lanatus ssp. cordophanus, C. lanatus ssp. vulgaris, and even South African Citrullus amarus, and evident introgression between the Libyan seed (UMB-6) and populations of C. lanatus. An unexpected new insight is that Citrullus appears to have initially been collected or cultivated for its seeds, not its flesh, consistent with seed damage patterns induced by human teeth in the oldest Libyan material.

Large variance in inbreeding within the Iberian wolf population.

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300-350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here we measured changes on genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared to other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

Molecular Clocks and Archeogenomics of a Late Period Egyptian Date Palm Leaf Reveal Introgression from Wild Relatives and Add Timestamps on the Domestication.

The date palm, Phoenix dactylifera, has been a cornerstone of Middle Eastern and North African agriculture for millennia. It was first domesticated in the Persian Gulf, and its evolution appears to have been influenced by gene flow from two wild relatives, P. theophrasti, currently restricted to Crete and Turkey, and P. sylvestris, widespread from Bangladesh to the West Himalayas. Genomes of ancient date palm seeds show that gene flow from P. theophrasti to P. dactylifera may have occurred by ∼2,200 years ago, but traces of P. sylvestris could not be detected. We here integrate archeogenomics of a ∼2,100-year-old P. dactylifera leaf from Saqqara (Egypt), molecular-clock dating, and coalescence approaches with population genomic tests, to probe the hybridization between the date palm and its two closest relatives and provide minimum and maximum timestamps for its reticulated evolution. The Saqqara date palm shares a close genetic affinity with North African date palm populations, and we find clear genomic admixture from both P. theophrasti, and P. sylvestris, indicating that both had contributed to the date palm genome by 2,100 years ago. Molecular-clocks placed the divergence of P. theophrasti from P. dactylifera/P. sylvestris and that of P. dactylifera from P. sylvestris in the Upper Miocene, but strongly supported, conflicting topologies point to older gene flow between P. theophrasti and P. dactylifera, and P. sylvestris and P. dactylifera. Our work highlights the ancient hybrid origin of the date palms, and prompts the investigation of the functional significance of genetic material introgressed from both close relatives, which in turn could prove useful for modern date palm breeding.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Dataset from RNAseq analysis of differential gene expression among developmental stages of two non-marine ostracodes.

We contribute transcriptomic data for two species of Ostracoda, an early-diverged group of small-sized pancrustaceans. Data include new reference transcriptomes for two asexual non-marine species (Dolerocypris sinensis and Heterocypris aff. salina), as well as single-specimen transcriptomic data that served to analyse gene expression across four developmental stages in D. sinensis. Data are evaluated by computing gene expression profiles of the different developmental stages which consistently placed eggs and small larvae (at the stage of instar A-8) similar to each other, and apart from adults which were distinct from all other developmental stages but closest to large larvae (instar A-4). We further evaluated the transcriptomic data with two newly sequenced low-coverage genomes of the target species. The new data thus document the feasibility of obtaining reliable transcriptomic data from single specimens - even eggs - of these small metazoans.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size.

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), 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,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, 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 Ne 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.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines.

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

Paleogenomes Reveal a Complex Evolutionary History of Late Pleistocene Bison in Northeastern China.

Steppe bison are a typical representative of the Mid-Late Pleistocene steppes of the northern hemisphere. Despite the abundance of fossil remains, many questions related to their genetic diversity, population structure and dispersal route are still elusive. Here, we present both near-complete and partial mitochondrial genomes, as well as a partial nuclear genome from fossil bison samples excavated from Late Pleistocene strata in northeastern China. Maximum-likelihood and Bayesian trees both suggest the bison clade are divided into three maternal haplogroups (A, B and C), and Chinese individuals fall in two of them. Bayesian analysis shows that the split between haplogroup C and the ancestor of haplogroups A and B dates at 326 ky BP (95% HPD: 397-264 ky BP). In addition, our nuclear phylogenomic tree also supports a basal position for the individual carrying haplogroup C. Admixture analyses suggest that CADG467 (haplogroup C) has a similar genetic structure to steppe bison from Siberia (haplogroup B). Our new findings indicate that the genetic diversity of Pleistocene bison was probably even higher than previously thought and that northeastern Chinese populations of several mammalian species, including Pleistocene bison, were genetically distinct.

Mitochondrial DNA sequencing of a wet-collection syntype demonstrates the importance of type material as genetic resource for lantern shark taxonomy (Chondrichthyes: Etmopteridae).

After initial detection of target archival DNA of a 116-year-old syntype specimen of the smooth lantern shark, Etmopterus pusillus, in a single-stranded DNA library, we shotgun-sequenced additional 9 million reads from this same DNA library. Sequencing reads were used for extracting mitochondrial sequence information for analyses of mitochondrial DNA characteristics and reconstruction of the mitochondrial genome. The archival DNA is highly fragmented. A total of 4599 mitochondrial reads were available for the genome reconstruction using an iterative mapping approach. The resulting genome sequence has 12 times coverage and a length of 16 741 bp. All 37 vertebrate mitochondrial loci plus the control region were identified and annotated. The mitochondrial NADH2 gene was subsequently used to place the syntype haplotype in a network comprising multiple E. pusillus samples from various distant localities as well as sequences from a morphological similar species, the shortfin smooth lantern shark Etmopterus joungi. Results confirm the almost global distribution of E. pusillus and suggest E. joungi to be a junior synonym of E. pusillus. As mitochondrial DNA often represents the only available reference information in non-model organisms, this study illustrates the importance of mitochondrial DNA from an aged, wet collection type specimen for taxonomy.

Ancient DNA from the Asiatic Wild Dog (Cuon alpinus) from Europe.

The Asiatic wild dog (Cuon alpinus), restricted today largely to South and Southeast Asia, was widespread throughout Eurasia and even reached North America during the Pleistocene. Like many other species, it suffered from a huge range loss towards the end of the Pleistocene and went extinct in most of its former distribution. The fossil record of the dhole is scattered and the identification of fossils can be complicated by an overlap in size and a high morphological similarity between dholes and other canid species. We generated almost complete mitochondrial genomes for six putative dhole fossils from Europe. By using three lines of evidence, i.e., the number of reads mapping to various canid mitochondrial genomes, the evaluation and quantification of the mapping evenness along the reference genomes and phylogenetic analysis, we were able to identify two out of six samples as dhole, whereas four samples represent wolf fossils. This highlights the contribution genetic data can make when trying to identify the species affiliation of fossil specimens. The ancient dhole sequences are highly divergent when compared to modern dhole sequences, but the scarcity of dhole data for comparison impedes a more extensive analysis.

Ancient Mitogenomes Provide New Insights into the Origin and Early Introduction of Chinese Domestic Donkeys.

Both molecular data and archaeological evidence strongly support an African origin for the domestic donkey. Recent genetic studies further suggest that there were two distinct maternal lineages involved in its initial domestication. However, the exact introduction time and the dispersal process of domestic donkeys into ancient China are still unresolved. To address these questions, we retrieved three near-complete mitochondrial genomes from donkey specimens excavated from Gaoling County, Shaanxi Province, and Linxia Basin, Gansu Province, China, dated at 2,349-2,301, 469-311, and 2,160-2,004 cal. BP, respectively. Maximum-likelihood and Bayesian phylogenetic analyses reveal that the two older samples fall into the two different main lineages (i.e., clade Ⅰ and clade Ⅱ) of the domestic donkey, suggesting that the two donkey maternal lineages had been introduced into Midwestern China at least at the opening of Silk Road (approximately the first century BC). Bayesian analysis shows that the split of the two donkey maternal lineages is dated at 0.323 Ma (95% CI: 0.583-0.191 Ma) using root-tip dating calibrations based on near-complete mitogenomes, supporting the hypothesis that modern domestic donkeys go back to at least two independent domestication events. Moreover, Bayesian skyline plot analyses indicate an apparent female population increase between 5,000 and 2,500 years ago for clade I followed by a stable population size to the present day. In contrast, clade II keeps a relatively stable population size over the past 5,000 years. Overall, our study provides new insights into the early domestication history of Chinese domestic donkeys.

Successful application of ancient DNA extraction and library construction protocols to museum wet collection specimens.

Millions of scientific specimens are housed in museum collections, a large part of which are fluid preserved. The use of formaldehyde as fixative and subsequent storage in ethanol is especially common in ichthyology and herpetology. This type of preservation damages DNA and reduces the chance of successful retrieval of genetic data. We applied ancient DNA extraction and single stranded library construction protocols to a variety of vertebrate samples obtained from wet collections and of different ages. Our results show that almost all samples tested yielded endogenous DNA. Archival DNA extraction was successful across different tissue types as well as using small amounts of tissue. Conversion of archival DNA fragments into single-stranded libraries resulted in usable data even for samples with initially undetectable DNA amounts. Subsequent target capture approaches for mitochondrial DNA using homemade baits on a subset of 30 samples resulted in almost complete mitochondrial genome sequences in several instances. Thus, application of ancient DNA methodology makes wet collection specimens, including type material as well as rare, old or extinct species, accessible for genetic and genomic analyses. Our results, accompanied by detailed step-by-step protocols, are a large step forward to open the DNA archive of museum wet collections for scientific studies.

Annotated genome sequences of the carnivorous plant Roridula gorgonias and a non-carnivorous relative, Clethra arborea.

OBJECTIVE: Plant carnivory is distributed across the tree of life and has evolved at least six times independently, but sequenced and annotated nuclear genomes of carnivorous plants are currently lacking. We have sequenced and structurally annotated the nuclear genome of the carnivorous Roridula gorgonias and that of a non-carnivorous relative, Madeira's lily-of-the-valley-tree, Clethra arborea, both within the Ericales. This data adds an important resource to study the evolutionary genetics of plant carnivory across angiosperm lineages and also for functional and systematic aspects of plants within the Ericales. RESULTS: Our assemblies have total lengths of 284 Mbp (R. gorgonias) and 511 Mbp (C. arborea) and show high BUSCO scores of 84.2% and 89.5%, respectively. We used their predicted genes together with publicly available data from other Ericales' genomes and transcriptomes to assemble a phylogenomic data set for the inference of a species tree. However, groups of orthologs showed a marked absence of species represented by a transcriptome. We discuss possible reasons and caution against combining predicted genes from genome- and transriptome-based assemblies.

'Barcode fishing' for archival DNA from historical type material overcomes taxonomic hurdles, enabling the description of a new frog species.

Taxonomic progress is often hindered by intrinsic factors, such as morphologically cryptic species that require a broad suite of methods to distinguish, and extrinsic factors, such as uncertainties in the allocation of scientific names to species. These uncertainties can be due to a wide variety of factors, including old and poorly preserved type specimens (which contain only heavily degraded DNA or have lost important diagnostic characters), inappropriately chosen type specimens (e.g. juveniles without diagnostic characters) or poorly documented type specimens (with unprecise, incorrect, or missing locality data). Thanks to modern sequencing technologies it is now possible to overcome many such extrinsic factors by sequencing DNA from name-bearing type specimens of uncertain assignment and assigning these to known genetic lineages. Here, we apply this approach to frogs of the Mantidactylus ambreensis complex, which was recently shown to consist of two genetic lineages supported by concordant differentiation in mitochondrial and nuclear genes. These lineages co-occur on the Montagne d'Ambre Massif in northern Madagascar but appear to have diverged in allopatry. We use a recently published bait set based on three mitochondrial markers from all known Malagasy frog lineages to capture DNA sequences from the 127-year-old holotype of Mantidactylus ambreensis Mocquard, 1895. With the obtained sequences we are able to assign the name M. ambreensis to the lowland lineage, which is rather widespread in the rainforests of northern Madagascar, leaving the microendemic high-elevation lineage on Montagne d'Ambre in north Madagascar in need of description. We describe this species as Mantidactylus ambony sp. nov., differing from M. ambreensis in call parameters and a smaller body size. Thus, using target enrichment to obtain DNA sequence data from this old specimen, we were able to resolve the extrinsic (nomenclatural) hindrances to taxonomic resolution of this complex. We discuss the broad-scale versatility of this 'barcode fishing' approach, which can draw on the enormous success of global DNA barcoding initiatives to quickly and efficiently assign type specimens to lineages.

Hyena paleogenomes reveal a complex evolutionary history of cross-continental gene flow between spotted and cave hyena.

The genus Crocuta (African spotted and Eurasian cave hyenas) includes several closely related extinct and extant lineages. The relationships among these lineages, however, are contentious. Through the generation of population-level paleogenomes from late Pleistocene Eurasian cave hyena and genomes from modern African spotted hyena, we reveal the cross-continental evolutionary relationships between these enigmatic hyena lineages. We find a deep divergence (~2.5 Ma) between African and Eurasian Crocuta populations, suggesting that ancestral Crocuta left Africa around the same time as early Homo. Moreover, we find discordance between nuclear and mitochondrial phylogenies and evidence for bidirectional gene flow between African and Eurasian Crocuta after the lineages split, which may have complicated prior taxonomic classifications. Last, we find a number of introgressed loci that attained high frequencies within the recipient lineage, suggesting some level of adaptive advantage from admixture.

Molecular identification of late and terminal Pleistocene Equus ovodovi from northeastern China.

The extant diversity of horses (family Equidae) represents a small fraction of that occurring over their evolutionary history. One such lost lineage is the subgenus Sussemionus, which is thought to have become extinct during the Middle Pleistocene. However, recent molecular studies and morphological analysis have revealed that one of their representatives, E. ovodovi, did exist in Siberia during the Late Pleistocene. Fossil materials of E. ovodovi have thus far only been found in Russia. In this study, we extracted DNA from three equid fossil specimens excavated from northeastern China dated at 12,770-12,596, 29,525-28,887 and 40,201-38,848 cal. yBP, respectively, and retrieved three near-complete mitochondrial genomes from the specimens. Phylogenetic analyses cluster the Chinese haplotypes together with previously published Russian E. ovodovi, strongly supporting the assignment of these samples to this taxon. The molecular identification of E. ovodovi in northeastern China extends the known geographical range of this fossil species by several thousand kilometers to the east. The estimated coalescence time of all E. ovodovi haplotypes is approximately 199 Kya, with the Chinese haplotypes coalescing approximately 130 Kya. With a radiocarbon age of 12,770-12,596 cal. yBP, the youngest sample in this study represents the first E. ovodovi sample dating to the terminal Pleistocene, moving the extinction date of this species forwards considerably compared to previously documented fossils. Overall, comparison of our three mitochondrial genomes with the two published ones suggests a genetic diversity similar to several extant species of the genus Equus.

Paleogenome Reveals Genetic Contribution of Extinct Giant Panda to Extant Populations.

Historically, the giant panda was widely distributed from northern China to southwestern Asia [1]. As a result of range contraction and fragmentation, extant individuals are currently restricted to fragmented mountain ranges on the eastern margin of the Qinghai-Tibet plateau, where they are distributed among three major population clusters [2]. However, little is known about the genetic consequences of this dramatic range contraction. For example, were regions where giant pandas previously existed occupied by ancestors of present-day populations, or were these regions occupied by genetically distinct populations that are now extinct? If so, is there any contribution of these extinct populations to the genomes of giant pandas living today? To investigate these questions, we sequenced the nuclear genome of an ∼5,000-year-old giant panda from Jiangdongshan, Tengchong County in Yunnan Province, China. We find that this individual represents a genetically distinct population that diverged prior to the diversification of modern giant panda populations. We find evidence of differential admixture with this ancient population among modern individuals originating from different populations as well as within the same population. We also find evidence for directional gene flow, which transferred alleles from the ancient population into the modern giant panda lineages. A variable proportion of the genomes of extant individuals is therefore likely derived from the ancient population represented by our sequenced individual. Although extant giant panda populations retain reasonable genetic diversity, our results suggest that this represents only part of the genetic diversity this species harbored prior to its recent range contractions.

Optimized DNA sampling of ancient bones using Computed Tomography scans.

The prevalence of contaminant microbial DNA in ancient bone samples represents the principal limiting factor for palaeogenomic studies, as it may comprise more than 99% of DNA molecules obtained. Efforts to exclude or reduce this contaminant fraction have been numerous but also variable in their success. Here, we present a simple but highly effective method to increase the relative proportion of endogenous molecules obtained from ancient bones. Using computed tomography (CT) scanning, we identify the densest region of a bone as optimal for sampling. This approach accurately identifies the densest internal regions of petrous bones, which are known to be a source of high-purity ancient DNA. For ancient long bones, CT scans reveal a high-density outermost layer, which has been routinely removed and discarded prior to DNA extraction. For almost all long bones investigated, we find that targeted sampling of this outermost layer provides an increase in endogenous DNA content over that obtained from softer, trabecular bone. This targeted sampling can produce as much as 50-fold increase in the proportion of endogenous DNA, providing a directly proportional reduction in sequencing costs for shotgun sequencing experiments. The observed increases in endogenous DNA proportion are not associated with any reduction in absolute endogenous molecule recovery. Although sampling the outermost layer can result in higher levels of human contamination, some bones were found to have more contamination associated with the internal bone structures. Our method is highly consistent, reproducible and applicable across a wide range of bone types, ages and species. We predict that this discovery will greatly extend the potential to study ancient populations and species in the genomics era.