Calibrating the tree of vipers under the fossilized birth-death model

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.

How many specimens make a sufficient training set for automated three-dimensional feature extraction?

Deep learning has emerged as a robust tool for automating feature extraction from three-dimensional images, offering an efficient alternative to labour-intensive and potentially biased manual image segmentation methods. However, there has been limited exploration into the optimal training set sizes, including assessing whether artficial expansion by data augmentation can achieve consistent results in less time and how consistent these benefits are across different types of traits. In this study, we manually segmented 50 planktonic foraminifera specimens from the genus Menardella to determine the minimum number of training images required to produce accurate volumetric and shape data from internal and external structures. The results reveal unsurprisingly that deep learning models improve with a larger number of training images with eight specimens being required to achieve 95% accuracy. Furthermore, data augmentation can enhance network accuracy by up to 8.0%. Notably, predicting both volumetric and shape measurements for the internal structure poses a greater challenge compared with the external structure, owing to low contrast differences between different materials and increased geometric complexity. These results provide novel insight into optimal training set sizes for precise image segmentation of diverse traits and highlight the potential of data augmentation for enhancing multivariate feature extraction from three-dimensional images.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

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.

The bigger the threat, the longer the gaze? A cross-cultural study of Somalis and Czechs

High fear reaction, preferential attention, or fast detection are only a few of the specific responses which snakes evoke in humans. Previous research has shown that these responses are shared amongst several distinct cultures suggesting the evolutionary origin of the response. However, populations from sub-Saharan Africa have been largely missing in experimental research focused on this issue. In this paper, we focus on the effect of snake threat display on human spontaneous attention. We performed an eye-tracking experiment with participants from Somaliland and the Czechia and investigated whether human attention is swayed towards snakes in a threatening posture. Seventy-one Somalis and 71 Czechs were tested; the samples were matched for gender and comparable in age structure and education level. We also investigated the effect of snake morphotype as snakes differ in their threat display. We found that snakes in a threatening posture were indeed gazed upon more than snakes in a relaxed (non-threatening) posture. Further, we found a large effect of snake morphotype as this was especially prominent in cobras, less in vipers, and mostly non-significant in other morphotypes. Finally, despite highly different cultural and environmental backgrounds, the overall pattern of reaction towards snakes was similar in Somalis and Czechs supporting the evolutionary origin of the phenomenon. We concluded that human attention is preferentially directed towards snakes, especially cobras and vipers, in threatening postures.

Estimating the Age of Poorly Dated Fossil Specimens and Deposits Using a Total-Evidence Approach and the Fossilized Birth-Death Process

Bayesian total-evidence approaches under the fossilized birth-death model enable biologists to combine fossil and extant data while accounting for uncertainty in the ages of fossil specimens, in an integrative phylogenetic analysis. Fossil age uncertainty is a key feature of the fossil record as many empirical data sets may contain a mix of precisely dated and poorly dated fossil specimens or deposits. In this study, we explore whether reliable age estimates for fossil specimens can be obtained from Bayesian total-evidence phylogenetic analyses under the fossilized birth-death model. Through simulations based on the example of the Baltic amber deposit, we show that estimates of fossil ages obtained through such an analysis are accurate, particularly when the proportion of poorly dated specimens remains low and the majority of fossil specimens have precise dates. We confirm our results using an empirical data set of living and fossil penguins by artificially increasing the age uncertainty around some fossil specimens and showing that the resulting age estimates overlap with the recorded age ranges. Our results are applicable to many empirical data sets where classical methods of establishing fossil ages have failed, such as the Baltic amber and the Gobi Desert deposits. [Bayesian phylogenetic inference; fossil age estimates; fossilized birth-death; Lagerstätte; total-evidence.].

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.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

Extensive Variation in the Activities of Pseudocerastes and Eristicophis Viper Venoms Suggests Divergent Envenoming Strategies Are Used for Prey Capture

