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

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.

Significance of pre-Quaternary climate change for montane species diversity: insights from Asian salamanders (Salamandridae: Pachytriton)

Despite extensive focus on the genetic legacy of Pleistocene glaciation, impacts of earlier climatic change on biodiversity are poorly understood. Because amphibians are highly sensitive to variations in precipitation and temperature, we use a genus of Chinese montane salamanders (Salamandridae: Pachytriton) to study paleoclimatic change in East Asia, which experienced intensification of its monsoon circulation in the late Miocene associated with subsequent Pliocene warming. Using both nuclear and mitochondrial DNA sequences, we reconstruct the species tree under a coalescent model and demonstrate that all major lineages originated before the Quaternary. Initial speciation within the genus occurred after the summer monsoon entered a stage of substantial intensification. Heavy summer precipitation established temporary water connectivity through overflows between adjacent stream systems, which may facilitate geographic range expansion by aquatic species such as Pachytriton. Species were formed in allopatry likely through vicariant isolation during or after range expansion. To evaluate the influence of Pliocene warming on these cold-adapted salamanders, we construct a novel temperature buffer-zone model, which suggests widespread physiological stress or even extinction during the warming period. A significant deceleration of species accumulation rate is consistent with Pliocene range contraction, which affected P. granulosus and P. archospotus the most because they lack large temperature buffer zones. In contrast, demographic growth occurred in species for which refugia persist. The buffer-zone model reveals the Huangshan Mountain as a potential climatic refugium, which is similar to that found for other East Asian organisms. Our approach can incorporate future climatic data to evaluate the potential impact of ongoing global warming on montane species (particularly amphibians) and to predict possible population declines.

Phylogeny and biogeography of the Anderson's crocodile newt, Echinotriton andersoni (Amphibia: Caudata), as revealed by mitochondrial DNA sequences

The Anderson's crocodile newt, Echinotriton andersoni, is considered a relic and endangered species distributed in the Central Ryukyus. To elucidate phylogenetic relationships and detailed genetic structures among populations, we analyzed variation in the mitochondrial cytochrome b gene. Results strongly support a primary dichotomy between populations from the Amami and Okinawa Island Groups with substantial genetic divergence, favoring a primary divergence between the two island groups. Within the latter, populations from the southern part of Okinawajima Island are shown to be more closely related to those from Tokashikijima Island than to those from the northern and central parts of Okinawajima. The prominent genetic divergence between the two island groups of the Central Ryukyus seems to have initiated in the Miocene, i.e., prior to formation of the strait that has consistently separated these island groups since the Pleistocene. The ancestor of the southern Okinawajima-Tokashikijima is estimated to have migrated from the northern and central parts of Okinawajima into southern Okinawajima at the Pleistocene, and dispersed into Tokashikijima subsequently.

Postglacial colonization of the Qinling Mountains: phylogeography of the swelled vent frog (Feirana quadranus)

BACKGROUND

The influence of Pleistocene climatic fluctuations on intraspecific diversification in the Qinling-Daba Mountains of East Asia remains poorly investigated. We tested hypotheses concerning refugia during the last glacial maximum (LGM) in this region by examining the phylogeography of the swelled vent frog (Feirana quadranus; Dicroglossidae, Anura, Amphibia).

METHODOLOGY/PRINCIPAL FINDINGS

We obtained complete mitochondrial ND2 gene sequences of 224 individuals from 34 populations of Feirana quadranus for phylogeographic analyses. Additionally, we obtained nuclear tyrosinase gene sequences of 68 F. quadranus, one F. kangxianensis and three F. taihangnica samples to test for mitochondrial introgression among them. Phylogenetic analyses based on all genes revealed no introgression among them. Phylogenetic analyses based on ND2 datasets revealed that F. quadranus was comprised of six lineages which were separated by deep valleys; the sole exception is that the Main Qinling and Micang-Western Qinling lineages overlap in distribution. Analyses of population structure indicated restricted gene flow among lineages. Coalescent simulations and divergence dating indicated that the basal diversification within F. quadranus may be associated with the dramatic uplifts of the Tibetan Plateau during the Pliocene. Coalescent simulations indicated that Wuling, Daba, and Western Qinling-Micang-Longmen Mountains were refugia for F. quadranus during the LGM. Demographic analyses indicated that the Daba lineage experienced population size increase prior to the LGM but the Main Qinling and the Micang-Western Qinling lineages expanded in population size and range after the LGM, and the other lineages almost have stable population size or slight slow population size decline.

