Ancient DNA reveals lack of postglacial habitat tracking in the arctic fox

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens)

BACKGROUND

There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity.

RESULTS

A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks.

CONCLUSIONS

Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Mitochondrial genetic diversity and phylogenetic relationships of Echinococcus multilocularis in Europe

The cestode Echinococcus multilocularis is the causative agent of alveolar echinococcosis, a fatal zoonotic parasitic disease of the northern hemisphere. Red foxes are the main reservoir hosts and, likely, the main drivers of the geographic spread of the disease in Europe. Knowledge of genetic relationships among E. multilocularis isolates at a European scale is key to understanding the dispersal characteristics of E. multilocularis. Hence, the present study aimed to describe the genetic diversity of E. multilocularis isolates obtained from different host species in 19 European countries. Based on the analysis of complete nucleotide sequences of the cob, atp6, nad2, nad1 and cox1 mitochondrial genes (4,968 bp), 43 haplotypes were inferred. Four haplotypes represented 62.56 % of the examined isolates (142/227), and one of these four haplotypes was found in each country investigated, except Svalbard, Norway. While the haplotypes from Svalbard were markedly different from all the others, mainland Europe appeared to be dominated by two main clusters, represented by most western, central and eastern European countries, and the Baltic countries and northeastern Poland, respectively. Moreover, one Asian-like haplotype was identified in Latvia and northeastern Poland. To better elucidate the presence of Asian genetic variants of E. multilocularis in Europe, and to obtain a more comprehensive Europe-wide coverage, further studies, including samples from endemic regions not investigated in the present study, especially some eastern European countries, are needed. Further, the present work proposes historical causes that may have contributed to shaping the current genetic variability of E. multilocularis in Europe.

The origins and diversification of Holarctic brown bear populations inferred from genomes of past and present populations

The brown bear (Ursus arctos) is one of the survivors of the Late Quaternary megafauna extinctions. However, despite being widely distributed across the Holarctic, brown bears have experienced extensive range reductions, and even extirpations in some geographical regions. Previous research efforts using genetic data have provided valuable insights into their evolutionary history. However, most studies have been limited to contemporary individuals or mitochondrial DNA, limiting insights into population processes that preceded the present. Here, we present genomic data from two Late Pleistocene brown bears from Honshu, Japan and eastern Siberia, and combine them with published contemporary and ancient genomes from across the Holarctic range of brown bears to investigate the evolutionary relationships among brown bear populations through time and space. By including genomic data from Late Pleistocene and Holocene individuals sampled outside the current distribution range, we uncover diversity not present in contemporary populations. Notably, although contemporary individuals display geographically structured populations most likely driven by isolation-by-distance, this pattern varies among the ancient samples across different regions. The inclusion of ancient brown bears in our analysis provides novel insights into the evolutionary history of brown bears and contributes to understanding the populations and diversity lost during the Late Quaternary.

In Search of the Elusive North: Evolutionary History of the Arctic Fox (Vulpes lagopus) in the Palearctic from the Late Pleistocene to the Recent Inferred from Mitogenomic Data

Despite the high level of interest, the population history of arctic foxes during the Late Pleistocene and Holocene remains poorly understood. Here we aimed to fill gaps in the demographic and colonization history of the arctic fox by analyzing new ancient DNA data from fossil specimens aged from 50 to 1 thousand years from the Northern and Polar Urals, historic DNA from museum specimens from the Novaya Zemlya Archipelago and the Taymyr Peninsula and supplementing these data by previously published sequences of recent and extinct arctic foxes from other regions. This dataset was used for reconstruction of a time-calibrated phylogeny and a temporal haplotype network covering four time intervals: Late Pleistocene (ranging from 30 to 13 thousand years bp), Holocene (ranging from 4 to 1 thousand years bp), historical (approximately 150 years), and modern. Our results revealed that Late Pleistocene specimens showed no genetic similarity to either modern or historical specimens, thus supporting the earlier hypothesis on local extinction rather than habitat tracking.

Pan‐European phylogeography of the European roe deer ( Capreolus capreolus )

