Mitochondrial DNA variation and the evolutionary history of chromosome races of collared lemmings (Dicrostonyx) in the Eurasian Arctic

Collared lemmings (Dicrostonyx) demonstrate extensive chromosome variation along their circumpolar distribution in the high Arctic. To reveal the history of this genus and the origin of chromosome races in the Palearctic, we studied the geographical pattern of mtDNA variation in lemmings from 13 localities by using eight tetranucleotide restriction enzymes. The main split in mtDNA phylogeny is at the Bering Strait and corresponds to the main chromosome division between the Beringian and the Eurasian groups of karyotypes. Nucleotide divergence estimate of 6.8% suggests that, despite the Bering Land Bridge, Palearctic and Nearctic forms have been separated since the mid-Pleistocene. Five distinct phylogenetic groups of mtDNA haplotypes, with average divergence of 1.5%, corresponding to geographical regions, were found along the Palearctic coast. Low nucleotide and haplotype diversity and a star-like phylogeny within phylogeographical groups of haplotypes suggest regional bottleneck events in the recent past, most probably due to warming events during the Holocene. There is congruence between phylogeographical pattern of mtDNA variation and geographical distribution of chromosome races; 69% of the total mtDNA variation is allocated among chromosome races. This congruence implies that historical events such as fragmentation and allopatric bottleneck events have been important for the origin of chromosome races. However, historical factors do not explain the fixed autosome fusions found to distinguish certain populations.     

Glacial survival or late glacial colonization? Phylogeography of the root vole (Microtus oeconomus) in north‐west Norway

Aim It has been proposed that the root vole subspecies, Microtus oeconomus finmarchicus, survived the last glacial period on islands on the north‐west coast of Norway. The Norwegian island of Andøya may have constituted the only site with permanent ice‐free conditions. Geological surveys and fossil finds from Andøya demonstrate that survival throughout the last glacial maximum was probably possible for some plants and animals. In this study we aim to infer the recent evolutionary history of Norwegian root vole populations and to evaluate the glacial survival hypothesis. Methods DNA sequence variation in the mitochondrial cytochrome b gene was studied in 46 root voles from 19 localities. Location Northern Fennoscandia and north‐west Russia with a focus on islands on the north‐west coast of Norway. Results The phylogeographical analyses revealed two North European phylogroups labelled ‘Andøya’ and ‘Fennoscandia’. The Andøya phylogroup contained root voles from the Norwegian islands of Andøya, Ringvassøya and Reinøya and two localities in north‐west Russia. The Fennoscandian phylogroup encompassed root voles from the three Norwegian islands of Kvaløya, Håkøya and Arnøya and the remaining specimens from Norway, northern Sweden and Finland. Nucleotide diversity within the Andøya and Fennoscandian phylogroups was similar, ranging from 0.5% to 0.7%. Main conclusions Both our genetic data and previously published morphological data are consistent with in situ glacial survival of root voles on Andøya during the last glacial maximum. However, the level of genetic diversity observed in the extant island populations, the past periods of severe climatic conditions on Andøya and the ecology of the root vole are somewhat difficult to reconcile with this model. A biogeographical scenario involving late glacial recolonization along the northern coasts of Russia and Norway therefore represents a viable alternative. Our results demonstrate that complex recolonization and extinction histories can generate intricate phylogeographical patterns and relatively high levels of genetic variation in northern populations.     

mtDNA evidence for a local northern latitude Pleistocene refugium for the root vole (Microtus oeconomus,Arvicolinae, Rodentia) from Eastern Poland

We analysed the genetic structure of 33 populations of the root vole (tundra vole, Microtus oeconomus, Pallas, 1776) inhabiting their typical habitats, located at different distances from the southern boundary of the species’ range (52°14′–53°56′ N) in eastern Poland. We determined its phylogeographic pattern as well as the possible occurrence of a small, local high‐latitude refugium of this species in southern Poland, previously suggested in palaeontological studies. 908 bp of cytochrome b sequences were analysed from 439 root voles, and 21 mtDNA cytb haplotypes belonging to the Central European (CE) phylogroup were found. Haplotype diversity in the examined populations varied between 0 and 0.872 (mean: 0.425 ± 0.332), while nucleotide diversity ranged between 0 and 0.62% (mean: 0.235% ± 0.217). Within the CE phylogroup of M. oeconomus, we identified with high bootstrap support a newly separated group of M. oeconomus that evolved from CE, denoted CE‐PL S. This group is located in the southern and central part of eastern Poland and most likely diverged from phylogroup CE in a small, cryptic refugium situated in southern Poland, in the Kraków‐Cz?stochowa Upland and/or the Holy Cross Mountains during the LGM and Younger Dryas.     