CONCLUSIONS/SIGNIFICANCE

The Qinling-Daba Mountains hosted three refugia for F. quadranus during the LGM. Populations that originated in the Daba Mountains colonized the Main Qinling Mountains after the LGM. Recent sharp expansion of the Micang-Western Qinling and Main Qinling lineages probably contribute to their present-day secondary contact.

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7–10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3–7 m) nor the wider nuclear cline (RAG-2; 34–42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2–3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.

Mitochondrial DNA data unveil highly divergent populations within the genus Hynobius (Caudata: Hynobiidae) in South Korea

Korean salamanders of the genus Hynobius are currently classified into 3 species, H. leechii, H. quelpaertensis, and H. yangi. To investigate the phylogenetic relationship of these species, we analyzed the partial sequence of mitochondrial cytochrome b gene (907 bp) of 197 specimens from 43 regions in South Korea. Of these specimens, 93 were additionally examined with 12S rRNA (799 bp). Based on the partial sequence of the mitochondrial cytochrome b gene and 12S rRNA, 89 and 36 haplotypes were defined, respectively, consisting of six subclades (H. leechii, H. quelpaertensis, H. yangi, HC1, HC2, and HC3). Among these subclades, the three subclades (HC1, HC2, and HC3) were clearly separated from the 3 previously reported species in the genus Hynobius. Pairwise sequence divergence between the six subclades ranged from 6.3 to 11.2% in cytochrome b gene and 2.0 to 4.3% in 12S rRNA. These results indicate there may be more divergent populations than the three currently described. Moreover, the estimation of divergence time revealed that the Hynobius species in South Korea diverged during the Miocene epoch, approximately 9 - 5 MYA. In addition, we confirmed the distribution of the three known species (H. leechii, H. quelpaertensis, and H. yangi) and determined the distributions of new, distinct groups (or subclades; HC1, HC1, and HC3). To more accurately establish the taxonomic status and population structure, further genetic, morphological, and ecological studies will be needed.

Widespread introgression does not leak into allotopy in a broad sympatric zone

Species that overlap over a large part of their range and habitat requirements are challenging for the study of speciation and hybridization. In this respect, the study of broadscale introgressive hybridization has raised recent interest. Here we studied hybridization between two closely related amphibians Lissotriton helveticus and Lissotriton vulgaris that reproduce over a wide sympatric zone. We used mitochondrial and microsatellite markers on 1272 individuals in 37 sites over Europe to detect hybrids at the individual-level and to analyse Hardy-Weinberg and linkage disequilibria at the population-level. Morphological traits showed a strong bimodal distribution. Consistently, hybrid frequency was low (1.7%). We found asymmetric introgression with five times more hybrids in L. vulgaris than in L. helveticus, a pattern probably explained by an unequal effective population size in a study part wherein L. helveticus numerically predominates. Strikingly, significant levels of introgression were detected in 73% of sites shared by both species. Our study showed that introgression is widespread but remains confined to the sites where the two species reproduce at the same time. This pattern may explain why these species remain genetically distinct over a broad sympatric zone.

Patterns of differential introgression in a salamander hybrid zone: inferences from genetic data and ecological niche modelling

Hybrid zones have yielded considerable insight into many evolutionary processes, including speciation and the maintenance of species boundaries. Presented here are analyses from a hybrid zone that occurs among three salamanders -Plethodon jordani, Plethodon metcalfi and Plethodon teyahalee- from the southern Appalachian Mountains. Using a novel statistical approach for analysis of non-clinal, multispecies hybrid zones, we examined spatial patterns of variation at four markers: single-nucleotide polymorphisms (SNPs) located in the mtDNA ND2 gene and the nuclear DNA ILF3 gene, and the morphological markers of red cheek pigmentation and white flecks. Concordance of the ILF3 marker and both morphological markers across four transects is observed. In three of the four transects, however, the pattern of mtDNA is discordant from all other markers, with a higher representation of P. metcalfi mtDNA in the northern and lower elevation localities than is expected given the ILF3 marker and morphology. To explore whether climate plays a role in the position of the hybrid zone, we created ecological niche models for P. jordani and P. metcalfi. Modelling results suggest that hybrid zone position is not determined by steep gradients in climatic suitability for either species. Instead, the hybrid zone lies in a climatically homogenous region that is broadly suitable for both P. jordani and P. metcalfi. We discuss various selective (natural selection associated with climate) and behavioural processes (sex-biased dispersal, asymmetric reproductive isolation) that might explain the discordance in the extent to which mtDNA and nuclear DNA and colour-pattern traits have moved across this hybrid zone.

Phylogeography and molecular adaptation of Siberian salamander Salamandrella keyserlingii based on mitochondrial DNA variation

We assessed the phylogeographic pattern of Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), which appear to be the most northern ectothermic, terrestrial vertebrate in Northern Eurasia, by sequence analysis of a 611-bp fragment of the mitochondrial cytochrome b gene in 159 specimens from different localities (Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka and others). The data revealed that cytochrome b lineages of S. keyserlingii are divided into haplogroups A, B and C. Haplogroup A and B sequences are widespread in the Far East region, whereas haplogroup C consisting of several phylogenetic clusters (C1, C2, C3) is present in the all range of S. keyserlingii. Among them, cluster C3 appears to be specific for Sakhalin; most likely, it has arisen in situ in this island, with the entry time of the founder mtDNA estimated at about 0.4 MY. Analysis of cytochrome b gene variation by using different neutrality tests (including those based on K(A)/K(S)-ratio) has shown that differences between haplogroups were statistically insignificant, thus suggesting selective neutrality. However, analysis of amino acid changes allowed us to detect a signature of molecular adaptation, which might have led to appearance of adaptive cytochrome b variants in haplogroup C, originating most likely at the end of Eopleistocene (about 0.64 MY based on the haplogroup C divergence level). It seems probable that this adaptive mechanism could promote subsequent populating of new regions.

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