To provide the most comprehensive picture of species phylogeny and phylogeography of European roe deer (Capreolus capreolus), we analyzed mtDNA control region (610 bp) of 1469 samples of roe deer from Central and Eastern Europe and included into the analyses additional 1541 mtDNA sequences from GenBank from other regions of the continent. We detected two mtDNA lineages of the species: European and Siberian (an introgression of C. pygargus mtDNA into C. capreolus). The Siberian lineage was most frequent in the eastern part of the continent and declined toward Central Europe. The European lineage contained three clades (Central, Eastern, and Western) composed of several haplogroups, many of which were separated in space. The Western clade appeared to have a discontinuous range from Portugal to Russia. Most of the haplogroups in the Central and the Eastern clades were under expansion during the Weichselian glacial period before the Last Glacial Maximum (LGM), while the expansion time of the Western clade overlapped with the Eemian interglacial. The high genetic diversity of extant roe deer is the result of their survival during the LGM probably in a large, contiguous range spanning from the Iberian Peninsula to the Caucasus Mts and in two northern refugia.

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability

Pleistocene environmental changes are generally assumed to have dramatically affected species’ demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120–14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome‐wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present‐day levels during the Last Glacial Maximum (LGM; 26–19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species’ demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species’ demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

Phylogeography and genetic diversity of the microbivalve Kidderia subquadrata, reveals new data from West Antarctic Peninsula

It is well established that Antarctic biodiversity has been strongly influenced by rapid climatic fluctuations during the Quaternary. Marine invertebrates from Antarctica constitute an interesting lens through which to study the impacts of the last glacial periods as glaciation impacted the distribution and intraspecific genetic variation of these animals. However, the impact on the spatial genetic distribution and historical demography of local processes in areas adjacent to the West Antarctic Peninsula (WAP) is less clear. Here we present new genetic information on the bivalve Kidderia subquadrata, a small mollusk that inhabits intertidal rocky island ecosystems throughout the WAP. Using a phylogeographical approach, we examined the spatial patterns of genetic diversity in this brooder species to test the hypothesis of strong genetic structure in incubating organisms and the hypothesis of glacial refugia in organisms with limited dispersion. We found evidence of strong genetic structure among populations of the WAP and a recent expansion in the South Shetland Islands. Our findings are concordant with the predictions that incubating organisms, abundant in Antarctica, present a strong genetic structure among their populations and also support the hypothesis of glacial refugia in organisms with limited dispersion. The effect of the coastal current pattern in the WAP is suggested as a driver to the local spatial dynamics of the genetic diversity distribution. Although genetic information about this microbivalve is still scarce, the knowledge reported here has increased our understanding of the evolutionary patterns of this organism that is endemic to the Southern Ocean.

Using paleo-archives to safeguard biodiversity under climate change

Strategies for 21st-century environmental management and conservation under global change require a strong understanding of the biological mechanisms that mediate responses to climate- and human-driven change to successfully mitigate range contractions, extinctions, and the degradation of ecosystem services. Biodiversity responses to past rapid warming events can be followed in situ and over extended periods, using cross-disciplinary approaches that provide cost-effective and scalable information for species' conservation and the maintenance of resilient ecosystems in many bioregions. Beyond the intrinsic knowledge gain such integrative research will increasingly provide the context, tools, and relevant case studies to assist in mitigating climate-driven biodiversity losses in the 21st century and beyond.

Were ancient foxes far more carnivorous than recent ones?-Carnassial morphological evidence

Crown shape variation of the first lower molar in the arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) was analyzed using five groups of morphotypes. Carnassial morphologies were compared between the species and between spatially and temporally distant populations: one Late Pleistocene (n = 45) and seven modern populations of the arctic fox (n = 259), and one Late Pleistocene (n = 35) and eight modern populations of the red fox (n = 606). The dentition of Holocene red foxes had larger morphotype variability than that of arctic foxes. The lower carnassials of the red fox kept have some primitive characters (additional cusps and stylids, complex shape of transverse cristid), whereas the first lower molars of the arctic fox have undergone crown shape simplification, with the occlusal part of the tooth undergoing a more pronounced adaptation to a more carnivorous diet. From the Late Pleistocene of Belgium to the present days, the arctic fox’s crown shape has been simplified and some primitive characters have disappeared. In the red fox chronological changes in the morphology of the lower carnassials were not clearly identified. The phyletic tree based on morphotype carnassial characteristics indicated the distinctiveness of both foxes: in the arctic fox line, the ancient population from Belgium and recent Greenland made separate branches, whereas in the red foxes the ancient population from Belgium was most similar to modern red foxes from Belgium and Italy.

