Phylogeography of the Middle Eastern tree frogs (Hyla, Hylidae, Amphibia) as inferred from nuclear and mitochondrial DNA variation, with a description of a new species

Comparative phylogeography of four lizard taxa within an oceanic island

Comparative phylogeography aims to detect common patterns of differentiation among taxa attributed to the same geological or environmental barriers and, thus, find shared drivers of genomic splits and landscape features that explain patterns of occurrence and genetic diversity. Following this approach, this work focuses on four endemic lizard taxa, namely Acanthodactylus schreiberi schreiberi, Ophisops elegans schlueteri, Phoenicolacerta troodica, and Laudakia cypriaca, that co-occur in Cyprus, an isolated island of the Mediterranean. In an attempt to discern possible effects of the geological and ecological history of the region on diversity patterns along the lineages under study, samples were collected from all parts of Cyprus, and phylogenetic reconstructions, species delimitation analyses, and biogeographic comparisons were conducted at both genetic and genomic levels (Sanger and ddRAD sequencing data, respectively) for all four taxa. The derived phylogenomic relationships provide evidence of a strong genetic structure within the island that agrees with the geographical origins of each population. A rare case of an island-to-mainland dispersal event is confirmed for one of the focal taxa, making it one of the few documented cases of this type. Estimated diversification times during the Miocene reflect the timing of the first arrival on the island, while further diversification seems to have been strongly affected by the Pleistocene oscillations.

The phylogeography of Middle Eastern tree frogs in Israel

Western Palearctic treefrogs of the genus Hyla provide an example of a morphologically and ecologically cryptic group. Up to three distinct Hyla species have been proposed as resident in Israel and this number has consistently been subject to taxonomical debates. Here, we analyzed 16S rRNA and COI gene fragments of 658 individuals sampled at 47 pools in nine regions across Israel and the West Bank in order to resolve the taxonomic status of Hyla frogs. We generated both Bayesian and Maximum Likelihood phylogenies, and constructed time-calibrated trees to provide an evolutionary and historical context of sequence variations. We further applied SAMOVA as well as Monmonier's maximum-difference algorithm to study the genetic structure among populations and to identify potential zones acting as barriers to gene flow across locations. Our results revealed two distinct haplogroups for each gene fragment, with 95% CI divergence times dated from 8.9-17.1 Mya (16S) and 7.1-23.6 Mya (COI), respectively. SAMOVA and barrier analyses partitioned the populations into three groups. Our results highlight that, while there are probably only two Hyla species in Israel, one population of one of the species might qualify as a separate evolutionarily significant unit. Our findings elucidate the taxonomic status of Hyla frogs in Israel and provide the basis for determining appropriate management and conservation priorities.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Mitogenomics provides new insights into the phylogenetic relationships and evolutionary history of deep-sea sea stars (Asteroidea)

The deep sea (> 200 m) is considered as the largest and most remote biome, which characterized by low temperatures, low oxygen level, scarce food, constant darkness, and high hydrostatic pressure. The sea stars (class Asteroidea) are ecologically important and diverse echinoderms in all of the world’s oceans, occurring from the intertidal to the abyssal zone (to about 6000 m). To date, the phylogeny of the sea stars and the relationships of deep-sea and shallow water groups have not yet been fully resolved. Here, we recovered five mitochondrial genomes of deep-sea asteroids. The A+T content of the mtDNA in deep-sea asteroids were significantly higher than that of the shallow-water groups. The gene orders of the five new mitogenomes were identical to that of other asteroids. The phylogenetic analysis showed that the orders Valvatida, Paxillosida, Forcipulatida are paraphyletic. Velatida was the sister order of all the others and then the cladeValvatida-Spinulosida-Paxillosida-Notomyotida versus Forcipulatida-Brisingida. Deep-sea asteroids were nested in different lineages, instead of a well-supported clade. The tropical Western Pacific was suggested as the original area of asteroids, and the temperate water was initially colonized with asteroids by the migration events from the tropical and cold water. The time-calibrated phylogeny showed that Asteroidea originated during Devonian-Carboniferous boundary and the major lineages of Asteroidea originated during Permian–Triassic boundary. The divergence between the deep-sea and shallow-water asteroids coincided approximately with the Triassic-Jurassic extinction. Total 29 positively selected sites were detected in fifteen mitochondrial genes of five deep-sea lineages, implying a link between deep-sea adaption and mitochondrial molecular biology in asteroids.

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.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Morphological diversification of Mediterranean anurans: the roles of evolutionary history and climate

