Large-scale molecular phylogeny, morphology, divergence-time estimation, and the fossil record of advanced caenophidian snakes (Squamata: Serpentes)

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.

Melanin-based color variation in response to changing climates in snakes

Melanism, the process of heavier melanin deposition, can interact with climate variation at both micro and macro scales, ultimately influencing color evolution in organisms. While the ecological processes regulating melanin production in relation to climate have been extensively studied, intraspecific variations of melanism are seldom considered. Such scientific gap hampers our understanding of how species adapt to rapidly changing climates. For example, dark coloration may lead to higher heat absorption and be advantageous in cool climates, but also in hot environments as a UV or antimicrobial protection mechanism. To disentangle such opposing predictions, here we examined the effect of climate on shaping melanism variation in 150 barred grass snakes (Natrix helvetica) and 383 green whip snakes (Hierophis viridiflavus) across Italy. By utilizing melanistic morphs (charcoal and picturata in N. helvetica, charcoal and abundistic in H. viridiflavus) and compiling observations from 2002 to 2021, we predicted that charcoal morphs in H. viridiflavus would optimize heat absorption in cold environments, while offering protection from excessive UV radiation in N. helvetica within warm habitats; whereas picturata and abundistic morphs would thrive in humid environments, which naturally have a denser vegetation and wetter substrates producing darker ambient light, thus providing concealment advantages. While picturata and abundistic morphs did not align with our initial humidity expectations, the charcoal morph in N. helvetica is associated with UV environments, suggesting protection mechanisms against damaging solar radiation. H. viridiflavus is associated with high precipitations, which might offer antimicrobial protection. Overall, our results provide insights into the correlations between melanin-based color morphs and climate variables in snake populations. While suggestive of potential adaptive responses, future research should delve deeper into the underlying mechanisms regulating this relationship.

Ecological drivers of ultraviolet colour evolution in snakes

Ultraviolet (UV) colour patterns invisible to humans are widespread in nature. However, research bias favouring species with conspicuous colours under sexual selection can limit our assessment of other ecological drivers of UV colour, like interactions between predators and prey. Here we demonstrate widespread UV colouration across Western Hemisphere snakes and find stronger support for a predator defence function than for reproduction. We find that UV colouration has evolved repeatedly in species with ecologies most sensitive to bird predation, with no sexual dichromatism at any life stage. By modelling visual systems of potential predators, we find that snake conspicuousness correlates with UV colouration and predator cone number, providing a plausible mechanism for selection. Our results suggest that UV reflectance should not be assumed absent in "cryptically coloured" animals, as signalling beyond human visual capacities may be a key outcome of species interactions in many taxa for which UV colour is likely underreported.

Embryonic development of the neotropical pit viper Bothrops atrox (Serpentes: Viperidae: Crotalinae), with emphasis on pit organ morphogenesis and its evolution in snakes

BACKGROUND

Bothrops atrox is a pit viper with a loreal pit organ, and its embryological development remains undescribed. Here, we provide a comprehensive description of the embryology of B. atrox, focusing on the loreal pit organ and cephalic scales.

RESULTS

We characterized 13 developmental stages of B. atrox based on external features consistent with the embryogenesis of previously described snake species. The loreal pit organ originates from the circumorbital region and migrates to its final position. In Crotalinae, the pit organ first becomes visible at stage 28, whereas in Pythonidae labial, pit organs appear at Stage 35. Pit organs evolved independently three times in Serpentes, encompassing Boidae, Pythonidae, and Crotalinae. Boidae lacks embryological information for pit organs. Furthermore, we observed that head scalation onset occurs at Stage 33 in B. atrox, with fusion of scales surrounding the loreal pit organ.

CONCLUSIONS

The embryology of pit organs in Pythonidae and Boidae species remains poorly understood. Our detailed embryological descriptions are critical for proposing developmental scenarios for pit organs and guiding future research on these structures.

Red-on-Yellow Queen: Bio-Layer Interferometry Reveals Functional Diversity Within Micrurus Venoms and Toxin Resistance in Prey Species

Snakes in the family Elapidae largely produce venoms rich in three-finger toxins (3FTx) that bind to the α 1 subunit of nicotinic acetylcholine receptors (nAChRs), impeding ion channel activity. These neurotoxins immobilize the prey by disrupting muscle contraction. Coral snakes of the genus Micrurus are specialist predators who produce many 3FTx, making them an interesting system for examining the coevolution of these toxins and their targets in prey animals. We used a bio-layer interferometry technique to measure the binding interaction between 15 Micrurus venoms and 12 taxon-specific mimotopes designed to resemble the orthosteric binding region of the muscular nAChR subunit. We found that Micrurus venoms vary greatly in their potency on this assay and that this variation follows phylogenetic patterns rather than previously reported patterns of venom composition. The long-tailed Micrurus tend to have greater binding to nAChR orthosteric sites than their short-tailed relatives and we conclude this is the likely ancestral state. The repeated loss of this activity may be due to the evolution of 3FTx that bind to other regions of the nAChR. We also observed variations in the potency of the venoms depending on the taxon of the target mimotope. Rather than a pattern of prey-specificity, we found that mimotopes modeled after snake nAChRs are less susceptible to Micrurus venoms and that this resistance is partly due to a characteristic tryptophan → serine mutation within the orthosteric site in all snake mimotopes. This resistance may be part of a Red Queen arms race between coral snakes and their prey.

