Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data

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.

Evolutionary history and systematics of European blind mole rats (Rodentia: Spalacidae: Nannospalax): Multilocus phylogeny and species delimitation in a puzzling group

Species delimitation is a powerful approach to assist taxonomic decisions in challenging taxa where species boundaries are hard to establish. European taxa of the blind mole rats (genus Nannospalax) display small morphological differences and complex chromosomal evolution at a shallow evolutionary divergence level. Previous analyses led to the recognition of 25 'forms' in their distribution area. We provide a comprehensive framework to improve knowledge on the evolutionary history and revise the taxonomy of European blind mole rats based on samples from all but three of the 25 forms. We sequenced two nuclear-encoded genetic regions and the whole mitochondrial cytochrome b gene for phylogenetic tree reconstructions using concatenation and coalescence-based species-tree estimations. The phylogenetic analyses confirmed that Aegean N. insularis belongs to N. superspecies xanthodon, and that it represents the second known species of this superspecies in Europe. Mainland taxa reached Europe from Asia Minor in two colonisation events corresponding to two superspecies-level taxa: N. superspecies monticola (taxon established herewith) reached Europe c. 2.1 million years ago (Mya) and was followed by N. superspecies leucodon (re-defined herewith) c. 1.5 Mya. Species delimitation allowed the clarification of the taxonomic contents of the above superspecies. N. superspecies monticola contains three species geographically confined to the western periphery of the distribution of blind mole rats, whereas N. superspecies leucodon is more speciose with six species and several additional subspecies. The observed geographic pattern hints at a robust peripatric speciation process and rapid chromosomal evolution. The present treatment is thus regarded as the minimum taxonomic content of each lineage, which can be further refined based on other sources of information such as karyological traits, crossbreeding experiments, etc. The species delimitation models also allowed the recognition of a hitherto unnamed blind mole rat taxon from Albania, described here as a new subspecies.

Are vipers prototypic fear-evoking snakes? A cross-cultural comparison of Somalis and Czechs

Snakes are known as highly fear-evoking animals, eliciting preferential attention and fast detection in humans. We examined the human fear response to snakes in the context of both current and evolutionary experiences, conducting our research in the cradle of humankind, the Horn of Africa. This region is characterized by the frequent occurrence of various snake species, including deadly venomous viperids (adders) and elapids (cobras and mambas). We conducted experiments in Somaliland and compared the results with data from Czech respondents to address the still unresolved questions: To which extent is human fear of snakes affected by evolutionary or current experience and local culture? Can people of both nationalities recognize venomous snakes as a category, or are they only afraid of certain species that are most dangerous in a given area? Are respondents of both nationalities equally afraid of deadly snakes from both families (Viperidae, Elapidae)? We employed a well-established picture-sorting approach, consisting of 48 snake species belonging to four distinct groups. Our results revealed significant agreement among Somali as well as Czech respondents. We found a highly significant effect of the stimulus on perceived fear in both populations. Vipers appeared to be the most salient stimuli in both populations, as they occupied the highest positions according to the reported level of subjectively perceived fear. The position of vipers strongly contrasts with the fear ranking of deadly venomous elapids, which were in lower positions. Fear scores of vipers were significantly higher in both populations, and their best predictor was the body width of the snake. The evolutionary, cultural, and cognitive aspects of this phenomenon are discussed.

Strong intraspecific phylogenetic and karyotypic diversification in Isophya modestior (Orthoptera: Tettigoniidae: Phaneropterinae)

Isophya modestior (Orthoptera: Tettigoniidae) is a species distributed in central and south-eastern Europe, where its distribution is largely separated by two large rivers (Sava and Danube). Since previous studies on the song and morphology of the stridulatory file across its complete geographic range showed that the species is separated into two main groups, we decided to use phylogenetic and karyological analyses in order to evaluate the status of the previously analysed populations.

Phylogenetic analyses showed the existence of two major clades within I. modestior with very high bootstrap values and posterior probabilities—Clade A: present on the Balkan Peninsula, Slovenia (Inner Carniola), Italy, Pannonian Serbia (Vršac Mts and Deronje) and Austria (Burgenland and Lower Austria); Clade B: present in Slovenia (Upper Carniola), Croatia and Austria (Carinthia), Pannonian Serbia (Fruška Gora Mt.) and Hungary.

