Genetic divergence and phylogeographic relationships among european perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization

We used the widely distributed freshwater fish, perch (Perca fluviatilis), to investigate the postglacial colonization routes of freshwater fishes in Europe. Genetic variability within and among drainages was assessed using mitochondrial DNA (mtDNA) D-loop sequencing and RAPD markers from 55 populations all over Europe as well as one Siberian population. High level of structuring for both markers was observed among drainages and regions, while little differentiation was seen within drainages and regions. Phylogeographic relationships among European perch were determined from the distribution of 35 mtDNA haplotypes detected in the samples. In addition to a distinct southern European group, which includes a Greek and a southern Danubian population, three major groups of perch are observed: the western European drainages, the eastern European drainages including the Siberian population, and Norwegian populations from northern Norway, and western side of Oslofjord. Our data suggest that present perch populations in western and northern Europe were colonized from three main refugia, located in southeastern, northeastern and western Europe. In support of this, nested cladistic analysis of mtDNA clade and nested clade distances suggested historical range expansion as the main factor determining geographical distribution of haplotypes. The Baltic Sea has been colonized from all three refugia, and northeastern Europe harbours descendants from both eastern European refugia. In the upper part of the Danube lineages from the western European and the southern European refugia meet. The southern European refugium probably did not contribute to the recolonization of other western and northern European drainages after the last glaciation. However, phylogenetic analyses suggest that the southern European mtDNA lineage is the most ancient, and therefore likely to be the founder of all present perch lineages. The colonization routes used by perch probably also apply to other freshwater species with similar distribution patterns.     

Unique mitochondrial DNA lineages in Irish stickleback populations: cryptic refugium or rapid recolonization?

Repeated recolonization of freshwater environments following Pleistocene glaciations has played a major role in the evolution and adaptation of anadromous taxa. Located at the western fringe of Europe, Ireland and Britain were likely recolonized rapidly by anadromous fishes from the North Atlantic following the last glacial maximum (LGM). While the presence of unique mitochondrial haplotypes in Ireland suggests that a cryptic northern refugium may have played a role in recolonization, no explicit test of this hypothesis has been conducted. The three‐spined stickleback is native and ubiquitous to aquatic ecosystems throughout Ireland, making it an excellent model species with which to examine the biogeographical history of anadromous fishes in the region. We used mitochondrial and microsatellite markers to examine the presence of divergent evolutionary lineages and to assess broad‐scale patterns of geographical clustering among postglacially isolated populations. Our results confirm that Ireland is a region of secondary contact for divergent mitochondrial lineages and that endemic haplotypes occur in populations in Central and Southern Ireland. To test whether a putative Irish lineage arose from a cryptic Irish refugium, we used approximate Bayesian computation (ABC). However, we found no support for this hypothesis. Instead, the Irish lineage likely diverged from the European lineage as a result of postglacial isolation of freshwater populations by rising sea levels. These findings emphasize the need to rigorously test biogeographical hypothesis and contribute further evidence that postglacial processes may have shaped genetic diversity in temperate fauna.     

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.     

Genetic diversity and phylogenetic origin of brown trout Salmo trutta populations in eastern Balkans

The study focuses on the phylogenetic origin and genetic diversity of brown trout in the eastern part of the Balkan Peninsula. It further aims to reveal the impact of human-mediated transfers and stocking with non-indigenous trout on the populations in this area. For these purposes, mtDNA control region and microsatellite variation of 204 individuals from 16 populations were analysed. The results indicate that mtDNA haplotypes from the lower Danube basin and southern Black Sea basins differ substantially from a subclade of the Danubian lineage consisting of haplotypes found so far in the most of the Danube basin and in the Caspian and Aral Sea basins. Considering also the results of demographic analyses, this study evidences a complex evolutionary history of brown trout in the southern and western parts of the Black Sea basin. In the Aegean Sea basin, a high frequency of the central haplotype of Adriatic mtDNA lineage has been found. The other Adriatic lineage haplotypes found in this basin differ from the central haplotype by one mutational step only, indicating a recent evolution of the Adriatic lineage in the Aegean Sea basin. Substantial genetic differentiation among populations and basins was revealed. The hybridization with Atlantic brown trout was indicated in both sea basins, but especially in the Danube basin. Compared to other European regions, it can be inferred that the introgression of exogenous brown trout in the eastern Balkan populations is rather low.     

Late glacial history of the cold-adapted freshwater fish Cottus gobio,revealed by microsatellites

