Phylogenomic insights into the diversity and evolution of Palearctic vipers.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of âˆ¼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

A new set of microsatellite markers for Phoxinus lumaireul senso lato, Phoxinus marsilii and Phoxinus krkae for population and molecular taxonomic studies.

Minnows of the genus Phoxinus are common and an often highly abundant fish species in Palearctic freshwater habitats. Phoxinus species have a complex evolutionary history, phylogenetic relationships are not well understood and there are a number of unresolved taxonomic problems. There are currently 23 different mitochondrial genetic lineages identified in the genus Phoxinus, 13 of which are recognized as valid species. The taxonomic status of these lineages requires resolution, including the degree to which they can interbreed. Suitable nuclear molecular markers for studies of population divergence and interbreeding between morphotypes and mitochondrial lineages are lacking for Phoxinus species. Therefore, the authors developed a set of microsatellite markers using genomic information from Phoxinus lumaireul and tested their suitability for this and two related species, Phoxinus krkae and Phoxinus marsilii. Out of 16 microsatellite candidate loci isolated, 12 were found to be in Hardy-Weinberg equilibrium when tested on two P. lumaireul senso lato populations. Seven loci amplified across the three species, enabling the study of intraspecific genetic diversity and population structure within P. marsilii and P. krkae. The markers were able to clearly resolve differences among the three tested species, including the recently described P. krkae, and are therefore suitable for the detection of introgression and hybridization among populations consisting of mixtures of two or more of P. lumaireul s. l., P. marsilii and P. krkae.

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders.

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Description of a new species of Phoxinus from the upper Krka River (Adriatic Basin) in Croatia (Actinopterygii: Leuciscidae), first discovered as a molecular clade.

Using morphometric, meristic and qualitative characters, we confirmed a specific status of a clade discovered based on molecular data only. The newly described species Phoxinus krkae has a very local distribution in the upper Krka River but is easily morphologically distinguished from geographically neighbouring species with much wider ranges: Phoxinus lumaireul, Phoxinus marsilii and Phoxinus csikii. Phoxinus krkae is well differentiated in statistical analyses and differs from other species of Phoxinus in northern Adriatic Sea drainages and right-bank tributaries of the middle and lower River Danube by a combination of characters, none of them unique: caudal peduncle length 1.9-2.2 times minimum depth of caudal peduncle; upper lip not projecting beyond lower lip; rostral fold absent; total number of scales in lateral series 74-84; lateral line incomplete and interrupted, commonly 17-49 (< 60% of total number of scales in lateral series) with no or few canal segments on posterior body and caudal peduncle; scales on belly extending forward to over middle of distance between pelvic-fin base and pectoral-fin origin; patches of breast scales commonly connected by scattered scales or forming 1-3 regular rows of scales; total vertebrae commonly 38-39 (22 + 17, 22 + 16 or 21 + 17). The study reveals unrecognised morphological diversity within the group and provides the groundwork for future taxonomic study in Phoxinus in the Adriatic region and the Danube Basin.

Phenotypic diversity in an endangered freshwater fish Squalius microlepis (Actinopterygii, Leuciscidae).

Squalius microlepis was examined from recent and historical collections within the known range of the species with special emphasis on intraspecific variability and variations, and compared to its closest relative species S. tenellus (in total, 193 specimens; 33 absolute and 52 proportional measurements and ratios, and 12 counts including vertebrae). Squalius tenellus was perfectly differentiated in all statistical analyses and can be diagnosed by 76-95 (vs. 64-80) scales in lateral series, 68-83 (vs. 58-77) lateral-line scales, (17)18-20 (vs. 13-16(17)) scales above lateral line, and (7)8-10 (vs. 4-7) scales below lateral line. Squalius microlepis was morphologically heterogeneous, with two phenotypes readily distinguishable (phenotype 1 corresponding to S. microlepis s. str. as defined by its lectotype) by a combination of many characters; those contributing most to the discrimination were number of gill rakers, length of lower jaw (% interorbital width), and head length (% SL). Only phenotype 1 was found in the Ricina-Prolosko Blato-Vrljika karst system; most of the specimens from the lower Matica and the Tihaljina-Trebizat karst system were identified as phenotype 2; the sample from karstic poljes near Vrgorac contained both phenotype 1 and 2, and individuals of intermediate morphology. As very limited molecular data exist on the two phenotypes of S. microlepis, we refrain from any taxonomic conclusions until new molecular approaches (and new markers) are used. We also report on a dramatic reduction of the area of distribution and abundance of S. microlepis in recent years.

