A review of phylogeographic analyses of animal taxa from the Aegean and surrounding regions

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

The Silk roads: phylogeography of Central Asian dice snakes (Serpentes: Natricidae) shaped by rivers in deserts and mountain valleys

Influenced by rapid changes in climate and landscape features since the Miocene, widely distributed species provide suitable models to study the environmental impact on their evolution and current genetic diversity. The dice snake Natrix tessellata, widely distributed in the Western Palearctic is one such species. We aimed to resolve a detailed phylogeography of N. tessellata with a focus on the Central Asian clade with 4 and the Anatolia clade with 3 mitochondrial lineages, trace their origin, and correlate the environmental changes that affected their distribution through time. The expected time of divergence of both clades began at 3.7 Mya in the Pliocene, reaching lineage differentiation approximately 1 million years later. The genetic diversity in both clades is rich, suggesting different ancestral areas, glacial refugia, demographic changes, and colonization routes. The Caspian lineage is the most widespread lineage in Central Asia, distributed around the Caspian Sea and reaching the foothills of the Hindu Kush Mountains in Afghanistan, and Eastern European lowlands in the west. Its distribution is limited by deserts, mountains, and cold steppe environments. Similarly, Kazakhstan and Uzbekistan lineages followed the Amu Darya and the Syr Darya water systems in Central Asia, with ranges delimited by the large Kyzylkum and Karakum deserts. On the western side, there are several lineages within the Anatolia clade that converged in the central part of the peninsula with 2 being endemic to Western Asia. The distribution of both main clades was affected by expansion from their Pleistocene glacial refugia around the Caspian Sea and in the valleys of Central Asia as well as by environmental changes, mostly through aridification.

Another new species of the genus Habrophlebia Eaton, 1881 (Ephemeroptera, Leptophlebiidae) from the Maghreb

A new species belonging to the genus Habrophlebia Eaton, 1881 is described at the nymphal stage from the Rif Mountains of Morocco. This species presents unique features, such as the chorionic arrangement of the egg and the ornamentation of the posterior margin of abdominal tergites. It is compared to all west European Habrophlebia species and a table with discriminating characters is given. A phylogenetic reconstruction based on COI sequences fully supports the hypothesis of a new species in the Rif Mountains, with possible further distribution in southern Spain.

Baetidae (Insecta: Ephemeroptera) of Aurès Mountains (Algeria): A New Species of the Baetis alpinus Species Group, with Notes on Baetis Laech, 1815 Biogeography within Maghreb

A new species, Baetis (Baetis) dihyaesp. nov., belonging to the Baetis alpinus species group, is described and illustrated based on larval material collected in the Aurès Mountains (northeastern Algeria) in 2020-2021. This new species is closely related to three European species, e.g., Baetis (B.) alpinus (Pictet, 1843); B. (B.) nubecularis Eaton, 1898; and B. (B.) pasquetorum Righetti & Thomas, 2002 by the combination of the following characteristics: (i) more than one short, stout bristle at the tip of segment II of the maxillary palp and (ii) a well-developed paracercus. However, the new species clearly differs from all congeners of the Baetis alpinus species group primarily by the (a) structure of mouthparts-with 14-18 long submarginal setae arranged in a single irregular row on the dorsal surface of the labrum; 2-6 short, stout bristles at the tip of segment II of the maxillary palp; and segment II of the labial palp without a considerably developed apico-internal lobe); (b) setation of abdominal terga, with a few triangular-shaped scales sparsely scattered near the posterior margin only; and (c) a well-developed paracercus, comprised of more than 50 segments. Primary data on the biology and distribution of this new species are provided, and molecular affinities are verified by the analysis of COI (barcode) sequences. Detailed notes on the distribution of mayfly species belonging to the Baetis alpinus species group common in Western Europe and the western part of North Africa are presented. The historical movement of Baetis representatives between Europe, North West Africa, and subsequently Algeria, with the land bridges 'Strait of Gibraltar' and 'Strait of Sicily' as colonization routes, is discussed in detail and identified in the present study as the Western Algeria colonization path and Eastern Algeria colonization path, respectively.

Phylogenomic and population genomic analyses reveal the spatial-temporal dynamics of diversification of the Nigella arvensis complex (Ranunculaceae) in the Aegean archipelago