Paedomorphosis and retention of juvenile diet lead speciation in a group of Neotropical snakes (Colubroides-Philodryadini)

Dipsadidae is one of the largest clades of extant reptiles, showing an impressive morphological and ecological diversity. Despite this fact, the developmental processes behind its diversity are still largely unknown. In this study, we used 3D reconstructions based on micro-CT data and geometric morphometrics to evaluate the skull morphology of Philodryas agassizii, a small, surface-dwelling dipsadid that consume spiders. Adult individuals of P. agassizii exhibit a cranial morphology frequently observed in juveniles of other surface-dwelling colubroideans, represented in our analysis by its close relative Philodryas patagoniensis. Large orbits, gibbous neurocranial roof and a relatively short jaw complex are features present in juveniles of the latter species. Furthermore, we performed an extensive survey about diet of P. patagoniensis in which we detected an ontogenetic dietary shift, indicating that arthropods are more frequently consumed by juveniles of this dietary generalist. Thus, we infer that P. agassizzii retained not only the ancestral juvenile skull morphology but also dietary preferences. This study reveals that morphological changes driven by heterochronic changes, specifically paedomorphosis, influenced the retention of ancestral life history traits in P. agassizii, and therefore promoted cladogenesis. In this way, we obtained first evidence that heterochronic processes lead speciation in the snake megadiverse clade Dipsadidae.

Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Comparative morphology of oral glands in snakes of the family Homalopsidae reveals substantial variation and additional independent origins of salt glands within Serpentes

Using diffusible iodine-based contrast-enhanced computed tomography (diceCT), we examined the morphology of the oral glands of 12 species of the family Homalopsidae. Snakes of this family exhibit substantial interspecific morphological variation in their oral glands. Particular variables are the venom glands, ranging from large (e.g., Subsessor bocourti) to small (e.g., Erpeton tentaculatum). The supra- and infralabial glands are more uniform in morphology, being the second most developed in almost all the sampled species. Premaxillary glands distinct from the supralabial glands were observed in five species (Myron richardsonii, Bitia hydroides, Cantoria violacea, Fordonia leucobalia, and Gerarda prevostiana), in addition to Cerberus rynchops, the only species in which this condition was previously documented associated with the excretion of salt. In the three species of the saltwater group of homalopsids (C. violacea, F. leucobalia, and G. prevostiana), the premaxillary glands also extend posteriorly, occupying a large area above the supralabial gland, a condition not observed in any other species of snake studied thus far. Character evolution analyses indicate that premaxillary glands differentiated from the supralabial gland and evolved independently three or four times in the family, always in lineages that invaded marine habitats. Our results suggest that the differentiated premaxillary glands are likely salt glands, as is the case in C. rynchops. If corroborated, this increases to six or seven the number of independent evolutionary origins of salt glands in snakes that have undergone an evolutionary transition to marine life.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

Foraging predicts the evolution of warning coloration and mimicry in snakes

Warning coloration and Batesian mimicry are classic examples of Darwinian evolution, but empirical evolutionary patterns are often paradoxical. We test whether foraging costs predict the evolution of striking coloration by integrating genetic and ecological data for aposematic and mimetic snakes (Elapidae and Dipsadidae). Our phylogenetic comparison on a total of 432 species demonstrated that dramatic changes in coloration were well predicted by foraging strategy. Multiple tests consistently indicated that warning coloration and conspicuous mimicry were more likely to evolve in species where foraging costs of conspicuous appearance were relaxed by poor vision of their prey, concealed habitat, or nocturnal activity. Reversion to crypsis was also well predicted by ecology for elapids but not for dipsadids. In contrast to a theoretical prediction and general trends, snakes' conspicuous coloration was correlated with secretive ecology, suggesting that a selection regime underlies evolutionary patterns. We also found evidence that mimicry of inconspicuous models (pitvipers) may have evolved in association with foraging demand for crypsis. These findings demonstrate that foraging is an important factor necessary to understand the evolution, persistence, and diversity of warning coloration and mimicry of snakes, highlighting the significance of additional selective factors in solving the warning coloration paradox.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae)

BACKGROUND

Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates.

RESULTS

Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development.

CONCLUSIONS

We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

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.

Ecological correlates of cranial evolution in the megaradiation of dipsadine snakes

BACKGROUND

