A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

Comparative electron microscopy study of spermatozoa in snakes (Lepidosauria, Squamata)

The ultrastructure of snake sperm has received substantial attention primarily because snakes exhibit considerable variability in reproductive characteristics between species, with a wide range of mating systems and reproductive behaviors. Variability of sperm morphology among snake species may be associated with the reproductive strategies of each taxon, such as competition or sperm storage. We provide a detailed description of the sperm ultrastructure of nine snake species (Anilius scytale, Tropidophis paucisquamis, Bothrops jararaca, Oxyrhopus guibei, Dipsas mikanii, Micrurus corallinus, Xenopholis scalaris, Acrochordus javanicus, and Cylindrophis ruffus) and compared this with sperm data from the literature for the following taxa: Liotyphlops beui, Amerotyphlops reticulatus, Trilepida koppesi, Anilios waitii, Anilios endoterus, Aspidites melanochephalus, Boa constrictor amarali, Corallus hortulana, Epicrates cenchria, Boa constrictor occidentalis, Eryx jayakari, Micrurus corallinus, Micrurus surinamensis, Micrurus frontalis, Micrurus altirostris, Oxyuranus microlepidotus, Bothrops alternatus, Bothrops diporus, Crotalus durissus, Agkistrodon contortrix, Vipera aspis, Boiga irregularis, Zamenis schrenckii, Zamenis scalaris, Stegonotus cuculatus, Nerodia sipedon, Liodytes pygaea, and Myrrophis chinensis. We found twelve polymorphic characters in the ultrastructure of sperm among the described snakes. Our work supports the importance of ultrastructural analysis of sperm morphology to understand snake reproduction, and provides sperm-derived morphological characters for phylogenetic analysis.

Pivotal temperature is not for everyone: Evidence for temperature-dependent sex determination in three gecko species

The prevalence of environmental sex determination (ESD) in squamate reptiles is often overestimated in the literature. This is surprising because we have reliable data demonstrating ESD in only a few species. The documentation of ESD in three species of geckos presented here has significantly increased our knowledge, given that satisfactory evidence for ESD existed in only eight other gecko species. For the first time, we document the occurrence of ESD in the family Sphaerodactylidae. Our finding of unexpected variability in the shapes of reaction norms among geckos highlights that traditional descriptions using parameters such as pivotal temperature, that is, temperature producing a 50:50 sex ratio, are unsatisfactory. For example, the gecko Pachydactylus tigrinus lacks any pivotal temperature and its sex ratios are strongly female-biased across the entire range of viable temperatures. We argue for the effective capture of the relationship between temperature and sex ratio using specific nonlinear models rather than using classical simplistic descriptions and classifications of reaction norms.

Differentiated sex chromosomes, karyotype evolution, and spontaneous triploidy in carphodactylid geckos

Geckos exhibit derived karyotypes without a clear distinction between macrochromosomes and microchromosomes and intriguing diversity in sex determination mechanisms. We conducted cytogenetic analyses in six species from the genera Nephrurus, Phyllurus, and Saltuarius of the gecko family Carphodactylidae. We confirmed the presence of a female heterogametic system with markedly differentiated and heteromorphic sex chromosomes in all examined species, typically with the W chromosome notably larger than the Z chromosome. One species, Nephrurus cinctus, possesses unusual multiple Z1Z1Z2Z2/Z1Z2W sex chromosomes. The morphology of the sex chromosomes, along with repetitive DNA content, suggests that the differentiation or emergence of sex chromosomes occurred independently in the genus Phyllurus. Furthermore, our study unveils a case of spontaneous triploidy in a fully grown individual of Saltuarius cornutus (3n = 57) and explores its implications for reproduction in carphodactylid geckos. We revealed that most carphodactylids retain the putative ancestral gekkotan karyotype of 2n = 38, characterized by predominantly acrocentric chromosomes that gradually decrease in size. If present, biarmed chromosomes emerge through pericentric inversions, maintaining the chromosome (and centromere) numbers. However, Phyllurus platurus is a notable exception, with a karyotype of 2n = 22 chromosomes. Its eight pairs of biarmed chromosomes were probably formed by Robertsonian fusions of acrocentric chromosomes. The family underscores a remarkable instance of evolutionary stability in chromosome numbers, followed by a profound transformation through parallel interchromosomal rearrangements. Our study highlights the need to continue generating cytogenetic data in order to test long-standing ideas about reproductive biology and the evolution of genome and sex determination.

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.

Comparative anatomy and evolution of the atlantoaxial complex in the fossorial lineage Amphisbaenia (Squamata: Lacertoidea)

The atlas and axis are the first two vertebrae from the cervical series; these two vertebrae are responsible for neck flexion, extension, and rotation movements, while providing insertion points for muscles and tendons. Amphisbaenia is a group of fossorial squamates known for having four distinctive head shapes, which are related to different excavation methods. However, little is known about the relationship between these different digging patterns and the anatomy and evolution of the atlantoaxial complex. In this study, we used computed microtomography data to describe in detail of the atlantoaxial complex for 15 species, belonging to all six current families of Amphisbaenia. Furthermore, we evaluate evolutionary scenarios of selected characters related to the atlantoaxial complex in the most recent phylogeny for Amphisbaenia, using the criteria of parsimony and maximum likelihood. Our results indicate that the evolutionary pattern of the atlantoaxial complex presents a diversification in its morphology that is not always correlated with the shape of the head. This analysis reinforces the hypothesis of remarkable morphological convergences in the evolutionary history of Amphisbaenia. Additionally, some of the characters studied may represent independent evolution through convergence in some cases (e.g., horizontal axis of the neural column) and parallelism in others (e.g., present or absent from the transverse process).

Placental ontogeny in the Yucca Night Lizard, Xantusia vigilis

