Ancient DNA and osteological analyses of a unique paleo-archive reveal Early Holocene faunal expansion into the Scandinavian Arctic.

Paleo-archives are essential for our understanding of species responses to climate warming, yet such archives are extremely rare in the Arctic. Here, we combine morphological analyses and bulk-bone metabarcoding to investigate a unique chronology of bone deposits sealed in the high-latitude Storsteinhola cave system (68°50' N 16°22' E) in Norway. This deposit dates to a period of climate warming from the end of the Late Glacial [~13 thousand calibrated years before the present (ka cal B.P.)] to the Holocene thermal maximum (~5.6 ka cal B.P.). Paleogenetic analyses allow us to exploit the 1000s of morphologically unidentifiable bone fragments resulting in a high-resolution sequence with 40 different taxa, including species not previously found here. Our record reveals borealization in both the marine and terrestrial environments above the Arctic Circle as a naturally recurring phenomenon in past periods of warming, providing fundamental insights into the ecosystem-wide responses that are ongoing today.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change.

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

Holocene deglaciation drove rapid genetic diversification of Atlantic walrus.

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

Sedimentary ancient DNA metabarcoding as a tool for assessing prehistoric plant use at the Upper Paleolithic cave site Aghitu-3, Armenia.

Current knowledge about Paleolithic human plant use is limited by the rare survival of identifiable plant remains as well as the availability of methods for plant detection and identification. By analyzing DNA preserved in cave sediments, we can identify organisms in the absence of any visible remains, opening up new ways to study details of past human behavior, including plant use. Aghitu-3 Cave contains a 15,000-yearlong record (from ∼39,000 to 24,000 cal BP) of Upper Paleolithic human settlement and environmental variability in the Armenian Highlands. Finds from this cave include stone artifacts, faunal remains, bone tools, shell beads, charcoal, and pollen, among others. We applied sedimentary ancient DNA (sedaDNA) metabarcoding to the Aghitu-3 sedimentary sequence and combined this with pollen data to obtain a temporal reconstruction of plant assemblages. Our results reveal a stratification of plant abundance and diversity where sedaDNA reflects periods of human occupation, showing higher diversity in layers with increased human activity. Low pollen concentrations combined with high sedaDNA abundance indicate plant remains may have been brought into the cave by animals or humans during the deposition of the lower two archaeological horizons. Most of the recovered plants are reported to be useful for food, flavor, medicine, and/or technical purposes, demonstrating the potential of the environment around Aghitu-3 Cave to support humans during the Upper Paleolithic. Moreover, we identified several specific plant taxa that strengthen previous findings about Upper Paleolithic plant use in this region (i.e., for medicine and the manufacturing and dyeing of textiles). This study represents the first application of plant sedaDNA analysis of cave sediments for the investigation of potential plant use by prehistoric humans.

Complex population structure of the Atlantic puffin revealed by whole genome analyses.

The factors underlying gene flow and genomic population structure in vagile seabirds are notoriously difficult to understand due to their complex ecology with diverse dispersal barriers and extensive periods at sea. Yet, such understanding is vital for conservation management of seabirds that are globally declining at alarming rates. Here, we elucidate the population structure of the Atlantic puffin (Fratercula arctica) by assembling its reference genome and analyzing genome-wide resequencing data of 72 individuals from 12 colonies. We identify four large, genetically distinct clusters, observe isolation-by-distance between colonies within these clusters, and obtain evidence for a secondary contact zone. These observations disagree with the current taxonomy, and show that a complex set of contemporary biotic factors impede gene flow over different spatial scales. Our results highlight the power of whole genome data to reveal unexpected population structure in vagile marine seabirds and its value for seabird taxonomy, evolution and conservation.

Metagenomics: A viable tool for reconstructing herbivore diet.

