Selection, drift, and introgression shape MHC polymorphism in lizards.

The major histocompatibility complex (MHC) has long served as a model for the evolution of adaptive genetic diversity in wild populations. Pathogen-mediated selection is thought to be a main driver of MHC diversity, but it remains elusive to what degree selection shapes MHC diversity in complex biogeographical scenarios where other evolutionary processes (e.g. genetic drift and introgression) may also be acting. Here we focus on two closely related green lizard species, Lacerta trilineata and L. viridis, to address the evolutionary forces acting on MHC diversity in populations with different biogeographic structure. We characterized MHC class I exon 2 and exon 3, and neutral diversity (microsatellites), to study the relative importance of selection, drift, and introgression in shaping MHC diversity. As expected, positive selection was a significant force shaping the high diversity of MHC genes in both species. Moreover, introgression significantly increased MHC diversity in mainland populations, with a primary direction of gene flow from L. viridis to L. trilineata. Finally, we found significantly fewer MHC alleles in island populations, but maintained MHC sequence and functional diversity, suggesting that positive selection counteracted the effect of drift. Overall, our data support that different evolutionary processes govern MHC diversity in different biogeographical scenarios: positive selection occurs broadly while introgression acts in sympatry and drift when the population sizes decrease.

Contribution to the study of Acanthodactylus (Sauria: Lacertidae) mtDNA diversity focusing on the A. boskianus species group.

The Acanthodactylus boskianus species group includes three species (A. boskianus, A. nilsoni, and A. schreiberi) of unclear phylogeny and phylogeographic history. By sequencing fragments of two mtDNA genes and performing phylogenetic, demographic, and chronophylogenetic analyses, we aimed at identifying their phylogenetic relationships while unravelling their biogeographic history. The analyses indicated that A. boskianus is a species complex, while A. s. schreiberi and A. s. ataturi show, both, low intraspecific genetic diversity. From a biogeographic point of view, the ancestor of A. s. schreiberi colonized Cyprus from the Middle East through overseas dispersal during the Pleistocene, whereas A. s. ataturi is considered to be a relict of a previously wider distribution.

Phylogenetic position, origin and biogeography of Palearctic and Socotran blind-snakes (Serpentes: Typhlopidae).

The majority of the family Typhlopidae occurs in the Neotropic, Australasian, Indo-Malayan and Afrotropic ecoregions. They show a restricted distribution in the western Palearctic, where they include few native species, i.e. Rhinotyphlops simoni, R. episcopus and Typhlops vermicularis. A unique species among typhlopids is T. socotranus, found in Socotra, one of the most endemic-rich archipelagoes. In this study we determine the phylogenetic position of the above mentioned species and discuss their systematics, origin and biogeography. For this purpose we use three protein-coding nuclear markers (AMEL-amelogenin, BDNF-brain-derived neurotrophic factor and NT3-neurotrophin 3) to construct a time-calibrated phylogeny of the family Typhlopidae. Our results show that T. socotranus is a sister-species to T. vermicularis, while R. simoni and R. episcopus are sister-species to each other and are found within the African clade of the family, although they are geographically distributed in west Asia. Additionally we discuss several hypotheses on their origin, as well as the occurence of typhlopids in Eurasia.

Molecular phylogenetics and historical biogeography of the west-palearctic common toads (Bufo bufo species complex).