Squamate placentas support physiological exchange between mothers and embryos. Uterine and embryonic epithelial cells provide sites for transporting mechanisms and extraembryonic membranes provide the scaffolding for embryonic epithelial cells and vascular systems. Diversity in placental structure involves variation in extraembryonic membrane development as well as epithelial cell specializations. Variation in placental ontogeny is known to occur and, although lineage specific patterns have been described, phylogenetic distribution of specific patterns is poorly understood. Xantusia vigilis is a viviparous lizard in a monophyletic clade, Xantusiidae, of predominantly viviparous species. Xantusiidae is one of two viviparous lineages within the clade Scincoidea that provides an important outgroup comparison for Scincidae, which includes the largest number of independent origins of viviparity among Squamata. Previous reports contain brief descriptions of placental structure of X vigilis but the developmental pattern is unknown including relevant details for comparison with skinks. We studied placental ontogeny in X. vigilis to address two hypotheses: (1) the pattern of development of placental architecture is similar to species of Scincidae and, (2) placental epithelial cell specializations are similar to species of Scincidae. The terminal placental stage of X. vigilis is similar to skinks in that it includes a chorioallantoic placenta and an omphaloplacenta. The chorioallantoic placenta is richly vascularized with thin, squamous epithelial cells separating the two vascular systems. This morphology differs from the elaborate epithelial cell specializations as occur in some skink species, but is similar to many species. Epithelial cells of the omphaloplacenta are enlarged, as they are in scincids, yet development of the omphaloplacenta includes a vascular pattern known to occur only in gerrhonotine lizards. Histochemical staining properties of the epithelium of the omphalopleure of the omphaloplacenta indicate the potential for protein transport, a function not previously reported for lizards.

Cytogenetic Analysis of Satellitome of Madagascar Leaf-Tailed Geckos

Satellite DNA (satDNA) consists of sequences of DNA that form tandem repetitions across the genome, and it is notorious for its diversity and fast evolutionary rate. Despite its importance, satDNA has been only sporadically studied in reptile lineages. Here, we sequenced genomic DNA and PCR-amplified microdissected W chromosomes on the Illumina platform in order to characterize the monomers of satDNA from the Henkel's leaf-tailed gecko U. henkeli and to compare their topology by in situ hybridization in the karyotypes of the closely related Günther's flat-tail gecko U. guentheri and gold dust day gecko P. laticauda. We identified seventeen different satDNAs; twelve of them seem to accumulate in centromeres, telomeres and/or the W chromosome. Notably, centromeric and telomeric regions seem to share similar types of satDNAs, and we found two that seem to accumulate at both edges of all chromosomes in all three species. We speculate that the long-term stability of all-acrocentric karyotypes in geckos might be explained from the presence of specific satDNAs at the centromeric regions that are strong meiotic drivers, a hypothesis that should be further tested.

Tissue damaging toxins in snake venoms: mechanisms of action, pathophysiology and treatment strategies

Snakebite envenoming is an important public health issue responsible for mortality and severe morbidity. Where mortality is mainly caused by venom toxins that induce cardiovascular disturbances, neurotoxicity, and acute kidney injury, morbidity is caused by toxins that directly or indirectly destroy cells and degrade the extracellular matrix. These are referred to as 'tissue-damaging toxins' and have previously been classified in various ways, most of which are based on the tissues being affected (e.g., cardiotoxins, myotoxins). This categorisation, however, is primarily phenomenological and not mechanistic. In this review, we propose an alternative way of classifying cytotoxins based on their mechanistic effects rather than using a description that is organ- or tissue-based. The mechanisms of toxin-induced tissue damage and their clinical implications are discussed. This review contributes to our understanding of fundamental biological processes associated with snakebite envenoming, which may pave the way for a knowledge-based search for novel therapeutic options.

Cranial anatomy of the "round-headed" Amphisbaenian Zygaspis quadrifrons (Squamata, Amphisbaenia) based on high-resolution x-ray computed tomography

Amphisbaenians are a poorly understood clade of fossorial lizards. Because of their derived anatomy and relative scarcity, the systematics of the clade and its placement within squamates has long been controversial. Traditional approaches grouped species into four assemblages according to burrowing behavior and cranial morphology, resulting in the recognition of "shovel-headed," "round-headed," "keel-headed," and "spade-headed" morphotypes. Recent phylogenetic analyses do not support the monophyly of the taxa that share those morphotypes. Detailed analyses of cranial osteology were previously accomplished using high-resolution x-ray computed tomography (HRXCT) for the "shovel-headed" Rhineura hatcherii (Rhineruidae) and the "spade-headed" Diplometopon zarudnyi (Trogonophidae). A detailed description of the "round-headed" Amphisbaena alba was previously completed based upon traditional "dry" skeletal specimens. Seven species of the "round-headed" Blanus (Blanidae) were also analyzed using HRXCT. The goal of that project was a comparative analysis of all extant species of Blanus rather than a detailed, bone-by-bone description of one species, but certainly is useful for comparison with another "round-headed" taxon. The "round-headed" morphotype is by far the most common among amphisbaenians and is much in need of further documentation. We use HRXCT imagery to provide additional data about the disparity in cranial morphology among amphisbaenians. Those data allow us to provide another detailed description of a "round-headed" amphisbaenian, the poorly known southern African species Zygaspis quadrifrons. HRXCT is ideal for this relatively rare and diminutive species. We are able to visualize and describe a detailed reconstruction of the entire skull as well as individual cranial elements. Comparisons with other species that were described in similar detail-D. zarudnyi, Spathorhynchus fossorium, R. hatcherii, and A. alba-and to a lesser degree with Blanus, reveal a complex mosaic of morphological features of the skull in Zygaspis. Preliminary data suggest that intraspecific variation is present within Z. quadrifrons, and interspecific variation among other species of Zygaspis may be sufficient for species-level recognition based on cranial osteology. Our description is, therefore, also intended to serve as a baseline for comparative analysis of other specimens of Z. quadrifrons and of other species within the genus.

Heteromorphic ZZ/ZW sex chromosomes sharing gene content with mammalian XX/XY are conserved in Madagascan chameleons of the genus Furcifer