The distribution of genetic diversity at 10 highly polymorphic microsatellite loci within the European freshwater fish, Cottus gobio, L. was examined. The sampling range comprised a large geographical scale including lineages known to be highly divergent at both mitochondrial DNA (mtDNA) and allozymes. An analysis of genetic variability within populations showed that expected heterozygosity and allelic richness could be explained largely by current effective population sizes. Evidence was found, however, that historical processes predating the last major glaciation affected allelic richness. In addition to confirming the large-scale patterns from earlier studies, the microsatellite data revealed new insights into recent processes by analysing genetic structure within ancient lineages defined by mtDNA data. Stepwise mutation model (SMM) and nonSMM-based methods demonstrated a clear genetic structuring within the Northwestern European lineage comprising populations from Britain and Belgium, and within the Central European lineage populations from the rivers Danube, Elbe and Main. Supported by an analysis of genetic variability within populations these results showed that the bullhead populations most probably persisted throughout the last major glaciation within the British Isles and within the drainages of the rivers Elbe and Main. Such observations provide the first genetic evidence for a glacial refugium in such close proximity to the European glacial margins.     

Glacial refuges for three‐spined stickleback in the Iberian Peninsula: mitochondrial DNA phylogeography

相似文献   

Mitochondrial DNA D‐loop sequences suggest a Southeast Asian and Indian origin of Zimbabwean village chickens

This study sought to assess mitochondrial DNA (mtDNA) diversity and phylogeographic structure of chickens from five agro‐ecological zones of Zimbabwe. Furthermore, chickens from Zimbabwe were compared with populations from other geographical regions (Malawi, Sudan and Germany) and other management systems (broiler and layer purebred lines). Finally, haplotypes of these animals were aligned to chicken sequences, taken from GenBank, that reflected populations of presumed centres of domestication. A 455‐bp fragment of the mtDNA D‐loop region was sequenced in 283 chickens of 14 populations. Thirty‐two variable sites that defined 34 haplotypes were observed. In Zimbabwean chickens, diversity within ecotypes accounted for 96.8% of the variation, indicating little differentiation between ecotypes. The 34 haplotypes clustered into three clades that corresponded to (i) Zimbabwean and Malawian chickens, (ii) broiler and layer purebred lines and Northwest European chickens, and (iii) a mixture of chickens from Zimbabwe, Sudan, Northwest Europe and the purebred lines. Diversity among clades explained more than 80% of the total variation. Results indicated the existence of two distinct maternal lineages evenly distributed among the five Zimbabwean chicken ecotypes. For one of these lineages, chickens from Zimbabwe and Malawi shared major haplotypes with chicken populations that have a Southeast Asian background. The second maternal lineage, probably from the Indian subcontinent, was common to the five Zimbabwean chicken ecotypes, Sudanese and Northwest European chickens as well as purebred broiler and layer chicken lines. A third maternal lineage excluded Zimbabwean and other African chickens and clustered with haplotypes presumably originating from South China.     

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.     

Evolutionary history of lamprey paired species Lampetra fluviatilis (L.) and Lampetra planeri (Bloch) as inferred from mitochondrial DNA variation

A remarkable trend in the evolution of lampreys is the occurrence in most genera of 'paired species', in which the parasitic anadromous lampreys are believed to have given rise to nonparasitic freshwater resident populations. The present work examines the phylogeography of the European paired species Lampetra fluviatilis and Lampetra planeri, in an attempt to elucidate species pair evolutionary history. We studied sequence variation in cytochrome b and ATPase 6, 8 mitochondrial genes in 63 individuals from 21 localities of the paired species throughout their distribution range. Results from the phylogenetic and nested clade analyses were largely consistent, suggesting the existence of three major evolutionary lineages: lineage I and possibly lineage II are widespread throughout Europe, while the most ancestral lineage III is apparently restricted to the Iberian Peninsula. The high genetic diversity observed in the Iberian Peninsula is probably the result of refugial persistence and subsequent accumulation of variation over several ice ages, whereas the low levels of genetic diversity observed in central and northern Europe should reflect a rapid postglacial colonization. Results suggest that L. planeri originated within at least two distinct evolutionary lineages, rejecting the single origin hypothesis. The observed lack of taxa monophyly within lineage I may be the result of ongoing gene flow if the two taxa are alternate life-history forms of a single species. However, structure within lineage I is also consistent with the hypothesis of divergence of taxa after postglacial dispersion (around 2000 generations ago) with incomplete lineage sorting. Further testing of the alternative hypotheses is warranted.     

