Mitochondrial sequence analysis of Salamandra taxa suggests old splits of major lineages and postglacial recolonizations of central Europe from distinct source populations of Salamandra salamandra

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.

Structure and evolution of the mitochondrial control region and flanking sequences in the European cave salamander Proteus anguinus

The European cave salamander Proteus anguinus Laurenti 1768 is one of the best-known subterranean animals, yet its evolutionary history and systematic relationships remain enigmatic. This is the first comprehensive study on molecular evolution within the taxon, using an mtDNA segment containing the control region (CR) and adjacent sequences. Two to seven tandem repeats of 24-32 bp were found in the intergenic spacer region (VNTR1), and three, four or six repeats, 59-77 bp each, in the 3' end of the CR (VNTR2). Different molecular mechanisms account for VNTR2 formation in different lineages of Proteus. The overall CR variation was lower than that of the spacer region, the 3' end of the cytb gene, or the tRNA genes. Individual genes and the concatenated non-repetitive sequences produced similar, well resolved maximum likelihood, Bayesian inference and parsimony trees. The numbers of repeat elements as well as the genealogy of the VNTR2 repeat units were mostly inconsistent with the groupings of the non-repetitive sequences. Different degrees of repeat array homogenization were detected in all major groups. Orthology was established for the first and the second VNTR2 elements of some populations. These two copies may therefore be used for analyses at the population level. The pattern of CR sequence variation points to strong genetic isolation of hydrographically separated populations. Genetic separation of the major groups of populations is incongruent with the current division into subspecies.

Distortion of symmetrical introgression in a hybrid zone: evidence for locus-specific selection and uni-directional range expansion

The fate of species integrity upon natural hybridization depends on the interaction between selection and dispersal. The relative significance of these processes may be studied in the initial phase of contact before selection and gene flow reach equilibrium. Here we study a hybrid zone of two salamander species, Lyciasalamandra antalyana and Lyciasalamandra billae, at the initial phase of hybridization. We quantify the degree and mode of introgression using nuclear and mtDNA markers. The hybrid zone can be characterized as an abrupt transition zone, the central hybrid zone being only c. 400 m, but introgressed genes were traced up to 3 km. Introgression was traced in both sexes but gene flow may be slightly male-biased. Indirect evidence suggests that hybrid males are less viable than females. Introgression occurred at two levels: (1) locus-specific selection led to different allelic introgression patterns independent of species, while (2) asymmetrical species-level introgression occurred predominately from L. antalyana to L. billae due to range expansion of the former. This indicates that foreign genes can be incorporated into novel genomic environments, which in turn may contribute to the great diversity of morphological variants in Lyciasalamandra.

Genetic differentiation of the marbled white butterfly, Melanargia galathea, accounts for glacial distribution patterns and postglacial range expansion in southeastern Europe

Isolation of Mediterranean species in the southern European peninsulas during the cold glacial phases often resulted in differentiation of several genetic lineages confined to the respective peninsulas. However, whilst there is good genetic evidence for multiple refugia in Iberia, there are only limited data available for the Balkans. Therefore, we wish to examine the hypothesis of a strong genetic structuring within southeastern Europe for the existence of multiple Balkan differentiation centres and/or several leading edges. As a model we use the marbled white butterfly, Melanargia galathea. We studied 18 allozyme loci of 564 individuals from 16 populations distributed over a large part of southeastern Europe. The single populations showed moderately high genetic diversity and no northward decline of genetic diversity was detected. The overall genetic differentiation between populations was considerable (F(ST) 7.0%). Cluster analysis discriminated three genetic groups: (i) a western flank in the former Yugoslavia, parts of eastern Austria and Hungary; (ii) an eastern flank with populations from Bulgaria and Romania (south of the southern Carpathians and eastern Carpathians); and (iii) the eastern Carpathian Basin. Hierarchical variance analysis distributed 53% of the variance among populations between these three groups. One sample from the Greek-Bulgarian border clustered within the eastern flank, but showed some tendency towards the eastern Carpathian Basin populations. Two populations from Carinthia clustered together with the eastern Carpathian Basin ones and a population from Styria showed an intermediate genetic composition between the three groups. Most probably, the eastern and the western flank groups are due to postglacial range expansion from the northeastern and the northwestern edges of the glacial differentiation centre (so-called leading edges). The eastern Carpathian Basin group may have resulted from postglacial expansion from northern Greece through valley systems of the central Balkan peninsula, maybe even expanding westwards north of the Balkan mountains reaching some parts of eastern Austria (e.g. Carinthia). Therefore, the Balkanic refugium of M. galathea may or may not have been continuous along the coastal areas of the Mediterranean, but must have been strongly genetically structured.

