Uwe Fritz

Mitochondrial phylogeography of the European pond turtle, Emys orbicularis (Linnaeus 1758)

The phylogeny and phylogeography of Emys orbicularis was inferred from mitochondrial nucleotide sequences of the cytochrome b gene analysed by DNA sequencing and RNA heteroduplex analysis. Within the family Emydidae the monotypic genus Emys is affiliated with the nearctic taxa Emydoidea blandingii and Clemmys marmorata. The analysis of 423 individuals of E. orbicularis, originating throughout its distribution range, revealed a remarkable intraspecific differentiation in 20 different haplotypes with distinct geographical ranges. Maximum parsimony analysis produced a star‐like phylogeny with seven main lineages which may reflect separations in the late Pliocene. The haplotype distribution examined by partial Mantel tests and analysis of molecular variance revealed a substantial effect of glacial periods. This historical perspective suggests the existence of multiple glacial refugia and considerable Holocene range expansion which was modulated by climatic traits. Further support is gained for the occurrence of long‐term parapatry in glacial refugia.

Impact of mountain chains, sea straits and peripheral populations on genetic and taxonomic structure of a freshwater turtle, Mauremys leprosa (Reptilia, Testudines, Geoemydidae)

Mauremys leprosa, distributed in Iberia and North‐west Africa, contains two major clades of mtDNA haplotypes. Clade A occurs in Portugal, Spain and Morocco north of the Atlas Mountains. Clade B occurs south of the Atlas Mountains in Morocco and north of the Atlas Mountains in eastern Algeria and Tunisia. However, we recorded a single individual containing a clade B haplotype in Morocco from north of the Atlas Mountains. This could indicate gene flow between both clades. The phylogenetically most distinct clade A haplotypes are confined to Morocco, suggesting both clades originated in North Africa. Extensive diversity within clade A in south‐western Iberia argues for a glacial refuge located there. Other regions of the Iberian Peninsula, displaying distinctly lower haplotype diversities, were recolonized from within south‐western Iberia. Most populations in Portugal, Spain and northern Morocco contain the most common clade A haplotype, indicating dispersal from the south‐western Iberian refuge, gene flow across the Strait of Gibraltar, and reinvasion of Morocco by terrapins originating in south‐western Iberia. This hypothesis is consistent with demographic analyses, suggesting rapid clade A population increase while clade B is represented by stationary, fragmented populations. We recommend the eight, morphologically weakly diagnosable, subspecies of M. leprosa be reduced to two, reflecting major mtDNA clades: Mauremys l. leprosa (Iberian Peninsula and northern Morocco) and M. l. saharica (southern Morocco, eastern Algeria and Tunisia). Peripheral populations could play an important role in evolution of M. leprosa because we found endemic haplotypes in populations along the northern and southern range borders. Previous investigations in another western Palearctic freshwater turtle (Emys orbicularis) discovered similar differentiation of peripheral populations, and phylogeographies of Emys orbicularis and Mauremys rivulata underline the barrier status of mountain chains, in contrast to sea straits, suggesting common patterns for western Palearctic freshwater turtles.

A new cryptic species of pond turtle from southern Italy, the hottest spot in the range of the genus Emys (Reptilia, Testudines, Emydidae)