Size of the lower carnassial in the arctic and the red fox from Late Pleistocene in Belgium compared to other ancient and extant populations

Lengths, widths, and size proportions (length to width) of the lower carnassial were measured in 45 teeth of the arctic fox and 35 teeth of the red fox from Belgium radiocarbon dated to 46 640–14 120 ka BP. Data the Late Pleistocene foxes from Belgium were compared to 20 ancient and extant populations form Europe, Asia, and North America. The Pleistocene arctic fox from Belgium showed larger carnassial than in all recent samples of this species, whereas the Belgian fossil red foxes were characterized by the carnassial size comparable to that of the recent Siberian red foxes. Both fox species from the Pleistocene of Belgium showed the highest index of the carnassials length to width, which means increase in carnivorous adaptation. We conclude that the higher level of carnivorous specialization reached by the Belgian arctic and red foxes at the end of the Late Pleistocene reflected their scavenging on kills of large carnivores and human hunters (remains of megafauna). Harsh environmental conditions of that period and specific composition of ecosystems led to adapting to a more carnivorous food niche in both foxes.

Consequences of past climate change and recent human persecution on mitogenomic diversity in the arctic fox

Ancient DNA provides a powerful means to investigate the timing, rate and extent of population declines caused by extrinsic factors, such as past climate change and human activities. One species probably affected by both these factors is the arctic fox, which had a large distribution during the last glaciation that subsequently contracted at the start of the Holocene. More recently, the arctic fox population in Scandinavia went through a demographic bottleneck owing to human persecution. To investigate the consequences of these processes, we generated mitogenome sequences from a temporal dataset comprising Pleistocene, historical and modern arctic fox samples. We found no evidence that Pleistocene populations in mid-latitude Europe or Russia contributed to the present-day gene pool of the Scandinavian population, suggesting that postglacial climate warming led to local population extinctions. Furthermore, during the twentieth-century bottleneck in Scandinavia, at least half of the mitogenome haplotypes were lost, consistent with a 20-fold reduction in female effective population size. In conclusion, these results suggest that the arctic fox in mainland Western Europe has lost genetic diversity as a result of both past climate change and human persecution. Consequently, it might be particularly vulnerable to the future challenges posed by climate change. This article is part of a discussion meeting issue 'The past is a foreign country: how much can the fossil record actually inform conservation?'

Unlocking the potential of ancient fish DNA in the genomic era

Fish are the most diverse group of vertebrates, fulfil important ecological functions and are of significant economic interest for aquaculture and wild fisheries. Advances in DNA extraction methods, sequencing technologies and bioinformatic applications have advanced genomic research for nonmodel organisms, allowing the field of fish ancient DNA (aDNA) to move into the genomics era. This move is enabling researchers to investigate a multitude of new questions in evolutionary ecology that could not, until now, be addressed. In many cases, these new fields of research have relevance to evolutionary applications, such as the sustainable management of fisheries resources and the conservation of aquatic animals. Here, we focus on the application of fish aDNA to (a) highlight new research questions, (b) outline methodological advances and current challenges, (c) discuss how our understanding of fish ecology and evolution can benefit from aDNA applications and (d) provide a future perspective on how the field will help answer key questions in conservation and management. We conclude that the power of fish aDNA will be unlocked through the application of continually improving genomic resources and methods to well-chosen taxonomic groups represented by well-dated archaeological samples that can provide temporally and/or spatially extensive data sets.

Recent dating of extinct Atlantic gray whale fossils, (Eschrichtius robustus), Georgia Bight and Florida, western Atlantic Ocean

The Atlantic gray whale (Eschrichtius robustus) presents an interesting case study of climate related dispersal and extinction. While (limited) fossil records confirm its presence in the Atlantic up until the 18th Century, its abundance and distribution within the Eastern and Western basins are still not well understood. The discovery of presumed gray whale fossil remains from the Georgia Bight and the Atlantic coast of Florida, from the mid-1980s to late-2000s, provides a new opportunity to recover additional data regarding their chronology within the Western basin. Here, we apply accelerator mass spectrometry radiocarbon techniques to six fossil whale finds, identifying dates within marine isotope stage 3 (59-24 ka) and the late Holocene, ∼2,000 yr BP. We additionally confirm the taxonomic identification of two fossil bone samples as E. robustus using collagen peptide mass fingerprinting (ZooMS). The obtained dates, when combined with a larger corpus of previously published Atlantic gray whale fossil dates, support the hypothesis for the decline of the Atlantic gray whale in the late Pleistocene and the late Holocene. These new data augment the findings of the Eastern Atlantic Basin and better incorporate the Western Atlantic Basin into a pan-ocean understanding for the species.