The continental-shelf islands of the Aegean Sea provide an ideal geographical setting for evolutionary-biogeographical studies but disentangling the relationships between palaeogeographical history and the times, orders of modes of taxon divergence is not straightforward. Here, we used phylogenomic and population genomic approaches, based on orthologous gene sequences and transcriptome-derived SNP data, to reconstruct the spatial-temporal evolution of the Aegean Nigella arvensis complex (Ranunculaceae; 11 out of 12 taxa). The group's early diversification in the Early/Mid-Pliocene (c. 3.77 Mya) resulted in three main lineages (Greek mainland vs. central Aegean + Turkish mainland/eastern Aegean islands), while all extant taxa are of Late Plio-/Early Pleistocene origin (c. 3.30-1.59 Mya). Demographic modelling of the outcrossing taxa uncovered disparate modes of (sub)speciation, including divergence with gene flow on the Greek mainland, para- or peripatric diversification across eastern Aegean islands, and a 'mixing-isolation-mixing (MIM)' mode of subspeciation in the Cyclades. The two selfing species (N. stricta, N. doerfleri) evolved independently from the outcrossers. Present-day island configurations are clearly insufficient to explain the spatial-temporal history of lineage diversification and modes of (sub)speciation in Aegean Nigella. Moreover, our identification of positively selected genes in almost all taxa calls into question that this plant group represents a case of 'non-adaptive' radiation. Our study revealed an episodic diversification history of the N. arvensis complex, giving new insight into the modes and drivers of island speciation and adaption across multiple spatiotemporal scales.

Genotyping the phenotypic diversity in Aegean Natrixnatrixmoreotica (Bedriaga, 1882) (Reptilia, Serpentes, Natricidae)

We examined the mitochondrial identity of Aegean Natrixnatrixmoreotica representing different morphotypes, with a focus on new material from Milos and Skyros. We found no correlation between distinct morphotypes and mitochondrial identity. Our results support that grass snake populations are polyphenetic and that southern subspecies, including island populations, show a higher variability than northern ones.

A review of Eupholidoptera (Orthoptera, Tettigoniidae) from Crete, Gavdos, Gavdopoula, and Andikithira

Being nocturnal, hiding in prickly bushes and shrubs during the day, Eupholidoptera species in Crete and its neighbouring islands are easily overlooked, and until now our knowledge about their distribution was based on some thirty sightings across 11 species. In this paper results are presented of a study of Eupholidoptera specimens collected between 1987 and 2020 by hand-catches and pitfall and fermenting traps on the Greek islands of Crete, Gavdos, Gavdopoula, and Andikithira. Diagnostic features of all known species are presented and illustrated with stacked images. An updated key to all species is provided. Eupholidoptera francisae Tilmans & Odé, sp. nov. from Andikithira and southwestern Crete and Eupholidoptera marietheresae Willemse & Kotitsa, sp. nov. from Mt. Dikti are described. Female E. cretica, E. gemellata, and E. mariannae are described, and the female of E. astyla is redescribed. Bioacoustics for E. francisae Tilmans & Odé, sp. nov., E. giuliae, and E. jacquelinae are presented for the first time. Eupholidoptera smyrnensis is reported for the first time from Crete. A substantial amount of new distribution data for Eupholidoptera species on Crete is presented. The current distribution pattern and first analyses of phylogeny based on molecular data of Eupholidoptera species on Crete are discussed in relation to paleogeographical events.

The evolutionary journey of freshwater crabs of the genus Potamon (Decapoda: Brachyura: Potamidae)

AIM

Freshwater ecosystems host a rich biodiversity, including freshwater crabs. The family Potamidae is a diverse group of freshwater crabs with a Palearctic and Oriental distribution. Specifically, the genus Potamon is found in the Middle East, Southern Europe, and Northern Africa. Potamon species are considered true freshwater crabs due to their total independence from the marine environment, which makes them idea organisms to study historical freshwater connectivity dynamics. In this study we aim to elucidate the phylogenetic relationships of the genus Potamon and discuss the historical freshwater connectivity in the Mediterranean region.

METHODS

Here we use up to eleven genetic markers to reconstruct a robust phylogenetic tree of the group and use Bayesian methods to time-calibrate the tree. We further use historical biogeography methods implemented in a Bayesian framework to assess the ancestral distribution ranges of the group.

RESULTS

In this study we generate the most complete dataset for the group, covering practically the whole distribution range of the genus. We obtained a robust phylogenetic hypothesis and evaluated the times of divergence of the group. The result of the historical biogeography shows the evolution of historical distribution ranges of species of Potamon.

MAIN CONCLUSIONS

The Mediterranean Sea is placed in a key intersection for the exchange of fauna and flora. Here we present evidence for the potential of freshwater fauna exchange, through the Mediterranean Sea after a period of desiccation and the loss of contact with the ocean. The origin of Potamon genus is found to be in Western Asia, probably Eastern Iran.

Stability in the South, Turbulence Toward the North: Evolutionary History of Aurinia saxatilis (Brassicaceae) Revealed by Phylogenomic and Climatic Modelling Data