In most pan-Eurasiatic species complexes, two phenomena have been traditionally considered key processes of their cladogenesis and biogeography. First, it is hypothesized that the origin and development of the Central Asian Deserts generated a biogeographic barrier that fragmented past continuous distributions in Eastern and Western domains. Second, Pleistocene glaciations have been proposed as the main process driving the regional diversification within each of these domains. The European common toad and its closest relatives provide an interesting opportunity to examine the relative contributions of these paleogeographic and paleoclimatic events to the phylogeny and biogeography of a widespread Eurasiatic group. We investigate this issue by applying a multiproxy approach combining information from molecular phylogenies, a multiple correspondence analysis of allozyme data and species distribution models. Our study includes 304 specimens from 164 populations, covering most of the distributional range of the Bufo bufo species complex in the Western Palearctic. The phylogenies (ML and Bayesian analyses) were based on a total of 1988 bp of mitochondrial DNA encompassing three genes (tRNAval, 16S and ND1). A dataset with 173 species of the family Bufonidae was assembled to estimate the separation of the two pan-Eurasiatic species complexes of Bufo and to date the main biogeographic events within the Bufo bufo species complex. The allozyme study included sixteen protein systems, corresponding to 21 presumptive loci. Finally, the distribution models were based on maximum entropy. Our distribution models show that Eastern and Western species complexes are greatly isolated by the Central Asian Deserts, and our dating estimates place this divergence during the Middle Miocene, a moment in which different sources of evidence document a major upturn of the aridification rate of Central Asia. This climate-driven process likely separated the Eastern and Western species. At the level of the Western Palearctic, our dating estimates place most of the deepest phylogenetic structure before the Pleistocene, indicating that Pleistocene glaciations did not have a major role in splitting the major lineages. At a shallow level, the glacial dynamics contributed unevenly to the genetic structuring of populations, with a strong influence in the European-Caucasian populations, and a more relaxed effect in the Iberian populations.

Neogene climatic oscillations shape the biogeography and evolutionary history of the Eurasian blindsnake.

Typhlops vermicularis is the only extant scolecophidian representative occurring in Europe. Its main distribution area, the eastern Mediterranean, has a complicated geological and climatic history that has left an imprint on the phylogenies and biogeography of many taxa, especially amphibians and reptiles. Since reptiles are sensitive indicators of palaeogeographical and palaeoclimatic events, we investigated the intraspecific genealogy of T. vermicularis in a phylogeographical framework. A total of 130 specimens were analyzed, while the use of formalin and ethanol as preservatives called for a special treatment of the samples. Partial sequences of two mitochondrial (12S and ND2) and one nuclear (PRLR) marker were targeted and the results of the phylogenetic analyses (NJ, ML and BI) and the parsimony-network revealed the existence of 10 evolutionary significant units within this species. In combination with the results of the dispersal-vicariance analysis, we may conclude that the Eurasian blindsnake has encountered a sequence of extinction events, followed by secondary expansion from refugia. Estimation of divergence times showed that severe climatic changes between significantly wetter and drier conditions in the Late Neogene have played a key role on the evolutionary and biogeographical history of T. vermicularis. Additionally, both markers (mtDNA and nDNA) distinguished a largely-differentiated evolutionary lineage (Jordan and south Syria), which could even be reckoned as a full species. Our study reveals the existence of cryptic evolutionary lineages within T. vermicularis, which calls for further attention both on the protection of intraspecific varieties and the respective geographic areas that hold them.

Phylogeography of the ocellated skink Chalcides ocellatus (Squamata, Scincidae), with the use of mtDNA sequences: a hitch-hiker's guide to the Mediterranean.

We analyze geographic genetic variation in C. ocellatus to evaluate the influences of major climatic, paleogeographic and anthropogenic factors in its biogeographic history. Ninety four specimens from 61 populations were collected across all of its geographical range and analyzed based on partial mitochondrial sequences (cyt b, 12S, and ND1). Our results demonstrate that an ancestral form of C. ocellatus, which expanded in northwestern Africa at the end of Miocene, diverged in at least three separate evolutionary lineages approximately 4.57Ma: C. humilis spread south of the Sahara, while the other two (C. ocellatus sensu stricto) were restricted in the coastal North African region. The complicated history of the ocellated skink is a result of multiple vicariant phenomena followed by multiple active or passive dispersals. The Messinian salinity crisis and the re-flooding of the Mediterranean basin, the climatic transition from Middle to Upper Pliocene, and the hyperarid phase of the Sahara, affected the distribution and diversification of C. ocellatus, while in historical times it was introduced in the central Mediterranean islands and eastern Mediterranean region from Tunisia and Cyrenaica, respectively.

Mitochondrial DNA reveals the genealogical history of the snake-eyed lizards (Ophisops elegans and O. occidentalis) (Sauria: Lacertidae).