Metagenomics can generate data on the diet of herbivores, without the need for primer selection and PCR enrichment steps as is necessary in metabarcoding. Metagenomic approaches to diet analysis have remained relatively unexplored, requiring validation of bioinformatic steps. Currently, no metagenomic herbivore diet studies have utilized both chloroplast and nuclear markers as reference sequences for plant identification, which would increase the number of reads that could be taxonomically informative. Here, we explore how in silico simulation of metagenomic data sets resembling sequences obtained from faecal samples can be used to validate taxonomic assignment. Using a known list of sequences to create simulated data sets, we derived reliable identification parameters for taxonomic assignments of sequences. We applied these parameters to characterize the diet of western capercaillies (Tetrao urogallus) located in Norway, and compared the results with metabarcoding trnL P6 loop data generated from the same samples. Both methods performed similarly in the number of plant taxa identified (metagenomics 42 taxa, metabarcoding 43 taxa), with no significant difference in species resolution (metagenomics 24%, metabarcoding 23%). We further observed that while metagenomics was strongly affected by the age of faecal samples, with fresh samples outperforming old samples, metabarcoding was not affected by sample age. On the other hand, metagenomics allowed us to simultaneously obtain the mitochondrial genome of the western capercaillies, thereby providing additional ecological information. Our study demonstrates the potential of utilizing metagenomics for diet reconstruction but also highlights key considerations as compared to metabarcoding for future utilization of this technique.

Ancient genomes reveal early Andean farmers selected common beans while preserving diversity.

All crops are the product of a domestication process that started less than 12,000 years ago from one or more wild populations1,2. Farmers selected desirable phenotypic traits (such as improved energy accumulation, palatability of seeds and reduced natural shattering3) while leading domesticated populations through several more or less gradual demographic contractions2,4. As a consequence, the erosion of wild genetic variation5 is typical of modern cultivars, making them highly susceptible to pathogens, pests and environmental change6,7. The loss of genetic diversity hampers further crop improvement programmes to increase food production in a changing world, posing serious threats to food security8,9. Using both ancient and modern seeds, we analysed the temporal dynamics of genetic variation and selection during the domestication process of the common bean (Phaseolus vulgaris) in the southern Andes. Here, we show that most domestic traits were selected for before 2,500 years ago, with no or only minor loss of whole-genome heterozygosity. In fact, most of the changes at coding genes and linked regions that differentiate wild and domestic genomes are already present in the ancient genomes analysed here, and all ancient domestic genomes dated between 600 and 2,500 years ago are highly variable (at least as variable as modern genomes from the wild). Single seeds from modern cultivars show reduced variation when compared with ancient seeds, indicating that intensive selection within cultivars in the past few centuries probably partitioned ancestral variation within different genetically homogenous cultivars. When cultivars from different Andean regions are pooled, the genomic variation of the pool is higher than that observed in the pool of ancient seeds from north and central western Argentina. Considering that most desirable phenotypic traits are probably controlled by multiple polymorphic genes10, a plausible explanation of this decoupling of selection and genetic erosion is that early farmers applied a relatively weak selection pressure2 by using many phenotypically similar but genetically diverse individuals as parents. Our results imply that selection strategies during the past few centuries, as compared with earlier times, more intensively reduced genetic variation within cultivars and produced further improvements by focusing on a few plants carrying the traits of interest, at the cost of marked genetic erosion within Andean landraces.

A molecular phylogeny of historical and contemporary specimens of an under-studied micro-invertebrate group.