Increased mitochondrial DNA diversity in ancient Columbia River basin Chinook salmon Oncorhynchus tshawytscha

The Columbia River and its tributaries provide essential spawning and rearing habitat for many salmonid species, including Chinook salmon (Oncorhynchus tshawytscha). Chinook salmon were historically abundant throughout the basin and Native Americans in the region relied heavily on these fish for thousands of years. Following the arrival of Europeans in the 1800s, salmon in the basin experienced broad declines linked to overfishing, water diversion projects, habitat destruction, connectivity reduction, introgression with hatchery-origin fish, and hydropower development. Despite historical abundance, many native salmonids are now at risk of extinction. Research and management related to Chinook salmon is usually explored under what are termed "the four H's": habitat, harvest, hatcheries, and hydropower; here we explore a fifth H, history. Patterns of prehistoric and contemporary mitochondrial DNA variation from Chinook salmon were analyzed to characterize and compare population genetic diversity prior to recent alterations and, thus, elucidate a deeper history for this species. A total of 346 ancient and 366 contemporary samples were processed during this study. Species was determined for 130 of the ancient samples and control region haplotypes of 84 of these were sequenced. Diversity estimates from these 84 ancient Chinook salmon were compared to 379 contemporary samples. Our analysis provides the first direct measure of reduced genetic diversity for Chinook salmon from the ancient to the contemporary period, as measured both in direct loss of mitochondrial haplotypes and reductions in haplotype and nucleotide diversity. However, these losses do not appear equal across the basin, with higher losses of diversity in the mid-Columbia than in the Snake subbasin. The results are unexpected, as the two groups were predicted to share a common history as parts of the larger Columbia River Basin, and instead indicate that Chinook salmon in these subbasins may have divergent demographic histories.

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum.

The Effects of Paleoclimatic Events on Mediterranean Trout: Preliminary Evidences from Ancient DNA

In this pilot study for the first time, ancient DNA has been extracted from bone remains of Salmo trutta. These samples were from a stratigraphic succession located in a coastal cave of Calabria (southern Italy) inhabited by humans from upper Palaeolithic to historical times. Seven pairs of primers were used to PCR-amplify and sequence from 128 to 410 bp of the mtDNA control region of eleven samples. Three haplotypes were observed: two (ADcs-1 and MEcs-1) already described in rivers from the Italian peninsula; one (ATcs-33) belonging to the southern Atlantic clade of the AT Salmo trutta mtDNA lineage (sensu Bernatchez). The prehistoric occurrence of this latter haplotype in the water courses of the Italian peninsula has been detected for the first time in this study. Finally, we observed a correspondence between frequency of trout remains and variation in haplotype diversity that we related with ecological and demographic changes resulting from a period of rapid cooling known as the Younger Dryas.

Paleodistribution modeling suggests glacial refugia in Scandinavia and out-of-Tibet range expansion of the Arctic fox

Quaternary glacial cycles have shaped the geographic distributions and evolution of numerous species in the Arctic. Ancient DNA suggests that the Arctic fox went extinct in Europe at the end of the Pleistocene and that Scandinavia was subsequently recolonized from Siberia, indicating inability to track its habitat through space as climate changed. Using ecological niche modeling, we found that climatically suitable conditions for Arctic fox were found in Scandinavia both during the last glacial maximum (LGM) and the mid-Holocene. Our results are supported by fossil occurrences from the last glacial. Furthermore, the model projection for the LGM, validated with fossil records, suggested an approximate distance of 2000 km between suitable Arctic conditions and the Tibetan Plateau well within the dispersal distance of the species, supporting the recently proposed hypothesis of range expansion from an origin on the Tibetan Plateau to the rest of Eurasia. The fact that the Arctic fox disappeared from Scandinavia despite suitable conditions suggests that extant populations may be more sensitive to climate change than previously thought.

Genes Left Behind: Climate Change Threatens Cryptic Genetic Diversity in the Canopy-Forming Seaweed Bifurcaria bifurcata

