Molecular phylogenetics and historical biogeography among salamandrids of the "true" salamander clade: rapid branching of numerous highly divergent lineages in Mertensiella luschani associated with the rise of Anatolia

A skin secretion metabolome analysis of the Greek Dodecanese Lycian salamanders: Preliminary evidence of dietary alkaloid sequestration in urodeles

Lyciasalamandra species, like most amphibians, secrete a wide array of compounds from their granular and mucous skin glands, including the internally synthesized samandarine alkaloids, making their skin a complex organ performing a variety of functions. Lyciasalamandra helverseni and L. luschani basoglui are insular endemics of the Dodecanese islands of SE Greece, bearing distinct isolated populations, with well-documented phylogenetic profiles. Here, we employ a metabolomics approach, utilizing UPLC-ESI-HRMS/MS data of the skin secretions sampled from a number of specimens found in the islands of Karpathos, Kasos and Kastellorizo, in an effort to reveal aspects of their chemistry and diversity across populations. The results indicated statistically significant variation between all taxa examined, based on various secreted compounds. The underlying factors of variation highlighted by the multivariate analysis were differences in samandarine and other alkaloid content as well as in animal size. Metabolite annotation, based on dereplication tools and most importantly HRMS and HRMS/MS spectra, yielded a number of known samandarine alkaloids, reported for the first time in the currently studied Lyciasalamandra species. We also present documentation for novel members of the samandarine alkaloid family, as well as preliminary evidence for a possible dietary alkaloid sequestration. This work can set the basis for further research of this often-neglected endemic species of the Salamandridae, as well as the structural investigation of the samandarine alkaloid group.

Phylogeography of a Typical Forest Heliothermic Lizard Reveals the Combined Influence of Rivers and Climate Dynamics on Diversification in Eastern Amazonia

The formation of the Amazon drainage basin has been considered an important driver of speciation of several taxa, promoting vicariant events or reinforcement of barriers that restrict gene flow between opposite river margins. Several recent studies reported a set of miscellaneous events involving climatic fluctuations, geomorphological changes, and dispersal mechanisms as propellers of diversification of Amazonian rainforest taxa. Here, we show the results of dated phylogenetic, biogeographic, and populational analyses to investigate which events could better explain the current distribution of a heliothermic, active foraging lizard in the central and eastern portions of the Amazonian rainforest (besides a disjunct distribution in part of the Atlantic Forest). We sampled Kentropyx calcarata from most of its area of occurrence in Amazonia and used mitochondrial and nuclear markers to evaluate if the genetic structure agrees with evolutionary scenarios previously proposed for Amazonia. We performed phylogenetic and populational analyses to better understand the dynamics of this species in the Amazonia rainforest over time. Phylogenetic inference recovered ten K. calcarata structured lineages in eastern Amazonia, some of them limited by the Amazon River and its southern tributaries (Tapajós, Xingu, and Tocantins), although we detected occasional haplotype sharing across some of the river banks. According to molecular dating, K. calcarata diversified since Miocene–Pliocene, and some of the lineages presented signs of demographic expansion during the Pleistocene, supposedly triggered by climatic dynamics. The putative ancestral lineage of K. calcarata was distributed on the Guiana Shield, later spreading south and southeastward by dispersion. Our results indicate that Amazonian rivers acted as barriers to the dispersal of Kentropyx calcarata, but they were not the sole drivers of diversification.

Differentiated historical demography and ecological niche forming present distribution and genetic structure in coexisting two salamanders (Amphibia, Urodela, Hynobiidae) in a small island, Japan

Background

The climatic oscillations in the Quaternary period considerably shaped the distribution and population genetic structure of organisms. Studies on the historical dynamics of distribution and demography not only reflect the current geographic distribution but also allow us to understand the adaption and genetic differentiation of species. However, the process and factors affecting the present distribution and genetic structure of many taxa are still poorly understood, especially for endemic organisms to small islands.

Methods

Here, we integrated population genetic and ecological niche modelling approaches to investigate the historical distribution and demographic dynamics of two co-existing salamanders on Tsushima Island, Japan: the true H. tsuensis (Group A), and Hynobius sp. (Group B). We also examined the hypothesis on the equivalency and similarity of niches of these groups by identity and background tests for ecological niche space.