Geographic variation in the mtDNA haplotypes (cytochrome b gene) of 127 European pond turtles from Italy was investigated. Thirty‐eight of the Italian samples were also studied by nuclear fingerprinting (ISSR PCR) and compared with samples from other parts of the range representing all nine currently known mtDNA lineages of Emys orbicularis. Our genetic findings were compared against morphological data sets (measurements, colour pattern) for 109 adult turtles from southern Italy. Italy is displaying on a small geographical scale the most complicated variation known over the entire distributional area of Emys (North Africa over Europe and Asia Minor to the Caspian and Aral Seas). The Tyrrhenic coast of the Apennine Peninsula, the Mt. Pollino area and Basilicata are inhabited by Emys orbicularis galloitalica, a subspecies harbouring a distinct mtDNA lineage. The same lineage is also found in Sardinia. Along the Adriatic coast of Italy and on the Salentine Peninsula (Apulia, southern Italy), another morphologically distinctive subspecies (Emys orbicularis hellenica) occurs, which also bears a different mtDNA lineage. A higher diversity of mtDNA haplotypes in the south of the Apennine Peninsula suggests that the glacial refugia of E. o. galloitalica and E. o. hellenica were located here. A further refuge of E. o. hellenica probably existed in the southern Balkans. The west coasts of the Balkans and Corfu have probably been colonized from Italy and not from the geographically closer southern Balkanic refuge. In Sicily, a third mtDNA lineage is distributed, which is sister to all other known lineages of Emys. Morphologically, Sicilian pond turtles resemble E. o. galloitalica. However, nuclear fingerprinting revealed a clear distinctiveness of the Sicilian taxon, whereas no significant divergence was detected between representatives of the other eight mtDNA lineages of Emys. Furthermore, nuclear fingerprinting provided no evidence for current or past gene flow between the Sicilian taxon and the mainland subspecies of E. orbicularis. Therefore, Sicilian pond turtles are described here as a species new to science. Some populations in Calabria and on the Salentine Peninsula comprise individuals of different mtDNA lineages. We interpret this as a natural contact. However, we cannot exclude that these syntopic occurrences are the result of human activity. For example, in other parts of Italy, the natural distribution pattern of Emys is obscured by allochthonous turtles. This could also be true for southern Italy. The discovery of the complex taxonomic differentiation in southern Italy requires reconsidering conservation strategies.

Phenotypic plasticity leads to incongruence between morphology-based taxonomy and genetic differentiation in western Palaearctic tortoises (Testudo graeca complex; Testudines, Testudinidae)

Tortoises of the Testudo graeca complex inhabit a patchy range that covers part of three continents (Africa, Europe, Asia). It extends approximately 6500 km in an east-west direction from eastern Iran to the Moroccan Atlantic coast and about 1600 km in a north-south direction from the Danube Delta to the Libyan Cyrenaica Peninsula. Recent years have seen a rapid increase of recognized taxa. Based on morphological investigations, it was suggested that this group consists of as many as 20 distinct species and is paraphyletic with respect to T. kleinmanni sensu lato and T. marginata . Based on samples from representative localities of the entire range, we sequenced the mitochondrial cytochrome b gene and conducted nuclear genomic fingerprinting with ISSR PCR. The T. graeca complex is monophyletic and sister to a taxon consisting of T. kleinmanni sensu lato and T. marginata . The T. graeca complex comprises six well-supported mtDNA clades (A-F). Highest diversity is found in the Caucasian Region, where four clades occur in close neighbourhood. This suggests, in agreement with the fossil record, the Caucasian Region as a radiation centre. Clade A corresponds to haplotypes from the East Caucasus. It is the sister group of another clade (B) from North Africa and western Mediterranean islands. Clade C includes haplotypes from western Asia Minor, the southeastern Balkans and the western and central Caucasus Region. Its sister group is a fourth, widely distributed clade (D) from southern and eastern Asia Minor and the Levantine Region (Near East). Two further clades are distributed in Iran (E, northwestern and central Iran; F, eastern Iran). Distinctness of these six clades and sister group relationships of (A + B) and (C + D) are well-supported; however, the phylogeny of the resulting four clades (A + B), (C + D), E and F is poorly resolved. While in a previous study (Fritz et al., 2005a) all traditionally recognized Testudo species were highly distinct using mtDNA sequences and ISSR fingerprints, we detected within the T. graeca complex no nuclear genomic differentiation paralleling mtDNA clades. We conclude that all studied populations of the T. graeca complex are conspecific under the Biological Species Concept. There is major incongruence between mtDNA clades and morphologically defined taxa. Morphologically well-defined taxa, like T. g. armeniaca or T. g. floweri , nest within clades comprising also geographically neighbouring, but morphologically distinctive populations of other taxa (clade A: T. g. armeniaca , T. g. ibera , T. g. pallasi ; clade D: T. g. anamurensis , T. g. antakyensis , T. g. floweri , T. g. ibera , T. g. terrestris ), while sequences of morphologically similar tortoises of the same subspecies ( T. g. ibera sensu stricto or T. g. ibera sensu lato) scatter over two or three genetically distinct clades (A, C or A, C, D, respectively). This implies that pronounced morphological plasticity, resulting in phenotypes shaped by environmental pressure, masks genetic differentiation. To achieve a more realistic taxonomic arrangement reflecting mtDNA clades, we propose reducing the number of T. graeca subspecies considerably and regard in the eastern part of the range five subspecies as valid ( T. g. armeniaca , T. g. buxtoni , T. g. ibera , T. g. terrestris , T. g. zarudnyi ). As not all North African taxa were included in the present study, we refrain from synonymizing North African taxa with T. g. graeca (mtDNA clade B) that represents a further valid subspecies.