Mitochondrial evidence for distinct phylogeographic units in the endangered Malagasy poison frog Mantella bernhardi

Mantella bernhardi is an endemic species of Malagasy poison frog threatened by loss and fragmentation of its natural habitat and collection for the pet trade. It is classified as threatened according to the International Union for Conservation of Nature and Natural Resources (IUCN) categories and included in Appendix II of the Convention on the International Trade of Endangered Species (CITES). A recent survey has increased the known distributional range of the species from one to eight populations across southeastern Madagascar, but little is known about its biology and genetic diversity. Here we estimate inter- and intrapopulation mitochondrial genetic variation of four populations. Populations from the northern and southern parts of the distributional range showed a high degree of divergence (maximum of 11.35% in cytochrome b) and were recovered as reciprocally monophyletic groups. Nine haplotypes were detected in the northern and 12 in the southern populations. The population from Ranomafana National Park showed the lowest number of haplotypes and nucleotide diversity, and shared its most common haplotype with the second northern population from Tolongoina. All the other detected haplotypes were unique to each of the four populations. This suggests the existence of important barriers to gene flow, pre-dating human colonization of Madagascar at about 2000 years ago, in distinct contrast to other Mantella species that show a high degree of haplotype sharing throughout their range. The continued habitat fragmentation within the distribution range of M. bernhardi prevents any connection between its populations. Our data indicate the existence of at least two different management units for conservation in this species, corresponding to the North and South of its distribution range, and highlight the existence of strong regional endemism in southeastern Madagascar.

Evolutionary history of coastal tiger beetles in Japan based on a comparative phylogeography of four species

To reveal the phylogeographical patterns of four species of coastal tiger beetles in Japan (Lophyridia angulata, Abroscelis anchoralis, Cicindela lewisii and Chaetodera laetescripta), we conducted phylogenetic and nested clade analysis (NCA) using the mitochondrial DNA sequences of two loci (COI and 16S rRNA), with specimens sampled from Japan and neighbouring countries. Abroscelis anchoralis and L. angulata have similar disjunct distributions in Japan. The NCA indicated past fragmentation involving three isolated areas of A. anchoralis. In contrast, local populations of L. angulata in Japan shared the same haplotype, indicating recent vicariance. Co-occurrence of haplotypes from several divergent clades in Japanese populations of Ch. laetescripta suggested ancient vicariance and subsequent intermixing of local populations. The tree topology of C. lewisii, with shallow branches and little geographical segregation of haplotypes between Japan and Korea or within Japan, suggested that the Japanese population was segregated from the Korean population only recently. Restricted gene flow, with isolation by distance, was inferred for various geographical associations of haplotypes for coastal tiger beetles in the NCA. Based on these phylogeographical patterns, coupled with a molecular clock approach, the evolutionary history of four species of coastal tiger beetles was deduced, with the additional consideration of the competitive relationships among those species. We also discuss the conservation of highly localized A. anchoralis populations in Japan, using the concept of evolutionarily significant units.

The structure of biodiversity - insights from molecular phylogeography

DNA techniques, analytical methods and palaeoclimatic studies are greatly advancing our knowledge of the global distribution of genetic diversity, and how it evolved. Such phylogeographic studies are reviewed from Arctic, Temperate and Tropical regions, seeking commonalities of cause in the resulting genetic patterns. The genetic diversity is differently patterned within and among regions and biomes, and is related to their histories of climatic changes. This has major implications for conservation science.

Genetic variation in an endemic salamander, Salamandra atra, using amplified fragment length polymorphism

The pattern of genetic differentiation of the endemic alpine salamander, Salamandra atra, has been studied using amplified fragment length polymorphism (AFLP) from 11 populations throughout the range of the two currently recognized subspecies, atra and aurorae. Five different primer combinations produced 706 bands and were analyzed by constructing a phylogenetic tree using NJ and principal component analysis. Significant genetic variation was revealed by AFLP between and within populations but, our results show a lack of genetic structure. AFLP markers seems to be unsuitable to investigate complex and recent diversification.

Genetic consequences of climatic oscillations in the Quaternary

