Staying out in the cold: glacial refugia and mitochondrial DNA phylogeography in ancient European brown bears

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records. This article is part of the theme issue 'Cross-disciplinary approaches to prehistoric demography'.

Genetic analysis of Bromus tectorum (Poaceae) in the Mediterranean region: biogeographical pattern of native populations

Genetic diversity within and among 42 native populations of Bromus tectorum (cheatgrass) was characterized within two regions, the eastern Mediterranean and the western Mediterranean. Two hypotheses were tested for the genetic diversity of these populations: (1) populations from the eastern Mediterranean are more genetically diverse compared with populations to the west, a potential consequence of the species' westward dispersal with the spread of agriculture, and (2) populations across the Mediterranean contain comparable genetic diversity but display high genetic differentiation, a potential consequence of both regions having served as refugia during glacial advances in the late Quaternary Period. Populations in the eastern Mediterranean possess 16 polymorphic loci and 37 multilocus genotypes. In contrast, populations from the western Mediterranean include a subset of these polymorphic loci (9) and fewer multilocus genotypes (19), consistent with the dispersal of B. tectorum with the east-west Holocene spread of agriculture. Among the 19 multilocus genotypes identified in populations from the western Mediterranean, 13 are undetected among eastern Mediterranean populations. Average genetic diversity within populations from the eastern Mediterranean is nonetheless comparable to the genetic diversity in populations from the Iberian Peninsula, whereas diversity is the lowest in the populations from southern France. Our results suggest a prominent role for agriculture in the grass's western spread, although glacial history and environmental heterogeneity also could have influenced the grass's genetic diversity. The exceptionally high level of self-pollination (>99%) in B. tectorum has contributed to preserving the genetic signature associated with the species' biogeographical history across the Mediterranean region.

Genetic turnovers and northern survival during the last glacial maximum in European brown bears

The current phylogeographic pattern of European brown bears (Ursus arctos) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid-Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter-specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought.

Ancient DNA from Giant Panda (Ailuropoda melanoleuca) of South-Western China Reveals Genetic Diversity Loss during the Holocene

The giant panda was widely distributed in China and south-eastern Asia during the middle to late Pleistocene, prior to its habitat becoming rapidly reduced in the Holocene. While conservation reserves have been established and population numbers of the giant panda have recently increased, the interpretation of its genetic diversity remains controversial. Previous analyses, surprisingly, have indicated relatively high levels of genetic diversity raising issues concerning the efficiency and usefulness of reintroducing individuals from captive populations. However, due to a lack of DNA data from fossil specimens, it is unknown whether genetic diversity was even higher prior to the most recent population decline. We amplified complete cytb and 12s rRNA, partial 16s rRNA and ND1, and control region sequences from the mitochondrial genomes of two Holocene panda specimens. We estimated genetic diversity and population demography by analyzing the ancient mitochondrial DNA sequences alongside those from modern giant pandas, as well as from other members of the bear family (Ursidae). Phylogenetic analyses show that one of the ancient haplotypes is sister to all sampled modern pandas and the second ancient individual is nested among the modern haplotypes, suggesting that genetic diversity may indeed have been higher earlier during the Holocene. Bayesian skyline plot analysis supports this view and indicates a slight decline in female effective population size starting around 6000 years B.P., followed by a recovery around 2000 years ago. Therefore, while the genetic diversity of the giant panda has been affected by recent habitat contraction, it still harbors substantial genetic diversity. Moreover, while its still low population numbers require continued conservation efforts, there seem to be no immediate threats from the perspective of genetic evolutionary potential.

Extra-Mediterranean glacial refugia in a Mediterranean faunal element: the phylogeography of the chalk-hill blue Polyommatus coridon (Lepidoptera, Lycaenidae)

Most warm-adapted Central European species are thought to have survived ice ages exclusively in Mediterranean refugia. During recent years, this point of view has been questioned. Therefore, we tested the hypothesis that extra-Mediterranean refugia also played a role in warm-adapted insect species and selected the chalk-hill blue, Polyommatus coridon. We sequenced two mitochondrial loci (COI, CR) in 150 individuals from 30 populations covering nearly the complete range. Minimum spanning networks and other statistical analyses concordantly revealed four genetic lineages with strong phylogeographic signal: a western group in Italy, France and western/central Germany, an eastern lineage in the Balkan Peninsula, the Carpathian Basin and eastern Central Europe, an Alpine group with populations in the Alps and southern Germany and a Pyrenean group. Our results are generally consistent with previous analyses for P. coridon based on allozymes and DNA sequences, but provide additional insights. We propose that these four lineages have evolved during allopatry in different glacial refugia, two in typical Mediterranean refugia (Apennines and Balkan Peninsulas), but two in extra-Mediterranean areas south of the Alps and Pyrenees. This supports survival of warm-adapted organisms in these regions in close geographic proximity to the refugia of high mountain species.

