Statistical guidelines for detecting past population shifts using ancient DNA

Ancient Sheep Genomes Reveal Four Millennia of North European Short-Tailed Sheep in the Baltic Sea Region

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

Evolutionary history of the extinct wolf population from France in the context of global phylogeographic changes throughout the Holocene

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

Influence of Holocene habitat availability on Pacific gray whale (Eschrichtius robustus) population dynamics as inferred from whole mitochondrial genome sequences and environmental niche modeling

Environmental changes since the Pleistocene and commercial whaling in the last few centuries have drastically reduced many whale populations, including gray whales in the North Pacific. Herein we use complete mitogenome sequences from 74 individuals to evaluate gray whale phylogeography and historical demography, then use environmental niche modeling to assess how habitat availability has changed through time for Pacific gray whales. We identify a large degree of haplotype sharing between gray whales sampled in Russian and Mexican waters, coupled with very limited matrilineal population structure. Confirming previous studies, our environmental niche models showed a decrease in available habitat during the Last Glacial Maximum, but we find no genetic signals of recent population declines in mitochondrial genomes despite both sustained habitat loss and a commercial whaling bottleneck. Our results illustrate the complex dynamics of baleen whale biogeography since the Holocene as well as the difficulty in detecting recent demographic bottlenecks from mitochondrial DNA sequences.

Natural history collections recapitulate 200 years of faunal change

Changing species assemblages represent major challenges to ecosystems around the world. Retracing these changes is limited by our knowledge of past biodiversity. Natural history collections represent archives of biodiversity and are therefore an unparalleled source to study biodiversity changes. In the present study, we tested the value of natural history collections for reconstructing changes in the abundance and presence of species over time. In total, we scrutinized 17 080 quality-checked records for 242 epibenthic invertebrate species from the North and Baltic Seas collected throughout the last 200 years. Our approaches identified eight previously reported species introductions, 10 range expansions, six of which are new to science, as well as the long-term decline of 51 marine invertebrate species. The cross-validation of our results with published accounts of endangered species and neozoa of the area confirmed the results for two of the approaches for 49 to 55% of the identified species, and contradicted our results for 9 to 10%. The results based on relative record trends were less validated. We conclude that, with the proper approaches, natural history collections are an unmatched resource for recovering early species introductions and declines.

A Multidisciplinary Review of the Inka Imperial Resettlement Policy and Implications for Future Investigations

The rulers of the Inka empire conquered approximately 2 million km2 of the South American Andes in just under 100 years from 1438-1533 CE. Inside the empire, the elite conducted a systematic resettlement of the many Indigenous peoples in the Andes that had been rapidly colonised. The nature of this resettlement phenomenon is recorded within the Spanish colonial ethnohistorical record. Here we have broadly characterised the resettlement policy, despite the often incomplete and conflicting details in the descriptions. We then review research from multiple disciplines that investigate the empirical reality of the Inka resettlement policy, including stable isotope analysis, intentional cranial deformation morphology, ceramic artefact chemical analyses and genetics. Further, we discuss the benefits and limitations of each discipline for investigating the resettlement policy and emphasise their collective value in an interdisciplinary characterisation of the resettlement policy.

Paleogenomics: reconstruction of plant evolutionary trajectories from modern and ancient DNA

How contemporary plant genomes originated and evolved is a fascinating question. One approach uses reference genomes from extant species to reconstruct the sequence and structure of their common ancestors over deep timescales. A second approach focuses on the direct identification of genomic changes at a shorter timescale by sequencing ancient DNA preserved in subfossil remains. Merged within the nascent field of paleogenomics, these complementary approaches provide insights into the evolutionary forces that shaped the organization and regulation of modern genomes and open novel perspectives in fostering genetic gain in breeding programs and establishing tools to predict future population changes in response to anthropogenic pressure and global warming.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

Human evolution: a tale from ancient genomes

The field of human ancient DNA (aDNA) has moved from mitochondrial sequencing that suffered from contamination and provided limited biological insights, to become a fully genomic discipline that is changing our conception of human history. Recent successes include the sequencing of extinct hominins, and true population genomic studies of Bronze Age populations. Among the emerging areas of aDNA research, the analysis of past epigenomes is set to provide more new insights into human adaptation and disease susceptibility through time. Starting as a mere curiosity, ancient human genetics has become a major player in the understanding of our evolutionary history.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Taming the Past: Ancient DNA and the Study of Animal Domestication

