Morten Rasmussen

Ancient human genome sequence of an extinct Palaeo-Eskimo

We report here the genome sequence of an ancient human. Obtained from ∼4,000-year-old permafrost-preserved hair, the genome represents a male individual from the first known culture to settle in Greenland. Sequenced to an average depth of 20×, we recover 79% of the diploid genome, an amount close to the practical limit of current sequencing technologies. We identify 353,151 high-confidence single-nucleotide polymorphisms (SNPs), of which 6.8% have not been reported previously. We estimate raw read contamination to be no higher than 0.8%. We use functional SNP assessment to assign possible phenotypic characteristics of the individual that belonged to a culture whose location has yielded only trace human remains. We compare the high-confidence SNPs to those of contemporary populations to find the populations most closely related to the individual. This provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of that giving rise to the modern Native Americans and Inuit.

Monitoring endangered freshwater biodiversity using environmental DNA

Freshwater ecosystems are among the most endangered habitats on Earth, with thousands of animal species known to be threatened or already extinct. Reliable monitoring of threatened organisms is crucial for data-driven conservation actions but remains a challenge owing to nonstandardized methods that depend on practical and taxonomic expertise, which is rapidly declining. Here, we show that a diversity of rare and threatened freshwater animals--representing amphibians, fish, mammals, insects and crustaceans--can be detected and quantified based on DNA obtained directly from small water samples of lakes, ponds and streams. We successfully validate our findings in a controlled mesocosm experiment and show that DNA becomes undetectable within 2u2003weeks after removal of animals, indicating that DNA traces are near contemporary with presence of the species. We further demonstrate that entire faunas of amphibians and fish can be detected by high-throughput sequencing of DNA extracted from pond water. Our findings underpin the ubiquitous nature of DNA traces in the environment and establish environmental DNA as a tool for monitoring rare and threatened species across a wide range of taxonomic groups.

Detection of a Diverse Marine Fish Fauna Using Environmental DNA from Seawater Samples

Marine ecosystems worldwide are under threat with many fish species and populations suffering from human over-exploitation. This is greatly impacting global biodiversity, economy and human health. Intriguingly, marine fish are largely surveyed using selective and invasive methods, which are mostly limited to commercial species, and restricted to particular areas with favourable conditions. Furthermore, misidentification of species represents a major problem. Here, we investigate the potential of using metabarcoding of environmental DNA (eDNA) obtained directly from seawater samples to account for marine fish biodiversity. This eDNA approach has recently been used successfully in freshwater environments, but never in marine settings. We isolate eDNA from ½-litre seawater samples collected in a temperate marine ecosystem in Denmark. Using next-generation DNA sequencing of PCR amplicons, we obtain eDNA from 15 different fish species, including both important consumption species, as well as species rarely or never recorded by conventional monitoring. We also detect eDNA from a rare vagrant species in the area; European pilchard (Sardina pilchardus). Additionally, we detect four bird species. Records in national databases confirmed the occurrence of all detected species. To investigate the efficiency of the eDNA approach, we compared its performance with 9 methods conventionally used in marine fish surveys. Promisingly, eDNA covered the fish diversity better than or equal to any of the applied conventional methods. Our study demonstrates that even small samples of seawater contain eDNA from a wide range of local fish species. Finally, in order to examine the potential dispersal of eDNA in oceans, we performed an experiment addressing eDNA degradation in seawater, which shows that even small (100-bp) eDNA fragments degrades beyond detectability within days. Although further studies are needed to validate the eDNA approach in varying environmental conditions, our findings provide a strong proof-of-concept with great perspectives for future monitoring of marine biodiversity and resources.

Glacial Survival of Boreal Trees in Northern Scandinavia

