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.     

Where are you from,stranger? The enigmatic biogeography of North African pond turtles (Emys orbicularis)

The European pond turtle (Emys orbicularis) is a Nearctic element in the African fauna and thought to have invaded North Africa from the Iberian Peninsula. All North African populations are currently identified with the subspecies E. o. occidentalis. However, a nearly range-wide sampling in North Africa used for analyses of mitochondrial and microsatellite DNA provides evidence that only Moroccan populations belong to this taxon, while eastern Algerian and Tunisian pond turtles represent an undescribed distinct subspecies. These two taxa are most closely related to E. o. galloitalica with a native distribution along the Mediterranean coast of northern Spain through southern France to western and southern Italy. This group is sister to a clade comprising several mitochondrial lineages and subspecies of E. orbicularis from Central and Eastern Europe plus Asia, and the successive sisters are E. o. hellenica and E. trinacris. Our results suggest that E. orbicularis has been present in North Africa longer than on the Iberian Peninsula and that after an initial invasion of North Africa by pond turtles from an unknown European source region, there was a phase of diversification in North Africa, followed by a later re-invasion of Europe by one of the African lineages. The differentiation of pond turtles in North Africa parallels a general phylogeographic paradigm in amphibians and reptiles, with deeply divergent lineages in the western and eastern Maghreb. Acknowledging their genetic similarity, we propose to synonymize the previously recognized Iberian subspecies E. o. fritzjuergenobsti with E. o. occidentalis sensu stricto. The seriously imperiled Moroccan populations of E. o. occidentalis represent two Management Units different in mitochondrial haplotypes and microsatellite markers. The conservation status of eastern Algerian pond turtles is unclear, while Tunisian populations are endangered. Considering that Algerian and Tunisian pond turtles represent an endemic taxon, their situation throughout the historical range should be surveyed to establish a basis for conservation measures.     

Gene flow across secondary contact zones of the Emys orbicularis complex in the Western Mediterranean and evidence for extinction and re‐introduction of pond turtles on Corsica and Sardinia (Testudines: Emydidae)

European pond turtles represent a phylogeographically deeply structured complex of distinct taxa. Here, we use mitochondrial DNA sequences (cytochrome b gene) and eight polymorphic microsatellite loci to investigate genetic differentiation and gene flow of Sicilian, Corsican and Sardinian pond turtles and of subspecies involved in two secondary contact zones in the Pyrenean region and Southern Italy. Mitochondrial and microsatellite differentiation is largely concordant in populations from the core regions of the distribution ranges of the studied taxa. Both marker systems provide no evidence for gene flow between Sicilian pond turtles (Emys trinacris) and Southern Italian subspecies of E. orbicularis. By contrast, in the contact zones limited gene flow occurs between distinct subspecies of E. orbicularis. Although the Southern Italian contact zone is significantly older than the Pyrenean contact zone of Holocene age, patterns of asymmetric introgression are similar. Introgressive hybridization leads to the exchange of mitochondria, but microsatellite data indicate only a few individuals with mixed ancestry. This suggests that incipient isolating mechanisms maintain largely discrete nuclear genomic gene pools. Furthermore, this implies that Southern Italy acted as a hotspot rather than as a melting pot of genetic diversity during the last glacial. Pond turtles from Corsica and Sardinia are not differentiated from continental populations of the subspecies E. o. galloitalica, neither in the mitochondrial nor in the quickly evolving microsatellite markers. As the fossil record argues for a continuous presence of pond turtles on both islands since the Middle Pleistocene, this suggests that the native island populations became extinct and the extant turtles were later introduced by prehistoric settlers. The lack of genetic differentiation of pond turtles from Corsica and Sardinia supports the view that the subspecies described from these islands are not valid.     

Some like it cold: the glossiphoniid parasites of the Sicilian endemic pond turtle Emys trinacris (Testudines,Emydidae), an example of ‘parasite inertia’?