Evolutionary relationships of the tree frogs from the Middle East and the demographic histories of their populations were studied using a combination of mitochondrial and nuclear genes. Hyla savignyi and neighboring populations of H. orientalis (former eastern populations of H. arborea) were the main focus taxa. Within H. savignyi, a deep phylogenetic divergence dated about 8.4Ma was discovered. Southern populations from Yemen, Jordan, southern Syria and extreme north-eastern Israel are hereby described as a new species, H. felixarabica sp. nov. Our study points to a biogeographic connection of the south-western Arabian Peninsula and southern Levant and to the importance of the Dead Sea Rift as a historical barrier geographically separating the new species from H. savignyi. Major genetic breaks revealed within species (H. felixarabica: Yemen vs. Jordan-Syria; H. savignyi sensu stricto: Levant vs. Turkey-Iran) are probably connected to climate changes during the Plio-Pleistocene boundary, while the finer phylogeographic structuring probably resulted from the Quaternary climate oscillations. The Cypriote population of H. savignyi originated from southern Anatolia relatively recently. Hyla orientalis from the southern Black Sea region seems to be genetically quite uniform, although two phylogeographic units with western Turkish and Caucasus-Caspian affinities might be detected. Hyla savignyi and H. orientalis carry signals of population expansions dated to the middle to late Pleistocene, while populations of H. felixarabica seem to have rather been constant in size, which might indicate more stable climatic conditions in the southern regions during the Quaternary.

Divergence in the face of gene flow: the case of two newts (amphibia: salamandridae)

Understanding the process of divergence requires the quantitative characterization of patterns of gene flow between diverging taxa. New and powerful coalescent-based methods give insight into these processes in unprecedented details by enabling the reconstruction of the temporal distribution of past gene flow. Here, we use sequence variation at eight nuclear markers and mitochondrial DNA (mtDNA) in multiple populations to study diversity, divergence, and gene flow between two subspecies of a salamander, the smooth newt (Lissotriton vulgaris kosswigi and Lissotriton vulgaris vulgaris) in Turkey. The ranges of both subspecies encompass mainly the areas of this important glacial refugial area. Populations in refugia where species have been present for a long time and differentiated in situ should better preserve the record of past gene flow than young populations in postglacial expansion areas. Sequence diversity in both subspecies was substantial (nuclear pi(sil) = 0.69% and 1.31%). We detected long-term demographic stability in these refugial populations with large effective population sizes (N(e)) of the order of 1.5-3 x 10(5) individuals. Gene trees and the isolation with migration (IM) analysis complemented by tests of nested IM models showed that despite deep, pre-Pleistocene divergence of the studied newts, asymmetric introgression from vulgaris to kosswigi has occurred, with signatures of recent gene flow in mtDNA and an anonymous nuclear marker, and evidence for more ancient introgression in nuclear introns. The distribution of migration times raises the intriguing possibility that even the initial divergence may have occurred in the face of gene flow.

Reduced genetic variation in the Japanese giant salamander, Andrias japonicus (Amphibia: Caudata)