Foraging mode, relative prey size and diet breadth: A phylogenetically explicit analysis of snake feeding ecology.

Foraging modes (ambush vs. active foraging) are often correlated with a suite of morphological, physiological, behavioural and ecological traits known as the "adaptive syndrome" or "syndrome hypothesis." In snakes, an ecological correlate often reported in the literature is that ambush-hunting snakes have a higher relative meal size compared to actively foraging snakes which feed on smaller prey items. This "large meal versus small meal" feeding hypothesis between ambush and active foragers has become a widely accepted paradigm of snake feeding ecology, despite the fact that no rigorous meta-analysis has been conducted to support this generalization. We conducted a phylogenetically explicit meta-analysis, which included ca. 100 species, to test this paradigm of snake feeding ecology. We gathered data on prey size by inducing regurgitation by palpation in free-ranging snakes and by examining the stomach contents of preserved museum specimens. When we found prey, we recorded both snake and prey mass to estimate relative prey mass (prey mass/snake mass). We also reviewed published studies of snake feeding ecology to gather similar information for other species. Ambush and active foragers did not differ in minimum or average meal size but the maximum meal sizes consumed by ambush-foraging snakes were larger than the maximum meal sizes eaten by active foragers. This results in ambush-foraging snakes consuming a significantly wider range of meal sizes, rather than being large meal specialists compared to active foragers. We argue that ambush foragers evolved to be more opportunistic predators because they encounter prey less frequently compared to active foragers. This hypothesis is further supported by the fact that ambush foragers also exhibited marginally wider diet breadths, consuming a broader range of prey types in comparison with active foragers. Our study challenges aspects of the foraging syndrome as it is currently conceived, and our results have important implications for our understanding of how foraging mode has shaped the behaviour and physiology of ambush-foraging snakes.

Population Genetic Analyses Using 10 New Polymorphic Microsatellite Loci Confirms Genetic Subdivision within the Olm, Proteus anguinus.

We provide a comparative population genetic study of the elusive amphibian, Proteus anguinus, by comparing the genetic diversity and divergence among 4 cave populations (96 individuals) sampled in the Dinaric Karst of Croatia. We developed 10 variable microsatellite markers using pyrosequencing and applied them to the 4 selected populations belonging to 4 different cave systems. The results showed strong genetic differentiation between the 4 caves corroborating with previous findings suggesting that Proteus might comprise several unrecognized taxa. Our results confirmed that gene flow should be high within the caves, whereas it is low between hydrographic systems since geological periods. Finally, we conclude that the high genetic subdivision suggests the necessity of treating the 4 studied Proteus populations as evolutionary significant units.

Distribution of endangered Italian gudgeon Romanogobio benacensis (Cypriniformes, Cyprinidae, Gobioninae) with remarks on distinguishing morphological characters.