Results

Our result showed that Group A is considered to have undergone a recent population expansion after the Last Glacial Maximum while it is unlikely to have occurred in Group B. The highest suitability was predicted for Group A in southern Tsushima Island, whereas the northern part of Tsushima Island was the potential distribution of Group B. The results also suggested a restricted range of both salamanders during the Last Interglacial and Last Glacial Maximum, and recent expansion in Mid-Holocene. The genetic landscape-shape interpolation analysis and historical suitable area of ecological niche modelling were consistent, and suggested refugia used during glacial ages in southern part for Group A, and in northern part of Tsushima Island for Group B. Additionally, we found evidence of nonequivalence for the ecological niche of the two groups of the salamanders, although our test could not show either niche divergence or conservatism based on the background tests. The environmental predictors affecting the potential distribution of each group also showed distinctiveness, leading to differences in selecting suitable areas. Finally, the combination of population genetics and ecological modeling has revealed the differential demographic/historical response between coexisting two salamanders on a small island.

Karyological Diversification in the Genus Lyciasalamandra (Urodela: Salamandridae)

Simple Summary

The Lycian salamanders of the genus Lyciasalamandra are characterized by a debated taxonomy and phylogenetic relationships. They have been the subject of various molecular and phylogenetic analyses, but their chromosomal diversity is completely unknown. We here present a comparative cytogenetic analysis on five out of the seven described species and seven subspecies of Lyciasalamandra, providing the first karyological assessment on the genus and comparing them to closely related representatives of the genus Salamandra. We analyzed the occurrence and distribution of different conserved (chromosome number and morphology) and highly variable karyological features. We found an impressive diversity in the configuration of nucleolus organizing regions (NORs), which alternatively occur either as heteromorphic or homomorphic loci on distinct regions of different chromosome pairs. We highlight that the observed peculiar taxon-specific pattern of chromosome markers supports the taxonomic validity of the different studied evolutionary lineages and is consistent with a scenario of synchronous evolution in the Lycian salamanders.

Abstract

We performed the first cytogenetic analysis on five out of the seven species of the genus Lyciasalamandra, including seven subspecies, and representatives of its sister genus Salamandra. All the studied species have a similar karyotype of 2n = 24, mostly composed of biarmed elements. C-bands were observed on all chromosomes, at centromeric, telomeric and interstitial position. We found a peculiar taxon-specific NOR configuration, including either heteromorphic and homomorphic NORs on distinct regions of different chromosomes. Lyciasalamandra a.antalyana and L. helverseni showed two homomorphic NORs (pairs 8 and 2, respectively), while heteromorphic NORs were found in L. billae (pairs 6, 12), L. flavimembris (pairs 2, 12), L. l. luschani (pairs 2, 12), L. l. basoglui (pairs 6, 12), L. l. finikensis (pairs 2, 6) and S. lanzai (pairs 8, 10). Homomorphic NORs with an additional supernumerary site were shown by S. s. salamandra (pairs 2, 8) and S. s. gigliolii (pairs 2, 10). This unexpected highly variable NOR configuration is probably derived from multiple independent NOR translocations and paracentric inversions and correlated to lineage divergence in Lyciasalamandra. These results support the taxonomic validity of the studied taxa and are consistent with a hypothesized scenario of synchronous evolution in the genus.

Discovery of Two Lineages of Hynobius tsuensis (Amphibia, Caudata) Endemic to Tsushima Island, Japan

The Tsushima salamander, Hynobius tsuensis, is a lotic-breeding species endemic to Tsushima Island, Japan. We investigated genetic variation in this species using samples covering all known localities on the island. Phylogenetic analyses based on mitochondrial cytochrome b gene sequences divided H. tsuensis populations into two genetic groups: true H. tsuensis and an unknown species (Hynobius sp.) showing a closer relationship to Hynobius nebulosus than to H. tsuensis. The genetic distance between the two groups was greater than the genetic distances between different species of Hynobius. Both groups were found to breed in mountain streams, sometimes in the same streams syntopically. Population structure analysis based on several nuclear microsatellite loci strongly indicated reproductive isolation between the two groups in the syntopic localities. This is the first report confirming the syntopic occurrence of two lotic salamanders of Hynobius on a small island. We propose that the lotic breeding habits of the two salamanders were independently acquired.