The snake-eyed lizards of the genus Ophisops (Lacertidae) have been through a series of taxonomical revisions, but still their phylogenetic relationships remain uncertain. In the present study we estimate the phylogeographic structure of O. elegans across its distributional range and we evaluate the relationships between O. elegans and the sympatric, in North Africa, species O. occidentalis, using partial mtDNA sequences (16S rRNA, COI, and cyt b). All phylogenetic analyses produced topologically identical trees where extant populations of O. elegans and O. occidentalis were found polyphyletic. Taking into account all the potential causes of polyphyly (introgressive hybridization, incomplete lineage sorting, and imperfect taxonomy) we suggest the inaccurate taxonomy as the most likely explanation for the observed pattern. Our results stress the need for re-evaluation of the current taxonomical status of these species and their subspecies. Furthermore, our biogeographic analyses and the estimated time of divergences suggest a late Miocene diversification within these species, where the present distribution of O. elegans and O. occidentalis was the result of several dispersal and vicariant events, which are associated with climatic oscillations (the late Miocene aridification of Asia and northern Africa) and paleogeographic barriers of late Miocene and Pliocene period.

Molecular phylogeny of three Mesalina (Reptilia: Lacertidae) species (M. guttulata, M. brevirostris and M. bahaeldini) from North Africa and the Middle East: another case of paraphyly?

Mesalina is a widespread lacertid genus occurring throughout the Saharo-Sindian region from North Africa to Pakistan. It has been through a series of taxonomic revisions, but the phylogenetic relationships among the species remain unclear. In this study we estimate the phylogeographic structure of M. guttulata across most of its distributional range and we evaluate the relationships between M. guttulata and the sympatric species M. brevirostris and M. bahaeldini using partial mitochondrial DNA (mtDNA) sequences (cyt b and 16S). M. guttulata and M. brevirostris represent species complexes, whereas M. bahaeldini considered before as M. guttulata is a recently described species with very restricted distribution. Here we present the first evidence that M. guttulata is a paraphyletic taxon with respect to M. bahaeldini, while M. brevirostris proves to be a polytypic species or even a species complex, confirming previous studies. Although mtDNA markers have several properties that make them suitable for phylogeographic studies, they are not free of difficulties. Phylogeographic inferences within and between closely related species can be mislead by introgression and retention of ancestral polymorphism (incomplete lineage sorting). However, the present distribution pattern, the estimated times of divergence and the significant variation in morphology within M. guttulata led us to accept that the paraphyletic pattern observed, is most likely due to inaccurate taxonomy. Our hypothesis is that what has hitherto been considered as intraspecific variation, actually reflects species-level variation. Furthermore, our biogeographic analyses and the estimated time of divergences suggest that the present distribution of M. guttulata was the result of several dispersal and vicariant events, which are associated with historical changes (climatic oscillations and paleogeographic barriers) of late Miocene and Pliocene period.

Molecular phylogeny of the Greek legless skink Ophiomorus punctatissimus (Squamata: Scincidae): the impact of the mid-Aegean trench in its phylogeography.

Sequence data derived from three mitochondrial markers (cytochrome b, 16S rRNA and 12S rRNA genes) were used to infer the evolutionary history of several insular and mainland populations of the Greek legless skink (Ophiomorus punctatissimus), covering most of its distributional range. All phylogenetic analyses produced topologically identical trees that revealed a well-resolved phylogeny. These trees support two O. punctatissimus clades, which are geographically separated (west and east of the mid-Aegean trench). The assumption of a clock-like evolution could not be rejected, and thus a local clock was calibrated for the O. punctatissimus lineages. The non-overlapping geographic distributions of the major clades suggest a spatial and temporal sequence of diversification that coincides with paleogeographic separations during the geological history of the Aegean region. It seems that O. punctatissimus is an old eastern Mediterranean species that has been differentiating in this region at least from middle Miocene. It is possible that the ancestral form of O. punctatissimus invaded the Aegean region from Anatolia before the complete formation of the mid-Aegean trench, when the Aegean was still a uniform landmass, while other vicariant events have led to its present distribution.

Mitochondrial phylogeography of Rana (Pelophylax) populations in the Eastern Mediterranean region.