The global redistribution of biodiversity will intensify in the coming decades of climate change, making projections of species range shifts and of associated genetic losses important components of conservation planning. Highly-structured marine species, notably brown seaweeds, often harbor unique genetic variation at warmer low-latitude rear edges and thus are of particular concern. Here, a combination of Ecological Niche Models (ENMs) and molecular data is used to forecast the potential near-future impacts of climate change for a warm-temperate, canopy forming seaweed, Bifurcaria bifurcata. ENMs for B. bifurcata were developed using marine and terrestrial climatic variables, and its range projected for 2040-50 and 2090-2100 under two greenhouse emission scenarios. Geographical patterns of genetic diversity were assessed by screening 18 populations spawning the entire distribution for two organelle genes and 6 microsatellite markers. The southern limit of B. bifurcata was predicted to shift northwards to central Morocco by the mid-century. By 2090-2100, depending on the emission scenario, it could either retreat further north to western Iberia or be relocated back to Western Sahara. At the opposing margin, B. bifurcata was predicted to expand its range to Scotland or even Norway. Microsatellite diversity and endemism were highest in Morocco, where a unique and very restricted lineage was also identified. Our results imply that B. bifurcata will maintain a relatively broad latitudinal distribution. Although its persistence is not threatened, the predicted extirpation of a unique southern lineage or even the entire Moroccan diversity hotspot will erase a rich evolutionary legacy and shrink global diversity to current (low) European levels. NW Africa and similarly understudied southern regions should receive added attention if expected range changes and diversity loss of warm-temperate species is not to occur unnoticed.

A 2.5-million-year perspective on coarse-filter strategies for conserving nature's stage

Climate change will require novel conservation strategies. One such tactic is a coarse-filter approach that focuses on conserving nature's stage (CNS) rather than the actors (individual species). However, there is a temporal mismatch between the long-term goals of conservation and the short-term nature of most ecological studies, which leaves many assumptions untested. Paleoecology provides a valuable perspective on coarse-filter strategies by marshaling the natural experiments of the past to contextualize extinction risk due to the emerging impacts of climate change and anthropogenic threats. We reviewed examples from the paleoecological record that highlight the strengths, opportunities, and caveats of a CNS approach. We focused on the near-time geological past of the Quaternary, during which species were subjected to widespread changes in climate and concomitant changes in the physical environment in general. Species experienced a range of individualistic responses to these changes, including community turnover and novel associations, extinction and speciation, range shifts, changes in local richness and evenness, and both equilibrium and disequilibrium responses. Due to the dynamic nature of species responses to Quaternary climate change, a coarse-filter strategy may be appropriate for many taxa because it can accommodate dynamic processes. However, conservationists should also consider that the persistence of landforms varies across space and time, which could have potential long-term consequences for geodiversity and thus biodiversity.

The long term response of birds to climate change: new results from a cold stage avifauna in northern England

The early MIS 3 (55–40 Kyr BP associated with Middle Palaeolithic archaeology) bird remains from Pin Hole, Creswell Crags, Derbyshire, England are analysed in the context of the new dating of the site’s stratigraphy. The analysis is restricted to the material from the early MIS 3 level of the cave because the upper fauna is now known to include Holocene material as well as that from the Late Glacial. The results of the analysis confirm the presence of the taxa, possibly unexpected for a Late Pleistocene glacial deposit including records such as Alpine swift, demoiselle crane and long-legged buzzard with southern and/or eastern distributions today. These taxa are accompanied by more expected ones such as willow ptarmigan /red grouse and rock ptarmigan living today in northern and montane areas. Finally, there are temperate taxa normally requiring trees for nesting such as wood pigeon and grey heron. Therefore, the result of the analysis is that the avifauna of early MIS 3 in England included taxa whose ranges today do not overlap making it a non-analogue community similar to the many steppe-tundra mammalian faunas of the time. The inclusion of more temperate and woodland taxa is discussed in the light that parts of northern Europe may have acted as cryptic northern refugia for some such taxa during the last glacial. These records showing former ranges of taxa are considered in the light of modern phylogeographic studies as these often assume former ranges without considering the fossil record of those taxa. In addition to the anomalous combination of taxa during MIS 3 living in Derbyshire, the individuals of a number of the taxa are different in size and shape to members of the species today probably due to the high carrying capacity of the steppe-tundra.

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.

Extinction and recolonization of coastal megafauna following human arrival in New Zealand