Investigation of the ecological and evolutionary mechanisms governing the origin and diversification of species requires integrative approaches that often have to accommodate strong discordance among datasets. A common source of conflict is the combination of morphological and molecular characters with different evolutionary rates. Resolution of these discordances is crucial to assess the relative roles of different processes in generating and maintaining biodiversity. Anuran amphibians provide many examples of morphologically similar, genetically divergent lineages, posing questions about the relative roles of phylogeny and ecological factors in phenotypic evolution. We focused on three circum-Mediterranean anuran genera (Hyla, Alytes and Discoglossus), characterizing morphological and environmental disparity and comparing diversity patterns across biological levels of organization. Using a comparative phylogenetic framework, we tested how shared ancestry and climatic factors come together to shape phenotypic diversity. We found higher morphological differentiation within Hyla and Alytes than in Discoglossus. Body size and limb morphology contributed most to inter- and intraspecific morphological variation in Hyla and Alytes, but there was no strong phylogenetic signal, indicating that shared ancestry does not predict patterns of phenotypic divergence. In contrast, we uncovered a significant association between morphology and climatic descriptors, supporting the hypothesis that morphological disparity between species results from adaptive evolution.

Pleistocene climate fluctuations as the major driver of genetic diversity and distribution patterns of the Caspian green lizard, Lacerta strigata Eichwald, 1831

Green lizards of the genus Lacerta have served as excellent models for studying the impact of Pleistocene climatic oscillations on genetic structures. The Caspian green lizard, Lacerta strigata, occupies various habitats across the Caucasus and the South Caspian Sea, with the Hyrcanian Forests and north of the Alborz Mountains forming the core of the range. This study aimed to re-examine the phylogenetic relationships of L. strigata with other congeneric members and to assess the genetic structure and historical demography of the species. Furthermore, Species Distribution Models (SDMs) were performed to infer the species' potential habitat suitability and were then projected on climate scenarios reflecting current and past (6 ky and 21 ky before present) conditions. A total of 39 individuals collected from most of the distribution range, together with additional lacertid species sequences from the GenBank database, were examined using mtDNA (Cyt b and 12S ribosomal RNA) and nuclear (C-mos and β-fibrinogen) sequence data. Based on the phylogenetic analyses, L. strigata was found to be a sister taxon to all other members of the genus. The species included two main clades (regional western and eastern) that diverged in a period between the Early and Middle Pleistocene. Based on the BBM and S-Diva analyses, both dispersal and vicariance events explained the phylogeographic structure of the species in the Hyrcanian Forests. The historical demographic analyses using Bayesian skyline plots showed a mild increase in the effective population size from about 120 Kya for the western regional clade. According to phylogeographic structures and SDMs evidence, as in other species within the region, it appears that the south of the Caspian Sea (Hyrcanian Forests), and the Alborz Mountains acted as multiple refugia during cold periods and promoted expansion outwards amid the warm periods. Overall, the results provided evidence that the genetic structure of the species has been influenced by the Pleistocene climatic fluctuations.

Phylogenetic Relationships and Adaptation in Deep-Sea Mussels: Insights from Mitochondrial Genomes

Mitochondrial genomes (mitogenomes) are an excellent source of information for phylogenetic and evolutionary studies, but their application in marine invertebrates is limited. In the present study, we utilized mitogenomes to elucidate the phylogeny and environmental adaptation in deep-sea mussels (Mytilidae: Bathymodiolinae). We sequenced and assembled seven bathymodioline mitogenomes. A phylogenetic analysis integrating the seven newly assembled and six previously reported bathymodioline mitogenomes revealed that these bathymodiolines are divided into three well-supported clades represented by five Gigantidas species, six Bathymodiolus species, and two "Bathymodiolus" species, respectively. A Common interval Rearrangement Explorer (CREx) analysis revealed a gene order rearrangement in bathymodiolines that is distinct from that in other shallow-water mytilids. The CREx analysis also suggested that reversal, transposition, and tandem duplications with subsequent random gene loss (TDRL) may have been responsible for the evolution of mitochondrial gene orders in bathymodiolines. Moreover, a comparison of the mitogenomes of shallow-water and deep-sea mussels revealed that the latter lineage has experienced relaxed purifying selection, but 16 residues of the atp6, nad4, nad2, cob, nad5, and cox2 genes have underwent positive selection. Overall, this study provides new insights into the phylogenetic relationships and mitogenomic adaptations of deep-sea mussels.

