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.

Introgressive hybridization is a widespread evolutionary phenomenon which may lead to increased allelic variation at selective neutral loci and to transfer of fitness-related traits to introgressed lineages. We inferred the population genetic structure of the European roe deer (Capreolus capreolus) in Poland from mitochondrial (CR and cyt b) and sex-linked markers (ZFX, SRY, DBY4 and DBY8). Analyses of CR mtDNA sequences from 452 individuals indicated widespread introgression of Siberian roe deer (C. pygargus) mtDNA in the European roe deer genome, 2000 km from the current distribution range of C. pygargus. Introgressed individuals constituted 16.6% of the deer studied. Nearly 75% of them possessed haplotypes belonging to the group which arose 23 kyr ago and have not been detected within the natural range of Siberian roe deer, indicating that majority of present introgression has ancient origin. Unlike the mtDNA results, sex-specific markers did not show signs of introgression. Species distribution modelling analyses suggested that C. pygargus could have extended its range as far west as Central Europe after last glacial maximum. The main hybridization event was probably associated with range expansion of the most abundant European roe deer lineage from western refugia and took place in Central Europe after the Younger Dryas (10.8-10.0 ka BP). Initially, introgressed mtDNA variants could have spread out on the wave of expansion through the mechanism of gene surfing, reaching high frequencies in European roe deer populations and leading to observed asymmetrical gene flow. Human-mediated introductions of C. pygargus had minimal effect on the extent of mtDNA introgression.

Large terrestrial carnivores are particularly prone to factors constraining levels of population genetic diversity because of their low densities and high spatial requirements. We studied the pattern of Eurasian lynx Lynx lynx population genetic variability in the westernmost part of its natural range from Scandinavia to the Carpathian Mountains (north-central Europe) based on 190 samples using 613 base pair-long sequences from the mitochondrial DNA control region (mtDNAcr). We examined whether the population history or contemporary habitat constraints of this large and mobile carnivore could have significantly affected its genetic structure. We recorded nine mtDNA haplotypes, including five not previously reported. Lynx from Latvia and Estonia had the highest variability with haplotype and nucleotide diversities of 0.81-0.88% and 0.44-0.47%, respectively. In contrast, there was no polymorphism present in peripheral populations from Norway and the Carpathian Mountains. Lynx populations were strongly differentiated [analysis of molecular variance (AMOVA): FST = 0.570, P < 0.001, FST = 0.464, P < 0.001)]. Spatial Analysis of Molecular Variance identified four separate groups of populations: (1) Norway, Finland and Estonia; (2) Latvia and North Eastern Poland; (3) the Białowieża Primeval Forest; (4) the Carpathians. The patterns of genetic diversity and differentiation suggested a number of discrete populations that are poorly connected by contemporary gene flow and could therefore be considered demographically independent. The peripheral location of these populations, habitat fragmentation and the strict territorial structure of lynx populations are factors likely contributing to the observed patterns. The study provides suggestions for active conservation/management decisions including translocations or reintroductions of lynx. bs_bs_banner Animal Conservation. Print ISSN 1367-9430 Animal Conservation 15 (2012) 603-612

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein-coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species-specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein-coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four-fold higher d N /d S ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well-preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.

The trajectories of postglacial range expansions, the occurrence of lineage patches and the formation and maintenance of secondary contact between lineages may mostly reflect neutral demographic processes, including density blocking, that may leave long-lasting genetic signatures. However, a few studies have recently shown that climate may also play a role. We used red deer, a large, mobile herbivore that is assumed to be sensitive to climate change, to test hypotheses of possible selection on the mitochondrial DNA cytochrome b gene (mtDNA cytb) and competitive and/or density-blocking (using mtDNA control region). We searched for a possible link between the phylogeographic structure and abiotic climatic variables. Finally, we tested for isolation by distance and isolation by environment and assessed the impact of human-mediated translocations on the genetic structure of red deer. Our analysis of 30 red deer populations in Poland using the mtDNA control region (N = 357) and cytochrome b (N = 50) markers not only confirmed the presence of the Western and South-Eastern lineages of the species but also indicated the presence of a previously unnoticed, rare relic haplotype that grouped together C. e. italicus from Italy (the Mesola deer). No significant signs of positive selection were detected for the mtDNA cytb gene in the studied red deer. However, a significant signal for purifying selection was found in our study that may explain the narrowness of the contact zone because gene flow between the Western and South-Eastern lineages should drive relatively strong mito-nuclear incompatibilities. MtDNA control region differentiation among red deer populations in Poland correlated with different abiotic climatic variables. Strikingly, the southernmost ice sheet limits during the Elsterian was the most important factor, and it explained the largest amount of variation. However, neither isolation by distance (IBD) nor isolation by environment (IBE) were recorded, and a very limited impact of human translocations was evident. The above-mentioned results suggest that in contemporary red deer populations in Poland, the phylogeographic pattern is well preserved, and long-term processes (density and/or competitive blocking) still play a major role.