Historic range dynamics in Kaiser's mountain newt (Neurergus kaiseri): Insights from phylogeographic analyses and species distribution modeling

Vulnerable Kaiser's mountain newt, Neurergus kaiseri, is endemic to highland streams, springs, and pools of the southwestern Zagros mountain, Iran. The present study aimed to use an integration of phylogeographical and species distribution modeling (SDM) approaches to provide new insights into the evolutionary history of the species throughout Quaternary climate oscillations. The phylogeographical analysis was followed by analyzing two mitochondrial DNA (mt-DNA) markers including 127 control region (D-loop) and 72 NADH dehydrogenase 2 (ND2) sequences from 15 populations in the entire species range that were obtained from GenBank. Potential recent and past distribution (the Last Glacial Maximum, LGM, 21 Kya and the Mid-Holocene, 6 Kya) reconstructed by ensemble SDM using nine algorithms with CCSM4, MIROC-ESM, and MPI-ESM-P models. N. kaiseri displayed two distinct lineages in the northern and southern regions that diverged in the Early-Pleistocene. The demographics analysis showed signs of a slight increase in effective population size for both northern and southern populations in the Mid-Pleistocene. Biogeography analysis showed that both vicariance and dispersal events played an important role in the formation of recent species distribution of N. kaiseri. Based on SDM projection onto paleoclimatic data, N. kaiseri displayed a scenario of past range expansion that followed by postglacial contraction. The models showed that the distribution range of the species may have shifted to a lower altitude during LGM while with amelioration of climatic during Mid-Holocene to recent conditions caused the species to shift to the higher altitude. The findings of the current study support the hypothesis that the Zagros mountains​ may be acting as climatic refugia and play an important role in the protection of isolated populations during climate oscillations.

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

Wide Distribution but Low Differentiation: Phylogenetic Relationships and Phylogeography of Hynobius nigrescens (Amphibia: Caudata)

We surveyed the genetic structure of Hynobius nigrescens Stejneger, 1907, a lentic breeding salamander widespread throughout montane and lowland regions of northeastern Japan. We performed a mitochondrial DNA analysis to explore intraspecific genetic variation and infer the evolutionary population history of H. nigrescens. Complete 1141 bp sequences of the mitochondrial cytochrome b gene were studied for 134 adult and larval individuals collected from 62 localities, encompassing the known range of the species. Hynobius nigrescens proved to be monophyletic, including two major clades (Clade II from southwestern Hokuriku and Clade I from all other localities). The latter clade comprises four well-supported and geographically structured subclades, which show genetic distances smaller than those seen in the widely sympatric species Hynobius lichenatus. Results of population statistical analyses indicated that Clade II of the westernmost range of H. nigrescens seems to have maintained a constant population size, while Clade I from most of the northeastern species range shows a tendency of recent population expansion, which is evident in Subclades I-A from the northernmost range and I-B from southern Tohoku to northern Kanto and eastern Chubu. In contrast, Subclades I-C from northeastern Chubu and Sado Is. and I-D from northwestern Chubu to Hokuriku seem to have been relatively stable in population size. Hynobius nigrescens differs greatly from other salamander species from northeastern Japan in its much more recent periods of genetic differentiation and its pattern thereof, and is suggested to be a young faunal element in this region.

Genetic diversity, phylogenetic position and morphometric analysis of Astacus colchicus (Decapoda, Astacidae): a new insight into Eastern European crayfish fauna