During the last decade, ancient DNA research has been revolutionized by the availability of increasingly powerful DNA sequencing and ancillary genomics technologies, giving rise to the new field of paleogenomics. In this review, we show how our understanding of the genetic basis of animal domestication and the origins and dispersal of livestock and companion animals during the Upper Paleolithic and Neolithic periods is being rapidly transformed through new scientific knowledge generated with paleogenomic methods. These techniques have been particularly informative in revealing high-resolution patterns of artificial and natural selection and evidence for significant admixture between early domestic animal populations and their wild congeners.

Native American Genomics and Population Histories

Studies of Native American genetic diversity and population history have been transformed over the last decade by important developments in anthropological genetics. During this time, researchers have adopted new DNA technologies and computational approaches for analyzing genomic data, and they have become increasingly sensitive to the views of research participants and communities. As new methods are applied to long-standing questions, and as more research is conducted in collaboration with indigenous communities, we are gaining new insights into the history and diversity of indigenous populations. This review discusses the recent methodological advances and genetic studies that have improved our understanding of Native American genomics and population histories. We synthesize current knowledge about Native American genomic variation and build a model of population history in the Americas. We also discuss the broader implications of this research for anthropology and related disciplines, and we highlight challenges and other considerations for future research.

Conservation archaeogenomics: ancient DNA and biodiversity in the Anthropocene

There is growing consensus that we have entered the Anthropocene, a geologic epoch characterized by human domination of the ecosystems of the Earth. With the future uncertain, we are faced with understanding how global biodiversity will respond to anthropogenic perturbations. The archaeological record provides perspective on human-environment relations through time and across space. Ancient DNA (aDNA) analyses of plant and animal remains from archaeological sites are particularly useful for understanding past human-environment interactions, which can help guide conservation decisions during the environmental changes of the Anthropocene. Here, we define the emerging field of conservation archaeogenomics, which integrates archaeological and genomic data to generate baselines or benchmarks for scientists, managers, and policy-makers by evaluating climatic and human impacts on past, present, and future biodiversity.

Ancient genomics

The past decade has witnessed a revolution in ancient DNA (aDNA) research. Although the field's focus was previously limited to mitochondrial DNA and a few nuclear markers, whole genome sequences from the deep past can now be retrieved. This breakthrough is tightly connected to the massive sequence throughput of next generation sequencing platforms and the ability to target short and degraded DNA molecules. Many ancient specimens previously unsuitable for DNA analyses because of extensive degradation can now successfully be used as source materials. Additionally, the analytical power obtained by increasing the number of sequence reads to billions effectively means that contamination issues that have haunted aDNA research for decades, particularly in human studies, can now be efficiently and confidently quantified. At present, whole genomes have been sequenced from ancient anatomically modern humans, archaic hominins, ancient pathogens and megafaunal species. Those have revealed important functional and phenotypic information, as well as unexpected adaptation, migration and admixture patterns. As such, the field of aDNA has entered the new era of genomics and has provided valuable information when testing specific hypotheses related to the past.

Implications of the circumpolar genetic structure of polar bears for their conservation in a rapidly warming Arctic

We provide an expansive analysis of polar bear (Ursus maritimus) circumpolar genetic variation during the last two decades of decline in their sea-ice habitat. We sought to evaluate whether their genetic diversity and structure have changed over this period of habitat decline, how their current genetic patterns compare with past patterns, and how genetic demography changed with ancient fluctuations in climate. Characterizing their circumpolar genetic structure using microsatellite data, we defined four clusters that largely correspond to current ecological and oceanographic factors: Eastern Polar Basin, Western Polar Basin, Canadian Archipelago and Southern Canada. We document evidence for recent (ca. last 1–3 generations) directional gene flow from Southern Canada and the Eastern Polar Basin towards the Canadian Archipelago, an area hypothesized to be a future refugium for polar bears as climate-induced habitat decline continues. Our data provide empirical evidence in support of this hypothesis. The direction of current gene flow differs from earlier patterns of gene flow in the Holocene. From analyses of mitochondrial DNA, the Canadian Archipelago cluster and the Barents Sea subpopulation within the Eastern Polar Basin cluster did not show signals of population expansion, suggesting these areas may have served also as past interglacial refugia. Mismatch analyses of mitochondrial DNA data from polar and the paraphyletic brown bear (U. arctos) uncovered offset signals in timing of population expansion between the two species, that are attributed to differential demographic responses to past climate cycling. Mitogenomic structure of polar bears was shallow and developed recently, in contrast to the multiple clades of brown bears. We found no genetic signatures of recent hybridization between the species in our large, circumpolar sample, suggesting that recently observed hybrids represent localized events. Documenting changes in subpopulation connectivity will allow polar nations to proactively adjust conservation actions to continuing decline in sea-ice habitat.