Extinctions can dramatically reshape biological communities. As a case in point, ancient mass extinction events apparently facilitated dramatic new evolutionary radiations of surviving lineages. However, scientists have yet to fully understand the consequences of more recent biological upheaval, such as the megafaunal extinctions that occurred globally over the past 50 kyr. New Zealand was the world's last large landmass to be colonized by humans, and its exceptional archaeological record documents a vast number of vertebrate extinctions in the immediate aftermath of Polynesian arrival approximately AD 1280. This recently colonized archipelago thus presents an outstanding opportunity to test for rapid biological responses to extinction. Here, we use ancient DNA (aDNA) analysis to show that extinction of an endemic sea lion lineage (Phocarctos spp.) apparently facilitated a subsequent northward range expansion of a previously subantarctic-limited lineage. This finding parallels a similar extinction-replacement event in penguins (Megadyptes spp.). In both cases, an endemic mainland clade was completely eliminated soon after human arrival, and then replaced by a genetically divergent clade from the remote subantarctic region, all within the space of a few centuries. These data suggest that ecological and demographic processes can play a role in constraining lineage distributions, even for highly dispersive species, and highlight the potential for dynamic biological responses to extinction.

Revealing the maternal demographic history of Panthera leo using ancient DNA and a spatially explicit genealogical analysis

Background

Understanding the demographic history of a population is critical to conservation and to our broader understanding of evolutionary processes. For many tropical large mammals, however, this aim is confounded by the absence of fossil material and by the misleading signal obtained from genetic data of recently fragmented and isolated populations. This is particularly true for the lion which as a consequence of millennia of human persecution, has large gaps in its natural distribution and several recently extinct populations.

Results

We sequenced mitochondrial DNA from museum-preserved individuals, including the extinct Barbary lion (Panthera leo leo) and Iranian lion (P. l. persica), as well as lions from West and Central Africa. We added these to a broader sample of lion sequences, resulting in a data set spanning the historical range of lions. Our Bayesian phylogeographical analyses provide evidence for highly supported, reciprocally monophyletic lion clades. Using a molecular clock, we estimated that recent lion lineages began to diverge in the Late Pleistocene. Expanding equatorial rainforest probably separated lions in South and East Africa from other populations. West African lions then expanded into Central Africa during periods of rainforest contraction. Lastly, we found evidence of two separate incursions into Asia from North Africa, first into India and later into the Middle East.

Conclusions

We have identified deep, well-supported splits within the mitochondrial phylogeny of African lions, arguing for recognition of some regional populations as worthy of independent conservation. More morphological and nuclear DNA data are now needed to test these subdivisions.

Persistence of within-species lineages: a neglected control of speciation rates

We present a framework distinguishing three principal controls of speciation rate: rate of splitting, level of persistence, and length of speciation duration. We contend that discussions on diversification become clearer in the light of this framework, because speciation rate variation could be attributed to any of these controls. In particular, we claim that the role of persistence of within-species lineages in controlling speciation rates has been greatly underappreciated. More emphasis on the persistence control would change expectations of the role of several biological traits and environmental factors, because they may drive speciation rate in one direction through the persistence control and in the opposite direction through the other two controls. Traits and environments have been little studied regarding their influence on speciation rate through the persistence control, with climatic fluctuations being a relatively well-studied exception. Considering the recent advances in genomic and phylogenetic analysis, we think that the time is ripe for applying the framework in empirical research. Variation among clades and areas (and thus among traits and environments) in the importance of the three rate controls could be addressed for example by dating splitting events, detecting within-species lineages, and scanning genomes for evidence of divergent selection.

Evaluating the significance of paleophylogeographic species distribution models in reconstructing quaternary range-shifts of nearctic chelonians