The phylogeny of European crayfish fauna, especially with respect to Eastern European species, is still far from being completely resolved. To fill this gap, we analyzed most of the European crayfish species focusing on the phylogenetic position of the endemic crayfish Astacus colchicus, inhabiting Georgia. Three mitochondrial and one nuclear marker were used to study evolutionary relationships among European crayfish species, resulting in the unique phylogenetic position of A. colchicus indicating independent species status to A. astacus. Phylogenetic analyses revealed a deep molecular divergence of A. colchicus in comparison to A. astacus (6.5-10.9% in mtDNA and 1.1% in nDNA) as well as to Pontastacus leptodactylus and P. pachypus (5.5-10.0% in mtDNA and 1.4-2.4% in nDNA). Absent ventral process on second male pleopod and abdominal somites II and III with pleura rounded lacking prominent spines clearly indicate taxonomic assignment to the genus Astacus; however, the species is distributed almost in the middle of Ponto-Caspian area typical by occurrence of the genus Pontastacus. Several morphological indices linked to head length, carapace, and total body length and width were found to demonstrate apparent differences between A. colchicus and A. astacus. Although this study provides a novel insight into European crayfish phylogeography, we also point out the gaps in comprehensive study of the P. leptodactylus species complex, which could reveal details about the potential species status of particular species and subspecies within this genus.

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

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

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow-spotted mountain newt, Neurergus derjugini (Caudata: Salamandridae) in the mid-Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Reconstructing reliable timescales for species evolution is an important and indispensable goal of modern biogeography. However, many factors influence the estimation of divergence times, and uncertainty in the inferred time trees remains a major issue that is often insufficiently acknowledged. We here focus on a fundamental problem of time tree analysis: the combination of slow-evolving (nuclear DNA) and fast-evolving (mitochondrial DNA) markers in a single time tree. Both markers differ in their suitability to infer divergences at different time scales (the 'genome-timescale-dilemma'). However, strategies to infer shallow and deep divergences in a single time tree have rarely been compared empirically. Using Mediterranean amphibians as model system that is exceptional in its geographic and taxonomic completeness of available genetic information, we analyze 202 lineages of western Palearctic amphibians across the entire Mediterranean region. We compiled data of four nuclear and five mitochondrial genes and used twelve fossil calibration points widely acknowledged for amphibian evolution. We reconstruct time trees for an extensive lineage-level data set and compare the performances of the different trees: the first tree is based on primary fossil calibration and mitochondrial DNA, while the second tree is based on a combination of primary fossil and on secondary calibrations taken from a nuclear tree using mitochondrial DNA (two-step protocol). Focusing on a set of nodes that are most likely explained by vicariance, we statistically compare the reconstructed alternative time trees by applying a biogeographical plausibility test. Our two-step protocol outperformed the alternative approach in terms of spatial and temporal plausibility. It allows us to infer scenarios for Mediterranean amphibian evolution in eight geographic provinces. We identified several tectonic and climatic events explaining the majority of Mediterranean amphibian divergences, with Plio-Pleistocene climatic fluctuations being the dominant driver for intrageneric evolution. However, often more than one event could be invoked for a specific split. We give recommendations for the use of secondary calibrations in future molecular clock analyses at the community level.

Phylogeographic investigation and ecological niche modelling of the endemic frog species Nanorana pleskei revealed multiple refugia in the eastern Tibetan Plateau

The largest plateau Tibetan Plateau supplied an excellent opportunity to investigate the influence of the Pleistocene events on the high-elevation species. To test for the alternative hypotheses of Pleistocene glacial refugia, we used partial sequences of two mitochondrial genes and one nuclear gene to examine the phylogeographic patterns of the endemic frog species Nanorana pleskei across its known range in the eastern Tibetan Plateau, and conducted species distribution modelling (SDM) to explore changes of its distribution range through current and paleo periods. In all data sets, the species was divided into lineage north occupying open plateau platform and lineage south colonizing the mountainous plateau. The divergence of two major clades was estimated at the early Pleistocene. In mtDNA, lineage north contained northeastern and northwestern sublineages, and lineage south had two overlapping-distributed sublineages. Different lineages possessed distinct demographic characteristics, i.e., subdivision in the northeastern sublineage, historical bottleneck effects and recent expansions in the northwestern sublineage and the southeastern sublineage. SDMs depicted that stable suitable habitats had existed in the upper-middle streams of the Yellow River, Dadu River, Jinsha River and Yalong River. These regions were also recognized as the ancestral areas of different lineages. In conclusion, Nanorana pleskei lineages have probably experienced long-term separations. Stable suitable habitats existing in upper-middle streams of major rivers on the eastern Tibetan Plateau and distinct demographic dynamics of different lineages indicated that the lineages possessed independent evolutionary processes in multiple glacial refugia. The findings verified the profound effects of Pleistocene climatic fluctuations on the plateau endemic species.

The geography and timing of genetic divergence in the lizard Phrynocephalus theobaldi on the Qinghai-Tibetan plateau