Chameleons are well-known lizards with unique morphology and physiology, but their sex determination has remained poorly studied. Madagascan chameleons of the genus Furcifer have cytogenetically distinct Z and W sex chromosomes and occasionally Z1Z1Z2Z2/Z1Z2W multiple neo-sex chromosomes. To identify the gene content of their sex chromosomes, we microdissected and sequenced the sex chromosomes of F. oustaleti (ZZ/ZW) and F. pardalis (Z1Z1Z2Z2/Z1Z2W). In addition, we sequenced the genomes of a male and a female of F. lateralis (ZZ/ZW) and F. pardalis and performed a comparative coverage analysis between the sexes. Despite the notable heteromorphy and distinctiveness in heterochromatin content, the Z and W sex chromosomes share approximately 90% of their gene content. This finding demonstrates poor correlation of the degree of differentiation of sex chromosomes at the cytogenetic and gene level. The test of homology based on the comparison of gene copy number variation revealed that female heterogamety with differentiated sex chromosomes remained stable in the genus Furcifer for at least 20 million years. These chameleons co-opted for the role of sex chromosomes the same genomic region as viviparous mammals, lacertids and geckos of the genus Paroedura, which makes these groups excellent model for studies of convergent and divergent evolution of sex chromosomes.

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.

Habu snakes (Protobothrops flavoviridis) show variation in thoracic aortic vasoreactivity between adjacent Japanese islands

Habu snakes (Protobothrops flavoviridis) are pit vipers found in the geographically adjacent but ecologically divergent islands of Tokunoshima and Amami-Oshima in southwestern Japan. Abiotic factors can cause variation in animal populations between the two islands, and Habu snakes may show such intraspecific physiological variation. We therefore evaluated the vasoreactivity in aortas isolated from the Habu of both islands. Tokunoshima Habu showed significantly greater contractile responses to angiotensin (Ang) II, acetylcholine (ACh) and noradrenaline, and significantly higher affinities (pEC50) for Ang II and ACh, than Amami-Oshima Habu. ACh caused contractions in aortas from both populations, a finding previously unreported in snakes. Our findings indicate that vasoreactivity may differ between Tokunoshima and Amami-Oshima Habu.

Exceptionally rapid tooth development and ontogenetic changes in the feeding apparatus of the Komodo dragon

Dental developmental and replacement patterns in extinct amniotes have attracted a lot of attention. Notable among these are Paleozoic predatory synapsids, but also Mesozoic theropod dinosaurs, well known for having true ziphodonty, strongly serrated carinae with dentine cores within an enamel cap. The Komodo dragon, Varanus komodoensis, is the only extant terrestrial vertebrate to exhibit true ziphodonty, making it an ideal model organism for gaining new insights into the life history and feeding behaviours of theropod dinosaurs and early synapsids. We undertook a comparative dental histological analysis of this extant apex predator in combination with computed tomography of intact skulls. This study allowed us to reconstruct the dental morphology, ontogeny, and replacement patterns in the largest living lizard with known feeding behaviour, and apply our findings to extinct taxa where the behaviour is largely unknown. We discovered through computed tomography that V. komodoensis maintains up to five replacement teeth per tooth position, while histological analysis showed an exceptionally rapid formation of new teeth, every 40 days. Additionally, a dramatic ontogenetic shift in the dental morphology of V. komodoensis was also discovered, likely related to changes in feeding preferences and habitat. The juveniles have fewer dental specializations, lack true ziphodonty, are arboreal and feed mostly on insects, whereas the adults have strongly developed ziphodonty and are terrestrial apex predators with defleshing feeding behaviour. In addition, we found evidence that the ziphodont teeth of V. komodoensis have true ampullae (interdental folds for strengthening the serrations), similar to those found only in theropod dinosaurs. Comparisons with other species of Varanus and successive outgroup taxa reveal a complex pattern of dental features and adaptations, including the evolution of snake-like tongue flicking used for foraging for prey. However, only the Komodo dragon exhibits this remarkable set of dental innovations and specializations among squamates.

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.

First Cytogenetic Analysis of Hemidactylus mercatorius Gray, 1842 Provides Insights on Interspecific Chromosomal Diversification in the Genus Hemidactylus (Squamata: Gekkonidae)