The climatic cycles of the Quaternary, during which global mean annual temperatures have regularly changed by 5-10°C, provide a special opportunity for studying the rate, magnitude, and effects of geographic responses to changing climates. During the Quaternary, high- and mid-latitude species were extirpated from regions that were covered by ice or otherwise became unsuitable, persisting in refugial retreats where the environment was compatible with their tolerances. In this study we combine modern geographic range data, phylogeny, Pleistocene paleoclimatic models, and isotopic records of changes in global mean annual temperature, to produce a temporally continuous model of geographic changes in potential habitat for 59 species of North American turtles over the past 320 Ka (three full glacial-interglacial cycles). These paleophylogeographic models indicate the areas where past climates were compatible with the modern ranges of the species and serve as hypotheses for how their geographic ranges would have changed in response to Quaternary climate cycles. We test these hypotheses against physiological, genetic, taxonomic and fossil evidence, and we then use them to measure the effects of Quaternary climate cycles on species distributions. Patterns of range expansion, contraction, and fragmentation in the models are strongly congruent with (i) phylogeographic differentiation; (ii) morphological variation; (iii) physiological tolerances; and (iv) intraspecific genetic variability. Modern species with significant interspecific differentiation have geographic ranges that strongly fluctuated and repeatedly fragmented throughout the Quaternary. Modern species with low genetic diversity have geographic distributions that were highly variable and at times exceedingly small in the past. Our results reveal the potential for paleophylogeographic models to (i) reconstruct past geographic range modifications, (ii) identify geographic processes that result in genetic bottlenecks; and (iii) predict threats due to anthropogenic climate change in the future.

The shaping of genetic variation in edge-of-range populations under past and future climate change

With rates of climate change exceeding the rate at which many species are able to shift their range or adapt, it is important to understand how future changes are likely to affect biodiversity at all levels of organisation. Understanding past responses and extent of niche conservatism in climatic tolerance can help predict future consequences. We use an integrated approach to determine the genetic consequences of past and future climate changes on a bat species, Plecotus austriacus. Glacial refugia predicted by palaeo-modelling match those identified from analyses of extant genetic diversity and model-based inference of demographic history. Former refugial populations currently contain disproportionately high genetic diversity, but niche conservatism, shifts in suitable areas and barriers to migration mean that these hotspots of genetic diversity are under threat from future climate change. Evidence of population decline despite recent northward migration highlights the need to conserve leading-edge populations for spearheading future range shifts.

Holarctic genetic structure and range dynamics in the woolly mammoth

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32–34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20–15 kyr BP, when a severe population size decline occurred.

Could brown bears (Ursus arctos) have survived in Ireland during the Last Glacial Maximum?

Brown bears are recorded from Ireland during both the Late Pleistocene and early-mid Holocene. Although most of the Irish landmass was covered by an ice sheet during the Last Glacial Maximum (LGM), Irish brown bears are known to have hybridized with polar bears during the Late Pleistocene, and it is suggested that the Irish brown bear population did not become extinct but instead persisted in situ through the LGM in a southwestern ice-free refugium. We use historical population modelling to demonstrate that brown bears are highly unlikely to have survived through the LGM in Ireland under any combination of life-history parameters shown by living bear populations, but instead would have rapidly become extinct following advance of the British-Irish ice sheet, and probably recolonized Ireland during the end-Pleistocene Woodgrange Interstadial from a closely related nearby source population. The time available for brown bear-polar bear hybridization was therefore restricted to narrow periods at the beginning or end of the LGM. Brown bears would have been extremely vulnerable to extinction in Quaternary habitat refugia and required areas substantially larger than southwestern Ireland to survive adverse glacial conditions.

The de novo assembly of mitochondrial genomes of the extinct passenger pigeon (Ectopistes migratorius) with next generation sequencing

The information from ancient DNA (aDNA) provides an unparalleled opportunity to infer phylogenetic relationships and population history of extinct species and to investigate genetic evolution directly. However, the degraded and fragmented nature of aDNA has posed technical challenges for studies based on conventional PCR amplification. In this study, we present an approach based on next generation sequencing to efficiently sequence the complete mitochondrial genome (mitogenome) of two extinct passenger pigeons (Ectopistes migratorius) using de novo assembly of massive short (90 bp), paired-end or single-end reads. Although varying levels of human contamination and low levels of postmortem nucleotide lesion were observed, they did not impact sequencing accuracy. Our results demonstrated that the de novo assembly of shotgun sequence reads could be a potent approach to sequence mitogenomes, and offered an efficient way to infer evolutionary history of extinct species.

Serial population extinctions in a small mammal indicate Late Pleistocene ecosystem instability