Mitochondrial phylogeography of European pond turtles ( Emys orbicularis , Emys trinacris ) – an update

Based on more than 1100 samples of Emys orbicularis and E. trinacris, data on mtDNA diversity and distribution of haplotypes are provided, including for the first time data for Armenia, Georgia, Iran, and the Volga, Ural and Turgay River Basins of Russia and Kazakhstan. Eight mitochondrial lineages comprising 51 individual haplotypes occur in E. orbicularis, a ninth lineage with five haplotypes corresponds to E. trinacris. A high diversity of distinct mtDNA lineages and haplotypes occurs in the south, in the regions where putative glacial refuges were located. More northerly parts of Europe and adjacent Asia, which were recolonized by E. orbicularis in the Holocene, display distinctly less variation; most refuges did not contribute to northern recolonizations. Also in certain southern European lineages a decrease of haplotype diversity is observed with increasing latitude, suggestive of Holocene range expansions on a smaller scale.

A rangewide phylogeography of Hermann's tortoise, Testudo hermanni (Reptilia: Testudines: Testudinidae): implications for taxonomy

Hermanns tortoise (Testudo hermanni), the best‐known western Palaearctic tortoise species, has a rare natural distribution pattern comprising the Mediterranean areas of the Iberian, Apennine, and Balkan Peninsulas, as well as Sicily, Corsica and Sardinia. The western part of this range is traditionally considered habitat for T. h. hermanni, while T. h. boettgeri occurs in the Balkans. Taxonomy of this tortoise has been challenged in recent years, with the two subspecies being considered full species and the central Dalmatian populations of T. h. boettgeri being considered a third species, T. hercegovinensis. Using an mtDNA fragment approximately 1150 bp long (cytochrome b gene and adjacent portion of tRNA‐Thr gene), we investigated mtDNA diversity with regard to contrasting concepts of two subspecies or three species. Seven closely related haplotypes were identified from the western Mediterranean and 15 different, in part much‐differentiated, haplotypes from the Balkans. Western Mediterranean haplotypes differ from Balkan haplotypes in 16–42 mutation steps. One to seven mutation steps occur within western Mediterranean populations. Balkan haplotypes, differing in 1−37 nucleotides, group in parsimony network analysis into three major assemblages that display, in part, a similar degree of differentiation to that of western Mediterranean haplotypes relative to Balkan haplotypes. Rates of sequence evolution are different in both regions, and low divergence, palaeogeography and the fossil record suggest a slower molecular clock in the western Mediterranean. While monophyly in western Mediterranean haplotypes is well‐supported, conflicting evidence is obtained for Balkan haplotypes; maximum parsimony supports monophyly of Balkan haplotypes, but other phylogenetic analyses (Bayesian, ML, ME) indicate Balkan haplotypes could be paraphyletic with respect to the western Mediterranean clade. These results imply a process of differentiation not yet complete despite allopatry in the western Mediterranean and the Balkans, and suggest all populations of T. hermanni are conspecific. In the western Mediterranean no clear geographical pattern in haplotype distribution is found. Distribution of Balkan haplotypes is more structured. One group of similar haplotypes occurs in the eastern Balkans (Bulgaria, Republic of Macedonia, Romania and the Greek regions Evvia, Macedonia, Peloponnese, Thessaly and Thrace). Two distinct haplotypes, differing in eight to nine mutation steps from the most common haplotype of the first group, are confined to the western slope of the Taygetos Mts. in the Peloponnese. Yet another group, connected over between four and 23 mutation steps with haplotypes of the eastern Balkan group, occurs along the western slope of the Dinarid and Pindos Mts. (Istria, Dalmatia, western Greece). Taygetos haplotypes are nested within other haplotypes in all phylogenetic analyses and support for monophyly of the other Balkan groups is at best weak. We conclude that using the traditional two subspecies model should be continued for T. hermanni. Phylogeographies of T. hermanni and Emys orbicularis, another codistributed chelonian, are markedly different, but share a few similarities. Both were forced to retreat to southern refuges during Pleistocene glaciations. With the advent of Holocene warming, E. orbicularis underwent rapid range expansion and temperate regions of Europe and adjacent Asia were recolonized from refuges in the Balkans and the northern Black Sea Region. By contrast, T. hermanni remained more or less confined to refuges and nearby regions, resulting in a much smaller range, and allopatric and parapatric distribution of haplotype groups and clades. MtDNA lineages are more diverse in E. orbicularis than they are in T. hermanni on southern European peninsulas, indicating several distinct glacial refuges in close proximity and extensive intergradation during Holocene range expansion for E. orbicularis. In T. hermanni it is likely that only on the Balkan Peninsula was more than one refuge located, corresponding to the parapatric ranges of haplotype groups currently there. On the old western Mediterranean islands Corsica and Sardinia no differentiated (E. orbicularis) or only weakly differentiated haplotypes (T. hermanni) occur, even though there is evidence for the presence of both species on Corsica since at least the Middle Pleistocene. High mountain chains constitute major barriers separating distinct mtDNA clades or groups in each species.