An appreciation of the scale and frequency of climatic oscillations in the past few million years is modifying our views on how evolution proceeds. Such major events caused extinction and repeated changes in the ranges of those taxa that survived. Their spatial effects depend on latitude and topography, with extensive extinction and recolonization in higher latitudes and altitudinal shifts and complex refugia nearer the tropics. The associated population dynamics varied with life history and geography, and the present genetic constitution of the populations and species carry attenuated signals of these past dynamics. Phylogeographic studies with DNA have burgeoned recently and studies are reviewed from the arctic, temperate and tropical regions, seeking commonalities of cause in the resulting genetic patterns. Arctic species show distinct shallow genetic clades with common geographical boundaries. Thus Beringia is distinct phylogeographically, but its role as a refugial source is complex. Arctic taxa do not show the common genetic pattern of southern richness and northern purity in north-temperate species. Temperate refugial regions in Europe and North America show relatively deep DNA divergence for many taxa, indicating their presence over several Ice Ages, and suggesting a mode of speciation by repeated allopatry. DNA evidence indicates temperate species in Europe had different patterns of postglacial colonization across the same area and different ones in previous oscillations, whereas the northwest region of North America was colonized from the north, east and south. Tropical montane regions contain deeply diverged lineages, often in a relatively small geographical area, suggesting their survival there from the Pliocene. Our poor understanding of refugial biodiversity would benefit from further combined fossil and genetic studies.

Adaptive divergence vs. environmental plasticity: tracing local genetic adaptation of metamorphosis traits in salamanders

In order to assess the significance of local adaptation relative to environmental plasticity on the evolution of life history traits, we analysed the possible genetic basis of differences between pond- and stream-breeding fire salamanders (Salamandra salamandra) in Germany. These salamanders typically deposit their larvae in small streams, where they grow until they are sufficiently large to metamorphose. However, some populations in Western Germany use ponds as larval habitat. Because habitat quality of streams differs from that of ponds one expects life history differences in the pond animals, which may result either from a plastic response or through genetic differentiation (i.e. local adaptation). Using a phylogeographical analysis of mitochondrial D-loop sequences, we show that both stream and pond populations in Western Germany are derived from a single lineage that recolonized following the last glaciation. This finding suggests that pond breeding originated very recently. Our studies of habitat quality and metamorphic behaviour of larvae in natural ponds and streams disclosed that pond larvae experience a significantly reduced food supply and greater risk of drying than do stream larvae. Pond larvae metamorphose earlier at the cost of reduced mass. Common-environment experiments with pond and stream larvae show that metamorphic behaviour of pond larvae under limited-food conditions is determined genetically and is not simply a plastic response to the differing habitat conditions. These results show that phenotypic plasticity is less important than local adaptation in explaining differences in ecological diversification within this species and suggests the possibility of rapid evolution of genetic adaptations when new habitats are exploited.

Historical isolation, range expansion, and secondary contact of two highly divergent mitochondrial lineages in spotted salamanders (Ambystoma maculatum)

The high species diversity of aquatic and terrestrial faunas in eastern North America has been attributed to range reductions and allopatric diversification resulting from historical climate change. The role these processes may have played in speciation is still a matter of considerable debate; however, their impacts on intraspecific genetic structure have been well documented. We use mitochondrial DNA sequences to reconstruct an intraspecific phylogeny of the widespread North American spotted salamander, Ambystoma maculatum, and test whether phylogenetic patterns conform to regional biogeographical hypotheses about the origins of diversity in eastern North America. Specifically, we address the number and locations of historical refugia, the extent and patterns of postglacial colonization by divergent lineages, and the origin and affinities of populations in the Interior Highland region. Despite apparent morphological uniformity, genetic discontinuities throughout the range of this species suggest that populations were historically fragmented in at least two refugia in the southern Appalachian Mountains. The ranges of these two highly divergent clades expanded northward, resulting in two widely distributed lineages that are sympatric in regions previously proposed as suture zones for other taxa. The evolutionary history of spotted salamander populations underscores the generality of biogeographical processes in eastern North America: despite differences in population size, glacial refugia, and vagility, similar signatures of differentiation are evident among and within widespread taxa.

Mitochondrial DNA phylogeography and population history of Meladema diving beetles on the Atlantic Islands and in the Mediterranean basin (Coleoptera, Dytiscidae)