A comparison of chromosomes of 51 specimens revealed discrete differences between their karyotypes. The physical characteristics of the karyotypes included chromosome number (2n), sex chromosome (X) morphology and C-banding patterns. The standard chromosome complement of 50 specimens from different localities is characterized by 2n = 30 + X0 in males. In one male collected in Fruška Gora Mt. (Andrevlje), the chromosome number was reduced to 2n = 28 + neo-XY. Therefore, further cytogenetic studies involving larger samples, especially from Fruška Gora Mt., are needed in order to gain a more comprehensive view of the chromosome evolution in this group of Isophya species.

Cytogenetic Analysis of the Bimodal Karyotype of the Common European Adder, Vipera berus (Viperidae)

Vipera berus is the species with the largest range of snakes on Earth and one of the largest among reptiles in general. It is also the only snake species found in the Arctic Circle. Vipera berus is the most involved species of the genus Vipera in the process of interspecific hybridization in nature. The taxonomy of the genus Vipera is based on molecular markers and morphology and requires clarification using SC-karyotyping. This work is a detailed comparative study of the somatic and meiotic karyotypes of V. berus, with special attention to DNA and protein markers associated with synaptonemal complexes. The karyotype of V. berus is a remarkable example of a bimodal karyotype containing both 16 large macrochromosomes and 20 microchromosomes. We traced the stages of the asynchronous assembly of both types of bivalents. The number of crossing-over sites per pachytene nucleus, the localization of the nucleolar organizer, and the unique heterochromatin block on the autosomal bivalent 6-an important marker-were determined. Our results show that the average number of crossing-over sites per pachytene nucleus is 49.5, and the number of MLH1 sites per bivalent 1 reached 11, which is comparable to several species of agamas.

Climate Change Is Driving Shifts in Dragonfly Species Richness across Europe via Differential Dynamics of Taxonomic and Biogeographic Groups

Understanding how changes in species richness pattern correlate with range changes in different taxonomic and biogeographic groups is important for conservation because it allows for generalizations about which species are at greatest risk. Here, we assessed whether changes in species richness patterns result from generalized range shifts across taxonomic and biogeographic groups or from changes in specific subsets of species. Using data from 1988 and from 2010, we studied changes in distributional range of European dragonfly species, using outline distribution maps for all dragonflies combined and separately for taxonomic suborders (Zygoptera and Anisoptera) and biogeographic groups (Boreo-alpine, Eurasian, Mediterranean, and Tropical). The results demonstrated differing range dynamics for Zygoptera and Anisoptera, with Anisoptera driving local turnover in species richness to a greater extent than Zygoptera. The distributional range of Tropical and Mediterranean species had expanded to a much greater extent than that of Eurasian and Boreo-alpine species. Large-scale changes in species richness arose from several divergent, group-specific processes. Overall, local diversity especially declined in parts of southern and south-eastern Europe, reflecting local losses in multiple species rather than major range contractions among Mediterranean or Eurasian species. In fact, among the biogeographic groups, overall range declines were most prominent among Boreo-alpine species, highlighting the particular threat from climate change to this group.

Biochemistry and toxicology of proteins and peptides purified from the venom of Vipera berus berus

The isolation and characterization of individual snake venom components is important for a deeper understanding of the pathophysiology of envenomation and for improving the therapeutic procedures of patients. It also opens possibilities for the discovery of novel toxins that might be useful as tools for understanding cellular and molecular processes. The variable venom composition, toxicological and immunological properties of the common vipers (Vipera berus berus) have been reviewed. The combination of venom gland transcriptomics, bottom-up and top-down proteomics enabled comparison of common viper venom proteomes from multiple individuals. V. b. berus venom contains proteins and peptides belonging to 10–15 toxin families: snake venom metalloproteinase, phospholipases A₂ (PLA₂), snake venom serine proteinase, aspartic protease, L-amino acid oxidase (LAAO), hyaluronidase, 5′-nucleotidase, glutaminyl-peptide cyclotransferase, disintegrin, C-type lectin (snaclec), nerve growth factor, Kunitz type serine protease inhibitor, snake venom vascular endothelial growth factor, cysteine-rich secretory protein, bradykinin potentiating peptide, natriuretic peptides. PLA₂ and LAAO from V. b. berus venom produce more pronounced cytotoxic effects in cancer cells than normal cells, via induction of apoptosis, cell cycle arrest and suppression of proliferation. Proteomic data of V. b. berus venoms from different parts of Russia and Slovakian Republic have been compared with analogous data for Vipera nikolskii venom. Proteomic studies demonstrated quantitative differences in the composition of V. b. berus venom from different geographical regions. Differences in the venom composition of V. berus were mainly driven by the age, sex, habitat and diet of the snakes. The venom variability of V. berus results in a loss of antivenom efficacy against snakebites. The effectiveness of antibodies is discussed. This review presents an overview with a special focus on different toxins that have been isolated and characterized from the venoms of V. b. berus. Their main biochemical properties and toxic actions are described.