Dipsadine snakes represent one of the most spectacular vertebrate radiations that have occurred in any continental setting, with over 800 species in South and Central America. Their species richness is paralleled by stunning ecological diversity, ranging from arboreal snail-eating and aquatic eel-eating specialists to terrestrial generalists. Despite the ecological importance of this clade, little is known about the extent to which ecological specialization shapes broader patterns of phenotypic diversity within the group. Here, we test how habitat use and diet have influenced morphological diversification in skull shape across 160 dipsadine species using micro-CT and 3-D geometric morphometrics, and we use a phylogenetic comparative approach to test the contributions of habitat use and diet composition to variation in skull shape among species.

RESULTS

We demonstrate that while both habitat use and diet are significant predictors of shape in many regions of the skull, habitat use significantly predicts shape in a greater number of skull regions when compared to diet. We also find that across ecological groupings, fossorial and aquatic behaviors result in the strongest deviations in morphospace for several skull regions. We use simulations to address the robustness of our results and describe statistical anomalies that can arise from the application of phylogenetic generalized least squares to complex shape data.

CONCLUSIONS

Both habitat and dietary ecology are significantly correlated with skull shape in dipsadines; the strongest relationships involved skull shape in snakes with aquatic and fossorial lifestyles. This association between skull morphology and multiple ecological axes is consistent with a classic model of adaptive radiation and suggests that ecological factors were an important component in driving morphological diversification in the dipsadine megaradiation.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Molecular phylogeny reveals cryptic diversity in Sibynophis from China (Serpentes: Sibynophiidae)

The elucidation of species diversity and distribution is critical within the fields of evolution, genetics, and conservation. The genus Sibynophis contains rare snakes that have historically received little attention. In this study, we conducted comprehensive sampling and use both mitochondrial and nuclear genetic markers to explore Sibynophis species diversity within China. Our findings revealed that S. c. miyiensis should be considered synonymous with S. c. grahami, and S. c. grahami should be gave a specific rank as S. grahami. In addition, we discovered S. triangularis was new to China and Myanmar. Based on the specimens and molecular phylogeny results, we redefined the species distribution boundaries of each Chinese species.

Unveiling underestimated species diversity within the Central American Coralsnake, a medically important complex of venomous taxa

Coralsnakes of the genus Micrurus are a diverse group of venomous snakes ranging from the southern United States to southern South America. Much uncertainty remains over the genus diversity, and understanding Micrurus systematics is of medical importance. In particular, the widespread Micrurus nigrocinctus spans from Mexico throughout Central America and into Colombia, with a number of described subspecies. This study provides new insights into the phylogenetic relationships within M. nigrocinctus by examining sequence data from a broad sampling of specimens from Mexico, Guatemala, Honduras, Nicaragua, Costa Rica, and Panama. The recovered phylogenetic relationships suggest that M. nigrocinctus is a species complex originating in the Pliocene and composed of at least three distinct species-level lineages. In addition, recovery of highly divergent clades supports the elevation of some currently recognized subspecies to the full species rank while others may require synonymization.

Identity of the holotype and type locality of Rhabdophis leonardi (Wall, 1923) (Colubridae: Natricinae), with notes on the morphology and natural history of the species in southwestern China

The original description of Natrix leonardi (currently Rhabdophis leonardi) by Frank Wall in 1923, based on a specimen from the "Upper Burma Hills," lacked important morphological details that have complicated the assignment of recently collected material. Furthermore, although the holotype was never lost, its location has been misreported in one important taxonomic reference, leading to further confusion. We report the correct repository of the holotype (Natural History Museum, London), together with its current catalog number. We also describe key features of that specimen that were omitted from the original description, and provide new details on the morphology of the species, including sexual dichromatism unusual for the genus, based upon specimens from southern Sichuan, China. Rhabdophis leonardi is distinguished from its congeners by the following characters: 15 or 17 DSR at midbody and 6 supralabials; distinct annulus around the neck, broad and red in males, and narrow and orange with a black border in females; dorsal ground color light green or olive; some lateral and dorsal scales possessing black edges, the frequency of black edges gradually increasing from anterior to posterior, forming irregular and ill-defined transverse black bands; eye with prominent green iris; black ventral spots with a red edge, most numerous at midbody but extending halfway down the length of the tail. In southwestern China, this species is frequently found at 1730-2230 m elevation. It has been documented to prey upon anuran amphibians, including toads. A recently published phylogenetic analysis showed this species to be deeply nested with the genus Rhabdophis, as a member of the R. nuchalis Group. That analysis also revealed the existence of two closely related but geographically distinct subclades in the molecular analysis, one of which may represent an unnamed taxon.

Ultraconserved elements-based phylogenomic systematics of the snake superfamily Elapoidea, with the description of a new Afro-Asian family

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes.

T. Lalremsanga EO, Bhosale H, Gowande G, Patel H, Idiatullina SS, Poyarkov NA. Systematics of Trimeresurus popeiorum Smith, 1937 with a revised molecular phylogeny of Asian pitvipers of the genus Trimeresurus Lacépède, 1804 sensu lato