The phylogeny and population history of Meladema diving beetles (Coleoptera, Dytiscidae) were examined using mitochondrial DNA sequence from 16S ribosomal RNA and cytochrome oxidase I genes in 51 individuals from 22 populations of the three extant species (M. imbricata endemic to the western Canary Islands, M. lanio endemic to Madeira and M. coriacea widespread in the Western Mediterranean and on the western Canaries), using a combination of phylogenetic and nested clade analyses. Four main lineages are observed within Meladema, representing the three recognized species plus Corsican populations of M. coriacea. Phylogenetic analyses demonstrate the sister relationship of the two Atlantic Island taxa, and suggest the possible paraphyly of M. coriacea. A molecular clock approach reveals that speciation within the genus occurred in the Early Pleistocene, indicating that the Atlantic Island endemics are not Tertiary relict taxa as had been proposed previously. Our results point to past population bottlenecks in all four lineages, with recent (Late-Middle Pleistocene) range expansion in non-Corsican M. coriacea and M. imbricata. Within the Canary Islands, M. imbricata seems to have independently colonized La Gomera and La Palma from Tenerife (although a colonization of La Palma from La Gomera cannot be discarded), and M. coriacea has independently colonized Tenerife and Gran Canaria from separate mainland lineages. In the Mediterranean basin this species apparently colonized Corsica on a single occasion, relatively early in its evolutionary history (Early Pleistocene), and has colonized Mallorca recently on multiple occasions. On the only island where M. coriacea and M. imbricata are broadly sympatric (Tenerife), we report evidence of bidirectional hybridization between the two species.

Dispersal of viviparity across contact zones in Iberian populations of fire salamanders (Salamandra) inferred from discordance of genetic and morphological traits

We used partial sequences of the cytochrome b mitochondrial DNA (mtDNA) gene, obtained from 76 individuals representing 45 populations of Iberian Salamandra salamandra plus 15 sequences of additional species of Salamandra and related genera, to investigate contact zones. These zones, identified by earlier allozymic and morphological analyses, are between populations of viviparous (S. s. bernardezi and S. s. fastuosa) and ovoviviparous (S. s. gallaica and S. s. terrestris) salamanders. The distribution of mtDNA and nuclear markers is mostly concordant at one contact zone (between S. s. gallaica and S. s. bernardezi), but at another (between S. s. fastuosa and S. s. terrestris) the markers are offset by about 250 km. The observed geographic variation fits a model of mtDNA capture. Among the potential mechanisms responsible for such discordance we favor a combination of range shifts due to climatic fluctuations and biased genetic admixture across moving contact zones. We apply our findings to the issue of possible homoplasy in the evolution of viviparity and conclude that viviparity likely arose only once within S. salamandra.

Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus, based on complete sequencing of the mitochondrial DNA control region

The complete mitochondrial DNA (mtDNA) control region (1043 base pairs) and 162-bp of flanking transfer RNA genes were sequenced in 316 European grayling, Thymallus thymallus, from 44 populations throughout the Western European range of the species. A total of 58 haplotypes were revealed with pairwise divergence ranging from 0.001 to 0.038. An inferred intraspecific phylogenetic tree revealed two well-supported clades within the Danube basin, one highly divergent clade in the Adriatic basin, and one large, diverse group representing most other populations. A deeply divergent haplotype fixed in the Loire basin in central France, more groups of haplotypes from distinct Danubian tributaries, and a relatively ancestral haplotype fixed in former tributaries of the Elbe in Denmark all suggest a complex pattern of interglacial and postglacial expansions originating from disjunct refugia throughout central Europe. Despite some evidence of human-mediated stock transfers, parsimony-network-based nested-clade analysis (NCA) supported specific inferences relating to corridors of postglacial expansion such as the lower Rhine (Moselle) and Elbe systems (Danish populations) serving as sources for expansion into the Baltic to the north as well as the upper Rhine and Danube to the south; and specific Rhine populations (Doller, Orbe and Reuss) serving as sources for colonization of the Rhone. The multiple divergent clades representing populations in the upper Danube, as well as the deeply divergent haplotypes found in the Adriatic and Loire basins (> 5% divergence from Asian outgroups) support the theory that European grayling have had a long history in Western Europe, pre-dating Pleistocene glacial cycles. The patterns of mtDNA divergence shown here support a perspective of rich inter- and intrabasin genetic diversity that should be protected from current trends to translocate brood stocks for rearing and release in response to declining populations, especially in southern European basins.

Molecular systematics of the western rattlesnake, Crotalus viridis (Viperidae), with comments on the utility of the D-loop in phylogenetic studies of snakes