The Qinghai-Tibetan Plateau (QTP) represents one of the earth’s most significant physical features and there is increasing interest in the historical generation of biodiversity within this region. We hypothesized that there should be clear geographically coherent genetic structuring within one of the world’s highest altitude lizards, Phrynocephalus theobaldi, due to considerable historical population fragmentation in this environment. This was tested using a major mitochondrial DNA (mtDNA) survey and sequencing of two nuclear markers (AME and RAG-1) from P. theobaldi, from across the southern QTP. A Bayesian method (BPEC) was used to detect four geographically structured mtDNA clusters. A Bayesian phylogenetic tree, together with associated dating analyses, supported four corresponding evolutionary lineages with a timing of 3.74–7.03 Ma for the most basal P. theobaldi split and Pliocene splits of 2.97–5.79 Ma and 2.40–5.39 Ma in the two daughter lineages. Himalayan uplift and changes in the Jilong basin may have contributed to these divergences, but uplift of the Gangdese mountains is rejected due to its timing. The nuclear markers appeared to be sorted between the four mtDNA groups, and species delimitation analyses supported the four phylogeographical groups as candidate species. The study contributes to our understanding of biodiversity on the QTP.

Historical Demography of an Endangered Salamander, Ranodon sibiricus (Amphibia, Urodela, Hynobiidae): A Reassessment

Chinese populations of the endangered Siberian salamander Ranodon sibiricus are reported to have diverged only about 120 years ago, and to have the lowest genetic diversity of any amphibian. However, these conclusions require verification, as the main range of the species is in Kazakhstan. Moreover, the generation time used for estimating divergence time has a weak ground. In order to clarify these problems, we investigated the molecular phylogenetic relationship and historical demography of the species covering its whole distribution range using the mitochondrial DNA region reported for Chinese population (1072 bp sequences of the control region), while conducting skeletochronological analysis to estimate accurate generation time. As a result, the range expansion was estimated at 88,000-50,000 YA, based on the generation time of 6-10 years. Degree of intraspecific genetic differentiation is actually very small, but, as a single species, is not so small as had been reported for Chinese population alone.

Divergence history of the Carpathian and smooth newts modelled in space and time

Information about demographic history is essential for the understanding of the processes of divergence and speciation. Patterns of genetic variation within and between closely related species provide insights into the history of their interactions. Here, we investigated historical demography and genetic exchange between the Carpathian (Lissotriton montandoni, Lm) and smooth (L. vulgaris, Lv) newts. We combine an extensive geographical sampling and multilocus nuclear sequence data with the approximate Bayesian computation framework to test alternative scenarios of divergence and reconstruct the temporal and spatial pattern of gene flow between species. A model of recent (last glacial period) interspecific gene flow was favoured over alternative models. Thus, despite the relatively old divergence (4-6 mya) and presumably long periods of isolation, the species have retained the ability to exchange genes. Nevertheless, the low migration rates (ca. 10(-6) per gene copy per generation) are consistent with strong reproductive isolation between the species. Models allowing demographic changes were favoured, suggesting that the effective population sizes of both species at least doubled as divergence reaching the current ca. 0.2 million in Lm and 1 million in Lv. We found asymmetry in rates of interspecific gene flow between Lm and one evolutionary lineage of Lv. We suggest that intraspecific polymorphism for hybrid incompatibilities segregating within Lv could explain this pattern and propose further tests to distinguish between alternative explanations. Our study highlights the importance of incorporating intraspecific genetic structure into the models investigating the history of divergence.

Seven at one blow: the origin of major lineages of the viviparous Lycian salamanders (Lyciasalamandra Veith and Steinfartz, 2004) was triggered by a single paleo-historic event

The number of tectonic and climatic events that are used to explain speciation processes in the eastern Mediterranean region is low compared to the western Mediterranean. Among them, the emergence of the mid-Aegean trench and the Messinian Salinity Crisis (MSC) often concurred with speciation time estimates that were inferred from molecular data. We here present a dated molecular phylogeny of Lyciasalamandra from Turkey and Greece based on ca. 4500 bp of the mitochondrial genome (3000 bp of three nuclear genes appeared to be completely inconclusive due to their extremely low degree of variation among taxa). Seven major lineages emerged simultaneously from a basal hard polytomy. A scenario that dates this polytomy to 12.3 and 10.2 million years ago, around the final emergence of the mid-Aegean trench, appears to be most plausible. The MSC can be made responsible for first intraspecific divergence events within L. luschani, L. fazilae and L. flavimembris. Further diversification can be explained by Pliocene and Pleistocene glaciations. Based on levels of molecular differentiation we suggest the recently described species L. arikani, L. irfani and L. yehudahi to be treated as subspecies of L. billae.