Phylogeography and cryptic variation within the Lacerta viridis complex (Lacertidae, Reptilia)

It is well known that the current genetic pattern of many European species has been highly influenced by climatic changes during the Pleistocene. While there are many well known vertebrate examples, knowledge about squamate reptiles is sparse. To obtain more data, a range‐wide sampling of Lacerta viridis was conducted and phylogenetic relations within the L. viridis complex were analysed using an mtDNA fragment encompassing part of cytochrome b, the adjacent tRNA genes and the noncoding control region. Most genetic divergence was found in the south of the distribution range. The Carpathian Basin and the regions north of the Carpathians and Alps are inhabited by the same mitochondrial lineage, corresponding to Lacerta viridis viridis. Three distinct lineages occurred in the south‐eastern Balkans — corresponding to L. v. viridis, L. v. meridionalis, L. v. guentherpetersi— as well as a fourth lineage for which no subspecies name is available. This distribution pattern suggests a rapid range expansion of L. v. viridis after the Holocene warming, leading to a colonization of the northern part of the species range. An unexpected finding was that a highly distinct genetic lineage occurs along the western Balkan coast. Phylogenetic analyses (Bayesian, maximum likelihood, maximum parsimony) suggested that this west Balkan lineage could represent the sister taxon of Lacerta bilineata. Due to the morphological similarity of taxa within the L. viridis complex this cryptic taxon was previously assigned to L. v. viridis. The distribution pattern of several parapatric, in part highly, distinct genetic lineages suggested the existence of several refuges in close proximity on the southern Balkans. Within L. bilineata sensu stricto a generally similar pattern emerged, with a high genetic diversity on the Apennine peninsula, arguing for two distinct refuges there, and a low genetic diversity in the northern part of the range. Close to the south‐eastern Alps, three distinct lineages (L. b. bilineata, L. v. viridis, west Balkan taxon) occurred within close proximity. We suggest that the west Balkan lineage represents an early offshoot of L. bilineata that was isolated during a previous Pleistocene glacial from the more western L. bilineata populations, which survived in refuges on the Apennine peninsula.