The freshwater turtles of the genus Emys and some leech species of the family Glossiphoniidae are the only Palaearctic representatives of primarily Nearctic taxa, which jointly colonized Eurasia and the Maghreb during the Miocene. The strict trophic relationships occurring between the glossiphoniid parasite leech Placobdella costata and its host, the emydid Emys orbicularis, make them a prime example of host–parasite cophylogenetic evolution. In the light of the discovery of the Sicilian cryptic endemic species Emys trinacris, which is the sister species to the widespread Palaearctic E. orbicularis, the possible cophylogenetic divergence of the turtle hosts and their leech parasites was investigated. In spite of the deep divergence scored between the two pond turtle species and of their allopatric distribution, their leech parasites proved to be conspecific and indistinguishable based on the implemented molecular marker. This unexpected decoupling might likely be ascribed to the different dispersal abilities of the two taxa and/or to the recent, human‐mediated introduction of the leech parasites in Sicily. If this last scenario is confirmed, the long‐term effects of the introduced leech parasite on the endemic Sicilian pond turtle Emys trinacris should be carefully monitored. In the frame of this study, representatives of the widely spread predatory leech Helobdella stagnalis were observed on E. trinacris. Molecular analyses of their stomach content allowed to rule out the possibility of the existence of a trophic relationships between these two taxa, in contrast to what was previously suspected, and suggest that H. stagnalis specimens were rather attached to the turtles for non‐nutritional reasons.     

Distribution of mitochondrial haplotypes (cytb) in Polish populations of Emys orbicularis (L., 1758)

The European pond turtle, Emys orbicularis, inhabits a wide distribution area in the western Palaearctic. Polish populations of pond turtle represent the nominotypical subspecies Emys orbicularis orbicularis. The mitochondrial DNA haplotype (cytb gene) variation among 131 turtles from 26 locations in five regions of Poland was investigated. Five haplotypes belonging to three distinct lineages were identified. Two clades (I and II) were represented by two haplotypes each, while the other clade (IV) was represented by one haplotype. Three haplotypes were reported for the first time in E. orbicularis. The eastern part of Poland is inhabited exclusively by turtles bearing haplotype Ia. The remaining four sequence variants were recorded in western Poland where only the IIb haplotype is considered endemic. The distribution of the other haplotypes in western Poland could thus reflect past introductions or accidental releases. The authors regarded the two locations (Drzeczkowo and Karpicko) that were first included in the western Poland populations as autochthonous catchment areas of haplotype Ia.     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

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.     

TAXONOMIC REVISIONS OF MORPHOLOGICALLY SIMILAR SPECIES FROM TWO EUGLENOID GENERA: EUGLENA (E. GRANULATA AND E. VELATA) AND EUGLENARIA (EU. ANABAENA,EU. CAUDATA,AND EU. CLAVATA)1

The establishment of epitypes (together with the emended diagnoses) for three species of Euglenaria Karnkowska, E. W. Linton et Kwiatowski [Eu. anabaena (Mainx) Karnkowska et E. W. Linton; Eu. caudata (Hübner) Karnkowska et E. W. Linton; and Eu. clavata (Skuja) Karnkowska et E. W. Linton] and two species of Euglena Ehrenberg [E. granulata (Klebs) Schmitz and E. velata Klebs] was achieved due to literature studies, verification of morphological diagnostic features (cell size, cell shape, number of chloroplasts, the presence of mucocysts), as well as molecular characters (SSU rDNA). Now all these species are easy to identify and distinguish, despite their high morphological similarity, that is, spindle‐shaped (or cylindrically spindle‐shaped) cells and parietal, lobed chloroplasts with a single pyrenoid, accompanied by bilateral paramylon caps located on both sides of the chloroplast. E. granulata is the only species in this group that has spherical mucocysts. E. velata is distinguished by the largest cells (90–150 μm) and has the highest number of chloroplasts (>30). Eu. anabaena has the fewest chloroplasts (usually 3–6), and its cells are always (whether the organism is swimming or not) spindle‐shaped or cylindrically spindle‐shaped, in contrast to the cells of Eu. clavata, which are club‐shaped (clavate) while swimming and only after stopping change to resemble the shape of a spindle or a cylindrical spindle; Eu. clavata has numerous chloroplasts (15–20). Eu. caudata is characterized by asymmetrical spindle‐shaped (fusiform) cells, that is, with an elongated rear section and a shorter front section; the number of chloroplasts normally ranges from 7 to 15.     