Resolution of relationships at lower taxonomic levels is crucial for answering many evolutionary questions, and as such, sufficiently varied species representation is vital. This latter goal is not always achievable with relatively fresh samples. To alleviate the difficulties in procuring rarer taxa, we have seen increasing utilization of historical specimens in building molecular phylogenies using high throughput sequencing. This effort, however, has mainly focused on large-bodied or well-studied groups, with small-bodied and under-studied taxa under-prioritized. Here, we utilize both historical and contemporary specimens, to increase the resolution of phylogenetic relationships among a group of under-studied and small-bodied metazoans, namely, cheilostome bryozoans. In this study, we pioneer the sequencing of air-dried cheilostomes, utilizing a recently developed library preparation method for low DNA input. We evaluate a de novo mitogenome assembly and two iterative methods, using the sequenced target specimen as a reference for mapping, for our sequences. In doing so, we present mitochondrial and ribosomal RNA sequences of 43 cheilostomes representing 37 species, including 14 from historical samples ranging from 50 to 149 years old. The inferred phylogenetic relationships of these samples, analyzed together with publicly available sequence data, are shown in a statistically well-supported 65 taxa and 17 genes cheilostome tree, which is also the most broadly sampled and largest to date. The robust phylogenetic placement of historical samples whose contemporary conspecifics and/or congenerics have been sequenced verifies the appropriateness of our workflow and gives confidence in the phylogenetic placement of those historical samples for which there are no close relatives sequenced. The success of our workflow is highlighted by the circularization of a total of 27 mitogenomes, seven from historical cheilostome samples. Our study highlights the potential of utilizing DNA from micro-invertebrate specimens stored in natural history collections for resolving phylogenetic relationships among species.

Not a limitless resource: ethics and guidelines for destructive sampling of archaeofaunal remains.

With the advent of ancient DNA, as well as other methods such as isotope analysis, destructive sampling of archaeofaunal remains has increased much faster than the effort to collect and curate them. While there has been considerable discussion regarding the ethics of destructive sampling and analysis of human remains, this dialogue has not extended to archaeofaunal material. Here we address this gap and discuss the ethical challenges surrounding destructive sampling of materials from archaeofaunal collections. We suggest ways of mitigating the negative aspects of destructive sampling and present step-by-step guidelines aimed at relevant stakeholders, including scientists, holding institutions and scientific journals. Our suggestions are in most cases easily implemented without significant increases in project costs, but with clear long-term benefits in the preservation and use of zooarchaeological material.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Ancient DNA reveals the chronology of walrus ivory trade from Norse Greenland.

The importance of the Atlantic walrus ivory trade for the colonization, peak, and collapse of the medieval Norse colonies on Greenland has been extensively debated. Nevertheless, no studies have directly traced medieval European ivory back to distinct Arctic populations of walrus. Analysing the entire mitogenomes of 37 archaeological specimens from Europe, Svalbard, and Greenland, we here discover that Atlantic walrus comprises two monophyletic mitochondrial (MT) clades, which diverged between 23 400 and 251 120 years ago. Our improved genomic resolution allows us to reinterpret the geographical distribution of partial MT data from 306 modern and nineteenth-century specimens, finding that one of these clades was exclusively accessible to Greenlanders. With this discovery, we ascertain the biological origin of 23 archaeological specimens from Europe (most dated between 900 and 1400 CE). These results reveal a significant shift in trade from an early, predominantly eastern source towards a near exclusive representation of Greenland ivory. Our study provides empirical evidence for how this remote Arctic resource was progressively integrated into a medieval pan-European trade network, contributing to both the resilience and vulnerability of Norse Greenland society.

A stainless-steel mortar, pestle and sleeve design for the efficient fragmentation of ancient bone.

Different types of milling equipment - such as oscillating ball mills, freezer mills, mortar and pestle - can be used to fragment ancient bone prior to DNA extraction. However, each of these tools is associated with practical drawbacks. Here, we present the design for a stainless-steel mortar and pestle, with a removable sleeve to contain bone material. The tool is easy to clean, practical and its simplicity allows university workshops equipped with a lathe, boring tools and a milling machine to make these components at local expense. This design allows for the efficient fragmentation of ancient bone and improves sample throughput. This design is recommended as a useful, economical addition to existing laboratory equipment for the handling of ancient bone.

Ancient DNA reveals the Arctic origin of Viking Age cod from Haithabu, Germany.