The Balkan Peninsula played an important role in the evolution of many Mediterranean plants and served as a major source for post-Pleistocene colonisation of central and northern Europe. Its complex geo-climatic history and environmental heterogeneity significantly influenced spatiotemporal diversification and resulted in intricate phylogeographic patterns. To explore the evolutionary dynamics and phylogeographic patterns within the widespread eastern Mediterranean and central European species Aurinia saxatilis, we used a combination of phylogenomic (restriction-site associated DNA sequencing, RADseq) and phylogenetic (sequences of the plastid marker ndhF) data as well as species distribution models generated for the present and the Last Glacial Maximum (LGM). The inferred phylogenies retrieved three main geographically distinct lineages. The southern lineage is restricted to the eastern Mediterranean, where it is distributed throughout the Aegean area, the southern Balkan Peninsula, and the southern Apennine Peninsula, and corresponds to the species main distribution area during the LGM. The eastern lineage extends from the eastern Balkan Peninsula over the Carpathians to central Europe, while the central lineage occupies the central Balkan Peninsula. Molecular dating places the divergence among all the three lineages to the early to middle Pleistocene, indicating their long-term independent evolutionary trajectories. Our data revealed an early divergence and stable in situ persistence of the southernmost, eastern Mediterranean lineage, whereas the mainland, south-east European lineages experienced more complex and turbulent evolutionary dynamics triggered by Pleistocene climatic oscillations. Our data also support the existence of multiple glacial refugia in southeast Europe and highlight the central Balkan Peninsula not only as a cradle of lineage diversifications but also as a source of lineage dispersal. Finally, the extant genetic variation within A. saxatilis is congruent with the taxonomic separation of peripatric A. saxatilis subsp. saxatilis and A. saxatilis subsp. orientalis, whereas the taxonomic status of A. saxatilis subsp. megalocarpa remains doubtful.

First occurrence of Psammophis (Serpentes) from Europe witnesses another Messinian herpetofaunal dispersal from Africa - biogeographic implications and a discussion of the vertebral morphology of psammophiid snakes

We here describe abundant new snake material from the late Miocene (MN 13) of Salobreña, Spain. Vertebral morphology suggests a referral of the specimens to the extant psammophiid Psammophis, documenting the first occurrence of this genus in Europe. The diversity and disparity across the vertebral morphology of different psammophiid genera is discussed. We identify vertebral features that could diagnose Psammophis and therefore enable the recognition of the genus in the fossil record. A comparison of the new Spanish form with other taxa is conducted. We provide a detailed review of the psammophiid fossil record. Material previously described from the middle Miocene of Beni Mellal, Morocco is here tentatively referred to as? Psammophis sp., an action that renders that occurrence as the oldest (probable) record of the genus and Psammophiidae as a whole, providing thus a potential calibration point. On the other hand, Eastern European Pliocene material that had been previously supposedly referred to Psammophis is here discarded as being rather fragmentary, not affording any more precise determination. The two psammophiid genera Psammophis and Malpolon appear almost simultaneously in the European fossil record (MN 13), with the former achieving only a short-lived and apparently geographically limited distribution in the continent, while the latter still exists in its modern herpetofauna. We assess biogeographic implications of the new find, suggesting a direct dispersal event from northwestern Africa to the Iberian Peninsula during the late Miocene, facilitated by the Messinian Salinity Crisis. This article is protected by copyright. All rights reserved.

The story of a rock-star: multilocus phylogeny and species delimitation in the starred or roughtail rock agama, Laudakia stellio (Reptilia: Agamidae)

Situated at the junction of three continents, Europe, Asia and Africa, the Eastern Mediterranean is an ideal region to study the effects of palaeogeography, ecology and long human presence on animal evolution. Laudakia stellio (Squamata: Agamidae) is found across this region and offers an excellent opportunity for such studies. The high morphological variation across their range suggests that these lizards might represent a species complex. This is the first study exploring their evolutionary history, using molecular markers and individuals from all described subspecies. We employed the latest phylogenetic and species-delimitation methods to identify all distinct evolutionary lineages, their genetic variation and divergence times. The phenotypical diversity of L. stellio matches its genetic differentiation: almost all subspecies correspond to well-supported retrieved subclades and additional distinct lineages representing intermediate morphs have been retrieved. ‘Laudakia stellio’ represents three distinct evolutionary entities that diverged during the Plio-Pleistocene transition, which we propose as distinct species. One includes Greek and Turkish populations, as well as cryptic Anatolian lineages. The second comprises all other Near East populations and the third is endemic to Cyprus. Our results indicate a role of humans in shaping present distribution patterns, and highlight the importance of the Aegean, Anatolia and the Levant as glacial refugia and diversity hotspots.

Hidden diversity revealed in the freshwater snails, Bythinella and Pseudamnicola, in the Island of Crete

Crete with its complex geomorphological history is the island with the highest number of endemism observed in animal and plant taxa throughout the Aegean archipelago. While other groups of organisms within Crete are well-studied, the freshwater gastropod fauna still remains poorly investigated. Bythinella and Pseudamnicola, are 2 genera of freshwater springsnails, both present on the island, inhabiting springs and other freshwater habitats. Here, we conduct a comprehensive study on the distribution of the different genetic lineages of the 2 gastropod genera in order to assess the mode of their differentiation on the island and infer the actual number of species present in the island. Towards these aims, sequence data from the mitochondrial gene were used and analyzed within a phylogenetic framework. For Bythinella, our results strongly support at least 5 delineated Bythinella spp. inhabiting Crete, which correspond to the already described species from previous studies with the addition of a new one. Bythinella analyses reveal an old time-frame of differentiation with vicariant phenomena being more likely the main drivers shaping the present-day distribution of the genus' genetic lineages. For Pseudamnicola, our data indicate the presence of at least 2 delineated Pseudamnicola spp. with a differentiation more consistent to an isolation-by-distance pattern of a relatively recent origin. Dispersion processes followed by isolation of the populations and/or recent speciation, seem to be the underlying process for the current distribution of Pseudamnicola lineages.

