Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Kamila Plis, Magdalena Niedziałkowska, Tomasz Borowik, Johannes Lang, Mike Heddergott, Juha Tiainen, Aleksey Bunevich, Nikica Šprem, Ladislav Paule, Aleksey Danilkin, Marina Kholodova, Elena Zvychaynaya, Nadezhda Kashinina, Boštjan Pokorny, Katarina Flajšman, Algimantas Paulauskas, Mihajla Djan, Zoran Ristić, Luboš Novák, Szilvia Kusza, Christine Miller, Dimitris Tsaparis, Stoyan Stoyanov, Maryna Shkvyria, Franz Suchentrunk, Miroslav Kutal, Vukan Lavadinović, Dragana Šnjegota, Ana-Maria Krapal, Gabriel Dănilă, Rauno Veeroja, Elżbieta Dulko & Bogumiła Jędrzejewska Title: Pan-European phylogeography of the European roe deer (Capreolus capreolus) Year: 2022 Version: Published version Copyright: The Author(s) 2022 Rights: CC BY 4.0 Rights url: http://creativecommons.org/licenses/by/4.0/ Please cite the original version: Plis, K., Niedziałkowska, M., Borowik, T., Lang, J., Heddergott, M., Tiainen, J., Bunevich, A., Šprem, N., Paule, L., Danilkin, A., Kholodova, M., Zvychaynaya, E., Kashinina, N., Pokorny, B., Flajšman, K., Paulauskas, A., Djan, M., Ristić, Z., Novák, L., … Jędrzejewska, B. (2022). Pan-European phylogeography of the European roe deer (Capreolus capreolus). Ecology and Evolution, 12, e8931. https://doi.org/10.1002/ece3.8931 Ecology and Evolution. 2022;12:e8931.  | 1 of 14 https://doi.org/10.1002/ece3.8931 www.ecolevol.org Received: 9 October 2021  | Revised: 13 April 2022  | Accepted: 25 April 2022 DOI: 10.1002/ece3.8931 R E S E A R C H A R T I C L E Pan- European phylogeography of the European roe deer (Capreolus capreolus) Kamila Plis1  | Magdalena Niedziałkowska1  | Tomasz Borowik1  | Johannes Lang2  | Mike Heddergott3  | Juha Tiainen4,5  | Aleksey Bunevich6 | Nikica Šprem7  | Ladislav Paule8  | Aleksey Danilkin9 | Marina Kholodova9 | Elena Zvychaynaya9 | Nadezhda Kashinina9 | Boštjan Pokorny10,11  | Katarina Flajšman11  | Algimantas Paulauskas12  | Mihajla Djan13  | Zoran Ristić14 | Luboš Novák15 | Szilvia Kusza16  | Christine Miller17 | Dimitris Tsaparis18  | Stoyan Stoyanov19  | Maryna Shkvyria20 | Franz Suchentrunk21 | Miroslav Kutal22 | Vukan Lavadinović23 | Dragana Šnjegota24  | Ana- Maria Krapal25 | Gabriel Dănilă26 | Rauno Veeroja27 | Elżbieta Dulko28,29 | Bogumiła Jędrzejewska1 1Mammal Research Institute, Polish Academy of Sciences, Białowieża, Poland 2Working Group for Wildlife Research at Clinic for Birds, Reptiles, Amphibians and Fish, Justus- Liebig- University Giessen, Gießen, Germany 3Department of Zoology, Musée National d’Histoire Naturelle, Luxembourg City, Luxembourg 4Lammi Biological Station, University of Helsinki, Lammi, Finland 5Natural Resources Institute Finland (Luke), Helsinki, Finland 6State National Park Belovezhskaya Pushcha, Kamenyuki, Republic of Belarus 7Department of Fisheries, Apiculture, Wildlife Management and Special Zoology, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia 8Department of Phytology, Technical University in Zvolen, Zvolen, Slovak Republic 9A.N. Severtsov Institute of Ecology and Evolution, Russian Academy of Sciences, Moscow, Russia 10Faculty of Environmental Protection, Velenje, Slovenia 11Department of Forest Ecology, Slovenian Forestry Institute, Ljubljana, Slovenia 12Department of Biology, Vytautas Magnus University, Kaunas, Lithuania 13Department of Biology and Ecology, Faculty of Sciences, University of Novi Sad, Novi Sad, Republic of Serbia 14Department of Geography, Tourism and Hotel Management, Faculty of Sciences, University of Novi Sad, Novi Sad, Serbia 15Department of Forest Protection and Wildlife Management, Mendel University in Brno, Brno, Czech Republic 16Centre for Agricultural Genomics and Biotechnology, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, Debrecen, Hungary 17Bureau of Wildlife Biology Bavaria, Rottach- Egern, Germany 18Institute of Marine Biology, Biotechnology and Aquaculture (IMBBC), Hellenic Centre for Marine Research, Heraklion, Crete, Greece 19Wildlife Management Department, University of Forestry, Sofia, Bulgaria 20Kyiv Zoological Park of National Importance, Kyiv, Ukraine 21Research Institute of Wildlife Ecology, University of Veterinary Medicine Vienna, Vienna, Austria 22Department of Forest Ecology, Faculty of Forestry and Wood Technology, Mendel University in Brno, Brno, Czech Republic 23Faculty of Forestry, University of Belgrade, Belgrade, Serbia 24Faculty of Natural Sciences and Mathematics, University of Banja Luka, Banja Luka, Bosnia and Herzegovina 25"Grigore Antipa" National Museum of Natural History, Bucharest, Romania 26Faculty of Forestry, Stefan cel Mare University of Suceava, Suceava, Romania This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. © 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 2 of 14  |     PLIS et aL. 1  |  INTRODUC TION Climate fluctuations during the Last Glacial Period between 115,000 and 11,500 years ago (Lokrantz & Sohlenius, 2006) played an im- portant role in shaping the current species composition, distribution, and genetic diversity of mammals in Europe. During glaciations, the ranges of cold sensitive species were limited to refugial areas located in the Balkan, the Iberian and the Apennine Peninsulas (Hewitt, 2004; Taberlet et al., 1998) and south- eastern part of the continent (the Black Sea region and the Caucasus Mts.; Markova & Puzachenko, 2019). Contribution of a given refugium into postglacial recoloniza- tion processes varied a lot among species such as, for example, red deer Cervus elaphus (Niedziałkowska, Doan, et al., 2021), wild boar Sus scrofa (Niedziałkowska, Tarnowska, et al., 2021), or common vole Microtus arvalis (Stojak et al., 2015). Cold- adapted species such as, for example, reindeer Rangifer tarandus (Sommer et al., 2014), saiga an- telope Saiga tatarica (Nadachowski et al., 2016), or arctic fox Alopex lagopus (Dalén et al., 2007), thrived in the glacial stages, expanded their ranges southward, and during the onset of interglacial they un- derwent massive, large- scale extinctions. Species of mammals that have very broad biogeographic niche (from the Mediterranean to the boreal zone) had more diverse response to geological time- scale pul- sation of climate, most likely with several glacial refugia located at both lower and higher latitudes, where subpopulations diverged into different lineages and could have developed adaptations to different climate, habitat, and food- related conditions. Examples of such spe- cies include the bank vole Clethrionomys glareolus (Tarnowska et al., 2019), weasel Mustela nivalis (McDevitt et al., 2012), common vole M. arvalis, and field vole Microtus agrestis (Baca et al., 2020; Stojak et al., 2019). Among European ungulates, the roe deer (Capreolus capreolus Linnaeus, 1758) (Figure 1) is the most widely distributed species (Lovari et al., 2016). It occurs throughout the continent (Figure 2), except for the northernmost Scandinavia and some islands. In south- ern Europe, roe deer lives in the Mediterranean zone, and the east- ern range of the species reaches the Caucasus Mountains, northern 27Department of Wildlife Monitoring, Estonian Environment Agency, Tallin, Estonia 28Department of Anesthesiology, University of Virginia Health System, Charlottesville, Virginia, USA 29Faculty of Biology, University of Warsaw, Warszawa, Poland Correspondence Kamila Plis, Mammal Research Institute, Polish Academy of Sciences, Stoczek 1, Białowieża 17- 230, Poland. Email: kplis@ibs.bialowieza.pl Funding information European Commission, Grant/Award Number: PIRSES- GA- 2009- 247652; Narodowe Centrum Nauki, Grant/Award Number: 2013/11/B/NZ8/00884 Abstract To provide the most comprehensive picture of species phylogeny and phylogeogra- phy 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 re- gions 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. K E Y W O R D S Capreolus capreolus, expansion, mitochondrial DNA, the Last Glacial Maximum refugia, the Quaternary history T A X O N O M Y C L A S S I F I C A T I O N Phylogenetics     |  3 of 14PLIS et aL. Iran, and Iraq (Danilkin, 1995). The eastern border of its range runs through western Russia. Roe deer is a habitat opportunist occupying both forests and open habitats (Apollonio et al., 2010; Geist, 1998). European roe deer is distinctly smaller that the Siberian roe deer (Capreolus pygargus (Pallas, 1771) occurring in Asia and in the south of the European part of Russia. The present ranges of the two spe- cies overlap in the Volga- Don region (Kashinina et al., 2018). In some areas, European roe deer is an important game species. In historical times, when the populations were locally exterminated or their numbers were low, there were attempts to strengthen them with introduced animals (Apollonio et al., 2014; Baker & Hoelzel, 2013). So far, phylogenetic studies on European roe deer covered a large part of their contemporary range and indicated three main mitochondrial DNA (mtDNA) clades— Central, Eastern, and Western (Baker & Hoelzel, 2014; Lorenzini et al., 2014; Randi et al., 2004; Tsaparis et al., 2019). The Central clade is the most common through- out Europe, the Eastern one is restricted mainly to the Balkans, and the Western clade occurs only in the Iberian Peninsula (Gentile et al., 2009; Lorenzini et al., 2014; Mucci et al., 2012; Randi et al., 2004). In addition, on the basis of internal structuring within clades, some researchers proposed the separation of the roe deer occurring in Italy as a distinct subspecies C. c. italicus (Gentile et al., 2009; Mucci et al., 2012). In Spain, there were also attempts to identify the Celtic- Iberian mtDNA clade (Baker & Hoelzel, 2014; Royo et al., 2007) and a subspecies C. c. garganta (Meunier, 1983). Recent studies showed the introgression of the Siberian roe deer mtDNA genes into the eastern populations of the European roe deer (Kashinina et al., 2018; Lorenzini et al., 2014; Markov et al., 2016; Matosiuk et al., 2014). This showed a new level of the genetic com- plexity in the species and opened a debate on the sources of that admixture (Kashinina et al., 2018; Lorenzini et al., 2014; Matosiuk et al., 2014; Świsłocka et al., 2019). However, the limited sampling (Lorenzini et al., 2014) did not allow to define whole spatial range at which the introgression has occurred. Despite multiple studies on European roe deer phylogeogra- phy, there still remain substantial gaps in the sampling coverage in Europe, which make the picture incomplete and thus inconclusive. There were no roe deer studies done in northern (Finland), eastern (Belarus, Estonia), central (Slovakia, Czech Republic), and south- eastern (Croatia, Slovenia) parts of the continent. The most compre- hensive phylogeographical analyses, done by Lorenzini et al. (2014) and Randi et al. (2004) focused mostly on identifying the origin of mtDNA clades and on disentangling the relationships among them. So far, there had been no attempts to determine the location of the contact zones between them. Moreover, the broad spatial spread of F I G U R E 1 Roe deer female with fawn. Photograph by Tibor Pataky F I G U R E 2 Roe deer (Capreolus capreolus) sampling locations in Europe: distribution of own and literature data in the species range 4 of 14  |     PLIS et aL. the clades suggests that, apart from separation of the C. c. italicus subspecies, there is some internal heterogeneity within the three main clades. In our study, we aimed at reconstructing the phylogeny of roe deer on the basis of intensive sampling in central and Eastern Europe, including many areas that have not been studied before. We pooled our new data on mtDNA sequences with already avail- able roe deer sequences in order to perform a more detailed and comprehensive analysis of roe deer phylogeography in the whole European range of the species (in 26 countries). The aims of our study, based on the analysis of a 610- bp fragment of mtDNA con- trol region, were to: (1) describe the genetic diversity of roe deer over the entire species range; (2) identify spatial patterns of the clade distribution; and (3) compare the demographic processes at the level of lineages, clades, and haplogroups. Finally, we discuss the possible LGM refugia and routes of postglacial recolonization of the continent by roe deer. 2  |  MATERIAL AND METHODS 2.1  |  Sample and data collection During 2008– 2017, we collected 1469 roe deer samples from 20 countries in central and Eastern Europe (Table 1, Figure 2). The study area ranged from Germany to Russia (6°35′– 43°23′E) and from Finland to Greece (67°42′– 38°44′N). Samples, which consisted of fresh fragments of skin or muscle from legally hunted animals, were in majority taken by hunters. Prior to DNA extraction, all samples were stored in 96% ethanol at – 20°C. We assigned the geographic coordinates of the samples on the basis of information on location provided by the hunters. Additionally, we included 1541 sequences of roe deer from western and southern Europe (see Table S1), based on data available in GenBank and their frequencies reported in pub- lications (Baker & Hoelzel, 2013; Biosa et al., 2015; Gentile et al., 2009; Lorenzini et al., 2014; Randi et al., 2004; Royo et al., 2007). 2.2  |  DNA extraction and sequencing We extracted total genomic DNA using the Qiagen DNeasy Blood and Tissue Kit following the manufacturer’s guidelines. A fragment of the mtDNA control region was amplified by PCR with the primers L- Pro and H- Phe (Randi et al., 1998). Cycling conditions were 95°C for 15 min; 35 cycles of 94°C for 15 s, 56°C for 15 s, and 72°C for 1 min; and a final step of 72°C for 10 min. PCR products were purified using Clean Up (A&A Biotechnology). Sequencing reactions were carried out in a 10 µl volume using the Big Dye sequencing kit v.3.1 (Applied Biosystems) with the forward primer. Products were purified with the Exterminator kit (A&A Biotechnology) and sequenced on an ABI 3130 xl Genetic Analyzer (Applied Biosystems). Sequencing results were analyzed with the ABI DNA Sequencing Analysis software and aligned in BioEdit v.7.0.5.3 (Hall, 1999). TA B L E 1 The number of roe deer Capreolus capreolus samples analyzed from the European countries Country Number of samples This study Literature Total Austria 16 6 22 Lorenzini et al. (2014) Belarus 135 – 135 – Bulgaria 42 – 42 – Croatia 25 – 25 – Czech Republic 58 – 58 – Denmark – 19 19 Lorenzini et al. (2014), Randi et al. (2004) Estonia 5 – 5 – Finland 106 – 106 – France – 24 24 Lorenzini et al. (2014), Randi et al. (2004) Germany 159 15 174 Randi et al. (2004) Greece 52 15 67 Lorenzini et al. (2014), Randi et al. (2004) Hungary 54 – 54 – Italy – 747 747 Biosa et al. (2015), Gentile et al. (2009), Lorenzini et al. (2014), Randi et al. (2004) Lithuania 16 13 29 Lorenzini et al. (2014) Norway 7 – 7 – Poland 382 8 390 Portugal – 23 23 Randi et al. (2004) Romania 15 10 25 Lorenzini et al. (2014) Russia 97 – 97 – Serbia 50 178 228 Randi et al. (2004) Slovakia 99 – 99 – Slovenia 65 – 65 – Spain – 156 156 Lorenzini et al. (2014), Randi et al. (2004) Sweden 17 11 28 Lorenzini et al. (2014), Randi et al. (2004) Ukraine 69 3 72 Lorenzini et al. (2014) United Kingdom – 313 313 Baker and Hoelzel (2013) Total 1469 1541 3010     |  5 of 14PLIS et aL. 2.3  |  Sequence and data analyses Sequencing resulted in good- quality mitochondrial control region fragments of 610 bp for all analyzed samples. Further, these frag- ments were aligned against the sequence of European roe deer from the north- eastern Poland (GenBank KM068161.1; Olano- Marin et al., 2014) and manually edited in BioEdit v.7.0.5.3 Hall (1999). Additional roe deer mtDNA data from GenBank were pooled with the obtained own fragments and shortened to keep the same length for all sequences. We assigned the obtained sequences to haplotypes using Arlequin 3.5.1.3 software (Excoffier & Lischer, 2010). Indels were considered as differences in haplotype definition. Each haplotype was assigned a code reflecting one of the two species: Cc— for the European roe deer (C. capreolus) or Cp— for the Siberian roe deer (C. pygargus). Internal structure of haplogroups was based on hap- lotype genealogy constructed in HapView (Salzburger et al., 2011). Italian haplogroup C. c. italicus (C7 in this study) has already been defined as such in the literature (Gentile et al., 2009; Mucci et al., 2012). To determine sequence evolution model, we used jModelTest2 (Darriba et al., 2012), using resources available on CIPRES Science Gateway (Miller et al., 2010). The best model chosen by jModelTest2 was the Hasegawa, Kishino, and Yano model (HKY; Hasegawa et al., 1985) with additional variation rate among sites and a proportion of invariable sites (+G, +I). Phylogenetic trees were constructed in Mega 7.0.14 (Kumar et al., 2016) and MrBayes 3 (Ronquist & Huelsenbeck, 2003) with the maximum likelihood method and 10,000 bootstrap replications, using the chosen model. Summary statistics were calculated in DnaSP 5.10.01 (Librado & Rozas, 2009) according to three classification levels: lineage, clade, and grouping based on HapView. The following statistics were computed: number of unique haplotypes (h), number of segregating (polymorphic) sites (S), haplotype diversity (Hd), nucleotide diversity (π), and average number of pairwise nucleotide differences (k). Additionally we in- cluded B index (Levins, 2020) to express the diversity of haplotypes, using the formula: where pi is the proportion of samples with haplotype i in a deme. The B index minimum value is 1 and its upper bound is equal to the number of haplotypes in the sample. To evaluate possible models of expansion, we performed two neutrality tests in DnaSP: Tajima's D (Tajima, 1989) and Fu's Fs B = 1 ∑ p2 i TA B L E 2 Estimates of genetic diversity of mtDNA control region (610 bp) in the roe deer (C. capreolus) lineages, clades, and haplogroups (see Figure 3). Number of samples is indicated by n (numbers in parentheses correspond to samples collected and analyzed in this study), h— number of haplotypes, S— number of polymorphic sites, Hd— haplotype diversity, π— nucleotide diversity, k— average number of pairwise differences. Haplotypes CcH53 and CcH55 (European lineage) were not assigned to any clade Lineage, clade, haplogroup n h S Hd π k Siberian lineage 266 (257) 28 41 0.831 0.010 5.804 European lineage 2744 (1212) 300 84 0.982 0.011 6.401 Eastern clade 566 (375) 76 46 0.954 0.006 3.783 E1 77 (67) 6 5 0.174 0.000 0.180 E2 160 (60) 14 15 0.816 0.004 2.370 E3 140 (92) 22 16 0.900 0.004 2.243 E4 187 (154) 33 30 0.905 0.006 3.602 Central clade 2030 (795) 206 68 0.972 0.008 4.740 C1 737 (374) 83 42 0.951 0.006 3.817 C2 350 (207) 34 22 0.847 0.004 2.442 C3 124 (0) 5 4 0.124 0.000 0.127 C4 388 (139) 46 29 0.871 0.007 3.978 C5 12 (9) 3 6 0.439 0.002 1.379 C6 40 (33) 12 11 0.892 0.005 2.915 C7 328 (0) 10 11 0.639 0.002 1.341 C8 49 (31) 12 10 0.844 0.004 2.459 Western clade 148 (42) 18 21 0.901 0.008 4.802 W1 73 (36) 9 14 0.804 0.008 4.846 W2 75 (6) 9 11 0.798 0.004 2.708 CcH53 2 (2) 1 – – – – CcH55 2 (2) 1 – – – – Total 3010 (1469) 328 95 0.984 0.015 9.049 6 of 14  |     PLIS et aL. (Fu, 1997). We also used mismatch distribution with the sudden expansion model and goodness- of- fit tests (sum of squared devia- tion; Harpending's raggedness index (R). Expansion time (T) was es- timated based on equation: T = τ /2µ, where τ is Tau estimated in Arlequin 3.5.1.3 and µ is the mutation rate described as units of sub- stitutions per locus per generation (Rogers & Harpending, 1992). We applied a mutation rate of 0.04– 0.08 substitutions per site per Myr (Lorenzini et al., 2014; Randi et al., 1998, 2004; Royo et al., 2007; Vernesi et al., 2002) and used 3 years as a generation time (Randi et al., 1998, 2004). Time since expansion was calculated using the Excel Spreadsheet provided by Schenekar and Weiss (2011). The Last Glacial Period had impact on large areas of the northern hemisphere and its names differ regionally. Throughout this manu- script we used the name “Weichselian” as reference, which is equiv- alent to the other terms such as “Vistulian” (van Loon et al., 2012), “Wurm” or “Devensian” (Kolfschoten, 2000). The reconstruction of the potential LGM refugial range of European roe deer during that period was based on fossil and subfossil records of the species pro- vided by Sommer et al. (2009) and Markova and Puzachenko (2019), supplemented with habitat availability data (Markova & Puzachenko, 2019). 3  |  RESULTS 3.1  |  Phylogenetics of roe deer Analysis of 3010 roe deer sequences of a fragment of the mtDNA control region (610 bp) revealed 328 haplotypes (Table 2) with 95 polymorphic sites (83 transitions, 18 transversions, and 4 indels). The numbers of samples representing each haplotype ranged from 1 to 177; 106 sequences were singletons (see Table S1). We iden- tified 200 new haplotypes, not described in literature so far. The obtained sequences were deposited in GenBank under accession numbers ON368373– ON368700. Phylogenetic analyses revealed two major mtDNA lineages assigned to two species: the European roe deer C. capreolus (Cc) and the Siberian roe deer C. pygargus (Cp) (Figure 3; Figures S1 and S2). The majority of haplotypes (297 out of 328; 90.5%) grouped with the European roe deer lineage (Table 2). We detected 73 transitions (80.2%), 15 transversions (16.5%), and 3 indels (3.3%) in this group. Fifty- seven mutation sites were polymor- phic in Cc, but monomorphic in Cp. Roe deer exhibited high values of mtDNA diversity in both lineages: haplotype diversity (Hd) was 0.982 in Cc and 0.831 in Cp, and nucleotide diversity (π) 0.11 and 0.10, respectively. The European lineage showed a clear division into three major mtDNA clades (Figure 3; see Figures S1– S3). Following the names previously defined by Randi et al. (2004), we classified them as Eastern, Western, and Central clades. We distinguished a sepa- rate haplogroup for the roe deer subspecies C. c. italicus within the Central clade. All of them showed a wide range of haplotype diver- sity (largely proportional to the number of analyzed samples), and a relatively low nucleotide diversity (Table 2). The Central clade was most numerous: 74% of analyzed specimens and 67% of identified haplotypes in the European lineage belonged to it (Table 2). The Western clade, in turn, was least numerous, but it had the same level of nucleotide diversity as the Central one (π = 0.008). The Western clade had the highest average number of pairwise differences among all clades (Table 2), which suggests its nonmonophyletic origin and reflects the distribution of the clade and haplogroups on phyloge- netic networks (Figure 3 and Figure S3). The Eastern clade had inter- mediate values in all estimates. Further analysis of the haplotype genealogy of European roe deer lineage indicated the presence of an internal structuring in the de- scribed clades (Figure 3; Figure S3). Only haplotypes CcH55 (2 samples from Czech Republic) and CcH53 (1 sample from Slovakia and 1 from Serbia) were not assigned to any specific group. Inside the Eastern clade, we detected four haplogroups E1– E4 (Figure 3). The most nu- merous was haplogroup E4 (n = 187), while E1 was the least numerous (n = 77). Out of all the haplotypes that formed the haplogroup E1, only one (CcH220) was previously described (see Table S1). F I G U R E 3 Internal structure of European roe deer (C. capreolus) clades defined based on mtDNA haplotype genealogy constructed in HapView program. Circles represent unique haplotypes, while their sizes correspond to the number of individuals with a given haplotype. Small undescribed circles denote missing haplotypes. Names of the clades (see Figures S1 and S2) were assigned according to the naming proposed by Randi et al. (2004). Details on haplotypes in Table S1     |  7 of 14PLIS et aL. The Central clade consisted of eight haplogroups (C1– C8), among which C1 had the highest number of individuals (n = 2030) and haplotypes (h = 206). The smallest recorded haplogroup was C5 (Table 2, Figure 3). Haplogroups C3 and C7 consisted of haplotypes found only in the published data. Those two haplogroups had the lowest average number of pairwise differences (0.127 and 1.341, respectively) and low nucleotide diversity (<0.001 and 0.002). All haplotypes grouped in haplogroup C7 were previously described in roe deer belonging to the subspecies C. c. italicus (comp. Randi et al., 2004). The Western clade consisted of two haplogroups (W1 and W2) represented by similar numbers of individuals (Table 2). The hap- logroup W1 had the highest average number of pairwise differences among all detected clades and haplogroups (k = 4.846) and the high- est value of nucleotide diversity among all haplogroups (π = 0.008). The Western clade had many missing haplotypes (Figure 3 and Figure S3). 3.2  |  Spatial distribution of lineages, clades, and haplogroups The spatial distribution of roe deer lineages in Europe showed a strong geographical pattern. Haplotypes of the Siberian line- age were recorded in western Russia, Finland, Estonia, Belarus, Lithuania, Poland, Ukraine, Slovakia, Romania, and Hungary (Figure S4). Their range fully overlapped with that of the European line- age. The highest proportion of the Siberian haplotypes was found in south- eastern and eastern parts of roe deer geographical range (populations 14 and 13; 57 and 53% of samples, respectively), and it declined westwards (Figure 4a). In the European lineage, we observed the spatial structuring of the three major clades (Figure S4). The Central clade was the most widespread and covered large part of the continent, with the high- est frequencies in western, northern, and central regions (100% in populations 2– 4 and 11, and 95% in population 5) (Figure 4a). This was the only roe deer clade which occurred in Great Britain and in the Scandinavian Peninsula. Some haplogroups of the Central clade were widely distributed (C1, C2, and especially C4), while others were geographically limited: C3— to Great Britain, C5— to the coastal parts of the Black Sea, C7— to the Apennine Peninsula (Figure 4b and Figure S5). The ranges of haplogroups C6 and C8 were restricted to east- central Europe. The Eastern clade was recorded in southern and east- central Europe between 8° and 39°E longitude and below 56°N latitude (Figure S4). In Greece and Bulgaria, 99% of roe deer belonged to it, and the share of that clade declined northward (Figure 4a). The haplogroup E1 occurred only in Belarus, Lithuania, and the adjoining borderlands of Poland, Russia, and Ukraine (Figure 4b and Figure S6). Haplogroups E2– E4 considerably overlapped across the study area, yet some differ- ences in their geographical ranges could still be observed (Figure 4b and Figure S6). For instance, individuals assigned to E3 were mainly located in the Balkan- Dinaric and the Carpathian regions. The Western clade was least numerous, but it was distributed over large areas of continental Europe from the Iberian Peninsula to western Russia (Figure S4 and Figure S7). None of the haplotypes from this clade was present north of 50°N. Except for Spain and Portugal (population 1), where it was recorded in 21% of roe deer samples, the Western clade contributed minimally (0.3– 6%) to roe deer populations in central and eastern Europe (Figure 4a). Despite the low number of samples, inconsistency between the genetic po- sition of the Western clade in roe deer phylogenies (Figure 3 and Figure S3), supported the division of this clade into two haplogroups (W1 and W2). Although haplogroups W1 and W2 showed substan- tial spatial overlap, their distributions were rather patchy (Figure 4b and Figure S4). Groups of related haplotypes in each of those hap- logroups were endemic to some regions of Europe (Figure 4b and Figure S7). We calculated the diversity index B of haplogroups belong- ing to the European lineage of roe deer in the 14 studied popu- lations (Figure 4b). B index could vary from 1 to 14 (total number of haplogroups in the Eastern, Central, and Western clades). The highest diversity of phylogenetic haplogroups was found in the F I G U R E 4 Haplotype frequencies of roe deer mtDNA groups in the defined 14 populations (see Figure 2) of European roe deer. Upper panel (a) depicts contributions of the Siberian lineage and the three main clades (Central, Eastern, and Western) in the European lineage in the 14 populations. Lower panel (b) shows the proportions of all haplogroups identified in the European lineage (see Figure 2). Numbers in italics are B indices indicating effective haplotype diversity 8 of 14  |     PLIS et aL. central part of roe deer range (populations 6, 9, 10, and 12; B from 3.7 to 5.9). All peripheral populations (except for population 1 in the Iberian Peninsula) showed low diversity (B from 1.3 to 2.8; Figure 4b). 3.3  |  Demographic history of roe deer The goodness- of- fit tests, which compared expansion model with the observed mismatch distribution among individuals, showed evidence for expansion processes in all clades and nearly all haplo- groups in the European lineage and in the Siberian lineage as well (Figure 5 and Figures S8 and S9). The observed distribution of pair- wise differences did not deviate from the distribution expected under the expansion model. Additionally, values of Tajima's D and Fu's FS tests confirmed the recent expansion of haplogroups C3 and E1 (see Table S2). The only haplogroups not showing expan- sion tendency were C8 and W2, while C4 had only a weak support for expansion (see Table S2). The Central and Eastern clades ex- hibited strong unimodal distributions of the pairwise differences, which indicated one main expansion event in each case (Figure 5). The Western clade (Figure 5) and the Siberian lineage (see Figure S8) had multimodal distributions. Expansion times were calculated for the spatial and demographic models based on two assumed mutation rates (see Methods). Therefore, the results— with very wide confidence intervals— can only give a broad view of the demographic history of European roe deer. Most of haplogroups of the Central and the Western clades showed expansion times predating the LGM (see Table S3). Possibly, only C2, C3, and C5 expanded during or after the LGM. The oldest expansion time, perhaps reaching as far as the Eemian interglacial (<100 ka BP), was found inhaplogroupW1. Three haplogroups in the Eastern clade (E1, E2, E3) showed the expansion time around or after the LGM (see Table S3). 4  |  DISCUSSION We presented the comprehensive picture of the genetic diversity of roe deer— the most numerous, widespread, and ecologically flex- ible ungulate species of Europe. We identified the phylogeny within the species, mapped the spatial distribution of lineages, clades, and haplogroups, and attempted to describe their demographic history. By pooling together a large set of new data from central and eastern Europe and available literature data, we showed a clear, continent- wide pattern of phylogeography of roe deer populations spanning from the Iberian Peninsula to western and south- western Russia. However, there are still two important regions where data are lack- ing. The first is nearly the whole France (except for its southern part), and the second are the easternmost south- eastern verges of roe deer range, including the Caucasus Mts., Turkey, and south- eastern regions of the European part of Russia, where the range of C. capreolus overlaps with that of C. pygargus. Future research, filling those gaps in data coverage and analysis, cannot only broaden the pattern we presented in this work, but also verify some of our con- clusions and interpretations. The general picture of roe deer phylogeny and phylogeography we obtained in our study fully confirmed earlier findings (Lorenzini et al., 2014; Olano- Marin et al., 2014; Randi et al., 2004; Royo et al., 2007; Sommer et al., 2009; Vernesi et al., 2002) showing two lin- eages (the European and the Siberian) of roe deer mtDNA in Europe, with the European lineage consisting of three clades (the Western, the Central, and the Eastern). Clear spatial phylogenetic patterns indicated, that translocations recorded in the historical times had rather minor impact and did not disturb the overall picture of the roe F I G U R E 5 Mismatch distribution. Graphs represent the mismatch distribution for the three main European roe deer clades (Central, Eastern, and Western). Dashed lines correspond to the frequency expected based on the sudden expansion model     |  9 of 14PLIS et aL. deer phylogeography in the continental scale. The major novelty of our study comes from the fact that we were able to divide the highly differentiated clades of the European lineage into well- defined hap- logroups, to map their ranges of occurrence, and quantify the spatial patterns of their frequencies in roe deer populations. Despite the sampling bias across the different populations and the inevitable inaccuracies in delimiting the phylogenetic haplogroups (marked as discrepancies in their branching pattern and haplotype membership in different phylogenetic trees and networks), the broad division of clades in the European lineage into smaller phylogenetic units ap- peared very informative for understanding the spatiotemporal his- tory of the species. 4.1  |  The Siberian lineage of roe deer mtDNA The introgression of mtDNA genes of the Siberian roe deer into Capreolus capreolus populations has already been detected in Hungary (Markov et al., 2016, 2017), Poland (Lorenzini et al., 2014; Matosiuk et al., 2014; Olano- Marin et al., 2014), Lithuania (Lorenzini et al., 2014), Belarus (A. Danilkin, personal communication), Russia (Kashinina et al., 2018; Plakhina et al., 2014), and Ukraine (Danilkin et al., 2017; Kashinina et al., 2018). We found that the spatial range of the introgression is much wider than previously reported (it occurred also among roe deer in Finland, Estonia, Slovakia, and Romania), and the share of the Siberian mtDNA lineage in local pop- ulations declined from east to west. Despite a vast range of occur- rence, the Siberian lineage was not very rich in haplotypes (28), yet it had a diversified phylogeny and most likely a complicated demo- graphic history, which requires further studies. The origin and sources of the introgression are not yet fully rec- ognized. Some authors stated that it had been caused by an ancient hybridization during the postglacial range shifts of the European and the Siberian roe deer with subsequent range overlapping (Lorenzini et al., 2014; Matosiuk et al., 2014; Świsłocka et al., 2019). Conversely, Danilkin et al. (2017) and Kashinina et al. (2018) argued that the sources of hybridization were numerous relocations of the Siberian roe deer into the range of C. capreolus in the 19th and 20th centuries (Apollonio et al., 2014). Those two sources of introgression are most likely nonexclusive, as the recent human- caused hybridization could have supplemented the ancient natural process. The problem can be elucidated by a future comparative study on the modern Siberian lineage of European roe deer and the modern and ancient phylog- eny and phylogeography of the Siberian roe deer from their whole geographic range. 4.2  |  The European lineage of roe deer mtDNA— changes in time and space The European roe deer was already present in Europe at least 600 ka BP and the large number of fossil and subfossil records evidenced that the species occurred in both glacial and interglacial phases since then (Sommer et al., 2009). Radiocarbon (C14) dating of roe deer fos- sil bones (a method that reaches only as far as 50– 60 ka BP) showed that between 60 and 21 ka BP, roe deer occurred not only in the Mediterranean peninsulas, but also it repeatedly reached central Europe during milder interstadials (Sommer et al., 2009). Rapid cli- matic changes from colder to warmer periods were conducive to recurrent expansion– retreat cycles of roe deer and thus their ge- netic diversification. Indeed, most of haplogroups of the Central and Eastern clades and W2 of the Western clade could have originated in the Weichselian Pleniglacial, thus long before the Last Glacial Maximum. Based on records of fossil and subfossil bones collated by Sommer et al. (2009), supplemented with data from Markova and Puzachenko (2019) for the eastern part of the continent, we reconstructed the LGM range of the European roe deer (Figure 6a). According to the synthesis of palynological data, Markova et al. (2008) proposed that roe deer dwelled in the Mediterranean and the Caucasian forests, and the southern variant of the periglacial forest- steppe. Stefaniak (2015) reported the fossil bones of roe deer during the LGM from the Caucasus Mts. Thus, the species range during LGM spanned throughout southern Europe from the Iberian Peninsula and south- ern France to the eastern ranges of the Caucasus and embraced the southern part of the Carpathian region (Figure 6a). Further research on genetic diversity of roe deer in the Caucasus Mts. region is badly needed to elucidate the phylogeography of the species in the south- eastern verges of its range. Below, we discuss the contemporary phylogeographic structure of roe deer and reconstruct the most likely phylogeography of the species during the LGM and its post- glacial expansion. The Central clade is by far the most widespread, numerous, and diversified among roe deer mtDNA clades. It occurs from Portugal, Spain, and Great Britain in western Europe to the Scandinavian Peninsula in the north, and the eastern shores of the Black Sea in the east. However, the Central clade is rare in the Balkan Peninsula, where its single, southernmost records were found in Bulgaria. Among eight haplogroups, C4 had the widest range (see Figure S5) and the oldest expansion time, however, with a wide confidence interval. Thus, we suggest that it might be a relic older than the Weichselian Glaciation. C4 might have then survived the LGM in southern France, northern Italy, and the north- eastern shores of the Black Sea (Figure 6a), the areas indicated as suitable for the species during the LGM (Markova & Puzachenko, 2019). Its expan- sion could have started from those regions after the glacial period. Haplogroups C1 and C2 now are common in western and central Europe; we propose that the refugial areas of C1 were in south- ern France, the Dinaric Mts., and the eastern Carpathian region (Moldova), while those of C2— in the Iberian peninsula and south- ern France. In addition to the earlier described Italian haplogroup, denoted as C7 in this paper, that occurred in the Apennine Peninsula only (Lorenzini et al., 2002; Mucci et al., 2012; Randi et al., 2004), we found three other haplogroups, each with a very restricted range: C5 re- corded over the northern side of the Black Sea, C6 found mainly in 10 of 14  |     PLIS et aL. Slovenia, Czech Republic and Slovakia (from Eastern Alps to Western Carpathians), and C3 found in Great Britain, only. The isolated hap- logroups C5, C6, and C7 most likely did not expand from their LGM refugia (Figure 6a). Our haplogroup C3 is equivalent to the group, which supposed to represent the remnant Scottish population de- scribed by Baker and Hoelzel (2014) from northern Great Britain. The LGM refuge of the haplogroup C3 could have been located in present- day France, near the land bridge existing in the past between north- ern France and southern England. During postglacial period, C3 roe deer colonized Great Britain, but— as proposed by Baker and Hoelzel (2014)— the subsequent human pressure led to the extinction of those ancient populations in the southern part of the British Isles and their replacement by introduced individuals. Today, the haplotypes belong- ing to haplogroup C3 from the primary postglacial colonization are still preserved in the northern populations (Baker & Hoelzel, 2014). The distribution of the haplogroup C8 is similar to the localities of the southern spatially isolated haplogroups: the core of its range is located in the northwest of the Balkan Peninsula (in the Dinaric Mts.), which might have been its LGM refugium (Figure 6a). The Eastern clade contained three haplogroups (E2– E4), which together constituted nearly 100% of roe deer samples in the south- eastern Balkan Peninsula (the sampling region no 8: Bulgaria and Greece) and their share in local roe deer populations declined to- ward north. Obviously, the Balkan Peninsula was the glacial refuge of haplogroups E2, E3, and E4. The situation of an isolated haplogroup E1 is different. It now occurs largely in regions that had been ice covered during the LGM (Lithuania, Belarus, western Russia), and its spatial expansion times suggested recent immigration. E1 consists of basically one haplotype with its few derivatives differing by one mutation step. We suggest that it survived the LGM in a northern refuge located at high latitude on the Belarussian- Russian border- land (Figure 6a), where boreal coniferous trees were documented to grow during the LGM (Svenning et al., 2008; Tzedakis et al., 2013). Heikkilä et al. (2009) documented the presence of Betula, Pinus, and Picea populations in the eastern Baltic region (Latvia) already in the Late Glacial. This suggested that the full- glacial locations of those trees were in the regions south and south- east of Latvia, possibly at 55– 60°N, making the area habitable even for large mammals. Interestingly, the pan- European study on molecular biogeography of wild boar S. scrofa (Niedziałkowska, Tarnowska, et al., 2021) discov- ered a rare, ancient clade of the species in the Russian part of the range of roe deer haplogroup E1. Also, the Europe- wide study of wolf Canis lupus population by SNP analysis (Stronen et al., 2013) found that the western and northern Belarus wolves showed the most divergent genotypes within north- central Europe. The Western clade, discovered for the first time in the western part of the continent, appeared “Western” no longer, as we found it also in central and eastern Europe, where it formed up to 5.3% of roe deer in the sampling region no 13 (Belarus, Lithuania, and western Russia). We decided to retain the name “Western” for this group, as it has already been well established in literature. The internal phylogeny of the Western clade (many missing, possibly extinct haplotypes), discontinuity of its occurrence over the con- tinent, spatial segregation of groups of haplotypes, and very old expansion times strongly suggest that the Western clade has been a relic of the Eemian interglacial (MIS 5) 130– 80 ka BP, when it probably had occurred throughout the continent. This could also explain why results from different reconstructions (see Figure 3 and Figure S3) do not seem coherent. Many missing haplotypes could have blurred the position of the Western clade in the phy- logeographic networks. The present occurrence and diversifica- tion of the Western clade suggest at least three glacial refugia: the Iberian Peninsula, southern Alps, and the Carpathian region sensu lato. In addition, the haplogroup W1 had survived the LGM in F I G U R E 6 Reconstruction of the LGM refugial range of the European roe deer (a) and the most probable postglacial recolonization routes by the major haplogroups in the Eastern (b) and the Central clade (c). The LGM range reconstructed after Sommer et al. (2009), Stefaniak (2015), Markova and Puzachenko (2019), with modification over the Carpathians according to deciduous and coniferous tree range in the LGM after Tzedakis et al. (2013). Symbols of haplogroups as in Table 2 and Figure 3. RD— roe deer refugial population of unknown genetic profile. Cp— the Siberian lineage expansion into the European roe deer. Solid lines in (b) and (c) show the present ranges of occurrence of the haplogroups (see Figures S5 and S7)     |  11 of 14PLIS et aL. the north- eastern refuge located in the present- day Belarussian- Russian borderland (Figure 6a). 4.3  |  Postglacial colonization of Europe by roe deer Our data suggest that different haplogroups of roe deer contrib- uted differently to the postglacial colonization of the European continent. Among 14 haplogroups belonging to three mtDNA clades, seven played the major roles in recolonization of the post- glacial habitats. These were, from the most successful: C2, C1, C4, E4, E2, E3, and C8 (Figure 6a,b). The other seven haplogroups showed rather weak expansion or most probably remained in their LGM refugia. Also, different regions of the large, glacial refugial area could have contributed unevenly to the modern range of roe deer. The southern France refugial population might have expanded to cover the whole western and north- western Europe, and possibly reached also central and eastern parts of the continent (Figure 6c). Then, the Balkan refugial area, especially its eastern part, gave rise to the vast range of central and eastern Europe being repopulated by the Eastern clade of roe deer (Figure 6b), whereas colonizers from the Dinaric Mts. region (mainly C8) dispersed northward in a fairly narrow belt, probably between the western and eastern expanding populations. The glacial populations of roe deer from the southernmost and easternmost refugial regions (the Iberian and the Apennine Peninsulas, the eastern Carpathians and the north- eastern Black Sea shore) did not seem to expand much in the postglacial pe- riod. While the Iberian and Apennine roe deer were located be- hind the expanding populations, the expansion from the eastern Carpathians and the Black Sea region might have been halted by the Siberian lineage of roe deer expanding from the East (Figure 6b,c). The study focused on the Siberian lineage is needed to understand the mechanisms of that expansion and its influence on the genetic profile of European roe deer populations in the eastern part of the species range. 5  |  CONCLUSIONS Our data suggest that the present phylogeographic pattern of roe deer in Europe might have been shaped by the Weichselian pleniglacial diversification of haplogroups, survival of the Eemian relic clade and other clades in multiple refugial areas, the postgla- cial recolonization, and hybridization with the Siberian roe deer. Modern relocations and exploitation by humans probably did not significantly affect the observed phylogeographical patterns and had an impact only in a local scale. The proposed recolonization patterns and hybridization with C. pygargus need further genetic studies, including also ancient samples of the roe deer, to reveal with more details the evolutionary history of those two sister spe- cies in Eurasia. AUTHOR CONTRIBUTIONS Kamila Plis: Data curation (lead); formal analysis (lead); visualization (equal); writing – original draft (lead). Magdalena Niedziałkowska: Conceptualization (equal); formal analysis (supporting); methodology (supporting); writing – original draft (supporting); writing – review and editing (lead). Tomasz Borowik: Data curation (supporting); writing – original draft (supporting). Johannes Lang: Data curation (supporting); Writing – original draft (supporting); writing – review and editing (supporting). Mike Heddergott: Data curation (equal); Writing – original draft (supporting); writing – review and edit- ing (equal). Juha Tiainen: Data curation (equal); writing – original draft (supporting); writing – review and editing (equal). Aleksey Bunevich: Data curation (equal). Nikica Šprem: Data curation (equal); writing – original draft (supporting); writing – review and editing (equal). Ladislav Paule: Data curation (equal); Supervision (equal); writing – original draft (equal); writing – review and edit- ing (equal). Aleksey Danilkin: Conceptualization (supporting); data curation (lead); writing – original draft (equal); writing – review and editing (equal). Marina Kholodova: Data curation (supporting); Writing – review and editing (equal). Elena Zvychaynaya: Data curation (supporting); writing – review and editing (supporting). Nadezhda Kashinina: Data curation (supporting); writing – review and editing (supporting). Boštjan Pokorny: Data curation (equal); writing – original draft (supporting); writing – review and editing (equal). Katarina Flajšman: Data curation (equal); writing – original draft (supporting); writing – review and editing (equal). Algimantas Paulauskas: Conceptualization (supporting); data curation (equal); writing – original draft (supporting); writing – review and edit- ing (equal). Mihajla Djan: Data curation (equal); writing – original draft (supporting); writing – review and editing (supporting). Zoran Ristić: Data curation (equal); writing – review and editing (sup- porting). Luboš Novák: Data curation (equal); writing – review and editing (supporting). Szilvia Kusza: Data curation (equal); writing – original draft (supporting); writing – review and editing (support- ing). Christine Miller: Data curation (equal); writing – original draft (supporting); writing – review and editing (supporting). Dimitris Tsaparis: Data curation (supporting); writing – original draft (equal); writing – review and editing (equal). Stoyan Stoyanov: Data curation (equal); writing – original draft (equal); writing – re- view and editing (equal). Maryna Shkvyria: Data curation (equal); writing – original draft (equal); writing – review and editing (equal). Franz Suchentrunk: Conceptualization (supporting); data curation (equal); writing – original draft (equal); writing – review and edit- ing (equal). Miroslav Kutal: Data curation (equal); writing – origi- nal draft (supporting); writing – review and editing (supporting). Vukan Lavadinović: Data curation (supporting); writing – review and editing (supporting). Dragana Šnjegota: Data curation (equal); formal analysis (supporting); writing – original draft (equal); writ- ing – review and editing (equal). Ana- Maria Krapal: Data curation (supporting); investigation (supporting); writing – original draft (supporting); writing – review and editing (supporting). Gabriel Dănilă: Data curation (supporting); writing – original draft (sup- porting); writing – review and editing (supporting). Rauno Veeroja: 12 of 14  |     PLIS et aL. Data curation (supporting); writing – original draft (supporting); writing – review and editing (supporting). Elżbieta Dulko: Data cu- ration (supporting); formal analysis (equal); writing – original draft (supporting); writing – review and editing (supporting). Bogumiła Jędrzejewska: Conceptualization (lead); funding acquisition (lead); investigation (equal); supervision (lead); visualization (equal); writ- ing – original draft (lead); writing – review and editing (lead). ACKNOWLEDG MENTS This study was financed by the National Science Centre in Poland under project No. 2013/11/B/NZ8/00884 for BJ. This study forms part of the PhD thesis of KP. The research was supported by the European Commission’s project BIOGEAST No. PIRSES- GA- 2009- 247652. The authors thank to Polish Hunting Association, J. Žák, K. Turek, M. Mati, V. Sidarovich, I. Dykiy, D. Vishnewsky, M. Kolesnikov, and V. Lobkov and numerous hunters for their help in collecting samples. They also thank the reviewers for their helpful comments, which have improved this manuscript. CONFLIC T OF INTERE S T The authors declare that they have no conflict of interest. DATA AVAIL ABILIT Y S TATEMENT The authors confirm that all data underlying the findings are fully available without restriction. The data has been submitted to dryad with the Data Identifier https://doi.org/10.5061/dryad.76hdr 7t01. Complete mtDNA sequences have been deposited in GenBank with accession numbers ON368373– ON368700. ORCID Kamila Plis https://orcid.org/0000-0002-7462-5324 Magdalena Niedziałkowska https://orcid. org/0000-0002-7746-4116 Tomasz Borowik https://orcid.org/0000-0003-0321-8200 Johannes Lang https://orcid.org/0000-0002-7387-795X Mike Heddergott https://orcid.org/0000-0003-4536-5508 Juha Tiainen https://orcid.org/0000-0001-5855-2594 Nikica Šprem https://orcid.org/0000-0002-3475-6653 Ladislav Paule https://orcid.org/0000-0002-5097-7472 Boštjan Pokorny https://orcid.org/0000-0002-2981-5516 Katarina Flajšman https://orcid.org/0000-0002-5583-653X Algimantas Paulauskas https://orcid.org/0000-0002-6823-3754 Mihajla Djan https://orcid.org/0000-0002-2427-0676 Szilvia Kusza https://orcid.org/0000-0002-5441-5303 Dimitris Tsaparis https://orcid.org/0000-0001-9626-5553 Stoyan Stoyanov https://orcid.org/0000-0003-2775-1800 Dragana Šnjegota https://orcid.org/0000-0003-2125-4635 Bogumiła Jędrzejewska https://orcid.org/0000-0002-9652-0953 R E FE R E N C E S Apollonio, M., Andersen, R., & Putman, R. (2010). European ungulates and their management in the 21st century. Cambridge University Press. ISBN: 978- 0- 521- 76061- 4. Apollonio, M., Scandura, M., & Šprem, N. (2014). Reintroductions as a management tool for European ungulates. In R. Putman & M. Apollonio (Eds.), Behaviour and management of European ungulates (pp. 46– 77). Whittles Publishing. Baca, M., Popović, D., Baca, K., Lemanik, A., Doan, K., Horáček, I., López- García, J. M., Bañuls- Cardona, S., Pazonyi, P., Desclaux, E., Crégut- Bonnoure, E., Berto, C., Lenardić, J. M., Miękina, B., Murelaga, X., Cuenca- Bescós, G., Krajcarz, M., Marković, Z., Petculescu, A., … Nadachowski, A. (2020). Diverse responses of common vole (Microtus arvalis) populations to Late Glacial and Early Holocene climate changes – Evidence from ancient DNA. Quaternary Science Reviews, 233, 106239. https://doi.org/10.1016/j.quasc irev.2020.106239 Baker, K. H., & Hoelzel, A. R. (2013). Evolution of population genetic structure of the British roe deer by natural and anthropogenic processes (Capreolus capreolus). Ecology and Evolution, 3, 89– 102. https://doi.org/10.1002/ece3.430 Baker, K. H., & Hoelzel, A. R. (2014). Influence of Holocene environmen- tal change and anthropogenic impact on the diversity and distribu- tion of roe deer. Heredity, 112, 607– 615. https://doi.org/10.1038/ hdy.2013.142 Biosa, D., Scandura, M., Tagliavini, J., Luccarini, S., Mattioli, L., & Apollonio, M. (2015). Patterns of genetic admixture between roe deer of different origin in central Italy. Journal of Mammalogy, 96, 827– 838. https://doi.org/10.1093/jmamm al/gyv098 Dalén, L., Nyström, V., Valdiosera, C., Germonpré, M., Sablin, M., Turner, E., Angerbjörn, A., Arsuaga, J. L., & Götherström, A. (2007). Ancient DNA reveals lack of postglacial habitat tracking in the arctic fox. Proceedings of the National Academy of Sciences of the United States of America, 104, 6726– 6729. https://doi.org/10.1073/pnas.07013 41104 Danilkin, A. (1995). Behavioural ecology of Siberian and European roe deer. Springer. Danilkin, A. A., Plakhina, D. A., Zvychaynaya, E. Y., Domnich, A. V., Kholodova, M. V., Sorokin, P. A., & Volokh, A. M. (2017). Siberian roe deer (Capreolus pygargus Pallas, 1771) in Ukraine: Analysis of the mitochondrial and nuclear DNA. Biology Bulletin, 44(6), 575– 583. https://doi.org/10.1134/S1062 35901 706005X Darriba, D., Taboada, G. L., Doallo, R., & Posada, D. (2012). jModelTest 2: More models, new heuristics and parallel computing. Nature Methods, 9, 772. https://doi.org/10.1038/nmeth.2109 Excoffier, L., & Lischer, H. E. L. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10, 564– 567. https:// doi.org/10.1111/j.1755- 0998.2010.02847.x Fu, Y. X. (1997). Statistical tests of neutrality of mutations against popu- lation growth, hitchhiking and background selection. Genetics, 147, 915– 925. https://doi.org/10.1093/genet ics/147.2.915 Geist, V. (1998). Deer of the world: Their evolution, behaviour and ecology. Stackpole Books. Gentile, G., Vernesi, C., Vicario, S., Pecchioli, E., Caccone, A., Bertorelle, G., & Sbordoni, V. (2009). Mitochondrial DNA variation in roe deer (Capreolus capreolus) from Italy: Evidence of admixture in one of the last C. c. italicus pure populations from central- southern Italy. Italian Journal of Zoology, 76, 16– 27. https://doi.org/10.1080/11250 00080 2018725 Hall, T. A. (1999). BioEdit: A user- friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series, 41, 95– 98. https://doi.org/10.14601/ Phyto pathol_Medit err- 14998 u1.29 Hasegawa, M., Kishino, H., & Yano, T. (1985). Dating of the human- ape splitting by a molecular clock of mitochondrial DNA. Journal of Molecular Evolution, 22, 160– 174. https://doi.org/10.1007/BF021 01694 Heikkilä, M., Fontana, S. L., & Seppä, H. (2009). Rapid Late glacial tree population dynamics and ecosystem changes in the eastern Baltic     |  13 of 14PLIS et aL. region. Journal of Quaternary Science, 24, 802– 815. https://doi. org/10.1002/jqs.1254 Hewitt, G. M. (2004). Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions of the Royal Society of London B: Biological Sciences, 359, 183– 195. https://doi.org/10.1098/ rstb.2003.1388 Kashinina, N. V., Danilkin, A. A., Zvychaynaya, E. Y., Kholodova, M. V., & Kiryakulov, V. M. (2018). On the gene pool of roe deer (Capreolus) of Eastern Europe: Analysis of the cyt b gene sequence variability. Russian Journal of Genetics, 54, 825– 831. https://doi.org/10.1134/ S1022 79541 8070049 Kumar, S., Stecher, G., & Tamura, K. (2016). MEGA7: Molecular evolu- tionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution, 33(7), 1870– 1874. https://doi.org/10.1093/ molbe v/msw054 Levins, R. (2020). Evolution in changing environments. https://doi. org/10.2307/j.ctvx5 wbbh Librado, P., & Rozas, J. (2009). DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25(11), 1451– 1452. https://doi.org/10.1093/bioin forma tics/btp187 Lorenzini, R., Garofalo, L., Qin, X., Voloshina, I., & Lovari, S. (2014). Global phylogeography of the genus Capreolus (Artiodactyla: Cervidae), a Palaearctic meso- mammal. Zoological Journal of the Linnean Society, 170, 209– 221. https://doi.org/10.1111/ zoj.12091 Lorenzini, R., Lovari, S., & Masseti, M. (2002). The rediscovery of the Italian roe deer: Genetic differentiation and management im- plications. Italian Journal of Zoology, 69, 367– 379. https://doi. org/10.1080/11250 00020 9356482 Lorenzini, R., Lovari, S., & Masseti, M. (2006). Ice marginal fluctuations during the Weichselian glaciation in Fennoscandia, a literature re- view. Swedish Nuclear Fuel and Waste Management Co. Lovari, S., Herrero, J., Masseti, M., Ambarli, H., Lorenzini, R., & Giannatos, G. (2016). Capreolus capreolus. The IUCN Red List of Threatened Species 2016. e.T42395A22161386. https://doi.org/10.2305/ IUCN.UK.2016- 1.RLTS.T4239 5A221 61386.en. Downloaded on 07 October 2021. Markov, G., Zvychaynaya, E., Danilkin, A., Kholodova, M., & Sugar, L. (2016). Genetic diversity and phylogeography of roe deer (Capreolus capreolus L.) in different biogeographical regions in Europe. Comptes Rendus De L Academie Bulgare Des Sciences, 69, 579– 584. Markov, G. G., Zvychaynaya, E. Y., Danilkin, A., Kholodova, M. V., Sugár, L., & Dimitrov, H. (2017). Indication for genetic diversity of European Roe Deer Capreolus capreolus (L.) in southeastern Europe revealed by mt DNA markers. Acta Zoologica Bulgarica, Suppl, 8, 61– 67. Markova, A., Kolfschoten, T., Bohnke, S., Kosinsev, P. A., Mol, J., Puzachenko, A., Simakova, A. N., Smirnov, N., Verpoorte, A., & Golovachev, I. V. (2008). Evolution of European ecosystems during pleistocene- Holocene Transition (24– 8 Kyr BP). KMK. Markova, A., & Puzachenko, A. (2019). Mammal assemblages during the last glacial maximum (LGM) (≤24 - ≥17 kyr BP). In A. Markova, T. Kolfschoten, & A. Puzachenko (Eds.), Evolution of European ecosys- tems during pleistocene - holocene transition (24– 8 kyr BP) (pp. 44– 56). GEOS Press. Matosiuk, M., Borkowska, A., Świsłocka, M., Mirski, P., Borowski, Z., Krysiuk, K., Danilkin, A. A., Zvychaynaya, E. Y., Saveljev, A. P., & Ratkiewicz, M. (2014). Unexpected population genetic structure of European roe deer in Poland: An invasion of the mtDNA genome from Siberian roe deer. Molecular Ecology, 23, 2559– 2572. https:// doi.org/10.1111/mec.12745 McDevitt, A. D., Zub, K., Kawałko, A., Oliver, M. K., Herman, J. S., & Wójcik, J. M. (2012). Climate and refugial origin influence the mito- chondrial lineage distribution of weasels (Mustela nivalis) in a phy- logeographic suture zone. Biological Journal of the Linnean Society, 106, 57– 69. https://doi.org/10.1111/j.1095- 8312.2012.01840.x Meunier, K. (1983). Das spanisches Reh. In R. R. Hoffman (Ed.), Wildbiologische Informationen für den Jäger (pp. 120– 131). Berlin. Miller, M. A., Pfeiffer, W., & Schwartz, T. (2010). Creating the CIPRES Science Gateway for inference of large phylogenetic trees. Gateway Computing Environments Workshop. https://doi.org/10.1109/ GCE.2010.5676129 Mucci, N., Mattucci, F., & Randi, E. (2012). Conservation of threat- ened local gene pools: Landscape genetics of the Italian roe deer (Capreolus c. italicus) populations. Evolutionary Ecology Research, 14, 897– 920. Nadachowski, A., Lipecki, G., Ratajczak, U., Stefaniak, K., & Wojtal, P. (2016). Dispersal events of the saiga antelope (Saiga tatarica) in Central Europe in response to the climatic fluctuations in MIS 2 and the early part of MIS 1. Quaternary International, 420, 357– 362. https://doi.org/10.1016/j.quaint.2015.11.068 Niedziałkowska, M., Doan, K., Górny, M., Sykut, M., Stefaniak, K., Piotrowska, N., Jędrzejewska, B., Ridush, B., Pawełczyk, S., Mackiewicz, P., Schmölcke, U., Kosintsev, P., Makowiecki, D., Charniauski, M., Krasnodębski, D., Rannamäe, E., Saarma, U., Arakelyan, M., Manaseryan, N., … Stanković, A. (2021). Winter temperature and forest cover have shaped red deer distribution in Europe and the Ural Mountains since the Late Pleistocene. Journal of Biogeography, 48(1), 147– 159. https://doi.org/10.1111/jbi.13989 Niedziałkowska, M., Tarnowska, E., Ligmanowska, J., Jędrzejewska, B., Podgórski, T., Radziszewska, A., Ratajczyk, I., Kusza, S., Bunevich, A. N., Danila, G., Shkvyria, M., Grzybowski, T., & Woźniak, M. (2021). Clear phylogeographic pattern and genetic structure of wild boar Sus scrofa population in Central and Eastern Europe. Scientific Reports, 11, 9680. https://doi.org/10.1038/s4159 8- 021- 88991 - 1 Olano- Marin, J., Plis, K., Sönnichsen, L., Borowik, T., Niedziałkowska, M., & Jędrzejewska, B. (2014). Weak population structure in European roe deer (Capreolus capreolus) and evidence of introgressive hy- bridization with Siberian roe deer (C. pygargus) in northeastern Poland. PLoS One, 9, e109147. https://doi.org/10.1371/journ al.pone.0109147 Plakhina, D. A., Zvychaynaya, E. Y., Kholodova, M. V., & Danilkin, A. A. (2014). Identification of European (Capreolus capreolus L.) and Siberian (C. pygargus Pall.) roe deer hybrids by microsatellite marker analysis. Russian Journal of Genetics, 50, 757– 762. https://doi. org/10.1134/S1022 79541 4070151 Randi, E., Alves, P. C., Carranza, J., Milošević- Zlatanović, S., Sfougaris, A., & Mucci, N. (2004). Phylogeography of roe deer (Capreolus ca- preolus) populations: the effects of historical genetic subdivisions and recent nonequilibrium dynamics. Molecular Ecology, 13, 3071– 3083. https://doi.org/10.1111/j.1365- 294x.2004.02279.x Randi, E., Pierpaoli, M., & Danilkin, A. (1998). Mitochondrial DNA polymorphism in populations of Siberian and European roe deer (Capreolus pygargus and C. capreolus). Heredity, 80, 429– 437. https://doi.org/10.1046/j.1365- 2540.1998.00318.x Rogers, A. R., & Harpending, H. (1992). Population growth makes waves in the distribution of pairwise genetic differences. Molecular Biology and Evolution, 9, 552– 569. https://doi.org/10.1093/oxfor djour nals. molbev.a040727 Ronquist, F., & Huelsenbeck, J. P. (2003). MrBayes 3: Bayesian phyloge- netic inference under mixed models. Bioinformatics, 19, 1572– 1574. https://doi.org/10.1093/bioin forma tics/btg180 Royo, L. J., Pajares, G., Álvarez, I., Fernández, I., & Goyache, F. (2007). Genetic variability and differentiation in Spanish roe deer (Capreolus capreolus): A phylogeographic reassessment within the European framework. Molecular Phylogenetics and Evolution, 42, 47– 61. https://doi.org/10.1016/j.ympev.2006.05.020 Salzburger, W., Ewing, G. B., & Haeseler, A. V. (2011). The performance of phylogenetic algorithms in estimating haplotype genealogies with migration. Molecular Ecology, 20, 1952– 1963. https://doi. org/10.1111/j.1365- 294x.2011.05066.x 14 of 14  |     PLIS et aL. Schenekar, T., & Weiss, S. (2011). High rate of calculation errors in mis- match distribution analysis results in numerous false inferences of biological importance. Heredity, 107, 511– 512. https://doi. org/10.1038/hdy.2011.48 Sommer, R. S., Fahlke, J. M., Schmölcke, U., Benecke, N., & Zachos, F. E. (2009). Quaternary history of the European roe deer Capreolus capreolus. Mammal Review, 39, 1– 16. https://doi. org/10.1111/j.1365- 2907.2008.00137.x Sommer, R. S., Kalbe, J., Ekström, J., Benecke, N., & Liljegren, R. (2014). Range dynamics of the reindeer in Europe during the last 25,000 years. Journal of Biogeography, 41, 298– 306. https://doi. org/10.1111/jbi.12193 Stefaniak, K. (2015). Neogene and quaternary cervidae from Poland. Habilitation thesis, Kraków, Poland: Institute of Systematics and Evolution of Animals Polish Academy of Sciences. Stojak, J., Borowik, T., Górny, M., McDevitt, A. D., & Wójcik, J. M. (2019). Climatic influences on the genetic structure and distribution of the common vole and field vole in Europe. Mammal Research, 64, 19– 29. https://doi.org/10.1007/s1336 4- 018- 0395- 8 Stojak, J., McDevitt, A. D., Herman, J. S., Searle, J. B., & Wójcik, J. M. (2015). Postglacial colonization of eastern Europe from the Carpathian refugium: Evidence from mitochondrial DNA of the common vole Microtus arvalis. Biological Journal of the Linnean Society, 115, 927– 939. Stronen, A. V., Jędrzejewska, B., Pertoldi, C., Demontis, D., Randi, E., Niedziałkowska, M., Pilot, M., Sidorovich, V. E., Dykyy, I., Kusak, J., Tsingarska, E., Kojola, I., Karamanlidis, A. A., Ornicans, A., Lobkov, V. A., Dumenko, V., & Czarnomska, S. D. (2013). North- south dif- ferentiation and a region of high diversity in European wolves (Canis lupus). PLoS One, 8, e76454. https://doi.org/10.1371/journ al.pone.0076454 Svenning, J.- C., Normand, S., & Kageyama, M. (2008). Glacial refu- gia of temperate trees in Europe: insights from species distri- bution modelling. Journal of Ecology, 96, 1117– 1127. https://doi. org/10.1111/j.1365- 2745.2008.01422.x Świsłocka, M., Czajkowska, M., Matosiuk, M., Saveljev, A. P., Ratkiewicz, M., & Borkowska, A. (2019). No evidence for recent introgressive hybridization between the European and Siberian roe deer in Poland. Mammalian Biology, 97, 59– 63. https://doi.org/10.1016/j. mambio.2019.04.008 Taberlet, P., Fumagalli, L., Wust- Saucy, A. G., & Cosson, J.- F. (1998). Comparative phylogeography and postglacial colonization routes in Europe. Molecular Ecology, 7, 453– 464. https://doi. org/10.1046/j.1365- 294x.1998.00289.x Tajima, F. (1989). The effect of change in population size on DNA poly- morphism. Genetics, 123, 597– 601. https://doi.org/10.1093/genet ics/123.3.597 Tarnowska, E., Niedziałkowska, M., & Jędrzejewska, B. (2019). Genetic structure of bank vole populations in the contact zone of two lin- eages in north- eastern Poland. Mammalian Biology, 96, 93– 101. https://doi.org/10.1016/j.mambio.2018.10.011 Tsaparis, D., Sotiropoulos, K., Legakis, A., Kotoulas, G., & Kasapidis, P. (2019). New phylogeographic insights support the distinctiveness and conservation value of the little- known Greek roe deer popu- lations. Mammalian Biology, 96, 23– 27. https://doi.org/10.1016/j. mambio.2019.03.010 Tzedakis, P. C., Emerson, B. C., & Hewitt, G. M. (2013). Cryptic or mystic? Glacial tree refugia in northern Europe. Trends in Ecology & Evolution, 28, 696– 704. https://doi.org/10.1016/j.tree.2013.09.001 van Kolfschoten, T. (2000). The Eemian mammal fauna of central Europe. Netherlands Journal of Geosciences, 79, 269– 281. https://doi. org/10.1017/S0016 77460 0021752 van Loon, A. J., Błaszkiewicz, M., & Degórski, M. (2012). The role of permafrost in shaping the Late Glacial relief of northern Poland. Netherlands Journal of Geosciences, 91, 223– 231. https://doi. org/10.1017/S0016 77460 000161X Vernesi, C., Pecchioli, E., Caramelli, D., Tiedemann, R., Randi, E., & Bertorelle, G. (2002). The genetic structure of natural and reintro- duced roe deer (Capreolus capreolus) populations in the Alps and central Italy, with reference to the mitochondrial DNA phyloge- ography of Europe. Molecular Ecology, 11, 1285– 1297. https://doi. org/10.1046/j.1365- 294X.2002.01534.x SUPPORTING INFORMATION Additional supporting information may be found in the online version of the article at the publisher’s website. How to cite this article: Plis, K., Niedziałkowska, M., Borowik, T., Lang, J., Heddergott, M., Tiainen, J., Bunevich, A., Šprem, N., Paule, L., Danilkin, A., Kholodova, M., Zvychaynaya, E., Kashinina, N., Pokorny, B., Flajšman, K., Paulauskas, A., Djan, M., Ristić, Z., Novák, L., … Jędrzejewska, B. (2022). Pan- European phylogeography of the European roe deer (Capreolus capreolus). Ecology and Evolution, 12, e8931. https://doi.org/10.1002/ece3.8931