Biogeography and evolution of a widespread Central American lizard species complex: Norops humilis, (Squamata: Dactyloidae)

Background

Caribbean anole lizards (Dactyloidae) have frequently been used as models to study questions regarding biogeography and adaptive radiations, but the evolutionary history of Central American anoles (particularly those of the genus Norops) has not been well studied. Previous work has hypothesized a north-to-south dispersal pattern of Central American Norops, but no studies have examined dispersal within any Norops lineages. Here we test two major hypotheses for the dispersal of the N. humilis/quaggulus complex (defined herein, forming a subset within Savage and Guyer’s N. humilis group).

Results

Specimens of the N. humilis group were collected in Central America, from eastern Mexico to the Canal Zone of Panama. Major nodes were dated for comparison to the geologic history of Central America, and ancestral ranges were estimated for the N. humilis/quaggulus complex to test hypothesized dispersal patterns. These lineages displayed a northward dispersal pattern. We also demonstrate that the N. humilis/quaggulus complex consists of a series of highly differentiated mitochondrial lineages, with more conserved nuclear evolution. The paraphyly of the N. humilis species group is confirmed. A spatial analysis of molecular variance suggests that current populations are genetically distinct from one another, with limited mitochondrial gene flow occurring among sites.

Conclusions

The observed south-to-north colonization route within the Norops humilis/quaggulus complex represents the first evidence of a Norops lineage colonizing in a south-to-north pattern, (opposite to the previously held hypothesis for mainland Norops). One previously described taxon (N. quaggulus) was nested within N. humilis, demonstrating the paraphyly of this species; while our analyses also reject the monophyly of the Norops humilis species group (sensu Savage and Guyer), with N. tropidonotus, N. uniformis, and N. marsupialis being distantly related to/highly divergent from the N. humilis/quaggulus complex. Our work sheds light on mainland anole biogeography and past dispersal events, providing a pattern to test against other groups of mainland anoles.

Electronic supplementary material

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

The Phylogeography and Population Demography of the Yunnan Caecilian (Ichthyophis bannanicus): Massive Rivers as Barriers to Gene Flow

Ichthyophis bannanicus is the only caecilian species in China. In this study, the phylogeography and population demography of I. bannanicus were explored, based on the mitochondrial DNA genes (cyt b and ND2) and 15 polymorphic microsatellite loci. Altogether 158 individuals were collected from five populations in Yunnan province, Guangxi province, Guangdong province, and Northern Vietnam. Phylogeographical and population structure analysis identified either two groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun-Deqing) or three groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun, and Deqing), indicating that the Red River and Pearl River systems may have acted as gene-flow barriers for I. bannanicus. Historical population expansion that happened 15–17 Ka ago was detected for mtDNA data and was possibly triggered by warmer weather after the Last Glacial Maximum. However, the Bayesian simulations of population history based on microsatellite data pinpointed population decline in all populations since 19,123 to 1,029 years ago, demonstrating a significant influence of anthropogenic habitat alteration on I. bannanicus.

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.

Population structure, historical biogeography and demographic history of the alpine toad Scutiger ningshanensis in the Tsinling Mountains of Central China