The colonization history of British water vole (Arvicola amphibius (Linnaeus, 1758)): origins and development of the Celtic fringe

The terminal Pleistocene and Early Holocene, a period from 15 000 to 18 000 Before Present (BP), was critical in establishing the current Holarctic fauna, with temperate-climate species largely replacing cold-adapted ones at mid-latitudes. However, the timing and nature of this process remain unclear for many taxa, a point that impacts on current and future management strategies. Here, we use an ancient DNA dataset to test more directly postglacial histories of the water vole (Arvicola amphibius, formerly A terrestris), a species that is both a conservation priority and a pest in different parts of its range. We specifically examine colonization of Britain, where a complex genetic structure can be observed today. Although we focus on population history at the limits of the species' range, the inclusion of additional European samples allows insights into European postglacial colonization events and provides a molecular perspective on water vole taxonomy.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

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.

Ancient population genomics and the study of evolution

Recently, the study of ancient DNA (aDNA) has been greatly enhanced by the development of second-generation DNA sequencing technologies and targeted enrichment strategies. These developments have allowed the recovery of several complete ancient genomes, a result that would have been considered virtually impossible only a decade ago. Prior to these developments, aDNA research was largely focused on the recovery of short DNA sequences and their use in the study of phylogenetic relationships, molecular rates, species identification and population structure. However, it is now possible to sequence a large number of modern and ancient complete genomes from a single species and thereby study the genomic patterns of evolutionary change over time. Such a study would herald the beginnings of ancient population genomics and its use in the study of evolution. Species that are amenable to such large-scale studies warrant increased research effort. We report here progress on a population genomic study of the Adélie penguin (Pygoscelis adeliae). This species is ideally suited to ancient population genomic research because both modern and ancient samples are abundant in the permafrost conditions of Antarctica. This species will enable us to directly address many of the fundamental questions in ecology and evolution.

Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes

Background

The genus Vulpes (true foxes) comprises numerous species that inhabit a wide range of habitats and climatic conditions, including one species, the Arctic fox (Vulpes lagopus) which is adapted to the arctic region. A close relative to the Arctic fox, the red fox (Vulpes vulpes), occurs in subarctic to subtropical habitats. To study the genetic basis of their adaptations to different environments, transcriptome sequences from two Arctic foxes and one red fox individual were generated and analyzed for signatures of positive selection. In addition, the data allowed for a phylogenetic analysis and divergence time estimate between the two fox species.

Results

The de novo assembly of reads resulted in more than 160,000 contigs/transcripts per individual. Approximately 17,000 homologous genes were identified using human and the non-redundant databases. Positive selection analyses revealed several genes involved in various metabolic and molecular processes such as energy metabolism, cardiac gene regulation, apoptosis and blood coagulation to be under positive selection in foxes. Branch site tests identified four genes to be under positive selection in the Arctic fox transcriptome, two of which are fat metabolism genes. In the red fox transcriptome eight genes are under positive selection, including molecular process genes, notably genes involved in ATP metabolism. Analysis of the three transcriptomes and five Sanger re-sequenced genes in additional individuals identified a lower genetic variability within Arctic foxes compared to red foxes, which is consistent with distribution range differences and demographic responses to past climatic fluctuations. A phylogenomic analysis estimated that the Arctic and red fox lineages diverged about three million years ago.

Conclusions

Transcriptome data are an economic way to generate genomic resources for evolutionary studies. Despite not representing an entire genome, this transcriptome analysis identified numerous genes that are relevant to arctic adaptation in foxes. Similar to polar bears, fat metabolism seems to play a central role in adaptation of Arctic foxes to the cold climate, as has been identified in the polar bear, another arctic specialist.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1724-9) contains supplementary material, which is available to authorized users.

