Phylogeography of the forest-dwelling European pine marten (Martes martes): new insights into cryptic northern glacial refugia

The European Pine Marten Martes martes (Linnaeus, 1758) Is Autochthonous in Sicily and Constitutes a Well-Characterised Major Phylogroup within the Species (Carnivora, Mustelidae)

Simple Summary

The faunal assemblages currently occurring on islands are often a mélange of native and non-native species, and the actual status of some of them is doubtful at present. Since different laws are enforced for native and non-native species, sound knowledge about their status is pivotal for both their management and our understanding of the natural history of the studied areas. In the frame of this work, the Sicilian population of the European pine marten is genetically characterized for the first time, based on mitochondrial DNA sequences. Our results prove that the European pine marten is native on the island, where it arrived during the Pleistocene, and is represented there by a well-differentiated and endemic evolutionary lineage.

Abstract

No molecular data are currently available for the Sicilian populations of the European pine marten Martes martes, thus preventing any sound inference about its native or non-native status on the island, as well as the local phylogeography of the species. In order to investigate these issues, we sequenced two mtDNA markers in road-killed specimens collected in Sicily. Both markers consistently demonstrated the existence of a well-characterised Sicilian clade of the species, which is endemic to the island and constitutes the sister group of a clade including the Mediterranean and Central–North European major phylogroups of the European pine marten. Such evidence supports the autochthony of Martes martes in Sicily and points to a natural Pleistocene colonisation of the island followed by isolation. The occurrence of a, to date undetected, major phylogroup of the species in Sicily calls for the dedicated monitoring of the Sicilian populations of the species in order to preserve this evolutionarily significant unit.

Not out of the woods yet: genetic insights related to the recovery of the pine marten (Martes martes) in Ireland

In this study, the history of the pine marten (Martes martes) in Ireland is reviewed, revealing that the population has undergone several retractions and expansions over the last few hundred years. Here, we consider the genetic legacy of this flux in fortunes and its likely impacts upon the conservation and future recovery of the species. Using nuclear DNA markers (microsatellites), we found that the genetic diversity present in Ireland today is like that of other Irish carnivores, but there is evidence of a genetic bottleneck and low effective population size that might result in further reductions of diversity in the future. There is a lack of genetic structure, showing that the population has not been fragmented genetically, despite the low percentage of woodland in Ireland. We also reviewed the mitochondrial DNA diversity present in the Irish population and showed that there is only one contemporary and one extinct haplotype present; a reduced diversity relative to other Irish carnivores. The Irish haplotypes, both extant and extinct, are shared or are genetically similar to haplotypes commonly present in southern Europe today. We discuss the possibility of reinforcing the Irish population with animals from these sources to help supplement and maintain genetic diversity for future generations.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

Patterns of mtDNA variation reveal complex evolutionary history of relict and endangered peat bog pine (Pinus uliginosa)

Estimates of genetic differentiation at intra- and interspecific level are often hindered by the lack of suitable molecular markers. Low phylogeographic resolution limits development of appropriate conservation strategies especially in case of endangered forest tree species with small and disjunct distribution. In this study, we assessed fine-scale genetic structure of relict and endangered peat bog pine (Pinus uliginosa) and two other closely related European pine species (Pinus mugo and Pinus uncinata) using a set of 15 newly developed maternally inherited and seed-mediated mitochondrial DNA (mtDNA) markers and two previously known polymorphic mtDNA regions (nad1, nad7). Three main groups, corresponding in general to three investigated species were revealed in the haplotype network analysis. However, only P. uncinata was clearly distinct at all levels of analysis, whereas great genetic similarity and haplotype sharing was observed between P. uliginosa and P. mugo. Strong phylogeographic structure was found in P. uliginosa that showed high differentiation at relatively short geographical distance among populations and the existence of mitochondrial lineages of different evolutionary history. Hybridization with other pine species has likely contributed to genetic differentiation of P. uliginosa as indicated by contemporary distribution of mtDNA haplotypes. The research emphasizes the importance of accurate assessments of genetic structure of endangered species with complex evolutionary history for development of efficient conservation strategies.

PSMC analysis of effective population sizes in molecular ecology and its application to black-and-white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (Ne ) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in Ne over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black-and-white flycatchers. Ne curves of Atlas, collared, pied and semicollared flycatcher converged 1-2 million years ago (Ma) at an Ne of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Ne trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100-200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in Ne differed among species, and for pied flycatcher, the temporal dynamics of Ne differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per-site filter of ≥10 reads and no more than 25% of missing data.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