Speciation patterns in the Forficula auricularia species complex: cryptic and not so cryptic taxa across the western Palaearctic region

Forficula auricularia (the European earwig) is possibly a complex of cryptic species. To test this hypothesis, we performed: (1) a phylogeographic study based on fragments of the mitochondrial COI and the nuclear ITS2 markers on a wide geographic sampling, (2) morphometric analyses of lineages present in Spain and (3) niche overlap analyses. We recovered five reciprocally monophyletic ancient phylogroups with unique historical patterns of distribution, climatic niches and diversification. External morphology was conserved and not correlated with speciation events, except in one case. Phylogenetic placement of the morphologically distinct taxon renders F. auricularia paraphyletic. Based on the congruence of the phylogenetic units defined by mtDNA and nuclear sequence data, we conclude that phylogroups have their own historical and future evolutionary trajectory and represent independent taxonomic units. Forficula auricularia is a complex of at least four species: the morphologically diagnosable Forficula aeolica González-Miguéns & García-París sp. nov., and the cryptic taxa: Forficula mediterranea González-Miguéns & García-París sp. nov., Forficula dentataFabricius, 1775stat. nov. and Forficula auriculariaLinnaeus, 1758s.s. We also provide new synonymy for F. dentata.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

A river runs through it: tree frog genomics supports the Dead Sea Rift as a rare phylogeographical break

Phylogeographical breaks can be viewed as regional hotspots of diversity where the genetic integrity of incipient species is put to the test. We focus on an understudied species transition from the Middle East, namely the Dead Sea Rift in the Levant region, which presumably divided the tree frogs Hyla savignyi and H. felixarabica. Combining multilocus genetic analyses (mitochondrial DNA and RAD-sequencing) with ecological niche modelling, we test whether the rift effectively acts as a biogeographical barrier preventing this pair from admixing and merging. The answer is yes: despite weak signs of introgression, all parapatric populations were assigned to either species without cyto-nuclear discordance. Yet, the projected distributions under present and glacial conditions largely overlapped in the area, meaning their current parapatric ranges do not represent an ecological transition. Instead, we hypothesize that H. savignyi and H. felixarabica are maintained apart by limited opportunities for dispersal across the barren Jordan Valley, combined with advanced reproductive isolation. Therefore, the Dead Sea Rift may represent a rare phylogeographical break, and we encourage international efforts to assess its contribution to the rich biodiversity of the Middle East.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

Complex hybridization patterns in European pond turtles (Emys orbicularis) in the Pyrenean Region

Hybrid zones are natural laboratories allowing insights in genetic processes like lineage diversification, speciation and introgression. Using large sampling, 15 microsatellite loci and a mitochondrial marker, we examined the Pyrenean contact zone of three pond turtle taxa (Emys orbicularis orbicularis, E. o. galloitalica, E. o. occidentalis). The Pyrenees are a biogeographically important region separating many lineages endemic to the Iberian Peninsula from their Western European counterparts. We discovered limited admixture, reflecting a complex biogeographic scenario. Simulations using Approximate Bayesian Computing supported that E. o. orbicularis invaded the Iberian Peninsula in the Holocene, circumventing the Pyrenees along the Mediterranean coast, and hybridized in the northern peninsula with the local coastal subspecies galloitalica, and to a lesser extent, with occidentalis. While E. o. occidentalis, and in particular E. o. orbicularis, expanded their ranges considerably during Holocene warming, E. o. galloitalica remained largely confined to its former Iberian refuge. Admixture among the three taxa is surprisingly low, and a future taxonomic investigation that includes the unstudied subspecies of E. orbicularis from North Africa, Eastern Europe and Asia has to determine whether their current status properly reflects their evolutionary divergence or whether certain taxa should be regarded as full species.

Historical explanation of genetic variation in the Mediterranean horseshoe bat Rhinolophus euryale (Chiroptera: Rhinolophidae) inferred from mitochondrial cytochrome-b and D-loop genes in Iran