Crotalus viridis, the western rattlesnake, ranges throughout western North America and has been divided into at least eight subspecies. However, the validity of and relationships among these subspecies and the monophyly of C. viridis as a whole are questionable. We used mitochondrial DNA sequence data from the D-loop region and ND2 gene to examine the relationships among 26 populations of C. viridis and to test the monophyly of this species. These data were analyzed separately and combined using maximum-likelihood and maximum-parsimony. The C. viridis group was monophyletic in all combined analyses, consisting of two strongly divergent clades. We recommend that these clades be recognized as two distinct evolutionary species: C. viridis and C. oreganus. Crotalus viridis should be restricted to the subspecies viridis and nuntius and the remaining subspecies be assigned to the species C. oreganus. Our data do not allow strong evaluation of the validity of the subspecies. We found that the ND2 gene had greater sequence divergences among closely related individuals than the D-loop region, but this relationship reversed at higher levels of divergence. This pattern is apparently due to: (1) ND2 third positions evolving faster than the D-loop but becoming saturated at higher levels of divergence, and (2) the D-loop evolving faster than ND2 second (and possibly first) positions. Our results suggest that the ND2 gene is preferable for examining intraspecific relationships and the D-loop may better resolve relationships between species of snakes. The latter result is contrary to the common perception of the phylogenetic utility of the D-loop. Another unusual result is that the 145 bp spacer region, adjacent to the 5' end of the light strand of the D-loop, provides greater phylogenetic resolution than the 1030 bp D-loop.

Phylogeography of the Alpine salamander, Salamandra atra (Salamandridae) and the influence of the Pleistocene climatic oscillations on population divergence

Fifty individuals of the endemic Alpine salamander, Salamandra atra, representing 13 populations throughout the range of the two currently recognized subspecies, atra and aurorae, were examined for sequence variation in a large portion (1050 bp) of the mitochondrial cytochrome b gene. We revealed a large number of mitochondrial DNA (mtDNA) haplotypes (10). Interpopulation sequence divergence was very low, ranging from 0 to 3.1%. The relationships among haplotypes were poorly resolved. The divergence time estimate between several mtDNA haplotypes suggested a pre-Pleistocene differentiation approximately 3 million years ago. Moreover, the impact of the Pleistocene glaciations on the phylogeographical patterns appears to have been secondary, although a somewhat reduced genetic variability was found in populations living in areas that were directly affected by the glaciation.

The persistence of Pliocene populations through the Pleistocene climatic cycles: evidence from the phylogeography of an Iberian lizard

Ancient climatic fluctuations have caused changes in the demography and distribution of many species. The genetic differentiation between populations of the same species and of sister species is often attributed largely to the more recent Pleistocene fluctuations. Recent interpretations, which implicate earlier episodes, have proved controversial. We address the timing of genetic divergence in the Iberian lizard Lacerta schreiberi by studying the phylogeography of the cytochrome b sequence. The species has a remarkable morphological uniformity, yet our evidence suggests that earlier events in the Pliocene initiated the main divergence between populations. This interpretation implies that the different populations survived through the Pleistocene in separate localities. This conclusion is robust to different molecular clock calibrations. The persistence of earlier differentiation through the Pleistocene has wide implications for our understanding of Pleistocene refugia in this species and, by extension, to the biogeography of the whole region.

Genetic and ecological divergence of a monophyletic cichlid species pair under fully sympatric conditions in Lake Ejagham, Cameroon

Although there is mounting evidence that speciation can occur under sympatric conditions, unambiguous examples from nature are rare and it is almost always possible to propose alternative allopatric or parapatric scenarios. To identify an unequivocal case of sympatric speciation it is, therefore, necessary to analyse natural settings where recent monophyletic species flocks have evolved within a small and confined spatial range. We have studied such a case with a cichlid species flock that comprises five Tilapia forms endemic to a tiny lake (Lake Ejagham with a surface area of approximately 0.49 km2) in Western Cameroon. Analysis of mitochondrial D-Loop sequences shows that the flock is very young (approximately 10(4) years) and has originated from an adjacent riverine founder population. We have focused our study on a particular pair of forms within the lake that currently appears to be in the process of speciation. This pair is characterized by an unique breeding colouration and specific morphological aspects, which can serve as synapomorphic characters to prove monophyly. It has differentiated into a large inshore and a small pelagic form, apparently as a response to differential utilization of food resources. Still, breeding and brood care occurs in overlapping areas, both in time and space. Analysis of nuclear gene flow on the basis of microsatellite polymorphisms shows a highly restricted gene flow between the forms, suggesting reproductive isolation between them. This reproductive isolation is apparently achieved by size assortative mating, although occasional mixed pairs can be observed. Our findings are congruent with recent theoretical models for sympatric speciation, which show that differential ecological adaptations in combination with assortative mating could easily lead to speciation in sympatry.