•

Vipera berus berus venom composition is variable among different populations.

•

Venom contains about 15 protein/peptide families.

•

It disturbs blood coagulation inducing pro- or anticoagulant effects.

•

Venom contains different types of blood factor X activators.

•

PLA₂ and L-amino acid oxidase produce cytotoxic effects in cancer cells.

Phylogeography and Wolbachia Infections Reveal Postglacial Recolonization Routes of the Parthenogenetic Plant Louse Cacopsylla myrtilli (W. Wagner 1947), (Hemiptera, Psylloidea)

To reveal the phylogeographic pattern of the parthenogenetic psyllid Cacopsylla myrtilli (W. Wagner 1947) (Hemiptera, Psylloidea), we sequenced a 638 bp fragment of the mitochondrial COI gene from 962 individuals. These insects originated from 46 sampling sites, which cover a significant part of the northern Palearctic distribution range of the species. The sequence data revealed 40 haplotypes, with three main (H1, H2, and H3) and 37 derived ones. The main haplotypes H1 or H2 or both were present at all sampling sites. The star-like shape of the haplotype networks indicated recent population expansion. In most cases, the derived haplotypes were specific for each country, suggesting that the main haplotypes H1 and H2 are of refugial origin, and the derived haplotypes have emerged after the postglacial recolonization process. Based on the haplotype sequences, we suggest H3 to represent the ancestral haplotype from which H1 and H2 have evolved. We suggest that the main haplotype H3 together with its derived haplotypes represents bisexual C. myrtilli, which shows a limited distribution on both sides of the border between Finland and Russia in northern Fennoscandia. The genetic diversity was the highest in Sjoa in southern Norway and also high in the White Sea region in northwest Russia. Higher diversity in Sjoa was attributed to both earlier recolonizations compared to that of the White Sea region and the absence of Wolbachia infection. We suggest that these sites were colonized from different Pleistocene refugia, i.e., from western and eastern refugia, respectively. From the White Sea region, recolonization continued eastwards to Ural Mountains and westwards to Finland and further north to Kola Peninsula. From northern Finland, recolonization continued to Finnmark, Norway, and further to Sweden and finally reached a secondary contact zone with colonizers from Norway in Central Sweden. The Caucasus and Siberian/Manchurian refugial regions have played an important role in the origin of C. myrtilli populations in Siberia and the Russian Far East.

Rare genetic admixture and unidirectional gene flow between Vipera aspis and Vipera berus at their contact zone in western France

Asp vipers (Vipera aspis) and adders (Vipera berus) coexist in the Loire-Atlantique department in France where the two species reach their respective range limits. This contact zone is of special interest since hybridization has been recently discovered there. We carried out extensive sampling to further investigate the hybrid status of morphologically typical individuals and to evaluate the proportion of genetically admixed individuals in this area. Using microsatellite markers, no hybrids with typical morphological traits of either asp vipers or adders were detected. All recently investigated individuals with intermediate morphological traits were shown to be hybrids. A rather low proportion of genetically admixed individuals (1.5-3%) and a very small number of supposed second-generation hybrids suggest reduced fertility of first-generation hybrids or low viability of their progeny. The investigation of mtDNA of newly sampled hybrids support the finding that hybridization only occurs between female V. aspis and male V. berus. Several possible explanations for the unidirectional hybridization are discussed and consequent future studies suggested.