Tree Refugia Ideas of how and when boreal plants spread to the formerly glaciated parts of the world following the retreat of the glaciers 9000 years ago are long debated. Models of the postglacial spread of boreal plants argue for dispersal from southern refugia; however, Parducci et al. (p. 1083) have shown that both spruce and pine were present in small ice-free regions of Scandinavia much earlier than thought. DNA haplotyping confirmed that a remnant mitochondrial type of spruce, once unique to Scandinavia, now lives alongside the more common spruce originating from Eastern Europe. Evidence from lake cores collected from central and northern Norway indicated the survival of conifers as early as 22,000 years before the present, when apart from ice-free pockets, most of Scandinavia was covered by ice. DNA from modern and ancient spruce and pine indicate that both survived in ice-free areas during the last glaciations. It is commonly believed that trees were absent in Scandinavia during the last glaciation and first recolonized the Scandinavian Peninsula with the retreat of its ice sheet some 9000 years ago. Here, we show the presence of a rare mitochondrial DNA haplotype of spruce that appears unique to Scandinavia and with its highest frequency to the west—an area believed to sustain ice-free refugia during most of the last ice age. We further show the survival of DNA from this haplotype in lake sediments and pollen of Trøndelag in central Norway dating back ~10,300 years and chloroplast DNA of pine and spruce in lake sediments adjacent to the ice-free Andøya refugium in northwestern Norway as early as ~22,000 and 17,700 years ago, respectively. Our findings imply that conifer trees survived in ice-free refugia of Scandinavia during the last glaciation, challenging current views on survival and spread of trees as a response to climate changes.

Pulling out the 1%: whole-Genome capture for the targeted enrichment of ancient dna sequencing libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062-147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217-73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.

The genetic prehistory of the New World Arctic

Introduction Humans first peopled the North American Arctic (northern Alaska, Canada, and Greenland) around 6000 years ago, leaving behind a complex archaeological record that consisted of different cultural units and distinct ways of life, including the Early Paleo-Eskimos (Pre-Dorset/Saqqaq), the Late Paleo-Eskimos (Early Dorset, Middle Dorset, and Late Dorset), and the Thule cultures. Genetic origins of Paleo-Eskimos and Neo-Eskimos. All Paleo-Eskimos represent a single migration pulse from Siberia into the Americas, independent of the Neo-Eskimo Thule people (ancestors of modern-day Inuit) and the related extinct Sadlermiut population. The Siberian Birnirk people were likely cultural and genetic ancestors of modern-day Inuit. We also show ancient admixture between the Paleo- and Neo-Eskimo lineages, occurring at least 4000 years ago. Rationale We addressed the genetic origins and relationships of the various New World Arctic cultures to each other and to modern-day populations in the region. We obtained 26 genome-wide sequences and 169 mitochondrial DNA sequences from ancient human bone, teeth, and hair samples from Arctic Siberia, Alaska, Canada, and Greenland, and high-coverage genomes of two present-day Greenlandic Inuit, two Siberian Nivkhs, one Aleutian Islander, and two Athabascan Native Americans. Twenty-seven ancient samples were radiocarbon dated for accurate cultural assignment, of which 25 were corrected for marine reservoir effect to account for the dominant marine component in these individuals’ diets. Results Nuclear and mitochondrial DNA data unequivocally show that the Paleo-Eskimos are closer to each other than to any other present-day population. The Thule culture represents a distinct people that are genetic and cultural ancestors of modern-day Inuit. We additionally find the Siberian Birnirk culture (6th to 7th century CE) as likely cultural and genetic ancestors of the Thule. The extinct Sadlermiut people from the Hudson Bay region (15th to 19th century CE), considered to be Dorset remnants, are genetically closely related to Thule/Inuit, rather than the Paleo-Eskimos. Moreover, there is no evidence of matrilineal gene flow between Dorset or Thule groups with neighboring Norse (Vikings) populations settling in the Arctic around 1000 years ago. However, we do detect gene flow between the Paleo-Eskimo and Neo-Eskimo lineages, dating back to at least 4000 years. Conclusion Our study has a number of important implications: Paleo-Eskimos likely represent a single migration pulse into the Americas from Siberia, separate from the ones giving rise to the Inuit and other Native Americans, including Athabascan speakers. Paleo-Eskimos, despite showing cultural differences across time and space, constituted a single population displaying genetic continuity for more than 4000 years. On the contrary, the Thule people, ancestors of contemporary Inuit, represent a population replacement of the Paleo-Eskimos that occurred less than 700 years ago. The long-term genetic continuity of the Paleo-Eskimo gene pool and lack of evidence of Native American admixture suggest that the Saqqaq and Dorset people were largely living in genetic isolation after entering the New World. Thus, the Paleo-Eskimo technological innovations and changes through time, as evident from the archaeological record, seem to have occurred solely by movement of ideas within a single resident population. This suggests that cultural similarities and differences are not solid proxies for population movements and migrations into new and dramatically different environments, as is often assumed. Arctic genetics comes in from the cold Despite a well-characterized archaeological record, the genetics of the people who inhabit the Arctic have been unexplored. Raghavan et al. sequenced ancient and modern genomes of individuals from the North American Arctic (see the Perspective by Park). Analyses of these genomes indicate that the Arctic was colonized 6000 years ago by a migration separate from the one that gave rise to other Native American populations. Furthermore, the original paleo-inhabitants of the Arctic appear to have been completely replaced approximately 700 years ago. Science, this issue 10.1126/science.1255832; see also p. 1004 Early Arctic humans differed from both present-day Inuit and Native Americans. [Also see Perspective by Park] The New World Arctic, the last region of the Americas to be populated by humans, has a relatively well-researched archaeology, but an understanding of its genetic history is lacking. We present genome-wide sequence data from ancient and present-day humans from Greenland, Arctic Canada, Alaska, Aleutian Islands, and Siberia. We show that Paleo-Eskimos (~3000 BCE to 1300 CE) represent a migration pulse into the Americas independent of both Native American and Inuit expansions. Furthermore, the genetic continuity characterizing the Paleo-Eskimo period was interrupted by the arrival of a new population, representing the ancestors of present-day Inuit, with evidence of past gene flow between these lineages. Despite periodic abandonment of major Arctic regions, a single Paleo-Eskimo metapopulation likely survived in near-isolation for more than 4000 years, only to vanish around 700 years ago.