Population genetic structure, historical biogeography and historical demography of the alpine toad Scutiger ningshanensis were studied using the combined data mtDNA cytochrome b (cyt b) and the mtDNA cytochrome c oxidase subunit I (COI) as the molecular markers. This species has high genetic variation. There was a significant genetic differentiation among most populations. Three lineages were detected. The phylogenetic relationship analyses and the SAMOVA (spatial analysis of molecular variance) results showed significant phylogeographic structure. 82.15% genetic variation occurred among populations whereas differentiation within populations only contributed 17.85% to the total. Mantel test results showed a significant correlation between the pairwise calculated genetic distance and pairwise calculated geographical distance of the populations (regression coefficient  = 0.001286, correlation coefficient  = 0.77051, p (rrand≥robs)  = 0.0185<0.05), indicating the existence of isolation-by-distance pattern of genetic divergence for cyt b + COI sequence, which suggests that the distribution of genetic variation is due to geographical separation rather than natural selection. The population expansion or contraction and genetic differentiation between populations or lineages could be explained by topography and the repetitive uplifts of the Tsinling Mountains and the climatic cycles during the late Pliocene and Pleistocene. S. ningshanensis experienced a rapid population expansion about 40,000 years before present. The current decline in population size was probably caused by anthropogenic disturbance. Current populations of S. ningshanensis are from different refugia though the location of these refugia could not be determined in our study. Topography, climatic changes and repetitive population expansion/contraction together led to the high level of genetic variation in S. ningshanensis. A total of three management units (MUs) was determined, which must be considered when conservation policy is made in the future.

A phylogeographic evaluation of the Amolops mantzorum species group: cryptic species and plateau uplift

Phylogeographic analysis brings organism phylogeny, regional geological/climatic history, and population demography together, and provides insights into species history and speciation. Using DNA sequence data from a mitochondrial marker (cytochrome b) and a nuclear marker (cmyc intron 2), we examined phylogeography of the Amolops mantzorum species group. We first constructed Bayesian and maximum-likelihood gene trees and medium-joint networks for the recovered haplotypes, and estimated divergence time for each major lineage. Species boundaries were then delineated using the general mixed Yule-coalescent model (GMYC) and a Bayesian species-delimitation method (BP&P). Potential gene flow between putative species was also estimated using the isolation-with-migration model. Furthermore, species-distribution modeling was used to probe linkage between climatic changes and speciation. Lastly, extended Bayesian skyline plotting was employed to reveal historical demography for each putative species. Our analyses clearly delimit nine species in the group, including four well-recognized species and five putative species, of which three are often listed as synonyms of A. mantzorum. The nominal species A. mantzorum may in fact include two cryptic species. Divergence-time estimates align the speciation events with the recent intense uplift of the Tibetan Plateau in the last 3.6 million years. The species-distribution modeling also reveals different habitat preferences among species that are potentially linked to climatic changes associated with the uplift. Furthermore, three species sustained continuous population growth through the last glaciation, while others expanded only after the last glaciation. The eastern escarpment of Tibet is a biodiversity hotspot; its rich species diversity, unique topography, and well-established geological history provide an excellent system for evolutionary studies.

Molecular phylogeography and population genetic structure of an endangered species Pachyhynobius shangchengensis (hynobiid Salamander) in a fragmented habitat of southeastern China

The salamander Pachyhynobius shangchengensis (Hynobiidae) is a vulnerable species restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes in the Mount Dabieshan region in southeastern China. However, molecular phylogeography and population genetic structure of P. shangchengensis remain poorly investigated. In this study, we explored the genetic structure and phylogeography of P. shangchengensis based on partial sequences of the mitochondrial DNA (mtDNA) cytochrome b and cytochrome c oxidase subunit I genes. Fifty-one haplotypes and four clades were found among 93 samples. Phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographic regions separated by complicated topography and long distances. The distinct geographic distributions of all lineages and the estimated divergence time suggest spatial and temporal separation coinciding with climatic changes during the Pleistocene. Analysis of molecular variance indicated that most of the observed genetic variation occurred among the four groups, implying long-term interruption of gene flow, and the possible separation of P. shangchengensis into four management units for conservation.

Mitochondrial cytochrome B phylogeny and historical biogeography of the Tohoku salamander, Hynobius lichenatus (Amphibia, Caudata)