Molecular phylogeography and species distribution modelling (SDM) suggest that late Quaternary glacial cycles have portrayed a significant role in structuring current population genetic structure and diversity. Based on phylogenetic relationships using Bayesian inference and maximum likelihood of 535 bp mtDNA (D-loop) and 745 bp mtDNA (Cytb) in 62 individuals of the Mediterranean Horseshoe Bat, Rhinolophus euryale, from 13 different localities in Iran we identified two subspecific populations with differing population genetic structure distributed in southern Zagros Mts. and northern Elburz Mts. Analysis of molecular variance (AMOVA) obtained from D-loop sequences indicates that 21.18% of sequence variation is distributed among populations and 10.84% within them. Moreover, a degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θCT = 0.68, p = .005). The positive and significant correlation between geographic and genetic distances (R² = 0.28, r = 0.529, p = .000) is obtained following controlling for environmental distance. Spatial distribution of haplotypes indicates that marginal population of the species in southern part of the species range have occupied this section as a glacial refugia. However, this genetic variation, in conjunction with results of the SDM shows a massive postglacial range expansion for R. euryale towards higher latitudes in Iran.

The complete mitochondrial genome of the alvinocaridid shrimp Shinkaicaris leurokolos (Decapoda, Caridea): Insight into the mitochondrial genetic basis of deep-sea hydrothermal vent adaptation in the shrimp

Deep-sea hydrothermal vent is one of the most extreme environments on Earth with low oxygen and high levels of toxins. Decapod species from the family Alvinocarididae have colonized and successfully adapted to this extremely harsh environment. Mitochondria plays a vital role in oxygen usage and energy metabolism, thus it may be under selection in the adaptive evolution of the hydrothermal vent shrimps. In this study, the mitochondrial genome (mitogenome) of alvinocaridid shrimp Shinkaicaris leurokolos (Kikuchi & Hashimoto, 2000) was determined through Illumina sequencing. The mitogenome of S. leurokolos was 15,903bp in length, containing 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of sequenced alvinocaridids. It has the longest concatenated sequences of protein-coding genes, tRNAs and shortest pooled rRNAs among the alvinocaridids. The control regions (CRs) of alvinocaridid were significantly longer (p<0.01) than those of the other caridaen. The alignment of the alvinocaridid CRs revealed two conserved sequence blocks (CSBs), and each of the CSBs included a noncanonical open reading frame (ORF), which may be involved in adjusting mitochondrial energy metabolism to adapt to the hydrothermal environment. Phylogenetic analysis supported that the deep-sea hydrothermal vent shrimps may have originated from those living in shallow area. Positive selection analysis reveals the evidence of adaptive change in the mitogenome of Alvinocarididae. Thirty potentially important adaptive residues were identified, which were located in atp6, cox1, cox3, cytb and nad1-5. This study explores the mitochondrial genetic basis of hydrothermal vent adaptation in alvinocaridid for the first time, and provides valuable clues regarding the adaptation.

Pre-Quaternary divergence and subsequent radiation explain longitudinal patterns of genetic and morphological variation in the striped skink, Heremites vittatus

Background

Many animal and plant species in the Middle East and northern Africa have a predominantly longitudinal distribution, extending from Iran and Turkey along the eastern Mediterranean coast into northern Africa. These species are potentially characterized by longitudinal patterns of biological diversity, but little is known about the underlying biogeographic mechanisms and evolutionary timescales. We examined these questions in the striped skink, Heremites vittatus, one such species with a roughly longitudinal distribution across the Middle East and northern Africa, by analyzing range-wide patterns of mitochondrial DNA (mtDNA) sequence and multi-trait morphological variation.

Results

The striped skink exhibits a basic longitudinal organization of mtDNA diversity, with three major mitochondrial lineages inhabiting northern Africa, the eastern Mediterranean coast, and Turkey/Iran. Remarkably, these lineages are of pre-Quaternary origin, and are characterized by p-distances of 9–10%. In addition, within each of these lineages a more recent Quaternary genetic diversification was observed, as evidenced by deep subclades and high haplotype diversity especially in the Turkish/Iranian and eastern Mediterranean lineages. Consistent with the genetic variation, our morphological analysis revealed that the majority of morphological traits show significant mean differences between specimens from northern Africa, the eastern Mediterranean coast, and Turkey/Iran, suggesting lineage-specific trait evolution. In addition, a subset of traits exhibits clinal variation along the eastern Mediterranean coast, potentially indicating selection gradients at the geographic transition from northern Africa to Anatolia. The existence of allopatric, morphologically and genetically divergent lineages suggests that Heremites vittatus might represent a complex with several taxa.

Conclusions