Spatial and Ecological Drivers of Genetic Structure in Greek Populations of Alkanna tinctoria (Boraginaceae), a Polyploid Medicinal Herb

Background and Aims: Quantifying genetic variation is fundamental to understand a species' demographic trajectory and its ability to adapt to future changes. In comparison with diploids, however, genetic variation and factors fostering genetic divergence remain poorly studied in polyploids due to analytical challenges. Here, by employing a ploidy-aware framework, we investigated the genetic structure and its determinants in polyploid Alkanna tinctoria (Boraginaceae), an ancient medicinal herb that is the source of bioactive compounds known as alkannin and shikonin (A/S). From a practical perspective, such investigation can inform biodiversity management strategies. Methods: We collected 14 populations of A. tinctoria within its main distribution range in Greece and genotyped them using restriction site-associated DNA sequencing. In addition, we included two populations of A. sieberi. By using a ploidy-aware genotype calling based on likelihoods, we generated a dataset of 16,107 high-quality SNPs. Classical and model-based analysis was done to characterize the genetic structure within and between the sampled populations, complemented by genome size measurements and chromosomal counts. Finally, to reveal the drivers of genetic structure, we searched for associations between allele frequencies and spatial and climatic variables. Key Results: We found support for a marked regional structure in A. tinctoria along a latitudinal gradient in line with phytogeographic divisions. Several analyses identified interspecific admixture affecting both mainland and island populations. Modeling of spatial and climatic variables further demonstrated a larger contribution of neutral processes and a lesser albeit significant role of selection in shaping the observed genetic structure in A. tinctoria. Conclusion: Current findings provide evidence of strong genetic structure in A. tinctoria mainly driven by neutral processes. The revealed natural genomic variation in Greek Alkanna can be used to further predict variation in A/S production, whereas our bioinformatics approach should prove useful for the study of other non-model polyploid species.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Past connections with the mainland structure patterns of insular species richness in a continental-shelf archipelago (Aegean Sea, Greece)

Recent research in island biogeography has highlighted the important role of late Quaternary sea-level fluctuations in shaping biogeographic patterns in insular systems but focused on oceanic systems. Through this study, we aim investigate how late Quaternary sea-level fluctuations shaped species richness patterns in continental-shelf island systems. Focusing on the Aegean archipelago, we first compiled maps of the area's geography using published data, under three sea-level stands: (a) current; (b) median sea-level over the last nine glacial-interglacial cycles (MSL); and (c) Last Glacial Maximum (LGM). We gathered taxon-island occurrences for multiple chorotypes of angiosperms, butterflies, centipedes, and reptiles. We investigated the impact of present-day and past geographic settings on chorological groups by analyzing island species-area relationships (ISARs) and using generalized linear mixed models (GLMMs) selection based on multiple metrics of goodness of fit. Our results confirm that the Aegean's geography has changed dramatically since the LGM, whereas the MSL only modestly differs from the present configuration. Apart for centipedes, paleogeographic changes affected both native and endemic species diversity through altering connections between land-bridge islands and the mainland. On land-bridge islands, we detected over-representation of native species and under-representation of endemics. Unlike oceanic islands, sea-level-driven increase of isolation and area contraction did not strongly shape patterns of species richness. Furthermore, the LGM configurations rather than the MSL configuration shaped patterns of endemic species richness. This suggests that even short episodes of increased connectivity with continental populations are sufficient to counteract the genetic differentiation of insular populations. On the other hand, the over-representation of native nonendemic species on land-bridge islands reflected MSL rather than LGM mainland connections. Our study shows that in terms of processes affecting species richness patterns, continental archipelagos differ fundamentally from oceanic systems because episodic connections with the mainland have profound effects on the biota of land-bridge islands.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Geographic isolation and human-assisted dispersal in land snails: a Mediterranean story of Helix borealis and its relatives (Gastropoda: Stylommatophora: Helicidae)