Mitochondrial DNA diversity of the alpine newt (Ichthyosaura alpestris) in the Carpathian Basin: evidence for multiple cryptic lineages associated with Pleistocene refugia

The phylogeography and molecular taxonomy of the Alpine newt, Ichthyosaura alpestris, has been intensively studied in the past. However, previous studies did not include a comprehensive sampling from the Carpathian Basin, possibly a key region in the evolution of the species. We used a 1251 bp long fragment of the mitochondrial genome to infer the species’ evolutionary history in central-eastern Europe by assigning isolated Carpathian Basin populations from 6 regions to previously defined mtDNA lineages. We also revised the morphology-based intraspecific taxonomy of the species in the light of new genetic data. Alpine newt populations from the Carpathian Basin represented two different mitochondrial lineages. The Mátra, Bükk and Zemplén Mts populations can be assigned to the Western lineage of the nominotypical subspecies. Bakony and Őrség populations showed high haplotype diversity and formed a separate clade within the Western lineage, suggesting that the Carpathian Basin might have provided cryptic refugia for Alpine newt populations in their cold-continental forest-steppe landscapes during the younger Pleistocene. Newts from Apuseni Mts were related to the Eastern lineage but formed a distinct clade within this lineage. Considering the morphological and genetic differentiation of the Bakony and Őrség populations, consistent with a long independent evolutionary history, we propose these populations be referred to as Ichthyosaura alpestris bakonyiensis (Dely, 1964). We provide a redescription of this poorly known subspecies.

Reconstruction of forest dynamics in the Western Palaearctic based on phylogeographic analysis of the ringlet butterfly Erebia aethiops

Glacial refugia are centers of high biodiversity. Therefore, knowledge on their locations and reactions of associated populations and landscapes to climatic changes is crucial for conservation management. We here investigated the biogeography of a butterfly species linked to open forest habitats. Using mitochondrial and nuclear markers in combination with Bayesian simulations, we analyzed the location and age of potential glacial refugia of the species. We identified five putative refugia in Europe. Considering the ecological needs of our study species, tree density within these refugial areas, in contrast to earlier assumptions, must have exceeded the level of individually scattered trees. Our results also provide evidence that especially the refuge areas in the Carpathians were previously underestimated regarding their age: the refugia in the Southern Carpathians presented suitable conditions throughout several glacial cycles, probably since the Mindel or Riss cycles. Additionally, our analyses provided support for a forest refugium near the Tatra Mountains persisting the last glacial maximum. Our results underline the usefulness of this and probably other butterfly species as indicators of forest refugia.

Water deprivation compromises maternal physiology and reproductive success in a cold and wet adapted snake Vipera berus

Droughts are becoming more intense and frequent with climate change. These extreme weather events can lead to mass mortality and reproduction failure, and therefore cause population declines. Understanding how the reproductive physiology of organisms is affected by water shortages will help clarify whether females can adjust their reproductive strategy to dry conditions or may fail to reproduce and survive. In this study, we investigated the consequences of a short period of water deprivation (2 weeks) during early pregnancy on the physiology and behaviour of a cold- and wet-adapted ectotherm (Vipera berus). We also examined water allocation to developing embryos and embryonic survival. Water-deprived females exhibited significant dehydration, physiological stress and loss of muscle mass. These effects of water deprivation on water balance and muscle loss were correlated with the number of developing embryos. While water-deprived females maintained water transfer to embryos at the expense of their own maintenance, water deprivation also led to embryonic mortality. Overall, water deprivation amplifies the reproductive costs of water allocation to support embryonic development. The deleterious impacts of water deprivation on female current reproductive performance and on potential survival and future reproduction could lead to severe population declines in this species.

Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi

Background

The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai–Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history.

Results

Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche.

Conclusions

The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.

Genetic and demographic vulnerability of adder populations: Results of a genetic study in mainland Britain