Rapid genetic assimilation of native wall lizard populations (Podarcis muralis) through extensive hybridization with introduced lineages

The Common Wall Lizard (Podarcis muralis) has established more than 150 non-native populations in Central Europe, stemming from eight geographically distinct evolutionary lineages. While the majority of these introduced populations are found outside the native range, some of these populations also exist at the northern range margin in southwestern Germany. To (i) infer the level of hybridization in contact zones of alien and native lineages; and (ii) compare the genetic diversity among purebred introduced, native and hybrid populations, we used a combination of maternally inherited markers (mtDNA: cytb) and Mendelian markers (microsatellites). Our results suggest a rapid genetic assimilation of native populations by strong introgression from introduced lineages. Discordant patterns of mtDNA and nDNA variation within hybrid populations may be explained by directed mate choice of females towards males of alien lineages. In contrast to previous studies, we found a nonlinear relationship between genetic diversity and admixture level. The genetic diversity of hybrid populations was substantially higher than in introduced and native populations belonging to a single lineage, but rapidly reaching a plateau of high genetic diversity at an admixture level of two. However, even introduced populations with low founder sizes and from one source population retained moderate levels of genetic diversity and no evidence for a genetic bottleneck was found. The extent of introgression and the dominance of alien haplotypes in mixed populations indicate that introductions of non-native lineages represent a serious threat to the genetic integrity of native populations due to the rapid creation of hybrid swarms.     

African DNA Lineages in the Mitochondrial Gene Pool of Europeans

Mitochondrial DNA (mtDNA) nucleotide sequences of African origin are found in various European populations at a low frequency (on average, less than 1%). Data on mtDNA variation in Eurasian and African populations have been analyzed, and African mtDNA lineages have been found in Europeans. It has been demonstrated that, despite the high diversity of mtDNA haplotypes of African origin in Europeans, few monophyletic clusters of African lineages are characterized by long-term diversity formed in Europe. Only two such mtDNA clusters (from haplogroups L1b and L3b) have been found, their evolutionary age not exceeding 6500 years. European and African populations have been compared with respect to the frequency distributions of the alleles of autosomal microsatellite loci found in Russian carriers of African mtDNA haplotypes. It has been demonstrated that alleles typical of Europeans are characteristic of the autosomal genotypes of these Russian individuals.     

Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L) as revealed by mitochondrial DNA

Bearded vulture populations in the Western Palearctic have experienced a severe decline during the last two centuries that has led to the near extinction of the species in Europe. In this study we analyse the sequence variation at the mitochondrial control region throughout the species range to infer its recent evolutionary history and to evaluate the current genetic status of the species. This study became possible through the extensive use of museum specimens to study populations now extinct. Phylogenetic analysis revealed the existence of two divergent mitochondrial lineages, lineage A occurring mainly in Western European populations and lineage B in African, Eastern European and Central Asian populations. The relative frequencies of haplotypes belonging to each lineage in the different populations show a steep East-West clinal distribution with maximal mixture of the two lineages in the Alps and Greece populations. A genealogical signature for population growth was found for lineage B, but not for lineage A; futhermore the Clade B haplotypes in western populations and clade A haplo-types in eastern populations are recently derived, as revealed by their peripheral location in median-joining haplotype networks. This phylogeographical pattern suggests allopatric differentiation of the two lineages in separate Mediterranean and African or Asian glacial refugia, followed by range expansion from the latter leading to two secondary contact suture zones in Central Europe and North Africa. High levels of among-population differentiation were observed, although these were not correlated with geographical distance. Due to the marked genetic structure, extinction of Central European populations in the last century re-sulted in the loss of a major portion of the genetic diversity of the species. We also found direct evidence for the effect of drift altering the genetic composition of the remnant Pyrenean population after the demographic bottleneck of the last century. Our results argue for the management of the species as a single population, given the apparent ecological exchangeability of extant stocks, and support the ongoing reintroduction of mixed ancestry birds in the Alps and planned reintroductions in Southern Spain.     

African DNA lineages in mitochondrial gene pool of Europeans