Speciation, hybrid zones and phylogeography - or seeing genes in space and time

The origins and development of the study of speciation, hybrid zones and phylogeography are outlined using evolutionary iconography. This traces the ideas in this field from Lamarck and Darwin through to the present as represented in diagrams and figures. A 'tree of trees' summarizes this growth and current vitality. The new facility to use various DNA sequences from nuclear, mitochondrial and chloroplast genomes to determine genetic variation throughout a species range is examined particularly. There is great genomic subdivision across species distributions, which can be interpreted in the light of the recent demonstrations of severe palaeoclimatic oscillations. Refugia and postglacial colonization routes are proposed for several organisms across Europe. The role of geography in speciation through the Pleistocene is considered. These emerging principles and analyses are applied to data available on a variety of organisms in other regions of the world, such as the Arctic, North America and the Tropics, and including the progress of Homo sapiens through the last ice age. Some suggestions are made for future research directions.

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

Phylogenetic relationships among salamandrids of the "true" salamander clade are investigated using 2019 aligned base positions (713 parsimony informative) of 20 mitochondrial DNA sequences from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and the tRNA(Cys) genes. Parsimony analysis produces a robust phylogenetic estimate for the relationships of the major groups of "true" salamanders. Strong support is provided for the sister taxon relationship of Chioglossa and Mertensiella caucasica and for the placement of Salamandra and Mertensiella luschani as sister taxa. These relationships suggest two vicariant events between Europe and Anatolia caused by the formation of seaways in the Mediterranean Basin. Molecular divergence indicates an Early Miocene separation of Chioglossa and M. caucasica and a Late Miocene separation of Salamandra and M. luschani. The traditional phylogenetic hypothesis of a monophyletic Mertensiella is statistically rejected, indicating that southwestern and northeastern Anatolian populations have separate historical biogeographic origins. Therefore, we recommend placement of M. luschani in the genus Salamandra. Within M. luschani, six highly divergent lineages showing 7.6 to 10.1% pairwise sequence divergence are identified. Tests using four-taxon subsamples suggest that these lineages diverged nearly simultaneously in the Late Miocene, approximately 6 to 8 million years ago, when extensive uplifting of Anatolia occurred in response to the Arabian collision.

Phylogeography of the bullhead Cottus gobio (Pisces: Teleostei: Cottidae) suggests a pre-pleistocene origin of the major central European populations

The bullhead Cottus gobio is a small, bottom-dwelling fish consisting of populations that have not been subject to transplantations or artificial stocking. It is therefore an ideal model species for studying the colonization history of central European freshwater systems, in particular with respect to the possible influences of the Pleistocene glaciation cycles. We sampled Cottus populations across most of its distribution range, with a special emphasis on southern Germany where the major European drainage systems are in closest contact. Mitochondrial D-loop sequencing of more than 400 specimens and phylogenetic network analysis allowed us to draw a detailed picture of the colonization of Europe by C. gobio. Moreover, the molecular distances between the haplotypes enabled us to infer an approximate time frame for the origin of the various populations. The founder population of C. gobio stems apparently from the Paratethys and invaded Europe in the Pliocene. From there, the first colonization into central Europe occurred via the ancient lower Danube, with a separate colonization of the eastern European territories. During the late Pliocene, one of the central European populations must have reached the North Sea in a second step after which it then started to colonize the Atlantic drainages via coastal lines. Accordingly, we found very distinct populations in the upper and lower Rhine, which can be explained by the fact that the lower Rhine was disconnected from the upper Rhine until approximately 1 million years ago (Ma). More closely related, but still distinct, populations were found in the Elbe, the Main and the upper Danube, all presumably of Pleistocene origin. Intriguingly, they have largely maintained their population identity, despite the strong disturbance caused by the glaciation cycles in these areas. On the other hand, a mixing of populations during postglacial recolonization could be detected in the lower Rhine and its tributaries. However, the general pattern that emerges from our analysis suggests that the glaciation cycles did not have a major impact on the general population structure of C. gobio in central Europe.