Genetic factors are often overlooked in conservation planning, despite their importance in small isolated populations. We used mitochondrial and microsatellite markers to investigate population genetics of the adder (Vipera berus) in southern Britain, where numbers are declining. We found no evidence for loss of heterozygosity in any of the populations studied. Genetic diversity was comparable across sites, in line with published levels for mainland Europe. However, further analysis revealed a striking level of relatedness. Genetic networks constructed from inferred first degree relationships suggested a high proportion of individuals to be related at a level equivalent to that of half-siblings, with rare inferred full-sib dyads. These patterns of relatedness can be attributed to the high philopatry and low vagility of adders, which creates high local relatedness, in combination with the polyandrous breeding system in the adder, which may offset the risk of inbreeding in closed populations. We suggest that reliance on standard genetic indicators of inbreeding and diversity may underestimate demographic and genetic factors that make adder populations vulnerable to extirpation. We stress the importance of an integrated genetic and demographic approach in the conservation of adders, and other taxa of similar ecology.

A reference DNA barcode library for Austrian amphibians and reptiles

In the last few years, DNA barcoding became an established method for species identification in biodiversity inventories and monitoring studies. Such studies depend on the access to a comprehensive reference data base, covering all relevant taxa. Here we present a comprehensive DNA barcode inventory of all amphibian and reptile species native to Austria, except for the putatively extinct Vipera ursinii rakosiensis and Lissotriton helveticus, which has been only recently reported for the very western edge of Austria. A total of 194 DNA barcodes were generated in the framework of the Austrian Barcode of Life (ABOL) initiative. Species identification via DNA barcodes was successful for most species, except for the hybridogenetic species complex of water frogs (Pelophylax spp.) and the crested newts (Triturus spp.), in areas of sympatry. However, DNA barcoding also proved powerful in detecting deep conspecific lineages, e.g. within Natrix natrix or the wall lizard (Podarcis muralis), resulting in more than one Barcode Index Number (BIN) per species. Moreover, DNA barcodes revealed the presence of Natrix helvetica, which has been elevated to species level only recently, and genetic signatures of the Italian water frog Pelophylax bergeri in Western Austria for the first time. Comparison to previously published DNA barcoding data of European amphibians and reptiles corroborated the results of the Austrian data but also revealed certain peculiarities, underlining the particular strengths and in the case of the genus Pelophylax also the limitations of DNA barcoding. Consequently, DNA barcoding is not only powerful for species identification of all life stages of most Austrian amphibian and reptile species, but also for the detection of new species, the monitoring of gene flow or the presence of alien populations and/or species. Thus, DNA barcoding and the data generated in this study may serve both scientific and national or even transnational conservation purposes.

Biogeographic barriers, Pleistocene refugia, and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species-specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.

Vipera berus berus Venom from Russia: Venomics, Bioactivities and Preclinical Assessment of Microgen Antivenom