The pit viper snake genus Trimeresurus Lacépède, 1804 sensu lato, is a diverse group of nocturnal serpents comprising over 61 species. The genus is morphologically heterogeneous and has been divided into several subgenera. We present an updated phylogeny of Asian pitvipers and propose a revised classification. Additionally, we revise the taxonomy of T. popeiorum Smith, 1937 and propose taxonomic changes with support from molecular and morphological data. We restrict T. popeiorumsensu stricto to northeastern India, Bangladesh, southern China, and northern Myanmar; populations beyond these areas require further assessment. We also synonymize T. yingjiangensisChen et al., 2019 with T. popeiorum based on overlapping morphological characters, molecular data, and distribution. The findings shed new light on the taxonomy of T. popeiorum, warranting the need for assessing the population of T. popeiorum from southeast Asia.

Ejaculate characteristics over seasons in five species of lancehead pitvipers (Bothrops spp) kept in captivity

Due to their major medical importance in Latin America, lancehead pitvipers are frequently kept and bred in captivity for venom extraction to the production of antivenom serums. Nevertheless, despite the great contribution given to captive breeding, much of the knowledge of Bothrops' reproductive biology derived from sporadic and insufficient data provided by zoological collections. Thus, we aimed to investigate seasonal changes in gonadosomatic index (GSI) and seminal parameters (e.g., volume, concentration, motility, viability, and acrosome integrity) of five species of lancehead pitvipers from different biomes and phylogenetic groups, maintained in the indoors serpentarium at Butantan Institute (Brazil). Patterns of variation in GSI and semen parameters differed from one species to another, suggesting that captive populations should perhaps be managed distinctly to maximize reproductive success. Furthermore, in none of the studied species did changes in GSI occur concomitantly with seminal variations. GSI remained unaltered year-round for Jararaca (Bothrops jararaca) and Brazilian lancehead (Bothrops moojeni), whereas it peaked in the autumn for Common lancehead (Bothrops atrox), Jararacussu (Bothrops jararacussu), and Whitetail lancehead (Bothrops leucurus). But surprisingly, the scenario was inverted when we estimated the total number of motile spermatozoa per season, as Jararaca and Brazilian lancehead displayed seasonal differences and the other species did not vary throughout the year. Potential ecological and evolutionary factors underlying these differences were also discussed in the present article. Together, these findings can help to better define breeding management strategies for each species in captivity, in addition to optimizing the future use of artificial insemination and semen cryopreservation.

Increasing taxon sampling suggests a complete taxonomic rearrangement in Echinantherini (Serpentes: Dipsadidae)

Although the recent advances on the relationship of its major groups, the systematics of the rich fauna of Neotropical snakes is far from being a consensus. In this sense, derived groups presenting continental distributions have represented a main challenge. The taxonomy of the snake tribe Echinantherini is one of the most contentious among the diverse family known as Dipsadidae. The tribe is poorly sampled in phylogenetic studies, resulting in conflicting hypotheses of relationships among its taxa. Moreover, several rare and micro endemic species of Echinantherini have never been evaluated within a comprehensive phylogenetic framework. Here, we assess for the first time the phylogenetic position of the rare Echinanthera amoena within Echinantherini. We based our analyses on a comprehensive multilocus dataset including 14 of the 16 species described for the tribe. Our results support the monophyly of Echinantherini and strongly indicate E. amoena as a unique lineage, phylogenetically positioned apart from all other congeners. From the three current genera (Echinanthera, Taeniophallus, and Sordellina) our results indicate that Echinanthera and Taeniophallus are paraphyletic, since the T. affinis species group is positioned as sister to Echinanthera (except E. amoena) clustering apart from the clade formed by the T. brevirostris and T. occipitalis groups. We describe new genera for the T. affinis and T. occipitalis species groups and an additional monospecific genus for E. amoena. Although we did not evaluate the phylogenetic position of T. nebularis, we described a new genus and removed it from Echinantherini since its morphology strikingly departs from all species now included in the tribe. Finally, we redefine the genera Echinanthera and Taeniophallus and we provide comments about further directions to study the biogeography and the evolution of morphological traits in Echinantherini.

Complete mitochondrial genomes of Boigakraepelini and Hebiuscraspedogaster (Reptilia, Squamata, Colubridae) and their phylogenetic implications

The complete sequence of the mitochondrial genome is a powerful tool for studying phylogenetic relationships and molecular evolution in various species. In this work, the mitogenomes of Boigakraepelini and Hebiuscraspedogaster were sequenced and characterized for the first time. The lengths of the B.kraepelini and H.craspedogaster mitogenomes were 17,124 bp and 17,120 bp, respectively, and both included 13 protein-coding genes, 22 tRNAs, two rRNAs and two control regions. The arrangements of these mitochondrial genes were the same in B.kraepelini and H.craspedogaster. In addition, both genome compositions showed A+T bias (59.03%, 60.93%) and had positive AT skews (0.179, 0.117) and negative GC skews (-0.397, -0.348). The phylogenetic results illustrated a close relationship between B.kraepelini and the genus Lycodon. Moreover, H.craspedogaster was clustered with other Hebius snakes and closely related to other Natricinae species. These results will provide references for further research on the phylogeny of Colubridae.