Mitochondrial DNA (mtDNA) nucleotide sequences of African origin have been found at low frequency (1%, in average) in different European populations. In the present study, data on mtDNA variability in populations of Eurasia and Africa are analyzed and search of African-specific lineages present in Europeans is conducted. The results of analysis indicate that, despite a high diversity of African mtDNA haplotypes found in Europeans, monophyletic clusters of African mtDNA lineages, arisen in Europe and characterized by long-term diversity, are nearly absent in Europe. Only two respective clusters (belonging to haplogroups L1b and L3b), which evolutionary age does not exceed 6.5 thousands years, were revealed. Comparative analysis of distribution of frequencies of autosomal microsatellite alleles found in Russian individuals, carrying the African-specific mitochondrial haplotypes, in populations of Europe and Africa has indicated that autosomal genotypes of those Russian individuals are characterized by the presence of alleles characteristic mostly for Europeans.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

Origin, genetic diversity, and population structure of Chinese domestic sheep

To characterize the origin, genetic diversity, and phylogeographic structure of Chinese domestic sheep, we here analyzed a 531-bp fragment of mtDNA control region of 449 Chinese autochthonous sheep from 19 breeds/populations from 13 geographic regions, together with previously reported 44 sequences from Chinese indigenous sheep. Phylogenetic analysis showed that all three previously defined lineages A, B, and C were found in all sampled Chinese sheep populations, except for the absence of lineage C in four populations. Network profiles revealed that the lineages B and C displayed a star-like phylogeny with the founder haplotype in the centre, and that two star-like subclades with two founder haplotypes were identified in lineage A. The pattern of genetic variation in lineage A, together with the divergence time between the two central founder haplotypes suggested that two independent domestication events have occurred in sheep lineage A. Considerable mitochondrial diversity was observed in Chinese sheep. Weak structuring was observed either among Chinese indigenous sheep populations or between Asian and European sheep and this can be attributable to long-term strong gene flow induced by historical human movements. The high levels of intra-population diversity in Chinese sheep and the weak phylogeographic structuring indicated three geographically independent domestication events have occurred and the domestication place was not only confined to the Near East, but also occurred in other regions.     

Pan-European phylogeography of the aquatic snail Theodoxus fluviatilis (Gastropoda: Neritidae)

Investigating the geographical distribution of genetic lineages within species is critical to our understanding of how species evolve. As many species inhabit large and complex ranges, it is important that phylogeographical research take into account the entire range of widespread species to clarify how myriad extrinsic variables have affected their evolutionary history. Using phylogenetic, nested clade, and mismatch distribution analyses on a portion of the mitochondrial COI gene, I demonstrate that the wide-ranging freshwater snail Theodoxus fluviatilis possesses in parallel many of the phylogeographical patterns seen in less widespread freshwater species of Europe. Fragmentary forces play a major part in structuring the range of this species, with 12 of 14 geographically structured nested clades displaying a distribution consistent with fragmentation or restricted dispersal. Certain regions of southern Europe harbour the majority of genetic diversity (total haplotype diversity, H = 0.87), particularly Italy (H = 0.87) and areas surrounding the Black Sea (H = 0.81). Post-Pleistocene range expansion is pronounced, with the majority of northern European populations (95% of sample sites) having arisen from northern Italian individuals that initially colonized northern Germany. Additionally, two highly divergent haplotype lineages present in northern Germany imply that there were at least two postglacial recolonization routes. Estuaries may also provide a means of dispersal given that no genetic differentiation was found between estuarine populations and neighbouring freshwater populations. Taken together, these data reveal a species with a complex genetic history resulting from the fragmentary effects of European geology as well as continuous and discrete range expansion related to their aquatic biology.     

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.     

DIFFERING LEVELS OF AMONG-POPULATION DIVERGENCE IN THE MITOCHONDRIAL DNA OF PERIODICAL CICADAS RELATED TO HISTORICAL BIOGEOGRAPHY

Mitochondrial DNA (mtDNA) haplotypes were determined for 118 individuals of 13-and 17-year periodical cicadas (genus Magicicada) collected from 16 localities throughout the Midwest and eastern United States. Two distinct mtDNA lineages, identified as A and B, differ by 2.5% based on analysis of fragment patterns and restriction maps. Observed levels of mtDNA diversity within each lineage are low compared to estimates for other taxa. The two lineages are regionally segregated, with the boundary line occurring at a latitude of approximately 33° North. The levels of mtDNA diversity and population genetic structure differ within the two lineages. There is a remarkably low level of mean mtDNA divergence and no genetic structure in lineage A, whereas lineage B exhibits an order of magnitude higher level of mtDNA diversity and significant genetic structure among sampled populations. The low level of mtDNA diversity in cicadas may be attributed to (1) a population bottleneck that most likely occurred during the Pleistocene, (2) recent colonization following the retreat of the glaciers and the expansion of deciduous forests, and/or (3) high among-family reproductive variance (as a consequence of large population size, high fecundity, aggregative behavior of adults, and clumping of eggs). The difference in mtDNA diversity and population genetic structure between the lineages suggests that they experienced different biogeographic histories; we relate this to Pleistocene changes.     