Coexistence of <Emphasis Type="Italic">Placobdella costata</Emphasis> (Fr. Müller, 1846) (Hirudinida: Glossiphoniidae) and mud turtle <Emphasis Type="Italic">Emys orbicularis</Emphasis>

Emys orbicularis is receding from Europe, mainly due to anthropogenic habitat changes. Its parasite, Placobdella costata, is widely distributed within both the former and the present distribution range of the host. Though closely associated with the mud turtle, it may have other hosts (birds, amphibians, reptiles). Its reproductive period coincides with that of its host’s migration to the breeding grounds, thus facilitating dispersal of the parasite. Based on literature data we have analyzed the geographic spread of P. costata and mud turtle to observe the possible overlap of their habitats. Observations on the population of mud turtle and the associated leech species (P. costata) were carried out in eastern Poland — Podlasie Lowland. The studies were conducted in spring and summer in 1986–1993. The leeches were collected from the turtles caught in the water and on land. Observations showed that most leeches were found on turtles inhabiting the lake or moving to a breeding area. The greatest intensity of invasion was observed in June and July and that most leeches were observed in female E. orbicularis characterized by greater length of the carapace and weight, compared with males and juvenile individuals.     

Big data can cause big mistakes: using the Societas Europaea Herpetologica atlas by Sillero et al. (2014), the distribution of <Emphasis Type="Italic">Emys orbicularis</Emphasis> will be misunderstood

Problems with the reliability of the distribution map for Emys orbicularis in the online atlas of European amphibians and reptiles (NA2RE) are highlighted. Two distinct species (E. orbicularis, E. trinacris) are lumped together there. For E. orbicularis, many regions are included into the native distribution range where the species does not occur naturally. On the contrary, records are missing for vast parts of the native range. Using data from the NA2RE atlas cannot be recommended presently, especially not for species distribution modelling.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Molecular phylogeny of Central and South American slider turtles: implications for biogeography and systematics (Testudines: Emydidae: Trachemys)

We analyse phylogeny, systematics and biogeography of slider turtles (Trachemys spp.) using sequence data of four mitochondrial genes (3242 bp) and five nuclear loci (3396 bp) of most South American and southern Central American taxa and representatives of northern Central American, West Indian and North American slider species (16 species and subspecies) and allied North American species (genera Chrysemys, Deirochelys, Graptemys, Malaclemys, Pseudemys). By applying maximum likelihood, relaxed molecular clock and ancestral range analyses, we provide evidence for two successive colonizations of South America by slider turtles. In addition, we show that the current species delineation of Central and South American slider turtles is incorrect. Our data suggest that Trachemys grayi is a distinct polytypic species that embraces, besides the nominotypical subspecies, T. g. emolli and T. g. panamensis. Trachemys ornata is also polytypic with the subspecies T. o. ornata, T. o. callirostris, T. o. cataspila, T. o. chichiriviche and T. o. venusta. Moreover, T. adiutrix should be regarded as a subspecies of T. dorbigni. All studied Trachemys species are inferred to have originated in the Late Miocene to Early Pliocene. The ancestor of the two subspecies of T. dorbigni colonized South America most probably prior to the establishment of the land bridge connecting Central and South America, whereas the two South American subspecies of T. ornata represent a younger independent immigration wave from Central America.     

(Not) far from home: No sex bias in dispersal,but limited genetic patch size,in an endangered species,the Spotted Turtle (Clemmys guttata)