Our work demonstrates that early divergence events in the Pliocene, likely driven by both climatic and geological factors, established the longitudinal patterns and distribution of Heremites vittatus. Subsequent radiation during the Pleistocene generated the genetic and morphological diversity observed today. Our study provides further evidence that longitudinal diversity patterns and species distributions in the Middle East and northern Africa were shaped by complex evolutionary processes, involving the region’s intricate geological history, climatic oscillations, and the presence of the Sahara.

Electronic supplementary material

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

The mitochondrial genome of a sea anemone Bolocera sp. exhibits novel genetic structures potentially involved in adaptation to the deep-sea environment

The deep sea is one of the most extensive ecosystems on earth. Organisms living there survive in an extremely harsh environment, and their mitochondrial energy metabolism might be a result of evolution. As one of the most important organelles, mitochondria generate energy through energy metabolism and play an important role in almost all biological activities. In this study, the mitogenome of a deep‐sea sea anemone (Bolocera sp.) was sequenced and characterized. Like other metazoans, it contained 13 energy pathway protein‐coding genes and two ribosomal RNAs. However, it also exhibited some unique features: just two transfer RNA genes, two group I introns, two transposon‐like noncanonical open reading frames (ORFs), and a control region‐like (CR‐like) element. All of the mitochondrial genes were coded by the same strand (the H‐strand). The genetic order and orientation were identical to those of most sequenced actiniarians. Phylogenetic analyses showed that this species was closely related to Bolocera tuediae. Positive selection analysis showed that three residues (31 L and 42 N in ATP6, 570 S in ND5) of Bolocera sp. were positively selected sites. By comparing these features with those of shallow sea anemone species, we deduced that these novel gene features may influence the activity of mitochondrial genes. This study may provide some clues regarding the adaptation of Bolocera sp. to the deep‐sea environment.

Neogene amphibians and reptiles (Caudata, Anura, Gekkota, Lacertilia, and Testudines) from the south of Western Siberia, Russia, and Northeastern Kazakhstan

Background

The present-day amphibian and reptile fauna of Western Siberia are the least diverse of the Palaearctic Realm, as a consequence of the unfavourable climatic conditions that predominate in this region. The origin and emergence of these herpetofaunal groups are poorly understood. Aside from the better-explored European Neogene localities yielding amphibian and reptile fossil remains, the Neogene herpetofauna of Western Asia is understudied. The few available data need critical reviews and new interpretations, taking into account the more recent records of the European herpetofauna. The comparison of this previous data with that of European fossil records would provide data on palaeobiogeographic affiliations of the region as well as on the origin and emergence of the present-day fauna of Western Siberia. An overview of the earliest occurrences of certain amphibian lineages is still needed. In addition, studies that address such knowledge gaps can be useful for molecular biologists in their calibration of molecular clocks.

Methods and Results

In this study, we considered critically reviewed available data from amphibian and reptile fauna from over 40 Western Siberian, Russian and Northeastern Kazakhstan localities, ranging from the Middle Miocene to Early Pleistocene. Herein, we provided new interpretations that arose from our assessment of the previously published and new data. More than 50 amphibians and reptile taxa were identified belonging to families Hynobiidae, Cryptobranchidae, Salamandridae, Palaeobatrachidae, Bombinatoridae, Pelobatidae, Hylidae, Bufonidae, Ranidae, Gekkonidae, Lacertidae, and Emydidae. Palaeobiogeographic analyses were performed for these groups and palaeoprecipitation values were estimated for 12 localities, using the bioclimatic analysis of herpetofaunal assemblages.

Conclusion

The Neogene assemblage of Western Siberia was found to be dominated by groups of European affinities, such as Palaeobatrachidae, Bombina, Hyla, Bufo bufo, and a small part of this assemblage included Eastern Palaearctic taxa (e.g. Salamandrella, Tylototriton, Bufotes viridis). For several taxa (e.g. Mioproteus, Hyla, Bombina, Rana temporaria), the Western Siberian occurrences represented their most eastern Eurasian records. The most diverse collection of fossil remains was found in the Middle Miocene. Less diversity has been registered towards the Early Pleistocene, potentially due to the progressive cooling of the climate in the Northern Hemisphere. The results of our study showed higher-amplitude changes of precipitation development in Western Siberia from the Early Miocene to the Pliocene, than previously assumed.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.

Genetic variations of Turkish bank vole, Myodes glareolus (Mammalia: Rodentia) inferred from mtDNA