Mitochondrial phylogeography of Testudo graeca in the Western Mediterranean: Old complex divergence in North Africa and recent arrival in Europe

We investigated the mitochondrial phylogeography of spur-thighed tortoises (Testudo graeca) in the Western Mediterranean. In North Africa, four major lineages (A-D) occur that together constitute a well-supported clade corresponding to one of the six major clades within T. graeca; the North African clade is sister to a Caucasian clade representing the subspecies T. g. armeniaca. Phylogenetic relationships between the North African lineages are badly resolved. Lineage A is distributed in Tunisia and adjacent Algeria, lineage B in Algeria and northern Morocco, lineage C in the Libyan Cyrenaica Peninsula, and lineage D north of the High Atlas Mts. and in the Souss Valley (southern Morocco). Lineage B is subdivided into two subgroups, B1 (eastern Morocco and Algeria) and B2 (north-western Morocco). Italian tortoises harbour haplotypes of lineage A, Spanish tortoises of subgroup B1. Based on a relaxed molecular clock calibrated with fossil evidence, the six major mtDNA clades of T. graeca are estimated to have diverged approximately 4.2-1.8 Ma ago; the split between the clades representing the eastern subspecies T. g. ibera and T. g. terrestris is younger than the split between Western Mediterranean tortoises and T. g. armeniaca. The Western Mediterranean lineages A-D were dated to have diverged at least 1.4-1.1 Ma ago; B1 and B2 split approximately 0.7 Ma ago. Our results suggest that Italian and Spanish tortoises were either introduced or originated from trans-oceanic dispersal in historic or prehistoric times. Spur-thighed tortoises invaded North Africa probably across Near Eastern landbridges that emerged in the Late Tertiary. Their diversification in North Africa seems to be correlated with habitat aridization cycles during the Pleistocene. The ranges of the Western Mediterranean lineages largely correspond to the distribution of morphologically defined subspecies in North Africa, with exception of T. g. graeca and T. g. whitei, and of T. g. lamberti and T. g. marokkensis, which are not differentiated. We propose to lump the first two subspecies under the name of T. g. graeca and the latter under the name of T. g. marokkensis. The complex differentiation of spur-thighed tortoises in North Africa implies that the model of a bipartite east-west differentiation, as proposed for other Maghrebian amphibians and reptiles, may be too simplistic, reflecting incomplete locality sampling rather than actual phylogeographic differentiation.

The freshwater turtle genus Mauremys (Testudines, Geoemydidae) — a textbook example of an east–west disjunction or a taxonomic misconcept?

Barth, D. Bernhard, D. Fritzsch, G. & Fritz, U. (2004): The freshwater turtle genus Mauremys— a textbook example of an east–west disjunction or a taxonomic misconcept? —Zoologica Scripta, 33, 213–221.

Go east: phylogeographies of Mauremys caspica and M. rivulata– discordance of morphology, mitochondrial and nuclear genomic markers and rare hybridization

In recent years many cases of hybridization and introgression became known for chelonians, requiring a better understanding of their speciation mechanisms. Phylogeographic investigations offer basic data for this challenge. We use the sister species Mauremys caspica and M. rivulata, the most abundant terrapins in the Near and Middle East and South‐east Europe, as model. Their phylogeographies provide evidence that speciation of chelonians fits the allopatric speciation model, with both species being in the parapatric phase of speciation, and that intrinsic isolation mechanisms are developed during speciation. Hybridization between M. caspica and M. rivulata is very rare, suggesting that the increasing numbers of hybrids in other species are caused by human impact on environment (breakdown of ecological isolation). Genetic differentiation within M. caspica and M. rivulata resembles the paradigm of southern genetic richness and northern purity of European biota. However, in west Asia this pattern is likely to reflect dispersal and vicariance events older than the Holocene. For M. caspica three distinct Pleistocene refuges are postulated (Central Anatolia, south coast of Caspian Sea, Gulf of Persia). Morphologically defined subspecies within M. caspica are not supported by genetic data. This is one of the few studies available about the phylogeography of west and central Asian species.