The phylogenetic relationships among 46 samples from 27 populations of the Japanese giant salamander, Andriasjaponicus and its congener, A. davidianus from China was investigated, using 3664 bp sequences of the mitochondrial genes NADH1, NADH3, cyt b and CR, partial NADH6 and intervening genes. In phylogenetic trees constructed by MP, ML, and Bayesian methods, the family Cryptobranchidae and the genus Andrias both form monophyletic groups. Japanese A. japonicus and Chinese A. davidianus are sister taxa and can be regarded as separate species despite a small degree of genetic differentiation. Andriasjaponicus is divided into central and western clades, but the phylogenetic relationships within the latter clade are unresolved. As previously reported from allozyme analyses, A. japonicus exhibits little genetic differentiation, in strong contrast to salamanders of the genus Hynobius with which their distributions overlap. This reduced genetic variability in A. japonicus is attributable to a unique mating system of polygyny, delayed sexual maturity, notable longevity, life in a stable aquatic environment, and gigantism, as well as bottleneck effects following habitat fragmentation and extinction of local populations during Quaternary glaciations. The species is thus susceptible to extinction by potential environmental fluctuations, and requires extensive conservation measures.

Phylogeography of Sylvirana latouchii (Anura, Ranidae) in Taiwan

Biogeographic studies are important for understanding the natural history of faunas. To comprehend the geographical patterns of genetic variation in anurans in Taiwan, we investigated the genetic structure of Sylvirana latouchii (Anura, Ranidae) from 31 populations by using mitochondrial DNA (mtDNA) cytochrome b sequences. A neighbor-joining tree of 38 haplotypes revealed three major divergences in Taiwanese S. latouchii: the northern, western, and eastern-and-southern clades. Each clade was restricted to a single geographical district and showed obvious differentiation. The patterns of geographical divergence in this species reflect common historical events experienced by other native animals distributed in Taiwan. The order of divergence times between clades was inferred using a molecular clock test. The population relationship of S. latouchii between Taiwan and mainland China is discussed. Further study employing more populations of S. latouchii from mainland China is necessary to clarify the original geographical patterns and migratory history of this species.

Origin and post-colonization evolution of the Chatham Islands skink (Oligosoma nigriplantare nigriplantare)

Island ecosystems provide an opportunity to examine a range of evolutionary and ecological processes. The Chatham Islands are an isolated archipelago situated approximately 800 km east of New Zealand. Geological evidence indicates that the Chatham Islands re-emerged within the last 1-4 million years, following a prolonged period of marine inundation, and therefore the resident flora and fauna is the result of long-distance overwater dispersal. We examine the origin and post-colonization evolution of the Chatham Islands skink, Oligosoma nigriplantare nigriplantare, the sole reptile species occurring on the archipelago. We sampled O. n. nigriplantare from across nine islands within the Chatham Islands group, and representative samples from across the range of its closest relative, the New Zealand mainland common skink (Oligosoma nigriplantare polychroma). Our mitochondrial sequence data indicate that O. n. nigriplantare diverged from O. n. polychroma 5.86-7.29 million years ago. This pre-dates the emergence date for the Chatham Islands, but indicates that O. n. nigriplantare colonized the Chatham Islands via overwater dispersal on a single occasion. Despite the substantial morphological variability evident in O. n. nigriplantare, only relatively shallow genetic divergences (maximum divergence approximately 2%) were found across the Chatham Islands. Our analyses (haplotypic diversity, Phi(ST), analysis of molecular variance, and nested clade phylogeographical analysis) indicated restricted gene flow in O. n. nigriplantare resulting in strong differentiation between islands. However, the restrictions to gene flow might have only arisen recently as there was also a significant pattern of isolation by distance, possibly from when the Chatham Islands were a single landmass during Pleistocene glacial maxima when sea levels were lower. The level of genetic and morphological divergence between O. n. nigriplantare and O. n. polychroma might warrant their recognition as distinct species.

Phylogenetic relationships and heterogeneous evolutionary processes among phrynosomatine sand lizards (Squamata, Iguanidae) revisited