The bank vole, Myodes glareolus, lives in deciduous forests throughout the Palearctic region. In Turkey, this species is distributed only in northern Anatolia (the Black Sea region) where these forests exist. This study reveals genetic differentiation among bank vole populations based on two regions of mitochondrial DNA (cytochrome b and D-loop). Populations in northern Anatolia are divided into two genetic lineages (the "eastern" and "western Black Sea" lineages) by the Kızılırmak Valley. While the western Black Sea lineage is close to the Balkan lineage, in accordance with their geographical proximities, surprisingly, the Uludag lineage, also situated in Western Turkey appears related to the eastern Black Sea population. The divergence time analyses suggest a separation between the Balkan and Turkish groups around 0.26 Mya, whereas the split between the eastern and western Black sea lineages appeared a little bit later (0.20 Mya). Our results suggest that regional refuges existed for this species in Turkey and that small-scale habitat fragmentations led to genetic differentiations between Myodes populations.

Historical biogeography resolves the origins of endemic Arabian toad lineages (Anura: Bufonidae): Evidence for ancient vicariance and dispersal events with the Horn of Africa and South Asia

Background

The Arabian Peninsula is home to a unique fauna that has assembled and evolved throughout the course of major geophysical events, including the separation of the Arabian Plate from Africa and subsequent collision with Eurasia. Opportunities for faunal exchanges with particular continents occurred in temporally distinct periods, and the presence of African, Western Eurasian, and South Asian derived taxa on the Arabian Peninsula signifies the complexity of these historical biogeographic events. The six true toad species (family Bufonidae) endemic to the Arabian Peninsula present a considerable taxonomic and biogeographic challenge because they are part of a global bufonid radiation, including several genera surrounding the Arabian Peninsula, and difficult to discriminate morphologically. As they could be derived from African, Western Eurasian, or South Asian toad groups, elucidating their evolutionary relationships has important implications for historical biogeography. Here, we analyze a global molecular data set of 243 bufonid lineages, with an emphasis on new sampling from the Horn of Africa, Western Eurasia, South Asia, and the Arabian Peninsula, to reconstruct the evolutionary relationships of the Arabian species. We produce a robust time-calibrated phylogeny to infer the biogeographic history of this group on and around the Arabian Peninsula.

Results

Our phylogenetic analyses indicate two of the endemic Arabian toad species, “Bufo” tihamicus and “Bufo” arabicus, evolved independently within the African genus Amietophrynus. We confirm the Arabian species Duttaphrynus dhufarensis is of South Asian origin, but do not find evidence for the Asian genus Duttaphrynus being present in the Horn of Africa, discrediting a previously proposed Asian bufonid dispersal event to Africa. We also do not find evidence of the African genus Amietophrynus occurring in South Asia, suggesting that unlike many other vertebrate taxa, toads have not used the Arabian Peninsula as a stepping-stone for trans-continental dispersal. Our divergence dating estimates strongly suggest the formation of the Red Sea drove simultaneous divergences between two of the Arabian species (A. tihamicus comb. nov. and A. arabicus comb. nov.) and their closest mainland African relatives in the Early Miocene. We estimate the divergence of D. dhufarensis with its closest South Asian relatives occurred in the mid to Late Miocene, suggesting the temporary or permanent land connections between the Arabian plate and Eurasia facilitated dispersal of this lineage to the Arabian Peninsula.

Conclusions

The Arabian bufonid assemblage, despite being comparatively depauperate with respect to surrounding continents, exemplifies the faunal pattern of the Arabian Peninsula, namely being a complex admixture of African, Western Eurasian, and South Asian elements. The historical biogeographic patterns exhibited by Arabian toads and their allies are concordant with studies of other vertebrate taxa, building support for the role of major geological events in driving simultaneous vicariance and dispersal events around the Arabian Peninsula. Although many taxa or groups exhibiting disjunct Afro-Arabian distributions appear to have dispersed more recently from the Horn of Africa via a southern land bridge or overwater dispersal, both Amietophrynus tihamicus and A. arabicus likely represent true African relicts resulting from vicariance associated with the Red Sea formation, a pattern that so far is rare among the vertebrate species investigated.

Electronic supplementary material

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

First assessment on the molecular phylogeny of Anatololacerta (Squamata, Lacertidae) distributed in Southern Anatolia: insights from mtDNA and nDNA markers