Consequences of a demographic bottleneck on genetic structure and variation in the Scandinavian brown bear

The Scandinavian brown bear went through a major decline in population size approximately 100 years ago, due to intense hunting. After being protected, the population subsequently recovered and today numbers in the thousands. The genetic diversity in the contemporary population has been investigated in considerable detail, and it has been shown that the population consists of several subpopulations that display relatively high levels of genetic variation. However, previous studies have been unable to resolve the degree to which the demographic bottleneck impacted the contemporary genetic structure and diversity. In this study, we used mitochondrial and microsatellite DNA markers from pre- and postbottleneck Scandinavian brown bear samples to investigate the effect of the bottleneck. Simulation and multivariate analysis suggested the same genetic structure for the historical and modern samples, which are clustered into three subpopulations in southern, central and northern Scandinavia. However, the southern subpopulation appears to have gone through a marked change in allele frequencies. When comparing the mitochondrial DNA diversity in the whole population, we found a major decline in haplotype numbers across the bottleneck. However, the loss of autosomal genetic diversity was less pronounced, although a significant decline in allelic richness was observed in the southern subpopulation. Approximate Bayesian computations provided clear support for a decline in effective population size during the bottleneck, in both the southern and northern subpopulations. These results have implications for the future management of the Scandinavian brown bear because they indicate a recent loss in genetic diversity and also that the current genetic structure may have been caused by historical ecological processes rather than recent anthropogenic persecution.

Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears

Polar bears are an arctic, marine adapted species that is closely related to brown bears. Genome analyses have shown that polar bears are distinct and genetically homogeneous in comparison to brown bears. However, these analyses have also revealed a remarkable episode of polar bear gene flow into the population of brown bears that colonized the Admiralty, Baranof and Chichagof islands (ABC islands) of Alaska. Here, we present an analysis of data from a large panel of polar bear and brown bear genomes that includes brown bears from the ABC islands, the Alaskan mainland and Europe. Our results provide clear evidence that gene flow between the two species had a geographically wide impact, with polar bear DNA found within the genomes of brown bears living both on the ABC islands and in the Alaskan mainland. Intriguingly, while brown bear genomes contain up to 8.8% polar bear ancestry, polar bear genomes appear to be devoid of brown bear ancestry, suggesting the presence of a barrier to gene flow in that direction.

Long-range gene flow and the effects of climatic and ecological factors on genetic structuring in a large, solitary carnivore: the Eurasian lynx

Due to their high mobility, large terrestrial predators are potentially capable of maintaining high connectivity, and therefore low genetic differentiation among populations. However, previous molecular studies have provided contradictory findings in relation to this. To elucidate patterns of genetic structure in large carnivores, we studied the genetic variability of the Eurasian lynx, Lynx lynx throughout north-eastern Europe using microsatellite, mitochondrial DNA control region and Y chromosome-linked markers. Using SAMOVA we found analogous patterns of genetic structure based on both mtDNA and microsatellites, which coincided with a relatively little evidence for male-biased dispersal. No polymorphism for the cytochrome b and ATP6 mtDNA genes and Y chromosome-linked markers were found. Lynx inhabiting a large area encompassing Finland, the Baltic countries and western Russia formed a single genetic unit, while some marginal populations were clearly divergent from others. The existence of a migration corridor was suggested to correspond with distribution of continuous forest cover. The lowest variability (in both markers) was found in lynx from Norway and Białowieża Primeval Forest (BPF), which coincided with a recent demographic bottleneck (Norway) or high habitat fragmentation (BPF). The Carpathian population, being monomorphic for the control region, showed relatively high microsatellite diversity, suggesting the effect of a past bottleneck (e.g. during Last Glacial Maximum) on its present genetic composition. Genetic structuring for the mtDNA control region was best explained by latitude and snow cover depth. Microsatellite structuring correlated with the lynx's main prey, especially the proportion of red deer (Cervus elaphus) in its diet. Eurasian lynx are capable of maintaining panmictic populations across eastern Europe unless they are severely limited by habitat continuity or a reduction in numbers. Different correlations of mtDNA and microsatellite population divergence patterns with climatic and ecological factors may suggest separate selective pressures acting on males and females in this solitary carnivore.

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.

Ancient DNA and population turnover in southern levantine pigs--signature of the sea peoples migration?