Phylogenetic analyses of DNA sequences were conducted to evaluate four alternative hypotheses of phrynosomatine sand lizard relationships. Sequences comprising 2871 aligned base pair positions representing the regions spanning ND1-COI and cyt b-tRNA(Thr) of the mitochondrial genome from all recognized sand lizard species were analyzed using unpartitioned parsimony and likelihood methods, likelihood methods with assumed partitions, Bayesian methods with assumed partitions, and Bayesian mixture models. The topology (Uma, (Callisaurus, (Cophosaurus, Holbrookia))) and thus monophyly of the "earless" taxa, Cophosaurus and Holbrookia, is supported by all analyses. Previously proposed topologies in which Uma and Callisaurus are sister taxa and those in which Holbrookia is the sister group to all other sand lizard taxa are rejected using both parsimony and likelihood-based significance tests with the combined, unparitioned data set. Bayesian hypothesis tests also reject those topologies using six assumed partitioning strategies, and the two partitioning strategies presumably associated with the most powerful tests also reject a third previously proposed topology, in which Callisaurus and Cophosaurus are sister taxa. For both maximum likelihood and Bayesian methods with assumed partitions, those partitions defined by codon position and tRNA stem and nonstems explained the data better than other strategies examined. Bayes factor estimates comparing results of assumed partitions versus mixture models suggest that mixture models perform better than assumed partitions when the latter were not based on functional characteristics of the data, such as codon position and tRNA stem and nonstems. However, assumed partitions performed better than mixture models when functional differences were incorporated. We reiterate the importance of accounting for heterogeneous evolutionary processes in the analysis of complex data sets and emphasize the importance of implementing mixed model likelihood methods.

Phylogeny and biogeography of the alpine newt Mesotriton alpestris (Salamandridae, Caudata), inferred from mtDNA sequences

In this paper, we performed phylogenetic analyses of Mesotriton alpestris populations from the entire range of species distribution, using fragments of two mtDNA genes, cytochrome b (309bp) and 16S rRNA ( approximately 500bp). Sequence diversity patterns and phylogenetic analyses reveal the existence of a relict lineage (Clade A) of late Miocene origin, comprising populations from south-eastern Serbia. This lineage is proposed to be ancestor to a western and an eastern lineage, which diverged during the middle Pliocene. The western lineage is further divided in two clades (Clades B, C) of middle Pliocene origin that represent populations from Italy (B) and populations from central Europe and Iberia (C). Further subdivision, dated back to the middle-late Pliocene, was found within the eastern lineage, representing southern (Clade D) and central-northern (Clade E) Balkan populations, respectively. Extensive sequence divergence, implying greater isolation in multiple refugia, is found within eastern clades, while the western clades seem to have been involved in the colonization of central, western and north-eastern Europe from a hypothetical refugium in central Europe. The extent of divergence does not support the current taxonomy indicating cryptic speciation in the Balkans, while paedomorphic lineages were found to have been evolved during early-middle Pleistocene probably as a response to the ongoing dramatic climatic oscillations.

Combining phylogeography with distribution modeling: multiple Pleistocene range expansions in a parthenogenetic gecko from the Australian arid zone

Phylogenetic and geographic evidence suggest that many parthenogenetic organisms have evolved recently and have spread rapidly. These patterns play a critical role in our understanding of the relative merits of sexual versus asexual reproductive modes, yet their interpretation is often hampered by a lack of detail. Here we present a detailed phylogeographic study of a vertebrate parthenogen, the Australian gecko Heteronotia binoei, in combination with statistical and biophysical modeling of its distribution during the last glacial maximum. Parthenogenetic H. binoei occur in the Australian arid zone and have the widest range of any known vertebrate parthenogen. They are broadly sympatric with their sexual counterparts, from which they arose via hybridization. We have applied nested clade phylogeographic, effective migration, and mismatch distribution analyses to mitochondrial DNA (mtDNA) sequences obtained for 319 individuals sampled throughout the known geographic ranges of two parthenogenetic mitochondrial lineages. These analyses provide strong evidence for past range expansion events from west to east across the arid zone, and for continuing eastward range expansion. Parthenogen formation and range expansion events date to the late Pleistocene, with one lineage expanding from the northwest of its present range around 240,000 years ago and the second lineage expanding from the far west around 70,000 years ago. Statistical and biophysical distribution models support these inferences of recent range expansion, with suitable climatic conditions during the last glacial maximum most likely limited to parts of the arid zone north and west of much of the current ranges of these lineages. Combination of phylogeographic analyses and distribution modeling allowed considerably stronger inferences of the history of this complex than either would in isolation, illustrating the power of combining complementary analytical approaches.

Multigenic and morphometric differentiation of ground squirrels (Spermophilus, Scuiridae, Rodentia) in Turkey, with a description of a new species