The genus Anatololacerta (Lacertidae) occurs mainly in Anatolia (western and southern Turkey) and on the Aegean islands Samos, Ikaria, and Rhodos. Although its taxonomy has long been debated and is currently nascent, three morphological species have been attributed to this genus: Anatololacerta anatolica, Anatololacerta oertzeni, and Anatololacerta danfordi. Here, we investigated the evolutionary history of A. oertzeni and Anatololacerta danfordi based on both mitochondrial and nuclear markers (16S rRNA and cmos). In total, 34 Anatololacerta specimens were analyzed using maximum likelihood (ML) and Bayesian inference (BI) methods. Our results supported the presence of four well-supported lineages: two belongs to A. oertzeni and two to A. danfordi. The temporal diversification of these lineages probably started with the divergence of the first A. oertzeni lineage from western Antalya at 7.9 Mya. The other two major splits may have occurred in early Pliocene (4.4 Mya: the divergence of the second A. oertzeni from A. danfordi) and in late Pliocene (2.7 Mya: the divergence of the two lineages of A. danfordi). The phylogeographical scenario suggests that the major diversification events (from late Miocene to late Pliocene) could be related with climatic oscillations (such as the late Miocene aridification and the Messinian Salinity Crisis) and tectonic movements (such as the uplift of the central Taurus mountain).

Cryptic diversity within the Anatololacerta species complex (Squamata: Lacertidae) in the Anatolian Peninsula: evidence from a multi-locus approach

The rapid development of innovative molecular tools for characterizing biodiversity is leading to an extensive and sometimes unexpected renovation of taxonomic classifications. Particularly, for species having allopatric or parapatric distributions or resulting from recent speciation processes, the absence of clear phenotypic differentiation may hinder the recognition of closely related taxa, while intraspecific polymorphism may be confused with the presence of more than one single species. In the present work, we apply different phylogenetic methods in order to infer relationships within the genus Anatololacerta, and to assess the taxonomy of this morphologically diversified group of lizards endemic to western and southern Anatolia and some neighboring Aegean islands. According to morphology, three species have been recognized (Anatololacerta anatolica, A. oertzeni and A. danfordi) as well as several subspecies, but small variation at immunological markers led some authors to join all the populations into one single taxon, A. danfordi. By selecting both mitochondrial and nuclear informative markers, we tested the effectiveness of classical "gene tree" (i.e. Bayesian Inference) vs. innovative (i.e. coalescent-based) "species tree" methods in resolving the Anatololacerta taxonomic enigma, as a case in point for similar studies on species complexes resulting from non-obvious and cryptic diversification patterns. According to our results, the gene tree method failed in resolving phylogenetic relationships among clades, whereas the multi-locus species tree approach, coupled with species delimitation methods, allowed the identification of four well distinct species. These species probably diversified in different allopatric refugia located in southern and western Anatolia, where isolated populations may have persisted during Pleistocene glacial cycles.

Evolutionary origins and genetic variation of the Seychelles treefrog, Tachycnemis seychellensis (Duméril and Bibron, 1841) (Amphibia: Anura: Hyperoliidae)

The hyperoliid frog Tachycnemis seychellensis, the only species of its genus, is endemic to the four largest granitic islands of the Seychelles archipelago and is reliant on freshwater bodies for reproduction. Its presence in the Seychelles is thought to be the product of a transoceanic dispersal, diverging from the genus Heterixalus, its closest living relative (currently endemic to Madagascar), between approximately 10-35Ma. A previous study documented substantial intraspecific morphological variation among island populations and also among populations within the largest island (Mahé). To assess intraspecific genetic variation and to infer the closest living relative(s) of T. seychellensis, DNA sequence data were generated for three mitochondrial and four nuclear markers. These data support a sister-group relationship between T. seychellensis and Heterixalus, with the divergence between the two occurring between approximately 11-19Ma based on cytb p-distances. Low levels of genetic variation were found among major mitochondrial haplotype clades of T. seychellensis (maximum 0.7% p-distance concatenated mtDNA), and samples from each of the islands (except La Digue) comprised multiple mitochondrial haplotype clades. Two nuclear genes (rag1 and tyr) showed no variation, and the other two (rho and pomc) lacked any notable geographic structuring, counter to patterns observed within presumably more vagile Seychelles taxa such as lizards. The low levels of genetic variation and phylogeographic structure support an interpretation that there is a single but morphologically highly variable species of Seychelles treefrog. The contrasting genetic and morphological intraspecific variation may be attributable to relatively recent admixture during low sea-level stands, ecophenotypic plasticity, local adaptation to different environmental conditions, and/or current and previously small population sizes. Low genetic phylogeographic structure but substantial morphological variation is unusual within anurans.