Ophidiomyces ophidiicola detection and infection: a global review on a potential threat to the world’s snake populations

Ophidiomyces ophidiicola (Oo) is one of the most relevant fungal pathogens for snakes. It is the etiological agent of ophidiomycosis, an emerging disease causing dysecdysis, skin abnormalities, crusting cutaneous lesions, and ulcerations. Despite this major tegumentary “tropism”, Oo infection can be systemic and it is capable of inducing visceral lesions. Moreover, ophidiomycosis may lead to abnormalities of reproductive physiology, hunting behavior, and thermoregulation, thus increasing the risks of sublethal effects and predation on affected snakes. Oo seems horizontally transmitted and can induce postnatal mortality. This article reviews published data on Oo detection and infection in all snake species in countries around the world and categorizes these data using new classification parameters. The presence of this fungus has been recorded in 11 states (considering the USA as a whole); however, in four states, the mycosis has only been reported in snakes held in captivity. Detection and/or infection of Oo has been ascertained in 62 snake species, divided into nine families. The taxa have been categorized with diagnostic criteria in order to report, for each species, the highest rank of categorization resulting from all cases. Therefore, 20 species have been included within the class “Ophidiomycosis and Oo shedder”, 11 within “Ophidiomycosis”, 16 in “Apparent ophidiomycosis”, and 15 within “Ophidiomyces ophidiicola present”. We also discuss the significance and limits of case classifications and Oo’s impact on wild populations, and we suggest methods for preliminary surveillance. Standardized methods, interdisciplinary studies, and cooperation between various research institutions may facilitate further Oo screening studies, elucidate the unclear aspects of the disease, and protect ophidiofauna from this emerging threat at the global level.

A new species of the genus Achalinus (Squamata: Xenodermidae) from Ningshan County, Shaanxi Province, China

A new species of odd-scaled snake in the genus Achalinus is described from Ningshan County, Shaanxi Province, China, based on seven female specimens. Morphologically, the new species can be distinguished from its congeners by the following combination of characters: (1) dorsum iridescent and uniformly dark brown in preservative, longitudinal vertebral line absent; (2) light brown beneath; (3) dotted black streak in the middle of the subcaudals absent; (4) tail length relatively short, TaL/ToL 12–16%, in females; (5) fewer subcaudals, 41–46 in females; (6) dorsal scales 23 rows throughout, strongly keeled; the outer-most rows on both sides are also keeled and slightly enlarged; (7) one loreal; (8) internasal not fused to prefrontal; (9) suture between internasals is similar in size when compared to the suture between prefrontals; (10) preocular and postocular absent; (11) 6 supralabials; (12) 5 infralabials, the first 3 (rarely 2) touching the first pair of chin shields; (13) 3 pairs of chin shields. A phylogenetic analysis using mitochondrial gene, cytochrome c oxidase subunit 1 (CO1), showed that this new species forms a monophyletic group with strong support. In addition, the uncorrected p-distances between the new species and other known congeners ranged from 3.4% to 13.0%. The recognition of the new species increases the number of described Achalinus species to 21.  

Taxonomical review of Prosymnaangolensis Boulenger, 1915 (Elapoidea, Prosymnidae) with the description of two new species

African Shovel-snout snakes (Prosymna Gray, 1849) are small, semi-fossorial snakes with a unique compressed and beak-like snout. Prosymna occur mainly in the savanna of sub-Saharan Africa. Of the 16 currently recognised species, four occur in Angola: Prosymnaambigua Bocage, 1873, P.angolensis Boulenger, 1915, P.frontalis (Peters, 1867), and P.visseri FitzSimons, 1959. The taxonomical status and evolutionary relationships of P.angolensis have never been assessed due to the lack of genetic material. This species is known to occur from western Angola southwards to Namibia, and eastwards to Zambia, Botswana and Zimbabwe. The species shows considerable variation in dorsal colouration across its range, and with the lower ventral scales count, an 'eastern race' was suggested. In recent years, Prosymna material from different parts of Angola has been collected, and with phylogenetic analysis and High Resolution X-ray Computed Tomography, the taxonomic status of these populations can be reviewed. Strong phylogenetic evidence was found to include the angolensis subgroup as part of the larger sundevalli group, and the existence of three phylogenetic lineages within the angolensis subgroup were identified, which each exhibit clear morphological and colouration differences. One of these lineages is assigned to the nominotypical P.angolensis and the other two described as new species, one of which corroborates the distinct eastern population previously detected. These results reinforce that a considerable part of Angolan herpetological diversity is still to be described and the need for further studies.

Evolution of the rostral scale and mimicry in the genus Xenodon Boie, 1826 (Serpentes: Dipsadidae: Xenodontinae)