Phylogenetic relationships in the western fraction of Rana (Pelophylax) have not been resolved till now, even though several works have been devoted to the subject. Here, we infer phylogenetic relationships among the species distributed in the area of the Eastern Mediterranean, comparing partial mitochondrial DNA sequences for the cytochrome b and 16S rRNA genes. The obtained molecular data clearly indicate that Western Palearctic water frogs underwent a basal radiation into at least 3 major lineages (the perezi, the lessonae, and the rindibunda/bedriagae lineages) advocating an upper Miocene speciation. Moreover, we consider that within the rindibunda/bedriagae lineage, R. (P.) ridibunda, R. (P.) epeirotica, R. (P.) cretensis, R. (P.) bedriagae, R. (P.) cerigensis and R. (P.) kurtmuelleri were differentiated from a common ancestor through a series of vicariant and dispersal events, during the last approximately 5Mya, even though the specific rank of some taxa may be questionable, such as R. cerigensis in respect with R. bedriagae and R. kurtmuelleri in respect with R. ridibunda.

Digestive performance in five Mediterranean lizard species: effects of temperature and insularity.

Temperature sensitivity of digestive processes has important ramifications for digestive performance in ectothermic vertebrates. We conducted a comparative analysis of temperature effects on digestive processes [gut passage times (GPTs) and apparent digestive efficiencies (ADEs)] in five lacertid lizards occurring in insular (Podarcis erhardii, P. gaigeae), and mainland (P. muralis, P. peloponnesiaca, Lacerta graeca) Mediterranean environments. GPTs were negatively correlated to temperature with mainland taxa having 10-20% longer GPTs than island taxa. In contrast to previous studies that estimate ADEs using bomb calorimetry, we compare ADEs by analyzing discrete efficiencies for lipids, sugars and proteins at three temperature regimes (20, 25, and 30 degrees C); each of these categories produces different results. ADEs for lipids and sugars showed a monotonic increase with temperature whereas ADEs for proteins decreased with temperature. Island taxa had consistently higher ADEs than their mainland counterparts for lipids and for proteins but not for sugars. They are characterized by superior energy acquisition abilities despite significantly shorter GPTs. Their increased digestive performance relative to the mainland species appears to allow them to maximize energy acquisition in unproductive island environments where food availability is spatially and seasonally clustered.

Phylogenetic relationships and biogeography of Podarcis species from the Balkan Peninsula, by bayesian and maximum likelihood analyses of mitochondrial DNA sequences.

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant nonavian reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data, the species of Podarcis falls into four main groups that have substantial geographic coherence (Western island group, southwestern group, Italian group, and Balkan Peninsula group). The Balkan Peninsula species are divided into two subgroups: the subgroup of P. taurica (P. taurica, P. milensis, P. gaigeae, and perhaps P. melisellensis), and the subgroup of P. erhardii (P. erhardii and P. peloponnesiaca). In the present study, the question of phylogenetic relationships among the species of Podarcis encountered in the Balkan Peninsula was addressed using partial mtDNA sequences for cytochrome b (cyt b) and 16S rRNA (16S). The data support the monophyly of Podarcis and suggest that there are three phylogenetic clades: the clade A (P. taurica, P. gaigeae, P. milensis, and P. melisellensis); the clade B (P. erhardii and P. peloponnesiaca), and the clade C (P. muralis and P. sicula). By examining intraspecific relationships it was found that extant populations of P. erhardii are paraphyletic. Furthermore, subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, suggesting that their status should be reconsidered. The distinct geographic distribution of the major clades of the phylogenetic tree and its topology suggest a spatial and temporal sequence of phylogenetic separations that coincide with some major paleogeographic separations during the geological history of the Aegean Sea. The results stress the need for a reconsideration of the evolutionary history of Balkan Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species and subspecies level.