Two newly recognized species of Hemidactylus (Squamata, Gekkonidae) from the Arabian Peninsula and Sinai, Egypt

A recent molecular phylogeny of the Arid clade of the genus Hemidactylus revealed that the recently described H. saba and two unnamed Hemidactylus species from Sinai, Saudi Arabia and Yemen form a well-supported monophyletic group within the Arabian radiation of the genus. The name ‘Hemidactylus saba species group’ is suggested for this clade. According to the results of morphological comparisons and the molecular analyses using two mitochondrial (12S and cytb) and four nuclear (cmos, mc1r, rag1, rag2) genes, the name Hemidactylus granosus Heyden, 1827 is resurrected from the synonymy of H. turcicus for the Sinai and Saudi Arabian species. The third species of this group from Yemen is described formally as a new species H. uliisp. n. The phylogenetic relationships of the members of ‘Hemidactylus saba species group’ are evaluated and the distribution and ecology of individual species are discussed.

Phylogeographic analyses strongly suggest cryptic speciation in the giant spiny frog (Dicroglossidae: Paa spinosa) and interspecies hybridization in Paa

Species identification is one of the most basic yet crucial issues in biology with potentially far-reaching implications for fields such as conservation, population ecology, and epidemiology. The widely distributed but threatened frog Paa spinosa has been speculated to represent a complex of multiple species. In this study, 254 individuals representing species of the genus Paa were investigated along the entire range of P. spinosa: 196 specimens of P. spinosa, 8 specimens of P. jiulongensis, 5 specimens of P. boulengeri, 20 specimens of P. exilispinosa, and 25 specimens of P. shini. Approximately 1333 bp of mtDNA sequence data (genes 12S rRNA and 16S rRNA) were used. Phylogenetic analyses were conducted using maximum parsimony, maximum likelihood and Bayesian inference. BEAST was used to estimate divergence dates of major clades. Results suggest that P. spinosa can be divided into three distinct major lineages. Each major lineage totally corresponds to geographical regions, revealing the presence of three candidate cryptic species. Isolation and differentiation among lineages are further supported by the great genetic distances between the lineages. The bifurcating phylogenetic pattern also suggests an east-west dispersal trend during historic cryptic speciation. Dating analysis estimates that P. spinosa from Western China split from the remaining P. spinosa populations in the Miocene and that P. spinosa from Eastern China diverged from Central China in the Pliocene. We also found that P. exilispinosa from Mainland China and Hong Kong might have a complex of multiple species. After identifying cryptic lineages, we then determine the discrepancy between the mtDNA and the morphotypes in several individuals. This discrepancy may have been caused by introgressive hybridization between P. spinosa and P. shini.

Out of Arabia: a complex biogeographic history of multiple vicariance and dispersal events in the gecko genus Hemidactylus (Reptilia: Gekkonidae)

The geological history of the Arabian Peninsula has played a crucial role in shaping current diversity and distribution patterns of many Arabian and African faunal elements. The gecko genus Hemidactylus is not an exception. In this study, we provide an insight into the phylogeny and systematics of 45 recognized species of the so-called Arid clade of the genus Hemidactylus from Arabia, the Horn of Africa, the Levant and Iran. The material comprises 358 specimens sequenced for up to two mitochondrial (12S rRNA, cytochrome b) and four nuclear (mc1r, cmos, rag1, rag2) genes with 4766 bp of the concatenated alignment length. A robust calibrated phylogeny and reconstruction of historical biogeography are inferred. We link the history of this genus with major geological events that occurred in the region within the last 30 million years. Two basal divergences correspond with the break-ups of the Arabian and African landmasses and subsequent separation of Socotra from the Arabian mainland, respectively, segregating the genus by means of vicariance. Formation of the Red Sea led to isolation and subsequent radiation in the Arabian Peninsula, which was followed by multiple independent expansions: 13.1 Ma to Iran; 9.8 Ma to NE Africa; 8.2 to Socotra Archipelago; 7–7.3 Ma two colonizations to the Near East; 5.9 Ma to NE Africa; and 4.1 to Socotra. Moreover, using multiple genetic markers we detected cryptic diversity within the genus, particularly in south-western Arabia and the Ethiopian highlands, and confirmed the existence of at least seven new species in the area. These findings highlight the role of Arabia and the Horn of Africa as an important Hemidactylus diversity hotspot.

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA(His) differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.