Multigene phylogeny, phylogeography and population structure of Podarcis cretensis species group in south Balkans

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.

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.

Dispersal patterns of Squamata in the Mediterranean: An evolutionary perspective

Mediterranean islands have a high diversity of squamates, although they are unevenly distributed. This variability in the composition of the reptile assemblages across islands may have been influenced by differences in the colonization abilities of these species. To evaluate the dispersal capacities of squamate species, we modeled their sea routes using cost surface models. We estimated the effects of some life-history traits and the phylogenetic signal in the characteristics of the modeled dispersal paths. We hypothesized that a significant phylogenetic signal should be present if the dispersal ability is enhanced by traits shared among evolutionarily related species. The results showed that no phylogenetic signal was present in the characteristics of the dispersal paths (i.e., in the distance traveled/bypassed sea depth). Thus, no superior island-colonizer lineages were detected in Mediterranean Squamata. However, our analyses also revealed that small-sized lizards were superior to other groups of squamates at dispersing over long distances on the sea.

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.

Morphological and genetic differentiation in the anguid lizard Pseudopus apodus supports the existence of an endemic subspecies in the Levant

The Levant represents one of the most important reptile diversity hotspots and centers of endemism in the Western Palearctic. The region harbored numerous taxa in glacial refugia during the Pleistocene climatic oscillations. Due to the hostile arid conditions in the warmer periods they were not always able to spread or come into contact with populations from more distant regions. One large and conspicuous member of the Levantine herpetofauna is the legless anguid lizardPseudopus apodus. This species is distributed from the Balkans to Central Asia with a portion of its range running along the eastern Mediterranean coast. Mitochondrial and nuclear DNA sequences, microsatellite genotypes, and morphology show that populations in this region differ from the two named subspecies and presumably had a long independent evolutionary history during the Quaternary. Here we describe the Levantine population as a new subspecies and present biogeographic scenarios for its origin and diversification. The new subspecies is genetically highly diverse, and it forms a sister lineage toPseudopusfrom the remaining parts of the range according to mtDNA. It is the largest-bodied of the three subspecies, but occupies the smallest range.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

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.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.