Revisiting the phylogeography and demography of European badgers (Meles meles) based on broad sampling, multiple markers and simulations

Although the phylogeography of European mammals has been extensively investigated since the 1990s, many studies were limited in terms of sampling distribution, the number of molecular markers used and the analytical techniques employed, frequently leading to incomplete postglacial recolonisation scenarios. The broad-scale genetic structure of the European badger (Meles meles) is of interest as it may result from historic restriction to glacial refugia and/or recent anthropogenic impact. However, previous studies were based mostly on samples from western Europe, making it difficult to draw robust conclusions about the location of refugia, patterns of postglacial expansion and recent demography. In the present study, continent-wide sampling and analyses with multiple markers provided evidence for two glacial refugia (Iberia and southeast Europe) that contributed to the genetic variation observed in badgers in Europe today. Approximate Bayesian computation provided support for a colonisation of Scandinavia from both Iberian and southeastern refugia. In the whole of Europe, we observed a decline in genetic diversity with increasing latitude, suggesting that the reduced diversity in the peripheral populations resulted from a postglacial expansion processes. Although MSVAR v.1.3 also provided evidence for recent genetic bottlenecks in some of these peripheral populations, the simulations performed to estimate the method's power to correctly infer the past demography of our empirical populations suggested that the timing and severity of bottlenecks could not be established with certainty. We urge caution against trying to relate demographic declines inferred using MSVAR with particular historic or climatological events.

Land-bridge calibration of molecular clocks and the post-glacial Colonization of Scandinavia by the Eurasian field vole Microtus agrestis

Phylogeography interprets molecular genetic variation in a spatial and temporal context. Molecular clocks are frequently used to calibrate phylogeographic analyses, however there is mounting evidence that molecular rates decay over the relevant timescales. It is therefore essential that an appropriate rate is determined, consistent with the temporal scale of the specific analysis. This can be achieved by using temporally spaced data such as ancient DNA or by relating the divergence of lineages directly to contemporaneous external events of known time. Here we calibrate a Eurasian field vole (Microtus agrestis) mitochondrial genealogy from the well-established series of post-glacial geophysical changes that led to the formation of the Baltic Sea and the separation of the Scandinavian peninsula from the central European mainland. The field vole exhibits the common phylogeographic pattern of Scandinavian colonization from both the north and the south, however the southernmost of the two relevant lineages appears to have originated in situ on the Scandinavian peninsula, or possibly in the adjacent island of Zealand, around the close of the Younger Dryas. The mitochondrial substitution rate and the timescale for the genealogy are closely consistent with those obtained with a previous calibration, based on the separation of the British Isles from mainland Europe. However the result here is arguably more certain, given the level of confidence that can be placed in one of the central assumptions of the calibration, that field voles could not survive the last glaciation of the southern part of the Scandinavian peninsula. Furthermore, the similarity between the molecular clock rate estimated here and those obtained by sampling heterochronous (ancient) DNA (including that of a congeneric species) suggest that there is little disparity between the measured genetic divergence and the population divergence that is implicit in our land-bridge calibration.

Large-scale genetic structuring of a widely distributed carnivore--the Eurasian lynx (Lynx lynx)

Over the last decades the phylogeography and genetic structure of a multitude of species inhabiting Europe and North America have been described. The flora and fauna of the vast landmasses of north-eastern Eurasia are still largely unexplored in this respect. The Eurasian lynx is a large felid that is relatively abundant over much of the Russian sub-continent and the adjoining countries. Analyzing 148 museum specimens collected throughout its range over the last 150 years we have described the large-scale genetic structuring in this highly mobile species. We have investigated the spatial genetic patterns using mitochondrial DNA sequences (D-loop and cytochrome b) and 11 microsatellite loci, and describe three phylogenetic clades and a clear structuring along an east-west gradient. The most likely scenario is that the contemporary Eurasian lynx populations originated in central Asia and that parts of Europe were inhabited by lynx during the Pleistocene. After the Last Glacial Maximum (LGM) range expansions lead to colonization of north-western Siberia and Scandinavia from the Caucasus and north-eastern Siberia from a refugium further east. No evidence of a Berinigan refugium could be detected in our data. We observed restricted gene flow and suggest that future studies of the Eurasian lynx explore to what extent the contemporary population structure may be explained by ecological variables.