Admixture and Genetic Diversity Distribution Patterns of Non-Recombining Lineages of Native American Ancestry in Colombian Populations

Genetic diversity of present American populations results from very complex demographic events involving different types and degrees of admixture. Through the analysis of lineage markers such as mtDNA and Y chromosome it is possible to recover the original Native American haplotypes, which remained identical since the admixture events due to the absence of recombination. However, the decrease in the effective population sizes and the consequent genetic drift effects suffered by these populations during the European colonization resulted in the loss or under-representation of a substantial fraction of the Native American lineages. In this study, we aim to clarify how the diversity and distribution of uniparental lineages vary with the different demographic characteristics (size, degree of isolation) and the different levels of admixture of extant Native groups in Colombia. We present new data resulting from the analyses of mtDNA whole control region, Y chromosome SNP haplogroups and STR haplotypes, and autosomal ancestry informative insertion-deletion polymorphisms in Colombian individuals from different ethnic and linguistic groups. The results demonstrate that populations presenting a high proportion of non-Native American ancestry have preserved nevertheless a substantial diversity of Native American lineages, for both mtDNA and Y chromosome. We suggest that, by maintaining the effective population sizes high, admixture allowed for a decrease in the effects of genetic drift due to Native population size reduction and thus resulting in an effective preservation of the Native American non-recombining lineages.     

Mitochondrial diversity of European pond turtles (Emys orbicularis) in Anatolia and the Ponto-Caspian Region: Multiple old refuges, hotspot of extant diversification and critically endangered endemics

The European pond turtle, Emys orbicularis (Linnaeus, 1758), is one of the world's most widely distributed chelonian species. We investigated its mitochondrial phylogeography and demography using ∼1300 cyt b sequences from the entire range, with a focus on the eastern part, in particular on Turkey, where the species currently suffers massive losses. Coloration data from >1450 turtles were compared with mtDNA differentiation to assess the validity of the currently accepted subspecies from Turkey, the Black Sea Region, the Transcaucasus, and Iran. Our study region harbors considerable part of the mtDNA diversity of Emys, including a newly discovered lineage and 16 new haplotypes. In this area corresponding to approximately one-third of the entire distribution range, six out of the ten mitochondrial lineages and about half of all 72 haplotypes occur. Two mitochondrial lineages (VIII, X) are confined to small ranges along the southern coast of Turkey, another lineage (I) occupies the remainder of Turkey, the entire Black Sea Region, and the north-eastern part of the species’ range. In the south-western corner of the Black Sea and in the Aegean Region, two lineages (II, IV) occur that have their main distribution areas farther west. In the Transcaucasus and northern Iran, another endemic lineage (VII) is found. Lineage I is the largest and most diverse of all lineages and has its greatest diversity in Anatolia. Phylogeographic and demographic data suggest Anatolia as an ancient glacial refuge for turtles harboring mitochondrial lineages I, VIII and X, and that Anatolia and the Black Sea coasts constitute a hotspot for a younger burst of diversification within lineage I. These two regions correspond to the glacial refuge from which lineage I turtles recolonized more northerly parts of the range in the Holocene; lineage II represents an off-shoot of lineage I that became isolated in a westward-located refuge in the south-eastern Balkans during a previous Pleistocene glacial. Our data on coloration indicate that such characters have only limited value for delineating evolutionarily significant units. We propose to reduce the number of subspecies using mtDNA lineages as arbiter, and to recognize three subspecies as valid in Turkey, the Black Sea Region, the Transcaucasus and Iran: Emys orbicularis orbicularis (mtDNA lineage I); E. o. eiselti Fritz et al., 1998 (X); and E. o. persica Eichwald, 1831 (VII). However, the southern Turkish lineage VIII most probably represents an additional undescribed subspecies. Both southern Turkish endemics are critically endangered, with only three surviving populations of fewer than 30 adults each. We recommend establishing sanctuaries for them, and including them in the IUCN Red List.