Sex-biased dispersal is common in many animals, with male-biased dispersal often found in studies of mammals and reptiles, including interpretations of spatial genetic structure, ostensibly as a result of male–male competition and a lack of male parental care. Few studies of sex-biased dispersal have been conducted in turtles, but a handful of studies, in saltwater turtles and in terrestrial turtles, have detected male-biased dispersal as expected. We tested for sex-biased dispersal in the endangered freshwater turtle, the spotted turtle (Clemmys guttata) by investigating fine-scale genetic spatial structure of males and females. We found significant spatial genetic structure in both sexes, but the patterns mimicked each other. Both males and females typically had higher than expected relatedness at distances <25 km, and in many distance classes greater than 25 km, less than expected relatedness. Similar patterns were apparent whether we used only loci in Hardy–Weinberg equilibrium (n = 7) or also included loci with potential null alleles (n = 5). We conclude that, contrary to expectations, sex-biased dispersal is not occurring in this species, possibly related to the reverse sexual dimorphism in this species, with females having brighter colors. We did, however, detect significant spatial genetic structure in males and females, separate and combined, showing philopatry within a genetic patch size of <25 km in C. guttata, which is concerning for an endangered species whose populations are often separated by distances greater than the genetic patch size.     

Demographic and genetic status of an isolated population of bog turtles (<Emphasis Type="Italic">Glyptemys muhlenbergii</Emphasis>): implications for managing small populations of long-lived animals

In this study, we sought to determine the population stability and genetic diversity of one isolated population of the federally-threatened bog turtle (Glyptemys muhlenbergii) in North Carolina. Using capture–recapture data, we estimated adult survival and population growth rate from 1992 to 2007. We found that the population decreased from an estimated 36 adult turtles in 1994 to approximately 11 adult turtles in 2007. We found a constant adult survival of 0.893 (SE = 0.018, 95% confidence interval, 0.853–0.924) between 1992 and 2007. Using 18 microsatellite markers, we compared the genetic status of this population with five other bog turtle populations. The target population displayed allelic richness (4.8 ± 0.5) and observed heterozygosity (0.619 ± 0.064) within the range of the other bog turtle populations. Coalescent analysis of population growth rate, effective population size, and timing of population structuring event also indicated the genetics of the target population were comparable to the other populations studied. Estimates of effective population size were a proportion of the census size in all populations except the target population, in which the effective population size was larger than the census size (30 turtles vs. 11 turtles). We attribute the high genetic diversity in the target population to the presence of multiple generations of old turtles. This study illustrates that the demographic status of populations of long-lived species may not be reflected genetically if a decline occurred recently. Consequently, the genetic integrity of populations of long-lived animals experiencing rapid demographic bottlenecks may be preserved through conservation efforts effective in addressing demographic problems.     

Turtle origins: insights from phylogenetic retrofitting and molecular scaffolds

Adding new taxa to morphological phylogenetic analyses without substantially revising the set of included characters is a common practice, with drawbacks (undersampling of relevant characters) and potential benefits (character selection is not biased by preconceptions over the affinities of the ‘retrofitted’ taxon). Retrofitting turtles (Testudines) and other taxa to recent reptile phylogenies consistently places turtles with anapsid‐grade parareptiles (especially Eunotosaurus and/or pareiasauromorphs), under both Bayesian and parsimony analyses. This morphological evidence for turtle–parareptile affinities appears to contradict the robust genomic evidence that extant (living) turtles are nested within diapsids as sister to extant archosaurs (birds and crocodilians). However, the morphological data are almost equally consistent with a turtle–archosaur clade: enforcing this molecular scaffold onto the morphological data does not greatly increase tree length (parsimony) or reduce likelihood (Bayesian inference). Moreover, under certain analytic conditions, Eunotosaurus groups with turtles and thus also falls within the turtle–archosaur clade. This result raises the possibility that turtles could simultaneously be most closely related to a taxon traditionally considered a parareptile (Eunotosaurus) and still have archosaurs as their closest extant sister group.