Snakes are a stimulating life form from an evolutionary perspective. Despite the basic morphological body shape (limbless, with a tubular body), these vertebrates are extremely diverse. The Neotropical region is one of the most diverse regions for snakes in the world, with >650 known species. Within this great diversity, the genus Xenodon includes 12 species with interesting adaptations to terrestrial and semi-fossorial habitats. Members of this genus are mostly diurnal and terrestrial, feed mainly on anurans and exhibit Batesian mimicry of venomous snakes of the genera Bothrops or Micrurus. Here, through phylogenetic analysis and ancestral state estimation, we explore the evolution of the rostral scale and mimicry within the genus Xenodon. Our results suggest that the ancestral lineage of Xenodon had a rounded rostral scale and exhibited Bothrops mimicry. The evolution of the rostral scale in Xenodon might be related to abiotic factors, as an adaptation for open and forested habitats, and mimicry is likely to be related to biotic factors, as a defensive strategy resembling those of venomous snakes.

Comprehensive taxon sampling and vetted fossils help clarify the time tree of shorebirds (Aves, Charadriiformes)

Shorebirds (Charadriiformes) are a globally distributed clade of modern birds and, due to their ecological and morphological disparity, a frequent subject of comparative studies. While molecular phylogenies have been key to establishing the suprafamilial backbone of the charadriiform tree, a number of relationships at both deep and shallow taxonomic levels remain poorly resolved. The timescale of shorebird evolution also remains uncertain as a result of extensive disagreements among the published divergence dating studies, stemming largely from different choices of fossil calibrations. Here, we present the most comprehensive non-supertree phylogeny of shorebirds to date, based on a total-evidence dataset comprising 353 ingroup taxa (90% of all extant or recently extinct species), 27 loci (15 mitochondrial and 12 nuclear), and 69 morphological characters. We further clarify the timeline of charadriiform evolution by time-scaling this phylogeny using a set of 14 up-to-date and thoroughly vetted fossil calibrations. In addition, we assemble a taxonomically restricted 100-locus dataset specifically designed to resolve outstanding problems in higher-level charadriiform phylogeny. In terms of tree topology, our results are largely congruent with previous studies but indicate that some of the conflicts among earlier analyses reflect a genuine signal of pervasive gene tree discordance. Monophyly of the plovers (Charadriidae), the position of the ibisbill (Ibidorhyncha), and the relationships among the five subfamilies of the gulls (Laridae) could not be resolved even with greatly increased locus and taxon sampling. Moreover, several localized regions of uncertainty persist in shallower parts of the tree, including the interrelationships of the true auks (Alcinae) and anarhynchine plovers. Our node-dating and macroevolutionary rate analyses find support for a Paleocene origin of crown-group shorebirds, as well as exceptionally rapid recent radiations of Old World oystercatchers (Haematopodidae) and select genera of gulls. Our study underscores the challenges involved in estimating a comprehensively sampled and carefully calibrated time tree for a diverse avian clade, and highlights areas in need of further research.

Simultaneously collecting coding and non-coding phylogenomic data using homemade full-length cDNA probes, tested by resolving the high-level relationships of Colubridae

Resolving intractable phylogenetic relationships often requires simultaneously analyzing a large number of coding and non-coding orthologous loci. To gather both coding and non-coding data, traditional sequence capture methods require custom-designed commercial probes. Here, we present a cost-effective sequence capture method based on homemade probes, to capture thousands of coding and non-coding orthologous loci simultaneously, suitable for all organisms. This approach, called “FLc-Capture,” synthesizes biotinylated full-length cDNAs from mRNA as capture probes, eliminates the need for costly commercial probe design and synthesis. To demonstrate the utility of FLc-Capture, we prepared full-length cDNA probes from mRNA extracted from a common colubrid snake. We performed capture experiments with these homemade cDNA probes and successfully obtained thousands of coding and non-coding genomic loci from 24 Colubridae species and 12 distantly related snake species of other families. The average capture specificity of FLc-Capture across all tested snake species is 35%, similar to the previously published EecSeq method. We constructed two phylogenomic data sets, one including 1,075 coding loci (∼817,000 bp) and the other including 1,948 non-coding loci (∼1,114,000 bp), to study the phylogeny of Colubridae. Both data sets yielded highly similar and well-resolved trees, with 85% of nodes having >95% bootstrap support. Our experimental tests show that FLc-Capture is a flexible, fast, and cost-effective sequence capture approach for simultaneously gathering coding and non-coding phylogenomic data sets to study intractable phylogenetic questions. We hope that this method will serve as a new data collection tool for evolutionary biologists working in the era of phylogenomics.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a “refugium” for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

VenomMaps: Updated species distribution maps and models for New World pitvipers (Viperidae: Crotalinae)