The Late Pleistocene global extinction of many terrestrial mammal species has been a subject of intensive scientific study for over a century, yet the relative contributions of environmental changes and the global expansion of humans remain unresolved. A defining component of these extinctions is a bias toward large species, with the majority of small-mammal taxa apparently surviving into the present. Here, we investigate the population-level history of a key tundra-specialist small mammal, the collared lemming (Dicrostonyx torquatus), to explore whether events during the Late Pleistocene had a discernible effect beyond the large mammal fauna. Using ancient DNA techniques to sample across three sites in North-West Europe, we observe a dramatic reduction in genetic diversity in this species over the last 50,000 y. We further identify a series of extinction-recolonization events, indicating a previously unrecognized instability in Late Pleistocene small-mammal populations, which we link with climatic fluctuations. Our results reveal climate-associated, repeated regional extinctions in a keystone prey species across the Late Pleistocene, a pattern likely to have had an impact on the wider steppe-tundra community, and one that is concordant with environmental change as a major force in structuring Late Pleistocene biodiversity.

Gene flow on ice: the role of sea ice and whaling in shaping Holarctic genetic diversity and population differentiation in bowhead whales (Balaena mysticetus)

Sea ice is believed to be a major factor shaping gene flow for polar marine organisms, but it remains unclear to what extent it represents a true barrier to dispersal for arctic cetaceans. Bowhead whales are highly adapted to polar sea ice and were targeted by commercial whalers throughout Arctic and subarctic seas for at least four centuries, resulting in severe reductions in most areas. Both changing ice conditions and reductions due to whaling may have affected geographic distribution and genetic diversity throughout their range, but little is known about range-wide genetic structure or whether it differed in the past. This study represents the first examination of genetic diversity and differentiation across all five putative stocks, including Baffin Bay-Davis Strait, Hudson Bay-Foxe Basin, Bering-Beaufort-Chukchi, Okhotsk, and Spitsbergen. We also utilized ancient specimens from Prince Regent Inlet (PRI) in the Canadian Arctic and compared them with modern stocks. Results from analysis of molecular variance and demographic simulations are consistent with recent and high gene flow between Atlantic and Pacific stocks in the recent past. Significant genetic differences between ancient and modern populations suggest PRI harbored unique maternal lineages in the past that have been recently lost, possibly due to loss of habitat during the Little Ice Age and/or whaling. Unexpectedly, samples from this location show a closer genetic relationship with modern Pacific stocks than Atlantic, supporting high gene flow between the central Canadian Arctic and Beaufort Sea over the past millennium despite extremely heavy ice cover over much of this period.

Phylogenetic Estimation of Timescales Using Ancient DNA: The Effects of Temporal Sampling Scheme and Uncertainty in Sample Ages

In recent years, ancient DNA has increasingly been used for estimating molecular timescales, particularly in studies of substitution rates and demographic histories. Molecular clocks can be calibrated using temporal information from ancient DNA sequences. This information comes from the ages of the ancient samples, which can be estimated by radiocarbon dating the source material or by dating the layers in which the material was deposited. Both methods involve sources of uncertainty. The performance of Bayesian phylogenetic inference depends on the information content of the data set, which includes variation in the DNA sequences and the structure of the sample ages. Various sources of estimation error can reduce our ability to estimate rates and timescales accurately and precisely. We investigated the impact of sample-dating uncertainties on the estimation of evolutionary timescale parameters using the software BEAST. Our analyses involved 11 published data sets and focused on estimates of substitution rate and root age. We show that, provided that samples have been accurately dated and have a broad temporal span, it might be unnecessary to account for sample-dating uncertainty in Bayesian phylogenetic analyses of ancient DNA. We also investigated the sample size and temporal span of the ancient DNA sequences needed to estimate phylogenetic timescales reliably. Our results show that the range of sample ages plays a crucial role in determining the quality of the results but that accurate and precise phylogenetic estimates of timescales can be made even with only a few ancient sequences. These findings have important practical consequences for studies of molecular rates, timescales, and population dynamics.

Ancient DNA studies: new perspectives on old samples

In spite of past controversies, the field of ancient DNA is now a reliable research area due to recent methodological improvements. A series of recent large-scale studies have revealed the true potential of ancient DNA samples to study the processes of evolution and to test models and assumptions commonly used to reconstruct patterns of evolution and to analyze population genetics and palaeoecological changes. Recent advances in DNA technologies, such as next-generation sequencing make it possible to recover DNA information from archaeological and paleontological remains allowing us to go back in time and study the genetic relationships between extinct organisms and their contemporary relatives. With the next-generation sequencing methodologies, DNA sequences can be retrieved even from samples (for example human remains) for which the technical pitfalls of classical methodologies required stringent criteria to guaranty the reliability of the results. In this paper, we review the methodologies applied to ancient DNA analysis and the perspectives that next-generation sequencing applications provide in this field.