Major transitions in human evolution revisited: a tribute to ancient DNA

The origin and diversification of modern humans have been characterized by major evolutionary transitions and demographic changes. Patterns of genetic variation within modern populations can help with reconstructing this ∼200 thousand year-long population history. However, by combining this information with genomic data from ancient remains, one can now directly access our evolutionary past and reveal our population history in much greater detail. This review outlines the main recent achievements in ancient DNA research and illustrates how the field recently moved from the polymerase chain reaction (PCR) amplification of short mitochondrial fragments to whole-genome sequencing and thereby revisited our own history. Ancient DNA research has revealed the routes that our ancestors took when colonizing the planet, whom they admixed with, how they domesticated plant and animal species, how they genetically responded to changes in lifestyle, and also, which pathogens decimated their populations. These approaches promise to soon solve many pending controversies about our own origins that are indecipherable from modern patterns of genetic variation alone, and therefore provide an extremely powerful toolkit for a new generation of molecular anthropologists.

Using Ancient DNA to Understand Evolutionary and Ecological Processes

Ancient DNA provides a unique means to record genetic change through time and directly observe evolutionary and ecological processes. Although mostly based on mitochondrial DNA, the increasing availability of genomic sequences is leading to unprecedented levels of resolution. Temporal studies of population genetics have revealed dynamic patterns of change in many large vertebrates, featuring localized extinctions, migrations, and population bottlenecks. The pronounced climate cycles of the Late Pleistocene have played a key role, reducing the taxonomic and genetic diversity of many taxa and shaping modern populations. Importantly, the complex series of events revealed by ancient DNA data is seldom reflected in current biogeographic patterns. DNA preserved in ancient sediments and coprolites has been used to characterize a range of paleoenvironments and reconstruct functional relationships in paleoecological systems. In the near future, genome-level surveys of ancient populations will play an increasingly important role in revealing, calibrating, and testing evolutionary processes.

Investigating population history using temporal genetic differentiation

The rapid advance of sequencing technology, coupled with improvements in molecular methods for obtaining genetic data from ancient sources, holds the promise of producing a wealth of genomic data from time-separated individuals. However, the population-genetic properties of time-structured samples have not been extensively explored. Here, we consider the implications of temporal sampling for analyses of genetic differentiation and use a temporal coalescent framework to show that complex historical events such as size reductions, population replacements, and transient genetic barriers between populations leave a footprint of genetic differentiation that can be traced through history using temporal samples. Our results emphasize explicit consideration of the temporal structure when making inferences and indicate that genomic data from ancient individuals will greatly increase our ability to reconstruct population history.

Demographic history of a recent invasion of house mice on the isolated Island of Gough

Island populations provide natural laboratories for studying key contributors to evolutionary change, including natural selection, population size and the colonization of new environments. The demographic histories of island populations can be reconstructed from patterns of genetic diversity. House mice (Mus musculus) inhabit islands throughout the globe, making them an attractive system for studying island colonization from a genetic perspective. Gough Island, in the central South Atlantic Ocean, is one of the remotest islands in the world. House mice were introduced to Gough Island by sealers during the 19th century and display unusual phenotypes, including exceptionally large body size and carnivorous feeding behaviour. We describe genetic variation in Gough Island mice using mitochondrial sequences, nuclear sequences and microsatellites. Phylogenetic analysis of mitochondrial sequences suggested that Gough Island mice belong to Mus musculus domesticus, with the maternal lineage possibly originating in England or France. Cluster analyses of microsatellites revealed genetic membership for Gough Island mice in multiple coastal populations in Western Europe, suggesting admixed ancestry. Gough Island mice showed substantial reductions in mitochondrial and nuclear sequence variation and weak reductions in microsatellite diversity compared with Western European populations, consistent with a population bottleneck. Approximate Bayesian computation (ABC) estimated that mice recently colonized Gough Island (~100 years ago) and experienced a 98% reduction in population size followed by a rapid expansion. Our results indicate that the unusual phenotypes of Gough Island mice evolved rapidly, positioning these mice as useful models for understanding rapid phenotypic evolution.