mapDamage: testing for damage patterns in ancient DNA sequences

SUMMARYnAncient DNA extracts consist of a mixture of contaminant DNA molecules, most often originating from environmental microbes, and endogenous fragments exhibiting substantial levels of DNA damage. The latter introduce specific nucleotide misincorporations and DNA fragmentation signatures in sequencing reads that could be advantageously used to argue for sequence validity. mapDamage is a Perl script that computes nucleotide misincorporation and fragmentation patterns using next-generation sequencing reads mapped against a reference genome. The Perl script outputs are further automatically processed in embedded R script in order to detect typical patterns of genuine ancient DNA sequences.nnnAVAILABILITY AND IMPLEMENTATIONnThe Perl script mapDamage is freely available with documentation and example files at http://geogenetics.ku.dk/all_literature/mapdamage/. The script requires prior installation of the SAMtools suite and R environment and has been validated on both GNU/Linux and MacOSX operating systems.

Molecular Diet Analysis of Two African Free-Tailed Bats (Molossidae) Using High Throughput Sequencing

Given the diversity of prey consumed by insectivorous bats, it is difficult to discern the composition of their diet using morphological or conventional PCR-based analyses of their faeces. We demonstrate the use of a powerful alternate tool, the use of the Roche FLX sequencing platform to deep-sequence uniquely 5′ tagged insect-generic barcode cytochrome c oxidase I (COI) fragments, that were PCR amplified from faecal pellets of two free-tailed bat species Chaerephon pumilus and Mops condylurus (family: Molossidae). Although the analyses were challenged by the paucity of southern African insect COI sequences in the GenBank and BOLD databases, similarity to existing collections allowed the preliminary identification of 25 prey families from six orders of insects within the diet of C. pumilus, and 24 families from seven orders within the diet of M. condylurus. Insects identified to families within the orders Lepidoptera and Diptera were widely present among the faecal samples analysed. The two families that were observed most frequently were Noctuidae and Nymphalidae (Lepidoptera). Species-level analysis of the data was accomplished using novel bioinformatics techniques for the identification of molecular operational taxonomic units (MOTU). Based on these analyses, our data provide little evidence of resource partitioning between sympatric M. condylurus and C. pumilus in the Simunye region of Swaziland at the time of year when the samples were collected, although as more complete databases against which to compare the sequences are generated this may have to be re-evaluated.