Distribution data on many freshwater fish species in Croatia are scarce and species identifications are difficult, requiring further detailed studies. This paper presents a report of the Italian gudgeon Romanogobio benacensis from the Mirna River in the Istra Peninsula in Croatia, in the south-east from its previously known distribution range. The identification of R. benacensis in Croatia was supported by a morphological comparison with R. benacensis from Italy and Slovenia, the common gudgeon Gobio gobio, and the Danubian gudgeon Gobio obtusirostris from geographically close locations. A combination of character states (number of scales between anus and anal-fin origin, branched pectoral-fin rays, lateral-line scales, total, abdominal, and caudal vertebrae, and the size and number of lateral blotches) distinguishes R. benacensis from both G. gobio and G. obtusirostris. The phylogenetic analyses using mitochondrial sequences of cytochrome b gene confirmed that specimens from the Mirna River belong to R. benacensis. Also, Reka River system (Adriatic Sea basin) in Slovenia is inhabited by a possibly introduced Danubian gudgeon, G. obtusirostris, and not by R. benacensis.

Monitoring health and reproductive status of olms (Proteus anguinus) by ultrasound.

The olm (Proteus anguinus) is a troglomorphic, neotenous amphibian with extraordinary life expectancy and unique adaptations that deserve further investigation. A low reproductive rate and habitat decline render it threatened by extinction. Establishing captive populations for maintenance and artificial breeding may one day become crucial to the species. Longitudinal, in-vivo assessment of inner organs is invaluable to our understanding of reproductive physiology, health, and behavior. Using ultrasound, we measured heart rate and assessed health and reproductive status of 13 captive olms at Zagreb Zoo. Heart rate averaged 42.9 ± 4.6 bpm (32-55 bpm), as determined via pulsed-wave Doppler at 4-12 MHz. By using frequencies of up to 70 MHz (ultrasound biomicroscopy), inner organs were visualized in detail. Assessment of the gastrointestinal tract provided insights into feeding status and digestive processes. Several subclinical pathologies were detected, including biliary sludge, subcutaneous edema, ascites, and skin lesions. Detection of skin lesions by ultrasound was more sensitive than visual adspection. Olms with ultrasonographically detected skin lesions tested positive for Saprolegnia and were treated. Three of the four affected individuals survived and subsequently tested negative for Saprolegnia. Sex was reliably determined; only one individual proved male. The reason for this extreme female-biased sex-ratio remains unknown. However, as most of the individuals were flushed from the caves by strong currents in spring, the sample may not be representative of natural populations. In female olms, different stages of ovarian follicular development were observed with diameters ranging between 0.1 and 1.1 mm. Results were confirmed by comparing ultrasound, necropsy, and histological findings of one dead specimen. In summary, ultrasound proved a valuable tool to support conservation and captive breeding programs by allowing non-invasive assessment of physiological parameters, clinical condition, and reproductive status in olms.

Surveying Europe's Only Cave-Dwelling Chordate Species (Proteus anguinus) Using Environmental DNA.

In surveillance of subterranean fauna, especially in the case of rare or elusive aquatic species, traditional techniques used for epigean species are often not feasible. We developed a non-invasive survey method based on environmental DNA (eDNA) to detect the presence of the red-listed cave-dwelling amphibian, Proteus anguinus, in the caves of the Dinaric Karst. We tested the method in fifteen caves in Croatia, from which the species was previously recorded or expected to occur. We successfully confirmed the presence of P. anguinus from ten caves and detected the species for the first time in five others. Using a hierarchical occupancy model we compared the availability and detection probability of eDNA of two water sampling methods, filtration and precipitation. The statistical analysis showed that both availability and detection probability depended on the method and estimates for both probabilities were higher using filter samples than for precipitation samples. Combining reliable field and laboratory methods with robust statistical modeling will give the best estimates of species occurrence.

Description of a new species of Alburnus Rafinesque, 1820 (Actinopterygii, Cyprinidae, Leuciscinae) from the Kolpa River in the Sava River system (upper Danube drainage), with remarks on the geographical distribution of shemayas in the Danube.