Snakes of the genera Pseudocerastes and Eristicophis (Viperidae: Viperinae) are known as the desert vipers due to their association with the arid environments of the Middle East. These species have received limited research attention and little is known about their venom or ecology. In this study, a comprehensive analysis of desert viper venoms was conducted by visualising the venom proteomes via gel electrophoresis and assessing the crude venoms for their cytotoxic, haemotoxic, and neurotoxic properties. Plasmas sourced from human, toad, and chicken were used as models to assess possible prey-linked venom activity. The venoms demonstrated substantial divergence in composition and bioactivity across all experiments. Pseudocerastes urarachnoides venom activated human coagulation factors X and prothrombin and demonstrated potent procoagulant activity in human, toad, and chicken plasmas, in stark contrast to the potent neurotoxic venom of P. fieldi. The venom of E. macmahonii also induced coagulation, though this did not appear to be via the activation of factor X or prothrombin. The coagulant properties of P. fieldi and P. persicus venoms varied among plasmas, demonstrating strong anticoagulant activity in the amphibian and human plasmas but no significant effect in that of bird. This is conjectured to reflect prey-specific toxin activity, though further ecological studies are required to confirm any dietary associations. This study reinforces the notion that phylogenetic relatedness of snakes cannot readily predict venom protein composition or function. The significant venom variation between these species raises serious concerns regarding antivenom paraspecificity. Future assessment of antivenom is crucial.

Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification

Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.

Widespread Evolution of Molecular Resistance to Snake Venom α-Neurotoxins in Vertebrates

Venomous snakes are important subjects of study in evolution, ecology, and biomedicine. Many venomous snakes have alpha-neurotoxins (α-neurotoxins) in their venom. These toxins bind the alpha-1 nicotinic acetylcholine receptor (nAChR) at the neuromuscular junction, causing paralysis and asphyxia. Several venomous snakes and their predators have evolved resistance to α-neurotoxins. The resistance is conferred by steric hindrance from N-glycosylated asparagines at amino acids 187 or 189, by an arginine at position 187 that has been hypothesized to either electrostatically repulse positively charged neurotoxins or sterically interfere with α-neurotoxin binding, or proline replacements at positions 194 or 197 of the nAChR ligand-binding domain to inhibit α-neurotoxin binding through structural changes in the receptor. Here, we analyzed this domain in 148 vertebrate species, and assessed its amino acid sequences for resistance-associated mutations. Of these sequences, 89 were sequenced de novo. We find widespread convergent evolution of the N-glycosylation form of resistance in several taxa including venomous snakes and their lizard prey, but not in the snake-eating birds studied. We also document new lineages with the arginine form of inhibition. Using an in vivo assay in four species, we provide further evidence that N-glycosylation mutations reduce the toxicity of cobra venom. The nAChR is of crucial importance for normal neuromuscular function and is highly conserved throughout the vertebrates as a result. Our research shows that the evolution of α-neurotoxins in snakes may well have prompted arms races and mutations to this ancient receptor across a wide range of sympatric vertebrates. These findings underscore the inter-connectedness of the biosphere and the ripple effects that one adaption can have across global ecosystems.

Origin, extinction and ancient DNA of a new fossil insular viper: molecular clues of overseas immigration

Viperinae is a subfamily of viperid snakes whose fossil record in the Mediterranean islands is, until now, restricted to 12 palaeontological deposits on seven islands. Revision of the material excavated 30 years ago from the Middle/Late Pleistocene–Holocene deposit of Es Pouàs [Eivissa (= Ibiza), Balearic Islands, western Mediterranean] revealed about 6000 bones of a small-sized viper across different stratigraphic levels. Its morphological characteristics are different enough to known species of Vipera to warrant the description of a new species, but the nearly complete mitochondrial genome obtained from this snake based on a sample dated to 16 130 ± 45 bp, suggested it belonged to a new insular population of Lataste’s viper (Vipera latastei), Vipera latastei ebusitana subsp. nov. Phylogenetic analysis indicates that the dispersal of the ancestors of V. l. ebusitana to Eivissa, most probably from a north-east Iberian population, occurred via overwater colonization < 1.5 Mya, well after the Messinian Salinity Crisis (5.97–5.32 Mya) when land bridges allowed terrestrial colonization of the Balearic Islands by mainland faunas. The morphological differences between V. l. ebusitana and the Iberian populations suggest that it is a new dwarf taxon resulting from insular evolutionary processes, becoming extinct shortly after the first human arrival to this island about 4000 years ago.