The common European adder, Vipera berus berus, is a medically relevant species, which is widely distributed in Russia and thus, is responsible for most snakebite accidents in Russia. We have investigated the toxic and enzymatic activities and have determined the proteomic composition of its venom. Phospholipases A₂ (PLA₂, 25.3% of the venom proteome), serine proteinases (SVSP, 16.2%), metalloproteinases (SVMP, 17.2%), vasoactive peptides (bradykinin-potentiating peptides (BPPs), 9.5% and C-type natriuretic peptides (C-NAP, 7.8%), cysteine-rich secretory protein (CRISP, 8%) and L-amino acid oxidase (LAO, 7.3%) represent the major toxin classes found in V. b. berus (Russia) venom. This study was also designed to assess the in vivo and in vitro preclinical efficacy of the Russian Microgen antivenom in neutralizing the main effects of V. b. berus venom. The results show that this antivenom is capable of neutralizing the lethal, hemorrhagic and PLA₂ activities. Third-generation antivenomics was applied to quantify the toxin-recognition landscape and the maximal binding capacity of the antivenom for each component of the venom. The antivenomics analysis revealed that 6.24% of the anti-V. b. berus F(ab’)₂ molecules fraction are toxin-binding antibodies, 60% of which represent clinically relevant antivenom molecules.

Tracing the maternal origin of the common wall lizard (Podarcis muralis) on the northern range margin in Central Europe

The maternal origin of isolated populations of the common wall lizard (Podracis muralis) in the Czech Republic, representing the north-eastern range border of the species, was addressed. We compared mitochondrial DNA sequences of the cytochrome b gene of samples from these populations with those from within the continuous range in Slovakia, the northern Balkan region, and those available from previous studies. We recorded five main haplogroups in the studied region, with all available Central European samples belonging to the same haplogroup. The star-like structure of this haplogroup suggests a scenario of relatively recent, post-glacial population expansion, which is further supported by a coalescent-based demographic analysis. The presence of unique haplotypes in two of the three isolated Czech populations together with close phylogenetic relationships to adjacent Slovak populations suggests either autochthonous origin or human-mediated introductions from geographically and genetically closest populations. We therefore support conservation programs for all three isolated Czech populations.

Extra-Mediterranean glacial refuges in barred and common grass snakes (Natrix helvetica, N. natrix)

Extra-Mediterranean glacial refugia of thermophilic biota, in particular in northern latitudes, are controversial. In the present study we provide genetic evidence for extra-Mediterranean refugia in two species of grass snake. The refuge of a widely distributed western European lineage of the barred grass snake (Natrix helvetica) was most likely located in southern France, outside the classical refuges in the southern European peninsulas. One genetic lineage of the common grass snake (N. natrix), distributed in Scandinavia, Central Europe and the Balkan Peninsula, had two distinct glacial refuges. We show that one was located in the southern Balkan Peninsula. However, Central Europe and Scandinavia were not colonized from there, but from a second refuge in Central Europe. This refuge was located in between the northern ice sheet and the Alpine glaciers of the last glaciation and most likely in a permafrost region. Another co-distributed genetic lineage of N. natrix, now massively hybridizing with the aforementioned lineage, survived the last glaciation in a structured refuge in the southern Balkan Peninsula, according to the idea of 'refugia-within-refugia'. It reached Central Europe only very recently. This study reports for the first time the glacial survival of a thermophilic egg-laying reptile species in Central Europe.

Chasing the phantom: biogeography and conservation of Vipera latastei-monticola in the Maghreb (North Africa)

The Maghreb region (North Africa) constitutes a major component of the Mediterranean Basin biodiversity hotspot. During the last centuries, a consistent human population growth has led to an unprecedented rate of habitat transformation and loss in the region and thus, threatening its biodiversity. The Western Mediterranean viper Vipera latastei-monticola inhabits humid and subhumid areas in the main mountain ranges of the Maghreb, facing such threatening factors; however, its elusive character and rarity hindered data collection for distinct biological purposes. Here, we study the biogeographical patterns and conservation status of the Maghrebian V. latastei-monticola resulting from recent sampling campaigns in Morocco and Tunisia. We update species distribution, and integrate phylogeographic and ecological niche modelling analyses at both species and lineage level to identify suitable areas, and to evaluate the impact of anthropogenic transformation and level of protection of their suitable space. We identified four highly divergent mitochondrial lineages, including a new lineage endemic to the Western High Atlas, with allopatric distributions and restricted to mountain ranges, supporting the role of mountains as past climatic refugia. Despite the remoteness of suitable areas, we report widespread habitat degradation and identify the low effectiveness of the current protected areas system in preserving the species and lineages range. Our study shows the urgent need to apply management actions for the long-term conservation of this vulnerable species and suggests a revaluation of the specific status of V. monticola, as these populations likely represent an ecotype of V. latastei.

Individual variability of venom from the European adder (Vipera berus berus) from one locality in Eastern Hungary

We have revealed intra-population variability among venom samples from several individual European adders (Vipera berus berus) within a defined population in Eastern Hungary. Individual differences in venom pattern were noticed, both gender-specific and age-related, by one-dimensional electrophoresis. Gelatin zymography demonstrated that these individual venoms have different degradation profiles indicating varying protease activity in the specimens from adders of different ages and genders. Some specimens shared a conserved region of substrate degradation, while others had lower or extremely low protease activity. Phospholipase A₂ activity of venoms was similar but not identical. Interspecimen diversity of the venom phospholipase A₂-spectra (based on the components' molecular masses) was detected by MALDI-TOF MS. The lethal toxicity of venoms (LD₅₀) also showed differences among individual snakes. Extracted venom samples had varying neuromuscular paralysing effect on chick biventer cervicis nerve-muscle preparations. The paralysing effect of venom was lost when calcium in the physiological salt solution was replaced by strontium; indicating that the block of twitch responses to nerve stimulation is associated with the activity of a phospholipase-dependent neurotoxin. In contrast to the studied V. b. berus venoms from different geographical regions so far, this is the first V. b. berus population discovered to have predominantly neurotoxic neuromuscular activity. The relevance of varying venom yields is also discussed. This study demonstrates that individual venom variation among V. b. berus living in particular area of Eastern Hungary might contribute to a wider range of clinical manifestations of V. b. berus envenoming than elsewhere in Europe.

Vipera berus (Linnaeus, 1758) remains from the Late Pleistocene of Slovakia

Completely preserved specimens of fossil snakes are extremely rare and ophidian palaeontologists are usually dependent only on disarticulated elements of a postcranial skeleton. Here we present an unusually well-preserved specimen of a small viperid snake from the Late Pleistocene firm travertine at the famous Gánovce-Hrádok Neanderthal mound in Slovakia. The complex study of both cranial and axial skeleton with well-preserved maxilla and basiparasphenoid confirms the presence of a viper from theVipera berusspecies group, and recent distribution ofV. berusspecies complex members supports identification of these preserved remains as belonging to the common adder,V. berus(Linnaeus, 1758). Associated faunal assemblages of the MFG-C and D mammalian fauna groups reported from the firm travertine indicate a humid climate in a predominantly woodland environment with typical forest species in the Gánovce-Hrádok vicinity throughout the Eemian optimum, and mixed forest and steppe environments in the late Eemian to early Weichselian stages. Occurrence ofV. berusdocuments the presence of open or semi-open biotopes with low vegetation. AlthoughV. berusoccurs in the Quaternary glacial/interglacial cycle and throughout the entire warm part, it never dominated herpetofaunal assemblages during the climatic optimum. Therefore, the presence ofV. berusmost likely indicates late Eemian or early Weichselian (interstadial) climatic conditions.

A new population and subspecies of the critically endangered Anatolian meadow viper Vipera anatolica Eiselt and Baran, 1970 in eastern Antalya province

We report on a new population of Vipera anatolica from the Geyik Mountain Range in eastern Antalya Province, Turkey. It represents only the second known location, and is situated in a valley about 200 km east from the terra typica at Kohu Dağ in western Antalya Province. We compare both populations and, based on marked differences in morphology, habitat, genetics, and its geographically isolated location, we describe the recently discovered population as a new subspecies. Aspects of ecology, threats, and conservation needs are discussed.

Retreat and extinction of the Late Pleistocene cave bear (Ursus spelaeus sensu lato)

The cave bear (Ursus spelaeus sensu lato) is a typical representative of Pleistocene megafauna which became extinct at the end of the Last Glacial. Detailed knowledge of cave bear extinction could explain this spectacular ecological transformation. The paper provides a report on the youngest remains of the cave bear dated to 20,930 ± 140 ¹⁴C years before present (BP). Ancient DNA analyses proved its affiliation to the Ursus ingressus haplotype. Using this record and 205 other dates, we determined, following eight approaches, the extinction time of this mammal at 26,100–24,300 cal. years BP. The time is only slightly earlier, i.e. 27,000–26,100 cal. years BP, when young dates without associated collagen data are excluded. The demise of cave bear falls within the coldest phase of the last glacial period, Greenland Stadial 3. This finding and the significant decrease in the cave bear records with cooling indicate that the drastic climatic changes were responsible for its extinction. Climate deterioration lowered vegetation productivity, on which the cave bear strongly depended as a strict herbivore. The distribution of the last cave bear records in Europe suggests that this animal was vanishing by fragmentation into subpopulations occupying small habitats. One of them was the Kraków-Częstochowa Upland in Poland, where we discovered the latest record of the cave bear and also two other, younger than 25,000 ¹⁴C years BP. The relatively long survival of this bear in karst regions may result from suitable microclimate and continuous access to water provided by deep aquifers, indicating a refugial role of such regions in the Pleistocene for many species.

The adder (Vipera berus) in Southern Altay Mountains: population characteristics, distribution, morphology and phylogenetic position

As the most widely distributed snake in Eurasia, the adder (Vipera berus) has been extensively investigated in Europe but poorly understood in Asia. The Southern Altay Mountains represent the adder's southern distribution limit in Central Asia, whereas its population status has never been assessed. We conducted, for the first time, field surveys for the adder at two areas of Southern Altay Mountains using a combination of line transects and random searches. We also described the morphological characteristics of the collected specimens and conducted analyses of external morphology and molecular phylogeny. The results showed that the adder distributed in both survey sites and we recorded a total of 34 sightings. In Kanas river valley, the estimated encounter rate over a total of 137 km transects was 0.15 ± 0.05 sightings/km. The occurrence of melanism was only 17%. The small size was typical for the adders in Southern Altay Mountains in contrast to other geographic populations of the nominate subspecies. A phylogenetic tree obtained by Bayesian Inference based on DNA sequences of the mitochondrial cytochrome b (1,023 bp) grouped them within the Northern clade of the species but failed to separate them from the subspecies V. b. sachalinensis. Our discovery extends the distribution range of V. berus and provides a basis for further researches. We discuss the hypothesis that the adder expands its distribution border to the southwest along the mountains' elevation gradient, but the population abundance declines gradually due to a drying climate.

Merging and comparing three mitochondrial markers for phylogenetic studies of Eurasian reindeer (Rangifer tarandus)

Phylogenetic analyses provide information that can be useful in the conservation of genetic variation by identifying intraspecific genetic structure. Reconstruction of phylogenetic relationships requires the use of markers with the appropriate amount of variation relative to the timeframe and purpose of the study. Here, genetic structure and clustering are inferred from comparative analyses of three widely used mitochondrial markers, the CR, cytb and the COI region, merged and separately, using Eurasian reindeer as a model. A Bayesian phylogeny and a MJ network, both based on the merged dataset, indicate several distinct maternal haplotype clusters within Eurasian reindeer. In addition to confirm previously described clusters, two new subclusters were found. When comparing the results from the merged dataset with the results from analyses of the three markers separately, similar clustering was found in the CR and COI phylogenies, whereas the cytb region showed poor resolution. Phylogenetic analyses of the merged dataset and the CR revealed congruent results, implying that single sequencing analysis of the CR is an applicable method for studying the haplotype structure in Eurasian reindeer.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central-marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central-marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold-adapted species likely used two isolated glacial refugia in southern France, in permafrost-free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.

Persistence and dispersal in a Southern Hemisphere glaciated landscape: the phylogeography of the spotted snow skink (Niveoscincus ocellatus) in Tasmania

Background

The aim of this research was to identify the effects of Pleistocene climate change on the distribution of fauna in Tasmania, and contrast this with biotic responses in other temperate regions in the Northern and Southern Hemisphere that experienced glacial activity during this epoch. This was achieved by examining the phylogeographic patterns in a widely distributed Tasmanian endemic reptile, Niveoscincus ocellatus. 204 individuals from 29 populations across the distributional range of N. ocellatus were surveyed for variation at two mitochondrial genes (ND2, ND4), and two nuclear genes (β-globin, RPS8). Phylogenetic relationships were reconstructed using a range of methods (maximum parsimony, Bayesian inference and haplotype networks), and the demographic histories of populations were assessed (AMOVA, Tajima’s D, Fu’s Fs, mismatch distributions, extended Bayesian skyline plots, and relaxed random walk analyses).

Results

There was a high degree of mitochondrial haplotype diversity (96 unique haplotypes) and phylogeographic structure, where spatially distinct groups were associated with Tasmania’s Northeast and a large area covering Southeast and Central Tasmania. Phylogeographic structure was also present within each major group, but the degree varied regionally, being highest in the Northeast. Only the Southeastern group had a signature of demographic expansion, occurring during the Pleistocene but post-dating the Last Glacial Maximum. In contrast, nuclear DNA had low levels of variation and a lack of phylogeographic structure, and further loci should be surveyed to corroborate the mitochondrial inferences.

Conclusions

The phylogeographic patterns of N. ocellatus indicate Pleistocene range and demographic expansion in N. ocellatus, particularly in the Southeast and Central areas of Tasmania. Expansion in Central and Southeastern areas appears to have been more recent in both demographic and spatial contexts, than in Northeast Tasmania, which is consistent with inferences for other taxa of greater stability and persistence in Northeast Tasmania during the Last Glacial Maximum. These phylogeographic patterns indicate contrasting demographic histories of populations in close proximity to areas directly affected by glaciers in the Southern Hemisphere during the LGM.

Electronic supplementary material

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