Knowledge of the range and chronology of historic trade and long-distance transport of natural resources is essential for determining the impacts of past human activities on marine environments. However, the specific biological sources of imported fauna are often difficult to identify, in particular if species have a wide spatial distribution and lack clear osteological or isotopic differentiation between populations. Here, we report that ancient fish-bone remains, despite being porous, brittle, and light, provide an excellent source of endogenous DNA (15-46%) of sufficient quality for whole-genome reconstruction. By comparing ancient sequence data to that of modern specimens, we determine the biological origin of 15 Viking Age (800-1066 CE) and subsequent medieval (1066-1280 CE) Atlantic cod (Gadus morhua) specimens from excavation sites in Germany, Norway, and the United Kingdom. Archaeological context indicates that one of these sites was a fishing settlement for the procurement of local catches, whereas the other localities were centers of trade. Fish from the trade sites show a mixed ancestry and are statistically differentiated from local fish populations. Moreover, Viking Age samples from Haithabu, Germany, are traced back to the North East Arctic Atlantic cod population that has supported the Lofoten fisheries of Norway for centuries. Our results resolve a long-standing controversial hypothesis and indicate that the marine resources of the North Atlantic Ocean were used to sustain an international demand for protein as far back as the Viking Age.

Combining bleach and mild predigestion improves ancient DNA recovery from bones.

The feasibility of genome-scale studies from archaeological material remains critically dependent on the ability to access endogenous, authentic DNA. In the majority of cases, this represents a few per cent of the DNA extract, at most. A number of specific pre-extraction protocols for bone powder aimed to improve ancient DNA recovery before library amplification have recently been developed. Here, we test the effects of combining two of such protocols, a bleach wash and a predigestion step, on 12 bone samples of Atlantic cod and domestic horse aged 750-1350 cal. years before present. Using high-throughput sequencing, we show that combined together, bleach wash and predigestion consistently yield DNA libraries with higher endogenous content than either of these methods alone. Additionally, the molecular complexity of these libraries is improved and endogenous DNA templates show larger size distributions. Other library characteristics, such as DNA damage profiles or the composition of microbial communities, are little affected by the pre-extraction protocols. Application of the combined protocol presented in this study will facilitate the genetic analysis of an increasing number of ancient remains and will reduce the cost of whole-genome sequencing.

DNA Metabarcoding Reveals Diet Overlap between the Endangered Walia Ibex and Domestic Goats - Implications for Conservation.

Human population expansion and associated degradation of the habitat of many wildlife species cause loss of biodiversity and species extinctions. The small Simen Mountains National Park in Ethiopia is one of the last strongholds for the preservation of a number of afro-alpine mammals, plants and birds, and it is home to the rare endemic Walia ibex, Capra walie. The narrow distribution range of this species as well as potential competition for resources with livestock, especially with domestic goat, Capra hircus, may compromise its future survival. Based on a curated afro-alpine taxonomic reference library constructed for plant taxon identification, we investigated the diet of the Walia ibex and addressed the dietary overlap with domestic goat using DNA metabarcoding of faecal samples. Faeces of both species were collected from different localities in the National Park. We show that both species are browsers, with forbs, shrubs and trees comprising the largest proportion of their diet, supplemented by grasses. There was a considerable overlap in dietary preferences. Several of the preferred diet items of the Walia ibex (Alchemilla sp., Hypericum revolutum, Erica arborea and Rumex sp.) were also among the most preferred diet items of the domestic goat. These results indicate that there is potential for competition between the two species, especially during the dry season, when resources are limited. Our findings, in combination with the expected increase in domestic herbivores, suggest that management plans should consider the potential threat posed by domestic goats to ensure future survival of the endangered Walia ibex.

Genomic signatures of the plateless phenotype in the threespine stickleback.