The Mediterranean basin is a major centre for land-snail diversity, with many localized endemics, but there are also species widely spread by humans. Both endemics and introductions can be found in the snail genus Helix, which comprises many large-bodied species used for human consumption in the past and present. The Mediterranean clade of Helix is currently distributed throughout this region, but the phylogenetic and biogeographic relationships among its forms from different parts of the basin remain enigmatic. The reasons include insufficient sampling, taxa with unclear taxonomy and a significant impact of human-assisted transport obscuring the natural distribution of phylogenetic lineages. We provide evidence that European and Anatolian populations of H. cincta and its relatives are not native to those regions, but originate from the northern Levant. These results have implications for taxonomy of the genus, but also for the understanding of its evolutionary history. We posit that the Mediterranean clade consists of four geographically separated groups, which diversified in Northern Africa, the Apennine Peninsula and Corsica, the Aegean and Greece, and the northern Levant. This geographic pattern has been subsequently blurred by multiple instances of human-assisted dispersal. However, revealing the founding populations with certainty requires thorough sampling in currently inaccessible countries.

Evolution of Poecilimon jonicus group (Orthoptera: Tettigoniidae): a history linked to the Aegean Neogene paleogeography

The Aegean archipelago is among the largest on Earth with astonishing biodiversity within Europe. Its territory underwent a massive geotectonic transformation in Neogene that resulted in multitude of changes in land-sea configuration and disintegrated the formerly united Aegean land to a complicated mainland-archipelago system. Therefore, it represents an excellent laboratory for studying evolution of terrestrial fauna. In the present study, we use a model group of flightless bush crickets with annual reproduction cycle—Poecilimon jonicus species group—to trace correlation of lineage diversification with the known paleogeographic events in the Aegean area. The group belongs to the hyperdiverse genus Poecilimon and has a disjunct distribution along the Hellenic arc from southwestern Anatolia through Crete to the western Balkans and the Apennines. To test our hypothesis, we inferred phylogenetic relationships of the P. jonicus group sensu lato using a nuclear fragment covering two spacers of the ribosomal cistron (ITS1 + ITS2). To study intra-group phylogeny, we compared mitochondrial phylogenies based on two matrices—(1) a concatenated ND2 and COI dataset of 1656 bp and (2) a 16S rRNA + 12S rRNA dataset of 1835 bp. As a second step, we estimated divergence times applying Bayesian approach with BEAST and a relative rate framework with RelTime on the mitochondrial matrices. We compare trees calibrated based on evolutionary rates and tectonic events and discuss radiation scenarios in concordance with known paleogeographic events in the Aegean area. Our results revealed robust phylogeny of the Poecilimon jonicus group and confirmed a strong link between its evolution and the Aegean paleogeography. The phylogenetic relationships of the group supported reconsideration of its systematics.

Molecular data suggest multiple origins and diversification times of freshwater gammarids on the Aegean archipelago

Our main aim was to investigate the diversity, origin and biogeographical affiliations of freshwater gammarids inhabiting the Aegean Islands by analysing their mtDNA and nDNA polymorphism, thereby providing the first insight into the phylogeography of the Aegean freshwater gammarid fauna. The study material was collected from Samothraki, Lesbos, Skyros, Evia, Andros, Tinos and Serifos islands as well as from mainland Greece. The DNA extracted was used for amplification of two mitochondrial (COI and 16S) and two nuclear markers (28S and EF1-alpha). The multimarker time-calibrated phylogeny supports multiple origins and different diversification times for the studied taxa. Three of the sampled insular populations most probably represent new, distinct species as supported by all the delimitation methods used in our study. Our results show that the evolution of freshwater taxa is associated with the geological history of the Aegean Basin. The biogeographic affiliations of the studied insular taxa indicate its continental origin, as well as the importance of the land fragmentation and the historical land connections of the islands. Based on the findings, we highlight the importance of studying insular freshwater biota to better understand diversification mechanisms in fresh waters as well as the origin of studied Aegean freshwater taxa.

Dictyocoela microsporidia diversity and co-diversification with their host, a gammarid species complex (Crustacea, Amphipoda) with an old history of divergence and high endemic diversity

BACKGROUND

Although the processes of co-evolution between parasites and their hosts are well known, evidence of co-speciation remains scarce. Microsporidian intracellular parasites, due to intimate relationships with their hosts and mixed mode of transmission (horizontal but also vertical, from mother to offspring), may represent an interesting biological model for investigating co-speciation. Amphipod crustaceans, especially gammarids, are regular hosts of microsporidian parasites, in particular the Dictyocoela spp., which have so far been found limited to these amphipods and are known to use a vertical mode of transmission. The amphipod genus Gammarus has a diversification history spanning the last 50-60 Mya and an extensive cryptic diversity in most of the nominal species. Here, we investigated the degree of co-diversification between Dictyocoela and Gammarus balcanicus, an amphipod with high degrees of ancient cryptic diversification and lineage endemism, by examining the genetic diversity of these parasites over the entire geographic range of the host. We hypothesised that the strong host diversification and vertical transmission of Dictyocoela would promote co-diversification.

RESULTS