Beyond providing critical information to biologists, species distributions are useful for naturalists, curious citizens, and applied disciplines including conservation planning and medical intervention. Venomous snakes are one group that highlight the importance of having accurate information given their cosmopolitan distribution and medical significance. Envenomation by snakebite is considered a neglected tropical disease by the World Health Organization and venomous snake distributions are used to assess vulnerability to snakebite based on species occurrence and antivenom/healthcare accessibility. However, recent studies highlighted the need for updated fine-scale distributions of venomous snakes. Pitvipers (Viperidae: Crotalinae) are responsible for >98% of snakebites in the New World. Therefore, to begin to address the need for updated fine-scale distributions, we created VenomMaps, a database and web application containing updated distribution maps and species distribution models for all species of New World pitvipers. With these distributions, biologists can better understand the biogeography and conservation status of this group, researchers can better assess vulnerability to snakebite, and medical professionals can easily discern species found in their area.

Measurement(s)	Species Distributions
Technology Type(s)	Geographic Information System • Species Distribution Model (MaxEnt/kuenm)
Factor Type(s)	Occurrence Records • Environmental Data
Sample Characteristic - Organism	Crotalinae
Sample Characteristic - Location	North America • South America

Australian Sea Snake Envenoming Causes Myotoxicity and Non-Specific Systemic Symptoms - Australian Snakebite Project (ASP-24)

Background: Sea snakes are venomous snakes found in the warm parts of the Indo-Pacific, including around Australia. Most sea snake envenoming causes myotoxicity, but previous Australian case reports describe neurotoxicity. We aimed to describe the epidemiology and clinical presentation of Australian sea snake envenoming and the effectiveness of antivenom.

Methods: Patients were recruited to the Australian Snakebite Project (ASP), an Australia-wide prospective observational study recruiting all patients with suspected or confirmed snakebite >2 years. Information about demographics, bite circumstances, species involved, clinical and laboratory features of envenoming, and treatment is collected and entered into a purpose-built database.

Results: Between January 2002 and August 2020, 13 patients with suspected sea snake bite were recruited to ASP, 11 were male; median age was 30 years. Bites occurred in Queensland and Western Australia. All patients were in or around, coastal waters at the time of bite. The species involved was identified in two cases (both Hydrophis zweifeli). Local effects occurred in 9 patients: pain (5), swelling (5), bleeding (2), bruising (1). Envenoming occurred in eight patients and was characterised by non-specific systemic features (6) and myotoxicity (2). Myotoxicity was severe (peak CK 28200 and 48100 U/L) and rapid in onset (time to peak CK 13.5 and 15.1 h) in these two patients. Non-specific systemic features included nausea (6), headache (6), abdominal pain (3), and diaphoresis (2). Leukocytosis, neutrophilia, and lymphopenia occurred in both patients with myotoxicity and was evident on the first blood test. No patients developed neurotoxicity or coagulopathy. Early Seqirus antivenom therapy was associated with a lower peak creatine kinase.

Conclusion: While relatively rare, sea snake envenoming is associated with significant morbidity and risk of mortality. Early antivenom appears to have a role in preventing severe myotoxicity and should be a goal of therapy.

Venom production and secretion in reptiles

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems.

First occurrence of Psammophis (Serpentes) from Europe witnesses another Messinian herpetofaunal dispersal from Africa - biogeographic implications and a discussion of the vertebral morphology of psammophiid snakes

We here describe abundant new snake material from the late Miocene (MN 13) of Salobreña, Spain. Vertebral morphology suggests a referral of the specimens to the extant psammophiid Psammophis, documenting the first occurrence of this genus in Europe. The diversity and disparity across the vertebral morphology of different psammophiid genera is discussed. We identify vertebral features that could diagnose Psammophis and therefore enable the recognition of the genus in the fossil record. A comparison of the new Spanish form with other taxa is conducted. We provide a detailed review of the psammophiid fossil record. Material previously described from the middle Miocene of Beni Mellal, Morocco is here tentatively referred to as? Psammophis sp., an action that renders that occurrence as the oldest (probable) record of the genus and Psammophiidae as a whole, providing thus a potential calibration point. On the other hand, Eastern European Pliocene material that had been previously supposedly referred to Psammophis is here discarded as being rather fragmentary, not affording any more precise determination. The two psammophiid genera Psammophis and Malpolon appear almost simultaneously in the European fossil record (MN 13), with the former achieving only a short-lived and apparently geographically limited distribution in the continent, while the latter still exists in its modern herpetofauna. We assess biogeographic implications of the new find, suggesting a direct dispersal event from northwestern Africa to the Iberian Peninsula during the late Miocene, facilitated by the Messinian Salinity Crisis. This article is protected by copyright. All rights reserved.

Isolated by dry lands: integrative analyses unveil the existence of a new species and a previously unknown evolutionary lineage of Brazilian Lanceheads (Serpentes: Viperidae: Bothrops) from a Caatinga moist-forest enclave