Gene conversion in the mitochondrial genome on interspecific hybridization in voles of the <Emphasis Type="Italic">Clethrionomys</Emphasis> genus

The phenomenon of interspecific hybridization accompanied by transfer of the mitochondrial genome from the northern red-backed vole (Clethrionomys rutilus) to the bank vole (Cl. glareolus) in northeastern Europe is well known already for 25 years. However, the possibility of recombination between homologous segments of maternal and paternal mtDNAs of the voles during fertilization was not previously studied. Analysis of data on variability of nucleotide sequences of the mitochondrial gene for cytochrome b in populations of red-backed and bank voles in the area of their sympatry has shown that as a result of interspecific hybridization, the mitochondrial gene pool of bank voles contains not only mtDNA haplotypes of red-backed vole females, but also mtDNA haplotypes of bank voles bearing short nucleotide tracts of red-backed vole mtDNA. This finding supports the hypothesis that an incomplete elimination of red-backed vole paternal mtDNA during the interspecific hybridization between bank vole females and red-backed vole males leads to the gene conversion of bank vole maternal mtDNA tracts by homologous ones of mtDNA of red-backed vole males.     

Genetic structure and colonization processes in European populations of the common vole, Microtus arvalis

The level of genetic differentiation within and between evolutionary lineages of the common vole (Microtus arvalis) in Europe was examined by analyzing mitochondrial sequences from the control region (mtDNA) and 12 nuclear microsatellite loci (nucDNA) for 338 voles from 18 populations. The distribution of evolutionary lineages and the affinity of populations to lineages were determined with additional sequence data from the mitochondrial cytochrome b gene. Our analyses demonstrated very high levels of differentiation between populations (overall FST: mtDNA 70%; nucDNA 17%). The affinity of populations to evolutionary lineages was strongly reflected in mtDNA but not in nucDNA variation. Patterns of genetic structure for both markers visualized in synthetic genetic maps suggest a postglacial range expansion of the species into the Alps, as well as a potentially more ancient colonization from the northeast to the southwest of Europe. This expansion is supported by estimates for the divergence times between evolutionary lineages and within the western European lineage, which predate the last glacial maximum (LGM). Furthermore, all measures of genetic diversity within populations increased significantly with longitude and showed a trend toward increase with latitude. We conclude that the detected patterns are difficult to explain only by range expansions from separate LGM refugia close to the Mediterranean. This suggests that some M. arvalis populations persisted during the LGM in suitable habitat further north and that the gradients in genetic diversity may represent traces of a more ancient colonization of Europe by the species.     

Phylogeography of field voles (Microtus agrestis) in Eurasia inferred from mitochondrial DNA sequences

In a distribution-wide phylogeographic survey of the field vole (Microtus agrestis), 75 specimens from 56 localities across Eurasia were examined for DNA sequence variation along the whole 1140 base pair (bp) mitochondrial (mt) cytochrome b gene. The species is subdivided into three main mtDNA phylogeographic groups - western, eastern and southern - with largely allopatric distributions. The western phylogeographical group is found in west and central Europe and spread most probably from a glacial refugium in the Carpathians. The eastern group covers a large range from Lithuania to central Asia, and probably originated from a southeast European source (e.g. the southern Urals or the Caucasus). The southern group occupies an area from Portugal to Hungary, with division into two distinct mtDNA sublineages that presumably derive from separate glacial refugia in the Iberian Peninsula. Molecular clock estimates suggest that the western and eastern field vole populations separated during the last glaciation, whereas the southern population dates back 0.5-0.9 Myr. High levels of mtDNA variation indicate relatively large population sizes and subdivisions within phylogeographic groups during the last glaciation. We report a possible new suture zone in east Europe.     

Colonization history of north European field voles (Microtus agrestis) revealed by mitochondrial DNA

The genetic structure of field vole (Microtus agrestis) populations from northern Europe was examined by restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) in 150 individuals from 67 localities. A total of 83 haplotypes was observed, most of which were rare and highly localized geographically. Overall nucleotide diversity was high (134%), but showed a tendency to decrease with higher latitude. Two major mtDNA lineages differing by 2% in nucleotide sequence were identified. A southern mtDNA lineage was observed in field voles from Britain, Denmark and southern and central Sweden, whereas voles from Finland and northern Sweden belonged to a northern lineage. The strict phylogeographic pattern suggests that the present population generic structure in field voles reflects glacial history: the two groups are derived from different glacial refugia, and recolonized Fennoscandia from two directions. A 150–200-km-wide secondary contact zone between the two mtDNA groups was found in northern Sweden. Distinct phylogeographic substructuring was observed within both major mtDNA groups.     