Historical DNA reveals the demographic history of Atlantic cod (Gadus morhua) in medieval and early modern Iceland

Atlantic cod (Gadus morhua) vertebrae from archaeological sites were used to study the history of the Icelandic Atlantic cod population in the time period of 1500-1990. Specifically, we used coalescence modelling to estimate population size and fluctuations from the sequence diversity at the cytochrome b (cytb) and Pantophysin I (PanI) loci. The models are consistent with an expanding population during the warm medieval period, large historical effective population size (NE), a marked bottleneck event at 1400-1500 and a decrease in NE in early modern times. The model results are corroborated by the reduction of haplotype and nucleotide variation over time and pairwise population distance as a significant portion of nucleotide variation partitioned across the 1550 time mark. The mean age of the historical fished stock is high in medieval times with a truncation in age in early modern times. The population size crash coincides with a period of known cooling in the North Atlantic, and we conclude that the collapse may be related to climate or climate-induced ecosystem change.

Inferring population size changes with sequence and SNP data: lessons from human bottlenecks

Reconstructing historical variation of population size from sequence and single-nucleotide polymorphism (SNP) data is valuable for understanding the evolutionary history of species. Changes in the population size of humans have been thoroughly investigated, and we review different methodologies of demographic reconstruction, specifically focusing on human bottlenecks. In addition to the classical approaches based on the site-frequency spectrum (SFS) or based on linkage disequilibrium, we also review more recent approaches that utilize atypical shared genomic fragments, such as identical by descent or homozygous segments between or within individuals. Compared with methods based on the SFS, these methods are well suited for detecting recent bottlenecks. In general, all these various methods suffer from bias and dependencies on confounding factors such as population structure or poor specification of the mutational and recombination processes, which can affect the demographic reconstruction. With the exception of SFS-based methods, the effects of confounding factors on the inference methods remain poorly investigated. We conclude that an important step when investigating population size changes rests on validating the demographic model by investigating to what extent the fitted demographic model can reproduce the main features of the polymorphism data.

Inferring recent historic abundance from current genetic diversity

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Ancient DNA reveals genetic stability despite demographic decline: 3,000 years of population history in the endemic Hawaiian petrel

In the Hawaiian Islands, human colonization, which began approximately 1,200 to 800 years ago, marks the beginning of a period in which nearly 75% of the endemic avifauna became extinct and the population size and range of many additional species declined. It remains unclear why some species persisted whereas others did not. The endemic Hawaiian petrel (Pterodroma sandwichensis) has escaped extinction, but colonies on two islands have been extirpated and populations on remaining islands have contracted. We obtained mitochondrial DNA sequences from 100 subfossil bones, 28 museum specimens, and 289 modern samples to investigate patterns of gene flow and temporal changes in the genetic diversity of this endangered species over the last 3,000 years, as Polynesians and then Europeans colonized the Hawaiian Islands. Genetic differentiation was found to be high between both modern and ancient petrel populations. However, gene flow was substantial between the extirpated colonies on Oahu and Molokai and modern birds from the island of Lanai. No significant reductions in genetic diversity occurred over this period, despite fears in the mid-1900s that this species may have been extinct. Simulations show that even a decline to a stable effective population size of 100 individuals would result in the loss of only 5% of the expected heterozygosity. Simulations also show that high levels of genetic diversity may be retained due to the long generation time of this species. Such decoupling between population size and genetic diversity in long-lived species can have important conservation implications. It appears that a pattern of dispersal from declining colonies, in addition to long generation time, may have allowed the Hawaiian petrel to escape a severe genetic bottleneck, and the associated extinction vortex, and persist despite a large population decline after human colonization.