Phylogeography of Balkan wall lizard (Podarcis taurica) and its relatives inferred from mitochondrial DNA sequences.

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data the species of Podarcis falls into four main groups that have substantial geographical conherence (western island group, southwestern group, Italian group and Balkan group). The Balkan species are divided in two subgroups: the subgroup of Podarcis taurica (P. taurica, P. milensis, P. gaigeae and perhaps P. melisellensis), and the subgroup of Podarcis erhardii (P. erhardii and P. peloponnesiaca). We addressed the question of phylogenetic relations among the species of the P. taurica subgroup encountered in Greece, as they can be inferred from partial mtDNA (cyt b and 16S) sequences. Our data support the monophyly of P. taurica subgroup and suggest that P. gaigeae, P. milensis and P. melisellensis form a clade, which thereinafter connects to P. taurica. Within the previous clade, P. gaigeae is more closely related to P. milensis than to P. melisellensis. However, the specimens of P. taurica were subdivided in two different groups. The first one includes the specimens from northeastern Greece, and the other group includes the specimens from the rest of continental Greece and Ionian islands. Because the molecular clock of the cyt b and 16 rRNA genes was not rejected in our model test, it is possible to estimate times of speciation events. Based on the splitting of the island of Crete from Peloponnisos [c. 5 million years ago (Ma)], the evolutionary rate for the cyt b is 1.55% per million years (Myr) and for the 16S rRNA is 0.46% per Myr. These results suggest that the evolutionary history of P. taurica in Greece is more complex than a single evolutionary invasion. The data analysed, stress the need for a reconsideration of the evolutionary history of Greek Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species level.

Phylogenetic relationships and evolutionary history of snake-eyed skink Ablepharus kitaibelii (Sauria: Scincidae).

Sequence data derived from two mitochondrial markers, 16S rRNA and cytochrome b genes, were used to infer the phylogenetic relationships of 38 populations of the snake-eyed skinks of the genus Ablepharus with emphasis on A. kitaibelii from Greece and Turkey. The partition-homogeneity tests indicated that the combined data set was homogeneous, and maximum-parsimony, maximum-likelihood, and Bayesian analyses produced topologically identical trees that revealed a well-resolved phylogeny. All species except A. kitaibelii form monophyletic units. The latter species appears paraphyletic with respect to A. budaki and A. chernovi with populations clustering into two distinct clades. A. chernovi and A. budaki, which have recently been raised to species status, were confirmed as genetically distinct forms. We used sequence divergence and paleogeographic history of the Aegean region to reconstruct a biogeographic evolutionary scenario for A. kitaibelii.

Molecular phylogeny and biogeography of the wall-lizard Podarcis erhardii (Squamata: Lacertidae).

Erhard's wall lizard, Podarcis erhardii (Sauria: Lacertidae), is highly diversified in Greece and especially in the southern Aegean region. Out of the 28 recognized subspecies, 27 are found in Greece from the North Sporades island-complex in the North Aegean (grossly south of the 39th parallel) to the island of Crete in the South. The species exhibits great morphological and ecological plasticity and inhabits many different habitats from rocky islets and sandy shores to mountaintops as high as 2000m. By examining intraspecific variability at a segment of the mitochondrial gene cytochrome b we have found that that extant populations of P. erhardii are paraphyletic. Furthermore, we have found that subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, so that their status should be reconsidered. The DNA based biogeographical and phylogenetic history of Podarcis in Southern Greece is congruent with available paleogeographic data of the region, which supports the view that DNA sequences may be a useful tool for the study of palaeogeography.

Cytochrome b sequences reveal Acomys minous (Rodentia, Muridae) paraphyly and answer the question about the ancestral karyotype of Acomys dimidiatus.

Sequences of the cytochrome b (cyt b) mitochondrial gene show that the spiny mouse Acomys from Crete, known as the endemic species A. minous, is composed of two distinct maternal lineages ("A" and "B"). Group "A" sequences cluster with A. nesiotes (Cyprus) and group "B" sequences cluster with A. cilicicus (Turkey), which is evidence of paraphyly of A. minous in regard to these two species. From cyt b sequences, the three taxa are very closely related to A. cahirinus (Egypt): the maximum divergence found among these sequences is 1.6%, which is equivalent to the intraspecific diversity observed in other Acomys species. Paleozoology evidenced that man unintentionally introduced Acomys into Crete and Cyprus during antiquity. The divergence time between Acomys cyt b sequences found in Crete was estimated at 0.4 Myr, which means that the diversity observed did not appear after the introduction but reflects a much more ancient polymorphism. Cytochrome b phylogeny and cytogenetic data together comprise evidence that, within the species A. dimidiatus (Saudi Arabia, Israel, Egypt), it is the karyotypic form with 36 chromosomes that derives from the form with 38 chromosomes, due to a single acrocentric fusion.