Eastern Beringian biogeography: historical and spatial genetic structure of singing voles in Alaska

Aim  Pleistocene climatic cycles have left marked signatures in the spatial and historical genetic structure of high‐latitude organisms. We examine the mitochondrial (cytochrome b) genetic structure of the singing vole, Microtus miurus (Rodentia: Cricetidae: Arvicolinae), a member of the Pleistocene Beringian fauna, and of the insular vole, Microtus abbreviatus, its putative sister species found only on the St Matthew Archipelago. We reconstruct the phylogenetic and phylogeographical structure of these taxa, characterize their geographical partitioning and date coalescent and cladogenetic events in these species. Finally, we compare the recovered results with the phylogenetic, coalescent and spatial genetic patterns of other eastern Beringian mammals and high‐latitude arvicoline rodents. Location  Continental Alaska (alpine and arctic tundra) and the St Matthew Archipelago (Bering Sea). Methods  We generated and analysed cytochrome b sequences of 97 singing and insular voles (M. miurus and M. abbreviatus) from Alaska. Deep evolutionary structure was inferred by phylogenetic analysis using parsimony, maximum likelihood and Bayesian approaches; the geographical structure of genetic diversity was assessed using analysis of molecular variance and network analysis; ages of cladogenetic and coalescent events were estimated using a relaxed molecular clock model with Bayesian approximation. Results  Regional nucleotide diversity in singing voles is higher than in other high‐latitude arvicoline species, but intra‐population diversity is within the observed range of values for arvicolines. Microtus abbreviatus specimens are phylogenetically nested within M. miurus. Molecular divergence date estimates indicate that current genetic diversity was formed in the last glacial (Wisconsinan) and previous interglacial (Sangamonian) periods, with the exception of a Middle Pleistocene split found between samples collected in the Wrangell Mountains region and all other singing vole samples. Main conclusions  High levels of phylogenetic and spatial structure are observed among analysed populations. This pattern is consistent with that expected for a taxon with a long history in Beringia. The spatial genetic structure of continental singing voles differs in its northern and southern ranges, possibly reflecting differences in habitat distribution between arctic and alpine tundra. Our phylogenetic results support the taxonomic inclusion of M. miurus in its senior synonym, M. abbreviatus.     

A new cytochrome b phylogroup of the common vole (Microtus arvalis) endemic to the Balkans and its implications for the evolutionary history of the species

The phylogeographic architecture of the common vole, Microtus arvalis, has been well‐studied using mitochondrial DNA and used to test hypotheses relating to glacial refugia. The distribution of the five described cytochrome b (cyt b) lineages in Europe west of Russia has been interpreted as a consequence of postglacial expansion from both southern and central European refugia. A recently proposed competing model suggests that the ‘cradle’ of the M. arvalis lineages is in western central Europe from where they dispersed in different directions after the Last Glacial Maximum. In the present study, we report a new cyt b lineage of the common vole from the Balkans that is not closely related to any other lineage and whose presence might help resolve these issues of glacial refugia. The Balkan phylogroup occurs along the southern distributional border of M. arvalis in central and eastern Bosnia and Herzegovina, Montenegro, and eastern Serbia. Further north and west in Slovenia, Bosnia and Herzegovina, and Serbia, common voles belong to the previously‐described Eastern lineage, whereas both lineages are sympatric in one site in Bosnia and Herzegovina. The Balkan phylogroup most reasonably occupied a glacial refugium already known for various Balkan endemic species, in contrast to the recently proposed model. South‐east Europe is an absolutely crucial area for understanding the postglacial colonization history of small mammals in Europe and the present study adds to the very few previous detailed phylogeographic studies of this region. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 788–796.     

Dynamic effects of predators on cyclic voles: field experimentation and model extrapolation