Using the parasite SSU rDNA as a molecular marker, analyzing 2225 host specimens from 88 sites covering whole host range, we found 31 haplogroups of Dictyocoela, 30 of which were novel, belonging to four Dictyocoela species already known to infect other Gammarus spp. The relationships between Dictyocoela and gammarids is therefore ancient, with the speciation in parasites preceding those of the hosts. Each novel haplogroup was nevertheless specific to G. balcanicus, leaving the possibility for subsequent co-diversification process during host diversification. A Procrustean Approach to Co-phylogeny (PACo) analysis revealed that diversification of Dictyocoela was not random with respect to that of the host. We found high degrees of congruence between the diversification of G. balcanicus and that of Dictyocoela roeselum and D. muelleri. However, we also found some incongruences between host and Dictyocoela phylogenies, e.g. in D. duebenum, probably due to host shifts between different G. balcanicus cryptic lineages.

CONCLUSION

The evolutionary history of Dictyocoela and Gammarus balcanicus represents an example of an overall host-parasite co-diversification, including cases of host shifts.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.

Dispersal versus vicariance in the Aegean: combining molecular and morphological phylogenies of eastern Mediterranean Dendarus (Coleoptera: Tenebrionidae) sheds new light on the phylogeography of the Aegean area

The Aegean archipelago, as an ‘evolutionary laboratory of nature’, is an ideal model for research in phylogeography. In this area, the darkling beetles of the genus Dendarus (distributed from Morocco to the Caucasus) exhibit a high level of diversity with 36 species, 27 of which are island endemics. However, their taxonomy is complex and unstable, having undergone continuous revision to address extensive morphological and ecological plasticity. Here, we examine the phylogenetic relationships of 23 species from Greece and Turkey, using mitochondrial and nuclear DNA sequences and 61 morphological characters, to unveil their phylogeny in the Aegean. This helps to clarify phylogeographic scenarios and historical processes that shaped the observed patterns. The analyses reveal 13 distinct lineages with several para- and polyphyletic cases that correspond to three major phylogroups [south/south-east Aegean (D. foraminosus complex, D. rhodius, D. sporadicus, D. wettsteini); central to north Aegean, Turkey and mainland Greece (D. crenulatus, D. moesiacus group, D. sinuatus complex, D. stygius) and mainland Greece (D. messenius, D. paganettii)], indicating the need for further taxonomic re-evaluation. Lineage topology and phylogeography suggest a spatial and temporal sequence of geographic isolation, following either a vicariant or a dispersal model coincident with major palaeogeographic separations in the Aegean.

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co-occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east-west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F1 individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.

Changes in Fish Taxonomy Affect Freshwater Biogeographical Regionalisations: Insights from Greece

Freshwater fishes are key indicators for delineating biogeographical maps worldwide. However, controversy in regional-scale ichthyogeographic boundaries still persists, especially in areas of high species endemicity, such as in Greece. One problem concerns the taxonomy of the fishes because there have been extensive changes, mainly due to an increased splitting of species in recent years in Europe. Here, we explore why ichthyogeographic boundary disagreements and uncertainties in region-scale biogeographical units persist. We compare cluster analyses of river basin fish fauna in Greece using two taxonomic datasets: the older fish taxonomy (from 1991) and the current taxonomy that now follows the phylogenetic species concept (PSC), which has become widely established in Europe after 2007. Cluster analyses using the older fish taxonomy depicts only two major biogeographical regional divisions, while the current taxonomy defines four major regional divisions in mainland Greece. Interestingly, some older maps from the pre-PSC taxonomy era also similarly show four ichthyogeographic divisions in Greece and we can assume that the older biogeographical work did not solely use numerical taxonomy but followed an expert-guided synthesis; the older regional definitions have persisted quite well despite radical changes in Europe’s fish taxonomy. Through the prism of biodiversity conservation planning, we hope this review may help identify ways to help standardize policy-relevant biogeographical mapping.

Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae)

The six-toothed bark beetle, Ips sexdentatus, is one of the most abundant scolytid species of the central and southern European countries. It mostly feeds on Pinus sp., whereas during population outbreaks it can also attack Picea sp. In spite of its broad distribution, its phylogeography has never been studied before. To do that, we employed an mtDNA marker on 489 individuals that covered most of its native range in Europe. Geographic distribution of the 86 haplotypes showed that at least three glacial refugia have played a significant role in shaping the currently observed pattern of genetic divergence in Europe, without excluding the contribution of minor refugial areas that acted in a similar manner. The revealed shallow structure can be considered an artifact of factors that reduced intraspecific diversity, at the same time favoring gene flow. As such, biological traits of the species itself (flying ability and host preference) and even human-mediated transport of wood seem to be the most prevailing and probable reasons that gave rise to the observed pattern.

Persistence of phylogeographic footprints helps to understand cryptic diversity detected in two marine amphipods widespread in the Mediterranean basin