Near Eastern wild boars possess a characteristic DNA signature. Unexpectedly, wild boars from Israel have the DNA sequences of European wild boars and domestic pigs. To understand how this anomaly evolved, we sequenced DNA from ancient and modern pigs from Israel. Pigs from Late Bronze Age (until ca. 1150 BCE) in Israel shared haplotypes of modern and ancient Near Eastern pigs. European haplotypes became dominant only during the Iron Age (ca. 900 BCE). This raises the possibility that European pigs were brought to the region by the Sea Peoples who migrated to the Levant at that time. Then, a complete genetic turnover took place, most likely because of repeated admixture between local and introduced European domestic pigs that went feral. Severe population bottlenecks likely accelerated this process. Introductions by humans have strongly affected the phylogeography of wild animals, and interpretations of phylogeography based on modern DNA alone should be taken with caution.

Persistence across Pleistocene ice ages in Mediterranean and extra-Mediterranean refugia: phylogeographic insights from the common wall lizard

BACKGROUND

Pleistocene climatic oscillations have played a major role in structuring present-day biodiversity. The southern Mediterranean peninsulas have long been recognized as major glacial refugia, from where Northern Europe was post-glacially colonized. However, recent studies have unravelled numerous additional refugia also in northern regions. We investigated the phylogeographic pattern of the widespread Western Palaearctic lizard Podarcis muralis, using a range-wide multilocus approach, to evaluate whether it is concordant with a recent expansion from southern glacial refugia or alternatively from a combination of Mediterranean and northern refugia.

RESULTS

We analyzed DNA sequences of two mitochondrial (cytb and nd4) and three nuclear (acm4, mc1r, and pdc) gene fragments in individuals from 52 localities across the species range, using phylogenetic and phylogeographic methods. The complex phylogeographic pattern observed, with 23 reciprocally monophyletic allo- parapatric lineages having a Pleistocene divergence, suggests a scenario of long-term isolation in multiple ice-age refugia across the species distribution range. Multiple lineages were identified within the three Mediterranean peninsulas - Iberia, Italy and the Balkans - where the highest genetic diversity was observed. Such an unprecedented phylogeographic pattern - here called "refugia within all refugia" - compasses the classical scenario of multiple southern refugia. However, unlike the southern refugia model, various distinct lineages were also found in northern regions, suggesting that additional refugia in France, Northern Italy, Eastern Alps and Central Balkans allowed the long-term persistence of this species throughout Pleistocene glaciations.

CONCLUSIONS

The phylogeography of Podarcis muralis provides a paradigm of temperate species survival in Mediterranean and extra-Mediterranean glacial refugia. Such refugia acted as independent biogeographic compartments for the long-term persistence of this species, for the differentiation of its genetic lineages, and for the short-distance post-glacial re-colonization of neighbouring areas. This finding echoes previous findings from recent phylogeographic studies on species from temperate ecoregions, thus suggesting the need for a reappraisal of the role of northern refugia for glacial persistence and post-glacial assembly of Holarctic biota.

A range-wide synthesis and timeline for phylogeographic events in the red fox (Vulpes vulpes)

BACKGROUND

Many boreo-temperate mammals have a Pleistocene fossil record throughout Eurasia and North America, but only few have a contemporary distribution that spans this large area. Examples of Holarctic-distributed carnivores are the brown bear, grey wolf, and red fox, all three ecological generalists with large dispersal capacity and a high adaptive flexibility. While the two former have been examined extensively across their ranges, no phylogeographic study of the red fox has been conducted across its entire Holarctic range. Moreover, no study included samples from central Asia, leaving a large sampling gap in the middle of the Eurasian landmass.

RESULTS

Here we provide the first mitochondrial DNA sequence data of red foxes from central Asia (Siberia), and new sequences from several European populations. In a range-wide synthesis of 729 red fox mitochondrial control region sequences, including 677 previously published and 52 newly obtained sequences, this manuscript describes the pattern and timing of major phylogeographic events in red foxes, using a Bayesian coalescence approach with multiple fossil tip and root calibration points. In a 335 bp alignment we found in total 175 unique haplotypes. All newly sequenced individuals belonged to the previously described Holarctic lineage. Our analyses confirmed the presence of three Nearctic- and two Japan-restricted lineages that were formed since the Mid/Late Pleistocene.

CONCLUSIONS