Here we describe a new species of Lancehead (Bothrops jabrensis sp. nov.) based on three individuals sampled from a previously unknown population from Pico do Jabre, an isolated and small Caatinga moist-forest enclave (CMFE) located in northeastern Brazil. Although this new species has an external morphology resembling those found in representatives of the Bothrops jararaca (Wied-Neuwied, 1824) species group, B. jabrensis can be diagnosed by a combination of meristic and color characters. Molecular phylogenetic analysis indicates the new species represents a unique and highly divergent lineage within Bothrops revealing the existence of a previously unknown phylogenetic lineage that has been evolving as an independent unit for more than 8 million years. Additionally, the estimated divergence time of this lineage conflicts with some proposed scenarios of historical processes associated with the evolution of CMFEs. Finally, the uniqueness of this species indicates its relevance for the maintenance of the phylogenetic diversity of Lanceheads in South America. Like other CMFEs, Pico do Jabre is consistently threatened by poaching, illegal fires, deforestation for agricultural purposes, and illegal logging. The restricted distribution of B. jabrensis, in a small and disturbed CMFE, strongly suggests that this species is critically endangered and is likely approaching extinction as a natural population.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

A “hotspot” within a hotspot: the reptiles of the Estação Ecológica and Área de Proteção Ambiental de Murici, Atlantic Forest of northeastern Brazil

Abstract: Currently the Atlantic Forest hotspot has less than 11% of its original coverage. However approximately 300 species of reptiles are known to inhabit this ecoregion, of which 34% are endemic. The creation of protected areas represents a strategy for preserving ecosystems and managing land use, and for attaining the proper management of these protected areas, information on local biodiversity is essential. Herein we provide the first list of reptile species for the Estação Ecológica and Área de Proteção Ambiental de Murici, two overlapping protected areas located in the state of Alagoas, one of the most important sets of forest remnants for the conservation of the Atlantic Forest in the northeast of Brazil. The species list was constructed based on expeditions and occasional encounters in the area between 1994 and 2022. A total of 89 reptile species were obtained during the 28 years of collection, being two species of Crocodylia, three species of Testudines and 84 species of Squamata. This richness is by far one of the greatest ever documented for the Atlantic Forest. Additionally, two species registered in the area are considered threatened according to the national list (Amerotyphlops paucisquamus and Bothrops muriciensis) and six are defined as data deficient for the assessment of their conservation status. Three species are recorded for the first time in the Atlantic Forest north of the São Francisco River: Dipsas indica, Trilepida salgueiroi and Cercophis auratus. We claim that the species list provided here will serve as a starting point for further studies in this rich reptile “hotspot” within the Atlantic Forest.

Questionnaire Background on the Hognose Snake (Heterodon nasicus) Breeders' Awareness of the Risk of Being Bitten

Simple Summary

Heterodon nasicus is a popular pet reptile. It represents dipsadine colubrids, non-front-fanged snakes (NFFS) having Duvernoy’s glands situated in the temporal region, which deliver toxic saliva via modified dentition during swallowing even though they are generally not viewed as venomous animals. As shown by previous descriptions, handlers were bitten while feeding their snakes with bare hands or when the smell of food persisted on their hands after feeding. The results of a survey for hognose snake breeders showed that 47.7% (n = 41) of the respondents admitted that they had had a situation in which H. nasicus tried to bite them; it occurred during routine approaches to daily snake handling. These situations are definitely preventable using proper feeding methods and by avoidance of routine risky habits.

Abstract

Several cases of Heterodon nasicus bites producing various symptoms have already been described. In this paper, a survey was conducted among Polish breeders of this species, where 47.7% of the respondents admitted that their snakes had tried to bite them, and 31.4% fed their snakes with their bare hands. When asked whether they feared being bitten by this snake during daily handling, an overwhelming majority of the respondents declared (74% of total negative responses, n = 63) that they were not afraid of being bitten. Based on the popularity of H. nasicus as a pet animal, it can only be assumed that the magnitude of accidental bites is much greater than that reported in published case reports. Therefore, it is important to try to educate pet breeders to report such cases and pay attention to the use of appropriate handling methods during contact with snake food and during feeding.

Multilocus phylogeny, natural history traits and classification of natricine snakes (Serpentes: Natricinae)

Natricine snakes are geographically widespread, species rich (with ~250 extant species) and both morphologically and ecologically diverse. We present a multilocus DNA sequence phylogeny for 249 natricine specimens representing 189 named species, including 69 specimens and 21 species not previously sampled. Our inferred Bayesian and maximum likelihood trees form the basis for evaluations of genus-level classification, historical biogeography, lineage diversification, and dietary, habit and reproductive-mode diversity and evolution, although several, mostly deeper, relationships remain poorly resolved. The optimal trees support natricine origins in Asia, with dispersals to Australo-Melanesia, sub-Saharan Africa (including Seychelles Archipelago, excluding Aldabra), Europe and North Africa and into North and Central America. Viviparity appears to have evolved independently three times in Natricinae but was not significantly associated with an aquatic habit. We found limited associations between habit and diet categories. We propose generic reallocations for four natricine species and highlight other points of uncertainty in natricine classification.