Amphipods of the genus Gammarus are a vital component of macrozoobenthic communities in European inland and coastal, marine and brackish waters of the Mediterranean and the Black Sea. Exceptional levels of cryptic diversity have been revealed for several widespread freshwater Gammarus species in Europe. No comprehensive assessment has yet been made for brackishwater counterparts, such as Gammarus aequicauda and G. insensibilis, which are among the most widely dispersed members of the so-called "G. locusta group" in the Mediterranean and in the Black Sea. Here we probe the diversity of these morphospecies examining the partitioning of mtDNA and nDNA across multiple populations along their distribution range and discuss it within the regional paleogeographic framework. We gathered molecular data from a collection of 166 individuals of G. aequicauda and G. insensibilis from 47 locations along their distribution range in the Mediterranean including the Black Sea. They were amplified for both mitochondrial COI and 16S rRNA as well as the nuclear 28S rRNA. All five MOTU delimitation methods (ABGD, BIN, bPTP, GMYC single and multiple threshold models) applied revealed deep divergence between Black Sea and Mediterranean populations in both G. aequicauda and G. insensibilis. There were eight distinct MOTUs delimited for G. aequicauda (6-18% K2P) and 4 MOTUs for G. insensibilis (4-14% K2P). No sympatric MOTUs were detected throughout their distribution range. Multimarker time-calibrated phylogeny indicated that divergence of both G. aequicauda and G. insensibilis species complexes started already in the late Oligocene/early Miocene with the split between clades inhabiting eastern and western part of the Mediterranean occurring in both species at the similar time. Our results indicate a high cryptic diversity within Mediterranean brackishwater Gammarus, similar to that observed for freshwater counterparts. Moreover, the phylogenetic history combined with the current geographic distribution indicate that the evolution of both studied Gammarus morphogroups has been strongly connected with the geological events in the Mediterranean Basin and it reflect the turbulent history of the area.

First endemic freshwater Gammarus from Crete and its evolutionary history-an integrative taxonomy approach

The Mediterranean islands are known as natural laboratories of evolution with a high level of endemic biodiversity. However, most biodiversity assessments have focused mainly on terrestrial and marine fauna, leaving the freshwater animals aside. Crete is one of the largest islands in the Mediterranean Basin, with a long history of isolation from the continental mainland. Gammarid amphipods are often dominant in macrozoobenthic communities in European inland waters. They are widely used in biomonitoring and exotoxicological studies. Herein, we describe Gammarus plaitisi sp. nov., endemic to Cretan streams, based on morphological characters and a set of molecular species delimitation methods using mitochondrial cytochrome oxidase subunit I and 16S rRNA genes as well as nuclear 28S rDNA, ITS1 and EF1-alpha genes. The divergence of the new species is strongly connected with the geological history of the island supporting its continental origin.

Diversification in continental island archipelagos: new evidence on the roles of fragmentation, colonization and gene flow on the genetic divergence of Aegean Nigella (Ranunculaceae)

Background and Aims

Disentangling the relative roles of past fragmentation (vicariance), colonization (dispersal) and post-divergence gene flow in the genetic divergence of continental island organisms remains a formidable challenge. Amplified fragment length polymorphisms (AFLPs) were used to (1) gain further insights into the biogeographical processes underlying the Pleistocene diversification of the Aegean Nigella arvensis complex; (2) evaluate the role of potential key factors driving patterns of population genetic variability (mating system, geographical isolation and historical contingencies); and (3) test the robustness of conclusions previously drawn from chloroplast (cp) DNA.

Methods

Genetic diversity was analysed for 235 AFLP markers from 48 populations (497 individuals) representing 11 taxa of the complex using population genetic methods and Bayesian assignment tests.

Key Results

Most designated taxa are identifiable as genetically distinct units. Both fragmentation and dispersal-driven diversification processes occurred at different geological time scales, from Early to Late Pleistocene, specifically (1) sea barrier-induced vicariant speciation in the Cyclades, the Western Cretan Strait and Ikaria; and (2) bi-regional colonizations of the 'Southern Aegean Island Arc' from the Western vs. Eastern Aegean mainland, followed by allopatric divergences in Crete vs. Rhodos and Karpathos/Kasos. Outcrossing island taxa experienced drift-related demographic processes that are magnified in the two insular selfing species. Population genetic differentiation on the mainland seems largely driven by dispersal limitation, while in the Central Aegean it may still be influenced by historical events (island fragmentation and sporadic long-distance colonization).

Conclusions

The biogeographical history of Aegean Nigella is more complex than expected for a strictly allopatric vicariant model of divergence. Nonetheless, the major phylogeographical boundaries of this radiation are largely congruent with the geography and history of islands, with little evidence for ongoing gene exchange between divergent taxa. The present results emphasize the need to investigate further biological and landscape features and contemporary vs. historical processes in driving population divergence and taxon diversification in Aegean plant radiations.

Neogene paleogeography provides context for understanding the origin and spatial distribution of cryptic diversity in a widespread Balkan freshwater amphipod

BACKGROUND