The phylogeographic history of red foxes is highly similar to that previously described for grey wolves and brown bears, indicating that climatic fluctuations and habitat changes since the Pleistocene had similar effects on these highly mobile generalist species. All three species originally diversified in Eurasia and later colonized North America and Japan. North American lineages persisted through the last glacial maximum south of the ice sheets, meeting more recent colonizers from Beringia during postglacial expansion into the northern Nearctic. Both brown bears and red foxes colonized Japan's northern island Hokkaido at least three times, all lineages being most closely related to different mainland lineages. Red foxes, grey wolves, and brown bears thus represent an interesting case where species that occupy similar ecological niches also exhibit similar phylogeographic histories.

The integration of multiple independent data reveals an unusual response to Pleistocene climatic changes in the hard tick Ixodes ricinus

In the last few years, improved analytical tools and the integration of genetic data with multiple sources of information have shown that temperate species exhibited more complex responses to ice ages than previously thought. In this study, we investigated how Pleistocene climatic changes affected the current distribution and genetic diversity of European populations of the tick Ixodes ricinus, an ectoparasite with high ecological plasticity. We first used mitochondrial and nuclear genetic markers to investigate the phylogeographic structure of the species and its Pleistocene history using coalescent-based methods; then we used species distribution modelling to infer the climatic niche of the species at last glacial maximum; finally, we reviewed the literature on the I. ricinus hosts to identify the locations of their glacial refugia. Our results support the scenario that during the last glacial phase, I. ricinus never experienced a prolonged allopatric divergence in separate glacial refugia, but persisted with interconnected populations across Southern and Central Europe. The generalist behaviour in host choice of I. ricinus would have played a major role in maintaining connections between its populations. Although most of the hosts persisted in separate refugia, from the point of view of I. ricinus, they represented a continuity of 'bridges' among populations. Our study highlights the importance of species-specific ecology in affecting responses to Pleistocene glacial-interglacial cycles. Together with other cases in Europe and elsewhere, it contributes to setting new hypotheses on how species with wide ecological plasticity coped with Pleistocene climatic changes.

Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles

Quaternary climatic fluctuations have had profound effects on the phylogeographic structure of many species. Classically, species were thought to have become isolated in peninsular refugia, but there is limited evidence that large, non-polar species survived outside traditional refugial areas. We examined the phylogeographic structure of the red fox (Vulpes vulpes), a species that shows high ecological adaptability in the western Palaearctic region. We compared mitochondrial DNA sequences (cytochrome b and control region) from 399 modern and 31 ancient individuals from across Europe. Our objective was to test whether red foxes colonised the British Isles from mainland Europe in the late Pleistocene, or whether there is evidence that they persisted in the region through the Last Glacial Maximum. We found red foxes to show a high degree of phylogeographic structuring across Europe and, consistent with palaeontological and ancient DNA evidence, confirmed via phylogenetic indicators that red foxes were persistent in areas outside peninsular refugia during the last ice age. Bayesian analyses and tests of neutrality indicated population expansion. We conclude that there is evidence that red foxes from the British Isles derived from central European populations that became isolated after the closure of the landbridge with Europe.

Glacial history of a modern invader: phylogeography and species distribution modelling of the Asian tiger mosquito Aedes albopictus

Background

The tiger mosquito, Aedes albopictus, is one of the 100 most invasive species in the world and a vector of human diseases. In the last 30 years, it has spread from its native range in East Asia to Africa, Europe, and the Americas. Although this modern invasion has been the focus of many studies, the history of the species’ native populations remains poorly understood. Here, we aimed to assess the role of Pleistocene climatic changes in shaping the current distribution of the species in its native range.

Methodology/Principal Findings

We investigated the phylogeography, historical demography, and species distribution of Ae. albopictus native populations at the Last Glacial Maximum (LGM). Individuals from 16 localities from East Asia were analyzed for sequence variation at two mitochondrial genes. No phylogeographic structure was observed across the study area. Demographic analyses showed a signature of population expansion that started roughly 70,000 years BP. The occurrence of a continuous and climatically suitable area comprising Southeast China, Indochinese Peninsula, and Sundaland during LGM was indicated by species distribution modelling.

Conclusions/Significance