Alburnus sava, new species, is described from the Kolpa River. The Kolpa is a tributary of the Sava, a major tributary of the Danube River, in the Black Sea basin. Alburnus sava is distinguished from its congeners in the Danube drainage, A. mento and A. sarmaticus, by having 23-27, usually 24-26, gill rakers; the ventral keel usually completely covered by scales (scaleless part maximum 15% of the keel length); 15-16, mode 15, branched pectoral-fin rays; the length of the gill raker at the junction of the arch limbs 65-70% of the length of the opposite outer gill filament; and a relatively long lower jaw (37-40% HL or 112-130% interorbital width). Alburnus sava is a large-sized potamadromous shemaya known to occur in the entire Sava drainage. The taxonomic status of A. mento and A. sarmaticus is confirmed. Alburnus danubicus is discussed and as there are no new arguments, it is kept as a valid species. New details on the distribution of shemayas in the Danube drainage are presented.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula.

BACKGROUND: Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges. RESULTS: We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity. CONCLUSIONS: Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Variation in the cranium shape of wall lizards (Podarcis spp.): effects of phylogenetic constraints, allometric constraints and ecology.

We used geometric morphometrics to explore the influence of phylogenetic and allometric constraints as well as ecology on variation in cranium shape in five species of monophyletic, morphologically similar Podarcis lizards (Podarcis erhardii, Podarcis melisellensis, Podarcis muralis, Podarcis sicula and Podarcis taurica). These species belong to different clades, they differ in their habitat preferences and can be classified into two distinct morphotypes: saxicolous and terrestrial. We found (i) no phylogenetic signal in cranium shape, (ii) diverging allometric slopes among species, and (iii) a significant effect of habitat on cranium shape. The saxicolous species (P. erhardii and P. muralis) had crania with elongated parietals, elongated cranium bases, shortened anterior parts of the dorsal cranium, reduced chambers of the jaw adductor muscles and larger subocular foramina. These cranial features are adaptations that compensate for a flattened cranium, dwelling on vertical surfaces and seeking refuge in crevices. The crania of the terrestrial species (P. melisellensis, P. sicula and P. taurica) tended to be more elongate and robust, with enlarged chambers of the jaw adductor muscle, reduced skull bases and shortened parietals. Terrestrial species exhibited more variation in cranium shape than saxicolous species. Our study suggests that shape variation in Podarcis sp. lizards is largely influenced by ecology, which likely affects species-specific patterns of static allometry.

First clinical experiences about the neurotoxic envenomings inflicted by lowland populations of the Balkan adder, Vipera berus bosniensis.

The first overall clinical description of envenomings by the lowland populations of the Balkan adder (Vipera berus bosniensis) is provided by this study. Fifty-four incidents have been collected retrospectively from the south-western Hungarian and the northern Croatian distribution area of the taxon. There were five (9%) asymptomatic, 24 (44%) mild, 12 (22%) moderate, 12 (22%) severe, and one fatal (2%) case according to the Poisoning Severity Score. The single death is a 60-year-old Hungarian case that was caused by V.b. bosniensis. Average hospitalisation was 2.75 days. The most common systemic symptoms were gastrointestinal disorders, ECG changes, persisting hypotension and neurological disorders. The initial phase of neurotoxic manifestations was always expressed in cranial nerve disturbances: ptosis, external ophthalmoplegia, diplopia, reduced focusing capability and blurred vision. Neuromuscular paralysis progressed to dyspnoea and lower limb paralysis in the most severe cases. Unusual symptoms were fluctuating arterial hypertension, drowsiness, and hypokalaemia. Laboratory results reveal leucocytosis, while deviation in the other laboratory values is not common. Envenomings by V.b. bosniensis significantly differ from those by the European adder (Vipera berus berus) in lower manifestation rate of extensive oedema, anaemia, CNS depression, and haematuria but the development of neuromuscular paralysis is high (20%). Their bites rather resulted in mild and moderate local symptoms in envenomed patients than those inflicted by the nominate form. This study presents the evidence of the frequent neurotoxic manifestations in Balkan adder-bitten patients for the first time, which strongly suggests that the venom of the lowland populations of V.b. bosniensis has neurotoxic activity.