The Balkans are a major worldwide biodiversity and endemism hotspot. Among the freshwater biota, amphipods are known for their high cryptic diversity. However, little is known about the temporal and paleogeographic aspects of their evolutionary history. We used paleogeography as a framework for understanding the onset of diversification in Gammarus roeselii: (1) we hypothesised that, given the high number of isolated waterbodies in the Balkans, the species is characterised by high level of cryptic diversity, even on a local scale; (2) the long geological history of the region might promote pre-Pleistocene divergence between lineages; (3) given that G. roeselii thrives both in lakes and rivers, its evolutionary history could be linked to the Balkan Neogene paleolake system; (4) we inspected whether the Pleistocene decline of hydrological networks could have any impact on the diversification of G. roeselii.

MATERIAL AND METHODS

DNA was extracted from 177 individuals collected from 26 sites all over Balkans. All individuals were amplified for ca. 650 bp long fragment of the mtDNA cytochrome oxidase subunit I (COI). After defining molecular operational taxonomic units (MOTU) based on COI, 50 individuals were amplified for ca. 900 bp long fragment of the nuclear 28S rDNA. Molecular diversity, divergence, differentiation and historical demography based on COI sequences were estimated for each MOTU. The relative frequency, geographic distribution and molecular divergence between COI haplotypes were presented as a median-joining network. COI was used also to reconstruct time-calibrated phylogeny with Bayesian inference. Probabilities of ancestors' occurrence in riverine or lacustrine habitats, as well their possible geographic locations, were estimated with the Bayesian method. A Neighbour Joining tree was constructed to illustrate the phylogenetic relationships between 28S rDNA haplotypes.

RESULTS

We revealed that G. roeselii includes at least 13 cryptic species or molecular operational taxonomic units (MOTUs), mostly of Miocene origin. A substantial Pleistocene diversification within-MOTUs was observed in several cases. We evidenced secondary contacts between very divergent MOTUs and introgression of nDNA. The Miocene ancestors could live in either lacustrine or riverine habitats yet their presumed geographic localisations overlapped with those of the Neogene lakes. Several extant riverine populations had Pleistocene lacustrine ancestors.

DISCUSSION

Neogene divergence of lineages resulting in substantial cryptic diversity may be a common phenomenon in extant freshwater benthic crustaceans occupying areas that were not glaciated during the Pleistocene. Evolution of G. roeselii could be associated with gradual deterioration of the paleolakes. The within-MOTU diversification might be driven by fragmentation of river systems during the Pleistocene. Extant ancient lakes could serve as local microrefugia during that time.

The Aegean archipelago: a natural laboratory of evolution, ecology and civilisations

The Aegean archipelago, comprising numerous islands and islets with great heterogeneity in topographic, geological, historical and environmental properties, offers an ideal natural laboratory for ecological and evolutionary research, and has been the stage for a very long interaction between human civilizations and local ecosystems. This work presents insights that have been gained from past and current relevant research in the area, highlighting also the importance of the Aegean archipelago as a useful model to address many major questions in biogeography, ecology and evolutionary processes. Among the most interesting findings from such studies concern the role of habitat heterogeneity as the most important determinant of species richness, the development of a new model (Choros) for the species–area–habitats relationship, the mechanistic aspects of the Small Island Effect, the very high rates of species turnover, the lack of a role for interspecific competition in shaping species co-occurrence patterns in most cases, the importance of non adaptive radiation in diversification of several taxa, the insights into the relative roles of vicariance and dispersal in speciation, the understanding of the interplay between human presence and the establishment of exotic species and extinction of indigenous biotas. Concluding, the Aegean archipelago is an ideal stage for research in evolution, ecology and biogeography, and has the potential to become a model study area at a global level, especially for land-bridge, continental islands.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.

Phylogeny, phylogeography, and evolution in the Mediterranean region: News from a freshwater mussel (Potomida, Unionida)

The Potomida genus (Bivalvia, Unionida) has a Circum-Mediterranean distribution and like other freshwater mussel species, its populations have suffered dramatic declines. Although this genus is currently considered as monotypic, it has a long history of taxonomic revisions and presently many aspects of its systematics and evolutionary history are unclear. We sampled a total of 323 individuals from 39 different sites across the Potomida genus distribution, and sequenced two mitochondrial (16S rDNA and Cytochrome c Oxidase Subunit I) and one nuclear (28S rDNA) genes to clarify its phylogeny and phylogeographic history. Our results show that the genus includes two well-supported clades, one comprising solely the western Mediterranean species Potomida littoralis, and the other including two eastern Mediterranean species, the Greek endemic P. acarnanica and the Anatolian and Middle Eastern P. semirugata. We suggest that Potomida started radiating during the upper Miocene, and that both vicariance and dispersal events shaped the diversification and distribution of the genus along the Mediterranean region. P. littoralis is further divided in two mitochondrial lineages, one restricted to Europe and the other occurring mostly in North Africa. Moreover, some European basins present both lineages in sympatry. The conservation status of the three recognized species should be reevaluated, particularly P. acarnanica, since it is restricted to two Greek river basins presenting a high risk of extinction. Overall, our results clarify some important gaps in knowledge concerning the phylogeny, phylogeography and evolution of the Potomida genus in the Mediterranean region with important taxonomical, ecological and conservational implications.