Our results suggest an evolutionary scenario in which, during the last glacial phase, Ae. albopictus did not experience a fragmentation phase but rather persisted in interconnected populations and experienced demographic growth. The wide ecological flexibility of the species probably played a crucial role in its response to glacial-induced environmental changes. Currently, there is little information on the impact of Pleistocene climatic changes on animal species in East Asia. Most of the studies focused on forest-associated species and suggested cycles of glacial fragmentation and post-glacial expansion. The case of Ae. albopictus, which exhibits a pattern not previously observed in the study area, adds an important piece to our understanding of the Pleistocene history of East Asian biota.

Extra-Mediterranean refugia: The rule and not the exception?

Some decades ago, biogeographers distinguished three major faunal types of high importance for Europe: (i) Mediterranean elements with exclusive glacial survival in the Mediterranean refugia, (ii) Siberian elements with glacial refugia in the eastern Palearctic and only postglacial expansion to Europe and (iii) arctic and/or alpine elements with large zonal distributions in the periglacial areas and postglacial retreat to the North and/or into the high mountain systems. Genetic analyses have unravelled numerous additional refugia both of continental and Mediterranean species, thus strongly modifying the biogeographical view of Europe. This modified notion is particularly true for the so-called Siberian species, which in many cases have not immigrated into Europe during the postglacial period, but most likely have survived the last, or even several glacial phases, in extra-Mediterranean refugia in some climatically favourable but geographically limited areas of southern Central and Eastern Europe. Recently, genetic analyses revealed that typical Mediterranean species have also survived the Last Glacial Maximum in cryptic northern refugia (e.g. in the Carpathians or even north of the Alps) in addition to their Mediterranean refuge areas.

50,000 years of genetic uniformity in the critically endangered Iberian lynx

Low genetic diversity in the endangered Iberian lynx, including lack of mitochondrial control region variation, is thought to result from historical or Pleistocene/Holocene population bottlenecks, and to indicate poor long-term viability. We find no variability in control region sequences from 19 Iberian lynx remains from across the Iberian Peninsula and spanning the last 50,000 years. This is best explained by continuously small female effective population size through time. We conclude that low genetic variability in the Iberian lynx is not in itself a threat to long-term viability, and so should not preclude conservation efforts.

Faunal histories from Holocene ancient DNA

Recent studies using ancient DNA have been instrumental in advancing understanding of the impact of Holocene climate change on biodiversity. Ancient DNA has been used to track demography, migration and diversity, and is providing new insights into the long-term dynamics of species and population distributions. The Holocene is key to understanding how the past has impacted on the present, as it bridges the gap between contemporary phylogeographic studies and those with inference on Pleistocene patterns, based on ancient DNA studies. Here, we examine the major patterns of Holocene faunal population dynamics and connectivity; highlighting the dynamic nature of species and population responses to Holocene climatic change, thereby providing an 'analogue' for understanding potential impacts of future change.

Modern and ancient red fox (Vulpes vulpes) in Europe show an unusual lack of geographical and temporal structuring, and differing responses within the carnivores to historical climatic change

Background

Despite phylogeographical patterns being well characterised in a large number of species, and generalised patterns emerging, the carnivores do not all appear to show consistent trends. While some species tend to fit with standard theoretical phylogeographic expectations (e.g. bears), others show little obvious modern phylogeographic structure (e.g. wolves). In this study we briefly review these studies, and present a new phylogeographical study of the red fox (Vulpes vulpes) throughout Europe, using a combination of ancient DNA sequences obtained from museum specimens, and modern sequences collated from GenBank. We used cytochrome b (250 bp) and the mitochondrial control region (268 bp) to elucidate both current and historical phylogeographical patterning.

Results

We found evidence for slight isolation by distance in modern populations, as well as differentiation associated with time, both of which can likely be attributed to random genetic drift. Despite high sequence diversity (11.2% cytochrome b, 16.4% control region), no evidence for spatial structure (from Bayesian trees) is found either in modern samples or ancient samples for either gene, and Bayesian skyline plots suggested little change in the effective population size over the past 40,000 years.

Conclusions

It is probable that the high dispersal ability and adaptability of the red fox has contributed to the lack of observable differentiation, which appears to have remained consistent over tens of thousands of years. Generalised patterns of how animals are thought to have responded to historical climatic change are not necessarily valid for all species, and so understanding the differences between species will be critical for predicting how species will be affected by future climatic change.

Ancient hybridization and an Irish origin for the modern polar bear matriline

BACKGROUND