This contribution provides the first karyotype description of Hemidactylus mercatorius and discusses the interspecific chromosome diversification in the genus. Chromosomal analysis was performed on samples from different Malagasy populations using standard karyotyping, Ag-NOR staining, and banding methods (sequential C-banding + Giemsa, + Chromomycin A3, +4',6-diamidino-2-phenylindole). Irrespective of sex or sampling locality, H. mercatorius shows a karyotype of 2n = 42 with metacentric (1, 18-21), submetacentric (4), subtelocentric (5, 11), and acrocentric pairs (all the remaining pairs). There was no heteromorphic chromosome pair and no clear distinction between macro- and microchromosomes. NORs were localised close to the centromeres of a medium acrocentric pair (14). Heterochromatic blocks were identified on the telomeric and centromeric regions of most chromosome pairs. A comparison with the karyotype of H. mabouia highlights that the different morphology of several chromosome pairs clearly distinguishes the two species, contrasting the previously proposed synonymy. The differences between the karyotypes of H. mercatorius and H. mabouia concern the number of biarmed and acrocentric elements, suggesting the occurrence of several chromosome inversions. Considering all the available karyotype data on Hemidactylus and its sister genus Cyrtodactylus, it is possible to advance an evolutionary hypothesis on their chromosomal evolution, starting from a common ancestor with 2n = 48 and all acrocentric elements. From this ancestral condition, the karyotype diversification in the two genera has been prevalently characterised by a progressive accumulation of fusions and inversions which have reduced the total chromosome count and increased the number of biarmed chromosomes.

Embryonic development of the pelvic girdle and hindlimb skeletal elements in Anilius scytale (Linnaeus, 1758) (Serpentes: Aniliidae)

Anilius scytale is the sister lineage of all other alethinophidian snakes. Morphology of the hind limb complex in adult A. scytale (Aniliidae) has been documented. We herein, for the first time, describe the embryology of the skeletal elements of its hind limb and pelvic girdle and contextualize the evolution of these structures. We identified pregnant females of A. scytale in the Herpetology Collection of the Museu Paraense Emílio Goeldi and separated 40 embryos. The embryos were sequentially staged using external and internal anatomy, collectively comprising a developmental series representing six stages. We cleared-stained one specimen of stages 31, 34, 36, and 37. Using the embryological information gleaned from A. scytale, we reinterpret evidence relating to the ossification of the pelvis and hindlimbs. In A. scytale hindlimb buds develop as transient structures that developed before Stage 30 and regresses in subsequent stages. There is no external or internal evidence of the forelimb or scapular girdle. From Stage 31 onwards the ischium, pubis, ilium, femur and zeugopodial cartilages are visible. Pubis and femur ossify towards the end of embryonic life, and cloacal spurs do not develop in the embryo. Skeletal elements of the hindlimb and pelvic girdle develop initially in the ventral zone of the cloaca-tail region. In subsequent stages the hindlimb and pelvic girdle elements migrate dorsally, with the pubis/ischium positioned medial to the ribs. A similar process may be associated with the achievement of the condition of the pelvic girdle in adults of scolecophidians, pythonids and boids.

Exceptionally preserved stomach contents of a young tyrannosaurid reveal an ontogenetic dietary shift in an iconic extinct predator

Tyrannosaurids were large carnivorous dinosaurs that underwent major changes in skull robusticity and body proportions as they grew, suggesting that they occupied different ecological niches during their life span. Although adults commonly fed on dinosaurian megaherbivores, the diet of juvenile tyrannosaurids is largely unknown. Here, we describe a remarkable specimen of a juvenile Gorgosaurus libratus that preserves the articulated hindlimbs of two yearling caenagnathid dinosaurs inside its abdominal cavity. The prey were selectively dismembered and consumed in two separate feeding events. This predator-prey association provides direct evidence of an ontogenetic dietary shift in tyrannosaurids. Juvenile individuals may have hunted small and young dinosaurs until they reached a size when, to satisfy energy requirements, they transitioned to feeding on dinosaurian megaherbivores. Tyrannosaurids occupied both mesopredator and apex predator roles during their life span, a factor that may have been key to their evolutionary success.

Turtle skull development unveils a molecular basis for amniote cranial diversity

Amniote skulls display diverse architectural patterns including remarkable variations in the number of temporal arches surrounding the upper and lower temporal fenestrae. However, the cellular and molecular basis underlying this diversification remains elusive. Turtles are a useful model to understand skull diversity due to the presence of secondarily closed temporal fenestrae and different extents of temporal emarginations (marginal reduction of dermal bones). Here, we analyzed embryos of three turtle species with varying degrees of temporal emargination and identified shared widespread coexpression of upstream osteogenic genes Msx2 and Runx2 and species-specific expression of more downstream osteogenic genes Sp7 and Sparc in the head. Further analysis of representative amniote embryos revealed differential expression patterns of osteogenic genes in the temporal region, suggesting that the spatiotemporal regulation of Msx2, Runx2, and Sp7 distinguishes the temporal skull morphology among amniotes. Moreover, the presence of Msx2- and/or Runx2-positive temporal mesenchyme with osteogenic potential may have contributed to their extremely diverse cranial morphology in reptiles.

A reference genome assembly for the continentally distributed ring-necked snake, Diadophis punctatus

Snakes in the family Colubridae include more than 2,000 currently recognized species, and comprise roughly 75% of the global snake species diversity on Earth. For such a spectacular radiation, colubrid snakes remain poorly understood ecologically and genetically. Two subfamilies, Colubrinae (788 species) and Dipsadinae (833 species), comprise the bulk of colubrid species richness. Dipsadines are a speciose and diverse group of snakes that largely inhabit Central and South America, with a handful of small-body-size genera that have invaded North America. Among them, the ring-necked snake, Diadophis punctatus, has an incredibly broad distribution with 14 subspecies. Given its continental distribution and high degree of variation in coloration, diet, feeding ecology, and behavior, the ring-necked snake is an excellent species for the study of genetic diversity and trait evolution. Within California, six subspecies form a continuously distributed "ring species" around the Central Valley, while a seventh, the regal ring-necked snake, Diadophis punctatus regalis is a disjunct outlier and Species of Special Concern in the state. Here, we report a new reference genome assembly for the San Diego ring-necked snake, D. p. similis, as part of the California Conservation Genomics Project. This assembly comprises a total of 444 scaffolds spanning 1,783 Mb and has a contig N50 of 8.0 Mb, scaffold N50 of 83 Mb, and BUSCO completeness score of 94.5%. This reference genome will be a valuable resource for studies of the taxonomy, conservation, and evolution of the ring-necked snake across its broad, continental distribution.

Investigating relationships among stress, reproduction, and immunity in three species of watersnake

Energy is a finite resource required for all physiological processes and must be allocated efficiently among essential activities to ensure fitness and survival. During the active season, adult organisms are expected to prioritize investment in reproduction over other energetically expensive processes, such as responding to immunological challenges. Furthermore, when encountering a stressor, the balance between reproduction and immunity might be disrupted in order to fuel the stress response. Because of the distinct differences in life histories across species, watersnakes provide a unique group of study in which to examine these tradeoffs. Over a two-year period, we captured three watersnake species throughout Northeast Arkansas. Animals were subjected to restraint stress and blood samples were collected throughout the acute stress response. Blood samples were used to assess innate immunity and steroid hormone concentrations. We found the peak in corticosterone concentration is season-specific, potentially because energetic reserves fluctuate with reproductive activities. We also found body condition was positively related to acute stress and negatively related to immunity. Watersnakes evidently prioritize reproduction over immunity, especially during the energetically intensive process of vitellogenesis. Energetic tradeoffs between reproduction, immunity, and the stress response are complex, and this study contributes to our understanding of energetic shifts in free-living organisms in the context of stress.

A landmarking protocol for geometric morphometric analysis of squamate endocasts

Landmark-based geometric morphometrics is widely used to study the morphology of the endocast, or internal mold of the braincase, and the diversity associated with this structure across vertebrates. Landmarks, as the basic unit of such methods, are intended to be points of correspondence, selected depending on the question at hand, whose proper definition is essential to guarantee robustness and reproducibility of results. In this study, 20 landmarks are defined to provide a framework to analyze the morphological variability in squamate endocasts. Ten species representing a cross-section of the diversity of Squamata from both phylogenetic and ecological (i.e., habitat) perspectives were considered, to select landmarks replicable throughout the entire clade, regardless of the degree of neuroanatomical resolution of the endocast. To assess the precision, accuracy, and repeatability of these newly defined landmarks, both intraobserver and interobserver error were investigated. Estimates of measurement error show that most of the landmarks established here are highly replicable, and preliminary results suggest that they capture aspects of endocast shape related to both phylogenetic and ecologic signals. This study provides a basis for further examinations of squamate endocast disparity using landmark-based geometric morphometrics.

Chromosome evolution in Iberolacerta, a genus that deviates from the standard karyotype formula of Lacertidae

This paper describes the preparation of flow-sorted chromosome paints from the Iberian Rock lizard Iberolacerta monticola, exemplifying their subsequent use in cross-species comparisons of chromosome painting. We carried out comparative analyses of chromosome evolution in the congeneric species I. galani and I. bonnali, as well as in two other species of Lacertini (Lacerta schreiberi and Timon lepidus) whose sex chromosomes were also studied through comparative genomic hybridization. Most species of Lacertini possess a diplod number of 2n = 38, with 36 acrocentric macrochromosomes and 2 microchromosomes. However, the nine species included in the genus Iberolacerta do not possess microchromosomes. Furthermore, very conspicuous differences from the standard Lacertini karyotype were observed in the three Pyrenean species of this genus, which included several biarmed metacentrics and a Z1Z2W multiple sex-chromosome system. With the possible exception of L. schreiberi, all the species of the family Lacertidae described to date appear to share homologous Z chromosomes, which date back to the last common ancestor of the whole group. We provide conclusive evidence that L. schreiberi should no longer be considered an exception to this rule, and demonstrate that the loss of microchromosomes in Iberolacerta was produced by their fusion to a middle-sized chromosome. Furthermore, we show that the multiple sex-chromosome system of the Pyrenean species of Iberolacerta originated from the fusion of the ancestral W chromosome with one of the shortest autosomes, and provide additional evidence of the fast evolution of DNA sequences linked to the W chromosome in Lacertini.

Endocast, brain, and bones: Correspondences and spatial relationships in squamates

Vertebrate endocasts are widely used in the fields of paleoneurology and comparative neuroanatomy. The validity of endocranial studies is dependent upon the extent to which an endocast reflects brain morphology. Due to the variable neuroanatomical resolution of vertebrate endocasts, direct information about the brain morphology can be sometimes difficult to assess and needs to be investigated across lineages. Here, we employ X-ray computed tomography (CT), including diffusible iodine-based contrast-enhanced CT, to qualitatively compare brains and endocasts in different species of squamates. The relative position of the squamate brain within the skull, as well as the variability that may exist in such spatial relationships, was examined to help clarify the neurological regions evidence on their endocasts. Our results indicate that squamate endocasts provide variable representation of the brain, depending on species and neuroanatomical regions. The olfactory bulbs and peduncles, cerebral hemispheres, as well as the medulla oblongata represent the most easily discernable brain regions from squamate endocasts. In contrast, the position of the optic lobes, the ventral diencephalon and the pituitary may be difficult to determine depending on species. Finally, squamate endocasts provide very limited or no information about the cerebellum. The spatial relationships revealed here between the brain and the surrounding bones may help to identify each of the endocranial region. However, as one-to-one correspondences between a bone and a specific region appear limited, the exact delimitation of these regions may remain challenging according to species. This study provides a basis for further examination and interpretation of squamate endocast disparity.

The affinities of the Late Triassic Cryptovaranoides and the age of crown squamates

Most living reptile diversity is concentrated in Squamata (lizards, including snakes), which have poorly known origins in space and time. Recently, †Cryptovaranoides microlanius from the Late Triassic of the United Kingdom was described as the oldest crown squamate. If true, this result would push back the origin of all major lizard clades by 30-65 Myr and suggest that divergence times for reptile clades estimated using genomic and morphological data are grossly inaccurate. Here, we use computed tomography scans and expanded phylogenetic datasets to re-evaluate the phylogenetic affinities of †Cryptovaranoides and other putative early squamates. We robustly reject the crown squamate affinities of †Cryptovaranoides, and instead resolve †Cryptovaranoides as a potential member of the bird and crocodylian total clade, Archosauromorpha. Bayesian total evidence dating supports a Jurassic origin of crown squamates, not Triassic as recently suggested. We highlight how features traditionally linked to lepidosaurs are in fact widespread across Triassic reptiles. Our study reaffirms the importance of critically choosing and constructing morphological datasets and appropriate taxon sampling to test the phylogenetic affinities of problematic fossils and calibrate the Tree of Life.

Functional Anatomy of the Thoracic Limb of the Komodo Dragon (Varanus komodoensis)

Since the Komodo dragon has been included on The International Union for Conservation of Nature (IUCN) Red List of Threatened Species, it is crucial to know in detail its biology as there is a limited availability of research material on these animals-mainly those who died in zoos or whose remains were found in the wild. Anatomy is essential for understanding physiology, identification of diseases, adaptations in the environment, and behavior. In this dissection study, the relationship of individual anatomical structures was analyzed, the anatomy of the active and passive movement system of the thoracic limb was described, photographs were taken, and a radiographic examination was conducted. This species has its own differences, even within closely related lizard species. Varanus komodoensis possesses triceps muscles with three heads, and the wrist is extended with additional bones for greater flexibility of the hand. The muscles of the forelimb are analogous to the hind limb; however, they differ in the mass of individual muscles, especially those predisposed to perform the most important antigravity and locomotive functions.

Description of the neuroanatomy of the brachial plexus in South American lizards. Phylogenetic implications

Few studies considered the anatomy of the nerve plexuses and musculature associated with them in ectothermic sauropsids. Based on differentiated Sudan Black B staining and conventional dissections, we describe the neuroanatomy of the brachial plexus, its main associated nerves, and muscles. For that, representatives of the genera Diplolaemus, Liolaemus, Phymaturus, and Tropidurus were selected. Based on this, potentially useful characters for phylogenetic analysis were described. Our results show that the brachial plexus can be formed by four, five, or six nerve branches. The brachial flexor trunk, circumflex, interosseous, median, radial, subscapulocoracoid, supracoracoid, and ulnar nerves were identified. Regarding the muscles innervated by the main nerves, the following muscles were identified: biceps brachii, deltoideus scapularis, latissimus dorsi, levator scapulae, pectoralis, serratus thoracis, trapezius, triceps longus caudalis, and triceps longus lateralis. Phylogenetic analyzes revealed 31 potential synapomorphies. There exists evidence that neuroanatomy studies in a phylogenetic context could provide useful information helping to elucidate the relationships between taxonomic groups.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

The Adrenal Gland of Squamata (Reptilia): A Comparative Overview

The adrenal gland is a complex endocrine organ composed of two components: a steroidogenic tissue, which produces steroid hormones, and a chromaffin tissue, which mainly produces norepinephrine and epinephrine. Through evolution, their relationships with each other changed. They begin as isolated chromaffin and steroidogenic cell aggregates, typical of fish, and end with the advanced compact gland, typical of mammals, which consists of an external steroidogenic cortical zone and an internal chromaffin medullary zone. The adrenal gland of reptiles is unique because, with few exceptions, it is near the gonads and genital ducts, and the chromaffin and steroidogenic tissues are closely associated. However, the degree of mixing is variable. For example, in Squamata, the mixing degree of chromaffin and steroidogenic tissues, their reciprocal position in the gland, and the relative quantities of norepinephrine and epinephrine secreted by the chromaffin cells are extremely variable. This variability could be related to the phylogenetic history of the species. After a brief discussion of the adrenal gland and its main functions in vertebrates, this overview will examine the general characteristics of the adrenal gland of squamates, the differences in morphology of the gland, and the possible relationships with the phylogeny of the different species.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

A New Species of Pareas (Squamata, Pareidae) from Guangxi Province, China

We described a new species of genus Pareas from Baise City, Guangxi Zhuang Autonomous Region, China, based on morphological and molecular evidence. Pareas baiseensis sp. nov. is distinguished from its congeners by the combination of (1) Yellowish-brown body colouration; (2) Frontal subhexagonal to diamond-shaped with its lateral sides converging posteriorly; (3) The anterior pair of chin shields is longer than it is broad; (4) Loreal not in contact with the eye, prefrontal in contact with the eye, two or three suboculars; (5) Rows of 15-15-15 dorsal scales, five rows of mid-dorsal scales keeled at the middle of the body, one vertebral scale row enlarged; (6) 187-191 ventrals, 89-97 subcaudals, all divided, cloacal plate single; (7) Two postocular stripes, the nuchal area forming a dark black four-pointed fork collar with the middle tines shorter than the outside tines. The genetic divergence (uncorrected p-distance) between the new species and other representatives of Pareas ranged from 13.9% to 24.4% for Cytochrome b (Cyt b) and 12.1% to 25.5% for NADH dehydrogenase subunit 4 (ND4). Phylogenetic analyses of mitochondrial DNA gene data recovered the new species from being the sister taxon to (P. boulengeri + P. chinensis) from China.

Large-scale snake genome analyses provide insights into vertebrate development

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Chromosome Diversity and Evolution of the Endemic Malagasy Velvet Geckos of the Genus Blaesodactylus (Reptilia, Gekkonidae)

We performed a molecular and phylogenetic analysis and a comparative cytogenetic study with standard karyotyping, silver staining (Ag-NOR) and sequential C-banding + Giemsa, + fluorochromes on several Blaesodactylus samples. The phylogenetic inference retrieved two main clades, the first comprises B. victori, B. microtuberculatus and B. boivini, while the second includes B. sakalava, B. antongilensis and B. ambonihazo. The available samples of B. sakalava form two different clades (here named B. sakalava clade A and clade B), which probably deserve a taxonomic re-evaluation. We found a karyological variability in Blaesodactylus in terms of chromosome number (2n = 40-42), morphology, location of NORs, and heterochromatin distribution pattern. Blaesodactylus antongilensis and B. sakalava clade A and B showed a karyotype of 2n = 40 mostly telocentric chromosomes. Pairs 1 and 6 were metacentric in B. sakalava clade A and B, while pair 1 was composed of subtelocentric/submetacentric elements in B. antongilensis. In contrast, B. boivini displayed a karyotype with 2n = 42 only telocentric chromosomes. NORs were on the first chromosome pair in B. boivini, and on the second pair in B. antongilensis. Adding our data to those available from the literature on evolutionarily related species, we highlight that the chromosome diversification in the genus probably proceeded towards a progressive reduction in the chromosome number and the formation of metacentric elements.

When lizards try out a more plant-based lifestyle: The macroevolution of mutualistic lizard-plant-interactions (Squamata: Sauria/Lacertilia)

Pollination and seed dispersal of plants by animals are key mutualistic processes for the conservation of plant diversity and ecosystem functioning. Although different animals frequently act as pollinators or seed dispersers, some species can provide both functions, so-called 'double mutualists', suggesting that the evolution of pollination and seed dispersal may be linked. Here, we assess the macroevolution of mutualistic behaviours in lizards (Lacertilia) by applying comparative methods to a phylogeny comprising 2,838 species. We found that both flower visitation (potential pollination) (recorded in 64 species [2.3% of total] across 9 families) and seed dispersal (recorded in 382 species [13,5% of total] across 26 families) have evolved repeatedly in Lacertilia. Furthermore, we found that seed dispersal activity pre-dated flower visitation and that the evolution of seed dispersal activity and flower visitation was correlated, illustrating a potential evolutionary mechanism behind the emergence of double mutualisms. Finally, we provide evidence that lineages with flower visitation or seed dispersal activity have higher diversification rates than lineages lacking these behaviours. Our study illustrates the repeated innovation of (double) mutualisms across Lacertilia and we argue that island settings may provide the ecological conditions under which (double) mutualisms persist over macroevolutionary timescales.

Hepatozoon Miller, 1908 parasites in the Colubridae snakes Clelia clelia (Daudin, 1803) and Drymarchon corais (Boie, 1827) from the Eastern Amazonia

Based on the genetic, morphological, and morphometric data of blood gamonts, we identified Hepatozoon parasites in colubrid snakes sampled in the Eastern Amazon region. Hepatozoon trigeminum was detected in the mussurana snake Clelia clelia and exhibited wide and elongated gamonts (mean dimensions: 14.25±0.65 × 4.31±0.43 μm) with an evident parasitophorous vacuole. Hepatozoon odwyerae sp. nov. was described in the indigo snake Drymarchon corais, whose gamonts have elongated and thin bodies (mean dimensions: 13.41±0.79 × 3.72±0.35 μm) with one end more tapered than the other. Phylogenetic analyses, based on the amplification of a 441 bp fragment of the 18S rRNA gene, revealed that the novel sequences of Hepatozoon spp. from our study were closely related to hemogregarine lineages found in lizards and snakes from Brazil, forming a well-supported monophyletic clade with them. The present study provides the first species description of Hepatoozon in D. corais and a new record of a host species for C. clelia using the integrated taxonomic data. We also highlight the importance of further investigations into the diversity of Hepatozoon in snakes, a rich but underestimated group of parasites, especially in the Amazonian biome.

Molecular characterization of Schellackia parasites in an urban population of sand lizards (Lacerta agilis) from Berlin, Germany

Lizards are hosts of several taxa of unicellular parasites of the phylum Apicomplexa, including Karyolysus, Schellackia, Lankesterella, and Hepatozoon. Parasite prevalence and the impact of infections on lizard biology remain largely unexplored. In this study, blood parasite infections were investigated in sand lizards (Lacerta agilis) from Berlin, Germany. Eighty-three individuals were investigated, and the detected blood parasites were identified as Schellackia sp. The combination of microscopic and molecular screening revealed a prevalence of 14.5%. Parasitemia values were low and most infections were subpatent. Phylogenetic analysis recovered a close relationship of the Schellackia parasites of this study with Schellackia sp. parasites of different Lacerta and Podarcis lizard species from Spain. Monitoring of Schellackia parasite infections in free-ranging lizards contributes to a better understanding of the distribution, diversity, and phylogenetic relationships of the neglected parasite taxon.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

The first description of dermal armour in snakes

Osteoderms, also called dermal armour, often play a role in predator defence. The presence of osteoderms is highly irregularly distributed across the squamate phylogeny and they have not been found in snakes. In this study, we searched for candidate snake species that would benefit from such armour to protect their body, focusing primarily on fossorial species with defensive tail displays. We examined the tail morphology of 27 snake species from different families using micro-computed tomography (µCT) and micro- radiography. We discovered dermal armour in four species of sand boas (Erycidae) that also feature enlarged and highly modified caudal vertebrae. This is the first description of dermal armour in snakes. Ancestral state reconstructions revealed that osteoderms likely evolved once or multiple times in Erycidae. We have not found osteoderms in any other examined snake species. Nevertheless, similar structures are known from unrelated squamate clades, such as gerrhosaurids and geckos. This supports the idea of underlying deep developmental homology. We propose the hypothesis that osteoderms protect sand boas like the "brigandine armour" of medieval warriors. We interpret it as another component of the sand boas' rich defence strategy.

Topological analysis of vertebrate phylogenetic trees utilizing Horton's first law

Phylogenetic trees visually represent evolution and diversification. While many studies have focused on the number and length of edges (branches), topological properties, such as edge connection patterns, are also important. In this study, the topological properties of phylogenetic trees were quantified, focusing on edge connection patterns. Horton's first law was applied to quantify the overall, rather than local, topological properties of phylogenetic trees. The topological properties of vertebrate phylogenetic trees for spiny-rayed fishes, Amphibians, turtles, Squamata, Aves, and placental mammals were analyzed. The topological features discussed herein include the number of first-order edges, maximum order, and bifurcation ratio. The average bifurcation ratio of all trees was approximately 3, suggesting that phylogenetic trees for different taxa have a common mechanism of evolution. Vertebrate phylogenetic trees were compared with artificial branching objects created from neutral stochastic branching model simulations. The topological properties of the actual vertebrate phylogenetic trees agreed with those of the artificial branching objects. Our study suggests that evolutionary events do not change the overall topological properties of actual phylogenetic trees, even if the number and length of the edges change. Specifically, non-neutral events (e.g., environmental changes and mass extinction) are not main factors associated with topological properties. The results instead demonstrate a relationship between the bifurcation ratio and symmetricity in the context of temporal changes of topological properties. When the number of first-order edges increased and the maximum order remained constant, the bifurcation ratio increased and symmetricity decreased. When the number of first-order edges increased and the maximum order increased by one, the bifurcation ratio decreased and symmetricity increased.

Material properties of skin in the flying snake Chrysopelea ornata

In snakes, the skin serves for protection, camouflage, visual signaling, locomotion, and its ability to stretch facilitates large prey ingestion. The flying snakes of the genus Chrysopelea are capable of jumping and gliding through the air, requiring additional functional demands: its skin must accommodate stretch in multiple directions during gliding and, perhaps more importantly, during high-speed, direct-impact landing. Is the skin of flying snakes specialized for gliding? Here, we characterized the material properties of the skin of Chrysopelea ornata and compared them with two nongliding species of colubrid snakes, Thamnophis sirtalis and Pantherophis guttatus, as well as with previously published values. The skin was examined using uniaxial tensile testing to measure stresses, and digital image correlation methods to determine strains, yielding metrics of strength, elastic modulus, strain energy, and extensibility. To test for loading orientation effects, specimens were tested from three orientations relative to the snake's long axis: lateral, circumferential, and ventral. Specimens were taken from two regions of the body, pre- and pos-tpyloric, to test for regional effects related to the ingestion of large prey. In comparison with T. sirtalis and P. guttatus, C. ornata exhibited higher post-pyloric and lower pre-pyloric extensibility in circumferential specimens. However, overall there were few differences in skin material properties of C. ornata compared to other species, both within and across studies, suggesting that the skin of flying snakes is not specialized for gliding locomotion. Surprisingly, circumferential specimens demonstrated lower strength and extensibility in pre-pyloric skin, suggesting less regional specialization related to large prey.

A real data-driven simulation strategy to select an imputation method for mixed-type trait data

Missing observations in trait datasets pose an obstacle for analyses in myriad biological disciplines. Considering the mixed results of imputation, the wide variety of available methods, and the varied structure of real trait datasets, a framework for selecting a suitable imputation method is advantageous. We invoked a real data-driven simulation strategy to select an imputation method for a given mixed-type (categorical, count, continuous) target dataset. Candidate methods included mean/mode imputation, k-nearest neighbour, random forests, and multivariate imputation by chained equations (MICE). Using a trait dataset of squamates (lizards and amphisbaenians; order: Squamata) as a target dataset, a complete-case dataset consisting of species with nearly complete information was formed for the imputation method selection. Missing data were induced by removing values from this dataset under different missingness mechanisms: missing completely at random (MCAR), missing at random (MAR), and missing not at random (MNAR). For each method, combinations with and without phylogenetic information from single gene (nuclear and mitochondrial) or multigene trees were used to impute the missing values for five numerical and two categorical traits. The performances of the methods were evaluated under each missing mechanism by determining the mean squared error and proportion falsely classified rates for numerical and categorical traits, respectively. A random forest method supplemented with a nuclear-derived phylogeny resulted in the lowest error rates for the majority of traits, and this method was used to impute missing values in the original dataset. Data with imputed values better reflected the characteristics and distributions of the original data compared to complete-case data. However, caution should be taken when imputing trait data as phylogeny did not always improve performance for every trait and in every scenario. Ultimately, these results support the use of a real data-driven simulation strategy for selecting a suitable imputation method for a given mixed-type trait dataset.

DNA barcoding of Chinese snakes reveals hidden diversity and conservation needs

DNA barcoding has greatly facilitated studies of taxonomy, biodiversity, biological conservation, and ecology. Here, we establish a reliable DNA barcoding library for Chinese snakes, unveiling hidden diversity with implications for taxonomy, and provide a standardized tool for conservation management. Our comprehensive study includes 1638 cytochrome c oxidase subunit I (COI) sequences from Chinese snakes that correspond to 17 families, 65 genera, 228 named species (80.6% of named species) and 36 candidate species. A barcode gap analysis reveals gaps, where all nearest neighbour distances exceed maximum intraspecific distances, in 217 named species and all candidate species. Three species-delimitation methods (ABGD, sGMYC, and sPTP) recover 320 operational taxonomic units (OTUs), of which 192 OTUs correspond to named and candidate species. Twenty-eight other named species share OTUs, such as Azemiops feae and A. kharini, Gloydius halys, G. shedaoensis, and G. intermedius, and Bungarus multicinctus and B. candidus, representing inconsistencies most probably caused by imperfect taxonomy, recent and rapid speciation, weak taxonomic signal, introgressive hybridization, and/or inadequate phylogenetic signal. In contrast, 43 species and candidate species assign to two or more OTUs due to having large intraspecific distances. If most OTUs detected in this study reflect valid species, including the 36 candidate species, then 30% more species would exist than are currently recognized. Several OTU divergences associate with known biogeographic barriers, such as the Taiwan Strait. In addition to facilitating future studies, this reliable and relatively comprehensive reference database will play an important role in the future monitoring, conservation, and management of Chinese snakes.

Proteomic analysis of the mandibular glands from the Chinese crocodile lizard, Shinisaurus crocodilurus - Another venomous lizard?

Based on its phylogenetic relationship to monitor lizards (Varanidae), Gila monsters (Heloderma spp.), and the earless monitor Lanthanotus borneesis, the Chinese crocodile lizard, Shinisaurus crocodilurus, has been assigned to the Toxicofera clade, which comprises venomous reptiles. However, no data about composition and biological activities of its oral secretion have been reported. In the present study, a proteomic analysis of the mandibular gland of S. crocodilurus and, for comparison, of the herbivorous Solomon Island skink Corucia zebrata, was performed. Scanning electron microscopy (SEM) of the teeth from S. crocodilurus revealed a sharp ridge on the anterior surface, but no grooves, whereas those of C. zebrata possess a flattened crown with a pointed cusp. Proteomic analysis of their gland extracts provided no evidence of venom-derived peptides or proteins, strongly supporting the non-venomous character of these lizards. Data are available via ProteomeXchange with identifier PXD039424.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

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.

Lower-level predictors and behavioral correlates of maximal aerobic capacity and sprint speed among individual lizards

The standard paradigm of organismal biology views lower-level traits (e.g. aspects of physiology) as determining organismal performance ability (e.g. maximal sprint speed), which in turn constrains behavior (e.g. social interactions). However, few studies have simultaneously examined all three levels of organization. We used focal observations to record movement behaviors and push-up displays in the field for adult male Sceloporus occidentalis lizards during the breeding season. We then captured animals, measured aspects of their physiology, morphology and performance, and counted ectoparasites and endoparasites as potential predictors of sprint speed and maximal oxygen consumption (V̇O2,max). Field behaviors were statistically repeatable, but not strongly so. Sprint speed and V̇O2,max were repeatable using residuals from regressions on body mass (speed: r=0.70; V̇O2,max: r=0.88). Both calf [standardized partial regression (path) coefficient B=0.53] and thigh [B=-0.37] muscle mass (as residuals from regressions on body mass) were significant predictors of sprint speed; hemoglobin concentration (B=0.42) was a predictor of V̇O2,max. In turn, V̇O2,max predicted the maximum number of four-legged push-ups per bout (B=0.39). In path analysis, log likelihood ratio tests indicated no direct paths from lower-level traits to behavior, supporting the idea that morphology, in the broad sense, only affects behavior indirectly through measures of performance. Our results show that inter-individual variation in field behaviors can be related to performance ability, which in turn reflect differences in morphology and physiology, although not parasite load. Given the low repeatability of field behaviors, some of the relationships between behavior and performance may be stronger than suggested by our results.