Understanding the genetic basis of traits involved in adaptive divergence and speciation is one of the most fundamental objectives in evolutionary biology. Toward that end, we look for signatures of extreme plate loss in the genome of freshwater threespine sticklebacks (Gasterosteus aculeatus). Plateless stickleback have been found in only a few lakes and streams across the world; they represent the far extreme of a phenotypic continuum (plate number) that has been studied for years, although plateless individuals have not yet been the subject of much investigation. We use a dense single nucleotide polymorphism dataset made using RADseq to study fish from three freshwater populations containing plateless and low plated individuals, as well as fish from full plated marine populations. Analyses were performed using FastStructure, sliding windows F ST, Bayescan and latent factor mixed models to search for genomic differences between the low plated and plateless phenotypes both within and among the three lakes. At least 18 genomic regions which may contribute to within-morph plate number variation were detected in our low plated stickleback populations. We see no evidence of a selective sweep between low and plateless fish; rather reduction of plate number within the low plated morph seems to be polygenic.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Use of ancient sedimentary DNA as a novel conservation tool for high-altitude tropical biodiversity.

Conservation of biodiversity may in the future increasingly depend upon the availability of scientific information to set suitable restoration targets. In traditional paleoecology, sediment-based pollen provides a means to define preanthropogenic impact conditions, but problems in establishing the exact provenance and ecologically meaningful levels of taxonomic resolution of the evidence are limiting. We explored the extent to which the use of sedimentary ancient DNA (sedaDNA) may complement pollen data in reconstructing past alpine environments in the tropics. We constructed a record of afro-alpine plants retrieved from DNA preserved in sediment cores from 2 volcanic crater sites in the Albertine Rift, eastern Africa. The record extended well beyond the onset of substantial anthropogenic effects on tropical mountains. To ensure high-quality taxonomic inference from the sedaDNA sequences, we built an extensive DNA reference library covering the majority of the afro-alpine flora, by sequencing DNA from taxonomically verified specimens. Comparisons with pollen records from the same sediment cores showed that plant diversity recovered with sedaDNA improved vegetation reconstructions based on pollen records by revealing both additional taxa and providing increased taxonomic resolution. Furthermore, combining the 2 measures assisted in distinguishing vegetation change at different geographic scales; sedaDNA almost exclusively reflects local vegetation, whereas pollen can potentially originate from a wide area that in highlands in particular can span several ecozones. Our results suggest that sedaDNA may provide information on restoration targets and the nature and magnitude of human-induced environmental changes, including in high conservation priority, biodiversity hotspots, where understanding of preanthropogenic impact (or reference) conditions is highly limited.

Fungal palaeodiversity revealed using high-throughput metabarcoding of ancient DNA from arctic permafrost.

The taxonomic and ecological diversity of ancient fungal communities was assessed by combining next generation sequencing and metabarcoding of DNA preserved in permafrost. Twenty-six sediment samples dated 16 000-32 000 radiocarbon years old from two localities in Siberia were analysed for fungal ITS. We detected 75 fungal OTUs from 21 orders representing three phyla, although rarefaction analyses suggested that the full diversity was not recovered despite generating an average of 6677 ± 3811 (mean ± SD) sequences per sample and that preservation bias likely has considerable effect on the recovered DNA. Most OTUs (75.4%) represented ascomycetes. Due to insufficient sequencing depth, DNA degradation and putative preservation biases in our samples, the recovered taxa probably do not represent the complete historic fungal community, and it is difficult to determine whether the fungal communities varied geographically or experienced a composition shift within the period of 16 000-32 000 bp. However, annotation of OTUs to functional ecological groups provided a wealth of information on the historic communities. About one-third of the OTUs are presumed plant-associates (pathogens, saprotrophs and endophytes) typical of graminoid- and forb-rich habitats. We also detected putative insect pathogens, coprophiles and keratinophiles likely associated with ancient insect and herbivore faunas. The detection of putative insect pathogens, mycoparasites, aquatic fungi and endophytes broadens our previous knowledge of the diversity of fungi present in Beringian palaeoecosystems. A large group of putatively psychrophilic/psychrotolerant fungi was also detected, most likely representing a modern, metabolically active fungal community.