Polar bears (Ursus maritimus) are among those species most susceptible to the rapidly changing arctic climate, and their survival is of global concern. Despite this, little is known about polar bear species history. Future conservation strategies would significantly benefit from an understanding of basic evolutionary information, such as the timing and conditions of their initial divergence from brown bears (U. arctos) or their response to previous environmental change.

RESULTS

We used a spatially explicit phylogeographic model to estimate the dynamics of 242 brown bear and polar bear matrilines sampled throughout the last 120,000 years and across their present and past geographic ranges. Our results show that the present distribution of these matrilines was shaped by a combination of regional stability and rapid, long-distance dispersal from ice-age refugia. In addition, hybridization between polar bears and brown bears may have occurred multiple times throughout the Late Pleistocene.

CONCLUSIONS

The reconstructed matrilineal history of brown and polar bears has two striking features. First, it is punctuated by dramatic and discrete climate-driven dispersal events. Second, opportunistic mating between these two species as their ranges overlapped has left a strong genetic imprint. In particular, a likely genetic exchange with extinct Irish brown bears forms the origin of the modern polar bear matriline. This suggests that interspecific hybridization not only may be more common than previously considered but may be a mechanism by which species deal with marginal habitats during periods of environmental deterioration.

Phylogeographic history of grey wolves in Europe

Background

While it is generally accepted that patterns of intra-specific genetic differentiation are substantially affected by glacial history, population genetic processes occurring during Pleistocene glaciations are still poorly understood. In this study, we address the question of the genetic consequences of Pleistocene glaciations for European grey wolves. Combining our data with data from published studies, we analysed phylogenetic relationships and geographic distribution of mitochondrial DNA haplotypes for 947 contemporary European wolves. We also compared the contemporary wolf sequences with published sequences of 24 ancient European wolves.

Results

We found that haplotypes representing two haplogroups, 1 and 2, overlap geographically, but substantially differ in frequency between populations from south-western and eastern Europe. A comparison between haplotypes from Europe and other continents showed that both haplogroups are spread throughout Eurasia, while only haplogroup 1 occurs in contemporary North American wolves. All ancient wolf samples from western Europe that dated from between 44,000 and 1,200 years B.P. belonged to haplogroup 2, suggesting the long-term predominance of this haplogroup in this region. Moreover, a comparison of current and past frequencies and distributions of the two haplogroups in Europe suggested that haplogroup 2 became outnumbered by haplogroup 1 during the last several thousand years.

Conclusions

Parallel haplogroup replacement, with haplogroup 2 being totally replaced by haplogroup 1, has been reported for North American grey wolves. Taking into account the similarity of diets reported for the late Pleistocene wolves from Europe and North America, the correspondence between these haplogroup frequency changes may suggest that they were associated with ecological changes occurring after the Last Glacial Maximum.

Population dynamic of the extinct European aurochs: genetic evidence of a north-south differentiation pattern and no evidence of post-glacial expansion

Background

The aurochs (Bos primigenius) was a large bovine that ranged over almost the entirety of the Eurasian continent and North Africa. It is the wild ancestor of the modern cattle (Bos taurus), and went extinct in 1627 probably as a consequence of human hunting and the progressive reduction of its habitat. To investigate in detail the genetic history of this species and to compare the population dynamics in different European areas, we analysed Bos primigenius remains from various sites across Italy.

Results

Fourteen samples provided ancient DNA fragments from the mitochondrial hypervariable region. Our data, jointly analysed with previously published sequences, support the view that Italian aurochsen were genetically similar to modern bovine breeds, but very different from northern/central European aurochsen. Bayesian analyses and coalescent simulations indicate that the genetic variation pattern in both Italian and northern/central European aurochsen is compatible with demographic stability after the last glaciation. We provide evidence that signatures of population expansion can erroneously arise in stable aurochsen populations when the different ages of the samples are not taken into account.

Conclusions

Distinct groups of aurochsen probably inhabited Italy and northern/central Europe after the last glaciation, respectively. On the contrary, Italian and Fertile Crescent aurochsen likely shared several mtDNA sequences, now common in modern breeds. We argue that a certain level of genetic homogeneity characterized aurochs populations in Southern Europe and the Middle East, and also that post-glacial recolonization of northern and central Europe advanced, without major demographic expansions, from eastern, and not southern, refugia.