High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species

Sympatric cryptic species, characterized by low morphological differentiation, pose a challenge to understanding the role of interspecific competition in structuring ecological communities. We used traditional (morphological) and novel molecular methods of diet analysis to study the diet of two cryptic bat species that are sympatric in southern England (Plecotus austriacus and P. auritus) (Fig. 1). Using Roche FLX 454 (Roche, Basel, CH) high-throughput sequencing (HTS) and uniquely tagged generic arthropod primers, we identified 142 prey Molecular Operational Taxonomic Units (MOTUs) in the diet of the cryptic bats, 60% of which were assigned to a likely species or genus. The findings from the molecular study supported the results of microscopic analyses in showing that the diets of both species were dominated by lepidopterans. However, HTS provided a sufficiently high resolution of prey identification to determine fine-scale differences in resource use. Although both bat species appeared to have a generalist diet, eared-moths from the family Noctuidae were the main prey consumed. Interspecific niche overlap was greater than expected by chance (Ojk = 0.72, P < 0.001) due to overlap in the consumption of the more common prey species. Yet, habitat associations of nongeneralist prey species found in the diets corresponded to those of their respective bat predator (grasslands for P. austriacus, and woodland for P. auritus). Overlap in common dietary resource use combined with differential specialist prey habitat associations suggests that habitat partitioning is the primary mechanism of coexistence. The performance of HTS is discussed in relation to previous methods of molecular and morphological diet analysis. By enabling species-level identification of dietary components, the application of DNA sequencing to diet analysis allows a more comprehensive comparison of the diet of sympatric cryptic species, and therefore can be an important tool for determining fine-scale mechanisms of coexistence.

A comparative study of ancient sedimentary DNA, pollen and macrofossils from permafrost sediments of northern Siberia reveals long-term vegetational stability

Although ancient DNA from sediments (sedaDNA) has been used to investigate past ecosystems, the approach has never been directly compared with the traditional methods of pollen and macrofossil analysis. We conducted a comparative survey of 18 ancient permafrost samples spanning the Late Pleistocene (46–12.5u2003thousand years ago), from the Taymyr Peninsula in northern Siberia. The results show that pollen, macrofossils and sedaDNA are complementary rather than overlapping and, in combination, reveal more detailed information on plant palaeocommunities than can be achieved by each individual approach. SedaDNA and macrofossils share greater overlap in plant identifications than with pollen, suggesting that sedaDNA is local in origin. These two proxies also permit identification to lower taxonomic levels than pollen, enabling investigation into temporal changes in species composition and the determination of indicator species to describe environmental changes. Combining data from all three proxies reveals an area continually dominated by a mosaic vegetation of tundra‐steppe, pioneer and wet‐indicator plants. Such vegetational stability is unexpected, given the severe climate changes taking place in the Northern Hemisphere during this time, with changes in average annual temperatures of >22u2003°C. This may explain the abundance of ice‐age mammals such as horse and bison in Taymyr Peninsula during the Pleistocene and why it acted as a refugium for the last mainland woolly mammoth. Our finding reveals the benefits of combining sedaDNA, pollen and macrofossil for palaeovegetational reconstruction and adds to the increasing evidence suggesting large areas of the Northern Hemisphere remained ecologically stable during the Late Pleistocene.