The Tohoku salamander, Hynobius lichenatus Boulenger, 1883, is a lentic breeding species widespread throughout montane regions of northeastern Japan. To explore intraspecific genetic variation and infer evolutionary history of H. lichenatus, we performed mitochondrial DNA analysis (complete 1141 bp sequences of the mitochondrial cytochrome b gene) using 215 adult and larval individuals collected from 75 localities, encompassing known distributional range of the species. Hynobius lichenatus proved to be monophyletic, including three well-supported and geographically structured clades (Clade I from northern Kanto, Clade II from southern Tohoku, and Clade III from northern Tohoku). These clades, respectively, comprise several subclades, and show genetic distances as large as those seen between different species of Hynobius. Results of population statistic analyses indicate that all clades and most subclades have maintained high genetic diversity and demographic stability over long periods. Molecular dating indicates divergence in H. lichenatus concords with topographic evolution of northeastern Japan from late Miocene to early Pleistocene, suggesting that paleogeographic events in this region, such as orogenesis, sea level change, and volcanic activity, have been crucial for shaping genetic patterns and diversity in this species. Hynobius lichenatus greatly differs from many other animal species from northeastern Japan in its much older periods and the pattern of genetic differentiation, and is suggested as an old faunal element in this region.

Phylogeny and genetic history of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870) inferred from complete mitochondrial genomes

We assessed phylogeny of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), the most northern ectothermic, terrestrial vertebrate in Eurasia, by sequence analysis of complete mitochondrial genomes in 26 specimens from different localities (China, Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka, Ural, European part of Russia). In addition, a complete mitochondrial genome of the Schrenck salamander, Salamandrella schrenckii, was determined for the first time. Bayesian phylogenetic analysis of the entire mtDNA genomes of S. keyserlingii demonstrates that two haplotype clades, AB and C, radiated about 1.4 million years ago (Mya). Bayesian skyline plots of population size change through time show an expansion around 250 thousand years ago (kya) and then a decline around the Last Glacial Maximum (25 kya) with subsequent restoration of population size. Climatic changes during the Quaternary period have dramatically affected the population genetic structure of the Siberian salamanders. In addition, complete mtDNA sequence analysis allowed us to recognize that the vast area of Northern Eurasia was colonized only by the Siberian salamander clade C1b during the last 150 kya. Meanwhile, we were unable to find evidence of molecular adaptation in this clade by analyzing the whole mitochondrial genomes of the Siberian salamanders.

Phylogeographic patterns of mtDNA variation revealed multiple glacial refugia for the frog species Feirana taihangnica endemic to the Qinling Mountains

Diversification patterns and demography of montane species are affected by Pleistocene climate fluctuations. Empirical cases from the Qinling Mountains (QM) region, which is a major biogeographic divider of East Asia, are few. We used DNA sequence data of the complete mitochondrial ND2 gene to detect effects of the Pleistocene glaciations on phylogeographic profiles of a frog species, Feirana taihangnica, which is endemic to the QM. Four distinct lineages consisting of seven sublineages were revealed. The strongest signal of biogeographical structure (F(ct) = 0.971, P < 0.01) was found when populations were grouped according to these seven sublineages. One narrow secondary contact zone was detected in the middle QM between the lineage from middle QM and the lineage from eastern QM. Coalescent simulations indicated that this species colonized the QM region by a stepping-stone model. Divergences among lineages had likely been influenced by the uplift of the Tibetan Plateau during the late Miocene-to-late Pleistocene, as well as by the Pleistocene climatic cycles. Coalescent simulations also suggested that F. taihangnica populations have persisted through the Pleistocene glacial periods in multiple refugia across the QM region. Demographic analyses indicated that all lineages, except the lineage in the Funiu Mountains, have been experienced postglacial expansion of population size and distribution range. In conclusion, Pleistocene climate fluctuations and tectonic changes during the late Miocene-late Pleistocene have profoundly influenced the phylogeography and historical demography of F. taihangnica.

Population structure of the salamander Hynobius retardatus inferred from a partial sequence of the mitochondrial DNA control region

We investigated population structure of the salamander Hynobius retardatus in Hokkaido, Japan using partial sequences of the mitochondrial DNA control region (490 bp) from 105 individuals. The salamanders were collected from 28 localities representing the entire regional distribution of this species. Twenty different haplotypes distributed across three haplotype groups were identified. Group 1 was widely distributed in central, northern, and eastern Hokkaido, except Erimo; Groups 2 and 3 appeared exclusively in Erimo and southern Hokkaido, respectively. The genetic distance between the three groups was not very large, but the distributions of the groups never overlapped spatially, indicating a hierarchical population structure comprising three regional groups, which was also supported by analysis of molecular variance. The results suggest that the present population structure is affected by current genetic barriers, as well as by historical transitions of climate and landscape.