Mechanisms generating the well-known 3-5 year cyclic fluctuations in densities of northern small rodents (voles and lemmings) have remained an ecological puzzle for decades. The hypothesis that these fluctuations are caused by delayed density-dependent impacts of predators was tested by replicated field experimentation in western Finland. We reduced densities of all main mammalian and avian predators through a 3 year vole cycle and compared vole abundances between four reduction and four control areas (each 2.5-3 km(2)). The reduction of predator densities increased the autumn density of voles fourfold in the low phase, accelerated the increase twofold, increased the autumn density of voles twofold in the peak phase, and retarded the initiation of decline of the vole cycle. Extrapolating these experimental results to their expected long-term dynamic effects through a demographic model produces changes from regular multiannual cycles to annual fluctuations with declining densities of specialist predators. This supports the findings of the field experiment and is in agreement with the predation hypothesis. We conclude that predators may indeed generate the cyclic population fluctuations of voles observed in northern Europe.     

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present‐day populations. In Europe, where both Mediterranean and extra‐Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra‐Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range‐wide colonization. The results provide clear evidence that both extra‐Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid‐latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%–47%) and the Carpathian (53%–72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end‐glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end‐glacial colonization history of this well‐studied organism and underscores the importance of genomic data in phylogeographic interpretation.     

Phylogeography of the common vole (Microtus arvalis) with particular emphasis on the colonization of the Orkney archipelago

To investigate the human introduction of the common vole Microtus arvalis onto the Orkney islands, the complete cytochrome b gene was sequenced in 41 specimens from both Orkney (four localities) and elsewhere in their range (26 localities). Orkney voles belonged to the same phylogenetic lineage, 'Western', as individuals from France and Spain indicating southwestern Europe as the most likely source area for the islands. This result is of interest with respect to the movement and trading links of the Neolithic people who likely transported the voles. As well as the Western lineage, our phylogenetic trees revealed three other purely European lineages: the 'Italian' (single specimen from N. Italy), the 'Central' (Germany, Netherlands, Denmark) and the 'Eastern' (Hungary, Slovakia, Poland, Ukraine, Finland, European Russia). Individuals from European Russia, W. Siberia, Georgia, Ukraine and Armenia formed a fifth distinct lineage coinciding with the distribution of the 'obscurus' chromosomal form of M. arvalis. These phylogeographical data suggest that M. arvalis occupied multiple refugia during the last glaciation.     

Transalpine colonisation and partial phylogeographic erosion by dispersal in the common vole (Microtus arvalis)

The colonisation history and genetic structure of the common vole ( Microtus arvalis ) was investigated in the region of the Alps by analysing the mitochondrial cytochrome b gene (mtDNA) and 19 microsatellite loci (nucDNA) for 137 voles from 52 localities. mtDNA data provided a much refined distribution of three highly divergent evolutionary lineages in the region compared to previous studies. Although high mountain ranges are widely accepted to be barriers for colonisation processes for many organisms and especially small terrestrial mammals, our phylogeographic analyses showed clear evidence of four transalpine colonisation events by the common vole. Individual-based phylogenetic analyses of nucDNA and two alternative Bayesian-clustering approaches revealed a deep genetic structure analogous to mtDNA. Incongruence between nucDNA and mtDNA at the individual level was restricted to the regions of contact between the lineages. mtDNA patterns and strong female philopatry in M. arvalis suggest that the crossings of the Alps occurred during the colonisation of the region when it was free from ice after the last glaciation. nucDNA patterns suggest that some of the transalpine elements of this phylogeographic pattern were subsequently eroded by male-biased gene flow. We conclude that the combination of phylogeography and landscape genetics at the individual level can provide very detailed insights into colonisation events and may even allow differentiation between historical and more recent processes.     

MITOCHONDRIAL DNA VARIATION IN THE FIELD VOLE (MICROTUS AGRESTIS): REGIONAL POPULATION STRUCTURE AND COLONIZATION HISTORY

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to examine the genetic structure among field voles (Microtus agrestis) from southern and central Sweden. A total of 57 haplotypes was identified in 158 voles from 60 localities. Overall mtDNA diversity was high, but both haplotype and nucleotide diversity exhibited pronounced geographic heterogeneity. Phylogenetic analyses revealed a shallow tree with seven primary mtDNA lineages separated by sequence divergences ranging from 0.6% to 1.0%. The geographic structure of mtDNA diversity and lineage distribution was complex but strongly structured and deviated significantly from an equilibrium situation. The extensive mtDNA diversity observed and the recent biogeographic history of the region suggests that the shallow mtDNA structure in the field vole cannot be explained solely by stochastic lineage sorting in situ or isolation by distance. Instead, the data suggest that the genetic imprints of historical demographic conditions and vicariant geographic events have been preserved and to a large extent determine the contemporary geographic distribution of mtDNA variation. A plausible historical scenario involves differentiation of mtDNA lineages in local populations in glacial refugia, a moving postglacial population structure, and bottlenecks and expansions of mtDNA lineages during the postglacial recolonization of Sweden. By combining the mtDNA data with an analysis of Y-chromosome variation, a specific population unit was identified in southwestern Sweden. This population, defined by a unique mtDNA lineage and fixation of a Y-chromosome variant, probably originated in a population bottleneck in southern Sweden about 12,000 to 13,000 calendar years ago.     

Cryptic variation in mole voles Ellobius (Arvicolinae,Rodentia) of Mongolia

The mole vole subgenus Ellobius is currently considered to include three species: Ellobius talpinus (distributed from SE Europe and Turkmenistan through Kazakhstan to SW Siberia), Ellobius alaicus (S Tianshan, Pamir-Alay) and Ellobius tancrei (East and West Central Asia, from the Amu-Darya to Mongolia and N China). A study focusing on the genetic variation in Ellobius from Mongolia was conducted using one mitochondrial and three nuclear markers. Two divergent allopatric lineages endemic to East Central Asia were revealed. The first lineage occurs from Dzungaria eastwards to central Mongolia and represents E. tancrei sensu stricto. The second lineage is found in East Gobi only and corresponds to a taxon described as Ellobius orientalis, which has been traditionally treated as a subspecies of E. tancrei. However, molecular and chromosome data indicate that orientalis is related not to E. tancrei but to E. talpinus, which is separated from the former by a distribution gap of ~2,000 km. The taxonomic status of the East Gobi mole vole is ambiguous, and its genetic distance from E. talpinus s. str. falls into the range characteristic for closely related vole species or semi-species. According to molecular estimates, the two taxa have been isolated since the late Middle Pleistocene. A similar divergence is observed between the East and West Central Asian lineages of E. tancrei. E. alaicus is placed as sister to the latter rendering E. tancrei sensu lato paraphyletic. The revealed phylogeographic pattern implies that East Central Asia was colonized by mole voles through multiple eastward dispersal events.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

Mitochondrial gene diversity in the common vole Microtus arvalis shaped by historical divergence and local adaptations

The phylogeography of the common vole (Microtus arvalis) was examined by analysing mitochondrial DNA (mtDNA) sequence variation in 1044 base pairs (bp) of the cytochrome b (cytb) gene and in 322 bp of the control region (ctr) among 106 individuals from 58 locations. The geographical distribution of four previously recognized cytb evolutionary lineages in Europe was refined and a new lineage was found in southern Germany. All lineages were distributed allopatrically, except in one sample that was probably located in a contact zone. The occurrence of several lineages in the Alps is in keeping with their recent recolonization from distinct sources. The translation of 84 cytb DNA sequences produced 33 distinct proteins with relationships that differed from those of the DNA haplotypes, suggesting that the mtDNA lineages did not diverge in response to selection. In comparison with M. agrestis, a neutrality test detected no overall evidence for selection in the cytb gene, but a closer examination of a structural model showed that evolutionarily conserved and functionally important positions were often affected. A new phylogeographical test of random accumulation of nonsynonymous mutations generated significant results in three lineages. We therefore conclude that the molecular diversity of cytb in M. arvalis is overall the result of the demographic history of the populations, but that there have been several episodes of local adaptation to peculiar environments.     

Hybridization between mtDNA-defined phylogeographic lineages of black ratsnakes (Pantherophis sp.)

Phylogeographic analyses using mitochondrial DNA (mtDNA) have revealed many examples of apparently deep historical subdivisions ('phylogroups') within many vertebrates. It remains unclear whether these phylogroups represent independently evolving, adaptively differentiated lineages or groups that show little functional differentiation and, hence, will merge on contact. Here, we use mtDNA sequence data to evaluate the phylogeographic relationships between two of the northernmost populations of black ratsnakes (Pantherophis obsoletus complex) in Ontario, Canada and previously analysed populations in the United States. We then use population-level analyses to evaluate the level of adaptive divergence between previously established mtDNA phylogroups. Phylogenetic analyses show that southern Ontario snakes have mtDNA haplotypes that fall within the Central mtDNA phylogroup, as designated by Burbrink et al. (2000). In contrast, snakes in eastern Ontario carry either Central or Eastern-specific haplotypes. Within the hybrid region, we found highly variable frequencies of mtDNA haplotypes among isolated sub-populations, no association between variation in cytonuclear (mtDNA) and nuclear (microsatellite DNA) markers, no difference in survival or reproductive success among snakes with different mtDNA haplotypes, and no effect of mate similarity in mtDNA on female clutch size. These results argue that the Eastern and Central phylogroups have merged in this region, likely due to a lack of adaptive differentiation between individuals in each lineage. Hence, in these snakes, phylogeographic structure in mtDNA is more a reflection of historical isolation rather than adaptive divergence. The observed reticulation between lineages and lack of evidence for hybrid disgenesis also bears on the classification of these lineages as distinct species.     

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine–Bair–Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.