This study reports on the molecular phylogenetics of ground squirrels, genus Spermophilus, in Turkey using cytochrome b (1140bp), part of the D-loop and flanking tRNAs (572bp), X chromosome (867-1051bp) and Y chromosome (983-989bp) DNA sequences. Individuals also were characterized by karyotype and with geometric morphometric analyses of mandibles and skulls. Two hundred fourteen individuals from 91 localities were studied. All the data support the recognition of a new species in SW Anatolia: the Taurus ground squirrel Spermophilus taurensis sp. nov. The new species has a small distribution in the Taurus Mountains in an area that is a hotspot for biodiversity. Molecular clock analysis suggests that the new species diverged from the European ground squirrel, Spermophilus citellus, about 2.5 million years ago and that the ancestor of these two species diverged from the widespread Anatolian ground squirrel, Spermophilus xanthoprymnus, about 5 million years ago. Morphometric differentiation in skull and mandible shape among the three species is incomplete, but statistically significant. S. xanthoprymnus is subdivided into five cytochrome b phylogroups and we use these data to infer the location of glacial refugia where the species lived during the last glacial maximum. This study illustrates the potential of combined molecular and morphometric studies to uncover new Anatolian species and to reconstruct their phylogeographic history. The new species is important for squirrel taxonomy and for understanding Eurasian mammal evolution.

Lineage divergence of a freshwater snail clade associated with post-Tethys marine basin development

The complex evolutionary history of the Eurasian gastropod lineage Theodoxus reflects the evolution of marine basins following the breakup of the Tethys Sea. Today, this clade inhabits the lakes, rivers, streams, and estuaries of Europe, southwestern Asia, and North Africa. Here we present the first phylogenetic hypothesis for this clade. Based upon extensive geographic and taxonomic sampling, portions of the mitochondrial genes for cytochrome c oxidase subunit I and 16S rRNA were sequenced and analysed using maximum likelihood and maximum parsimony methods. Results from bootstrap analyses, Bayesian analysis, and sensitivity analyses lend support to six deep phylogenetic subdivisions within Theodoxus. These major clades are geographically associated with the major post-Tethyan marine basins. Estimates of divergence times using a penalized likelihood approach indicate that divergence of these major lineages occurred during the Miocene, simultaneous with the breakup of the Mediterranean and Paratethys Seas. The resulting major subclades later diversified during the Pliocene, primarily within geographic regions associated with the eastern and western Mediterranean Sea, the Pannonian Basin, and the Black Sea, thus producing the extant species assemblages. Finally, these phylogenetic results imply that much of the current taxonomy is flawed, therefore we offer recommendations for revising the classification of Theodoxus species based on phylogenetic systematics.

A molecular assessment of phylogenetic relationships and lineage accumulation rates within the family Salamandridae (Amphibia, Caudata)

We examine phylogenetic relationships among salamanders of the family Salamandridae using approximately 2700 bases of new mtDNA sequence data (the tRNALeu, ND1, tRNAIle, tRNAGln, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and COI genes and the origin for light-strand replication) collected from 96 individuals representing 61 of the 66 recognized salamandrid species and outgroups. Phylogenetic analyses using maximum parsimony and Bayesian analysis are performed on the new data alone and combined with previously reported sequences from other parts of the mitochondrial genome. The basal phylogenetic split is a polytomy of lineages ancestral to (1) the Italian newt Salamandrina terdigitata, (2) a strongly supported clade comprising the "true" salamanders (genera Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), and (3) a strongly supported clade comprising all newts except S. terdigitata. Strongly supported clades within the true salamanders include monophyly of each genus and grouping Chioglossa and Mertensiella as the sister taxon to a clade comprising Lyciasalamandra and Salamandra. Among newts, genera Echinotriton, Pleurodeles, and Tylototriton form a strongly supported clade whose sister taxon comprises the genera Calotriton, Cynops, Euproctus, Neurergus, Notophthalmus, Pachytriton, Paramesotriton, Taricha, and Triturus. Our results strongly support monophyly of all polytypic newt genera except Paramesotriton and Triturus, which appear paraphyletic, and Calotriton, for which only one of the two species is sampled. Other well-supported clades within newts include (1) Asian genera Cynops, Pachytriton, and Paramesotriton, (2) North American genera Notophthalmus and Taricha, (3) the Triturus vulgaris species group, and (4) the Triturus cristatus species group; some additional groupings appear strong in Bayesian but not parsimony analyses. Rates of lineage accumulation through time are evaluated using this nearly comprehensive sampling of salamandrid species-level lineages. Rate of lineage accumulation appears constant throughout salamandrid evolutionary history with no obvious fluctuations associated with origins of morphological or ecological novelties.