Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species

Inside the head of Crotalus durissus LINNAEUS, 1758 (Serpentes, Viperidae, Crotalinae): Macroscopic description of the brain with ontogenetic insights

Neuroanatomy studies in vertebrates have garnered significant attention in recent years, particularly driven by advancements in computerized tomography imaging techniques. Nonetheless, these advancements remain largely constrained to specific vertebrate groups, notably mammals, birds, and fish, leaving studies in reptiles at an incipient stage. In this work, we aim to describe in detail the macroscopic morphology of the brain of Crotalus durissus based on a sample of four young and four adult individuals-three male and five female specimens, providing the first detailed description of the brain with a relatively modest sample available for reptiles. Our results show that the major macroscopic features identified in C. durissus suggest a brain structure typical of a multi-habitat and cathemeral/nocturnal alethinophidian species, thereby contributing significant data to the understanding of brain morphological evolution in snakes. Brain measurements showed distinct scaling patterns related to snout-vent length and head length, with variables such as brain length and cerebral hemisphere length decreasing with SVL, while others like olfactory bulb length and medulla oblongata width increased. Additional differences were observed comparing juveniles and adults, with adults generally exhibiting larger mean values for cerebellum and medulla oblongata measurements. However, the small sample size highlights the need for future studies with larger datasets to validate these findings and explore the developmental trajectories in greater detail.

Evolutionary History and Climatic Correlates of Hypermelanism in Viperidae

Body colorations have been investigated intensely concerning their adaptive significance from the ecological and evolutionary perspectives. Studies on melanism have caught growing interest thanks to its marked variability across space, time and taxon and, in ectotherms, it has been hypothesised to be driven by thermal advantages. Among reptiles, vipers show conspicuous inter- and intraspecific patterns of variation, making them excellent models to address evolutionary and adaptive patterns. We investigated the thermal melanism hypothesis across Viperidae by performing a phylogenetic comparative approach to assess whether its occurrence is phylogenetically driven or, alternatively, whether it is influenced by climate. Phylogenetic signal was detected and reconstructed the 'non-melanistic' form as the ancestral state at the root of their phylogeny, whereas a climatic effect was found so that melanism is more frequent in colder environments. With this work, we provide on a large geographical scale strong support for the putative advantages provided by melanism in colder climates; moreover, melanism appears to have evolved in multiple events throughout the diversification of vipers, but it has been rarely maintained over time. We hypothesise that it is maintained only when environmental conditions, such as low thermal regimes, render it favourable; nevertheless, experimental evidence is necessary to further support this hypothesis.

Iguanas rafted more than 8,000 km from North America to Fiji

Founder-event speciation can occur when one or more organisms colonize a distant, unoccupied area via long-distance dispersal, leading to the evolution of a new species lineage. Species radiations established by long-distance, and especially transoceanic, dispersal can cause substantial shifts in regional biodiversity. Here, we investigate the occurrence and timing of the greatest known long-distance oceanic dispersal event in the history of terrestrial vertebrates-the rafting of iguanas from North America to Fiji. Iguanas are large-bodied herbivores that are well-known overwater dispersers, including species that colonized the Caribbean and the Galápagos islands. However, the origin of Fijian iguanas had not been comprehensively tested. We estimated the phylogenetic relationships and evolutionary timescale of the iguanid lizard radiation using genome-wide exons and ultraconserved elements (UCEs). Those data indicate that the closest living relative of extant Fijian iguanas is the North American desert iguana and that the two taxa likely diverged during the late Paleogene near or after the onset of volcanism that produced the Fijian archipelago. Biogeographic models estimate North America as the most probable ancestral range of Fijian iguanas. Our analyses support the hypothesis that iguanas reached Fiji via an extraordinary oceanic dispersal event from western North America, and which spanned a fifth of the earth's circumference (>8,000 km). Overwater rafting of iguanas from North America to Fiji strengthens the importance of founder-event speciation in the diversification of iguanids and elucidates the scope of long-distance dispersal across terrestrial vertebrates.

Habitat openness and squamate color evolution over deep time

While the ecological roles of colored integument have been extensively studied, what regulates global patterns of color variation remains poorly understood. Here, using a global dataset of 1249 squamates, we evaluate whether and how six key eco-environmental variables and their interactions shaped the evolutionary history of their coloration. We show that only habitat openness consistently associates with brightness evolution, with brighter integuments favored in open habitats, possibly for enhanced heat reflection. Furthermore, brightness evolution rates likely track δ18O (a temperature proxy) changes and increase during global aridification phases, such as those in the Miocene and Pliocene. This trend may be due to the establishment of an arid climate that promoted habitat openness shifts, ultimately inducing adaption to new niches. Our findings suggest that a single environmental variable is associated with color variation in the largest extant tetrapod order.

Papillary and Callous Scales in the Integument of Agamid Lizards (Agamidae, Sauria) as a Phenomenon of Extraordinary Development of the Corneous Layers

Scaled integument of six species of the genus Acanthocercus and Laudakia nupta, family Agamidae was studied using light and scanning electron microscopy. Gross observation revealed the presence of two types of modified scales in the males. The enlarged scales covered with an extremely thick β-corneous layer were detected in the dorsal and ventral surfaces of the tail base and on the palmar and plantar limb surfaces of all species. After detachment of the β-layer, the surface of such scales was covered with high papillae ("papillary scales"). The callous scales were found in the precloacal region of Acanthocercus species and in both precloacal and mid-ventral regions of L. nupta. Modified scales were found in some females and subadut specimens, and absent in juveniles. A prominent papillary layer characterized the dermis of both scale types. It was assumed that well-developed dermal papillae in such scales expended the total surface area of stratum germinativum and created a pool of proliferated cells in the interpapillary loops to increase the production of differentiating keratinocytes. These processes were undoubtedly associated with the formation of a thick and resistant corneous layer that distinguished both types of scales. Functional role, a taxonomic value, and analogies with normal and pathological epidermis in birds and mammals are discussed.

A time-calibrated salamander phylogeny including 765 species and 503 genes

Recent time-calibrated amphibian phylogenies agree on the family-level relationships among extant salamanders but had disparate sampling regimes and inferred very different divergence times. For example, a recent phylogenomic study based on 220 nuclear loci had limited taxon sampling (41 species) and estimated relatively young divergence dates, whereas a more extensive supermatrix study based on 15 genes and 481 species estimated dates that were 22-45 million years older for major clades. Here, we combined phylogenomic and supermatrix approaches to estimate the largest salamander phylogeny to date based on molecular markers. Our matrix contained 765 salamander species and 503 genes (with 92.3% missing data overall). We included 284 more species than the previous largest salamander phylogeny (59% increase) and sampled approximately 93% of all currently described salamander species. Our dating analyses incorporated more than twice as many fossil calibration points within salamanders as previous studies. Maximum-likelihood estimates of tree topology yielded family-level relationships that were consistent with earlier studies. Nearly all species were placed in the expected genera, despite extensive missing data in many species. Bootstrap support was generally high across the tree but was poor in some clades where sampling of genes was limited (e.g., among some bolitoglossine salamanders). The dating analyses yielded age estimates for major clades that were generally intermediate between those from the previous phylogenomic and supermatrix analyses. We also provide a set of 200 time-calibrated trees for use in comparative analyses.

A straightforward workflow to explore species diversity using the Patagonian lizards of the Diplolaemus genus (Iguania: Leiosauridae) as a study case, with the description of a new species

Disputes over species descriptions, stemming from conceptual disparities and arbitrary species boundaries, are among the primary challenges of modern taxonomy. In this study, we introduce a straightforward workflow, grounded in evolutionary theory, designed to tackle these challenges. We exemplified this approach using Patagonian lizards from the Diplolaemus clade. This workflow involves assigning specimens to putative evolutionary lineages, conducting primary species delimitations, constructing a species tree, comparing lineages for evolutionary independence, and using post-hoc analyses to separate well-supported from ambiguous lineages. This approach aims to establish a reliable foundation for exploring the taxonomic and evolutionary diversity of challenging groups. Applying this workflow to the Diplolaemus clade, we used various analytical methods on genetic (mitochondrial and nuclear markers) and phenotypic data (meristic, linear, and geometric morphometrics). We identified ten lineages with varying degrees of evolutionary independence in a clade where only four species had been described. Among the newly identified lineages, two exhibited low support for evolutionary independence, three showed strong support but had non-conclusive information, and one was recognized and described as a new species. In summary, our hierarchical workflow not only facilitated comprehensive comparisons but also enabled us to draw robust conclusions.

Walking or hanging: the role of habitat use for body shape evolution in lacertid lizards

Differences in habitat use impose ecological constraints which in turn lead to functional and morphological differences through adaptation. In fact, a convergent evolutionary pattern is evident when species exhibit similar responses to similar environments. In this study, we examine how habitat use influences the evolution of body shape in lizards from the family Lacertidae. We divided our species into two categories: ground-dwellers and climbers, which encompasses the verticality and horizontality aspects of the habitat. We performed phylogenetic comparative analyses employing 186 species and seven linear morphological traits. Our results show contrasting patterns between head and limb shapes, which are considered distinct functional blocks. We observed differences in forelimb proportions, but not in hindlimb length, contrary to what was documented in other lizard groups, demonstrating a novel axis in the limb-locomotion-habitat relationship in this family. In addition, a clear effect of habitat use on head shape was detected. We observed that climbing species present on average flatter heads than ground-dwelling species, as well as different evolutionary trajectories. These findings suggest the complex interplay between habitat use and morphological evolution in lizards, highlighting how distinct selective pressures drive divergent adaptations in different functional traits.

Viviparity is associated with larger female size and higher sexual size dimorphism in a reproductively bimodal lizard

Squamate reptiles are central for studying phenotypic correlates of evolutionary transitions from oviparity to viviparity because these transitions are numerous, with many of them being recent. Several models of life-history theory predict that viviparity is associated with increased female size, and thus more female-biased sexual size dimorphism (SSD). Yet, the corresponding empirical evidence is overall weak and inconsistent. The lizard Zootoca vivipara, which occupies a major part of Northern Eurasia and includes four viviparous and two non-sister oviparous lineages, represents an excellent model for testing these predictions. We analysed how sex-specific body size and SSD is associated with parity mode, using body length data for nearly 14,000 adult individuals from 97 geographically distinct populations, which cover almost the entire species' range and represent all six lineages. Our analyses controlled for lineage identity, climatic seasonality (the strongest predictor of geographic body size variation in previous studies of this species) and several aspects of data heterogeneity. Parity mode, lineage and seasonality are significantly associated with female size and SSD; the first two predictors accounted for 14%-26% of the total variation each, while seasonality explained 5%-7%. Viviparous populations exhibited a larger female size than oviparous populations, with no concomitant differences in male size. The variation of male size was overall low and poorly explained by our predictors. Albeit fully expected from theory, the strong female bias of the body size differences between oviparous and viviparous populations found in Z. vivipara is not evident from available data on three other lizard systems of closely related lineages differing in parity mode. We confront this pattern with the data on female reproductive traits in the considered systems and the frequencies of evolutionary changes of parity mode in the corresponding lizard families and speculate why the life-history correlates of live-bearing in Z. vivipara are distinct. Comparing conspecific populations, our study provides the most direct evidence for the predicted effect of parity mode on adult body size but also demonstrates that the revealed pattern may not be general. This might explain why across squamates, viviparity is only weakly associated with larger size.

How life became colourful: colour vision, aposematism, sexual selection, flowers, and fruits

Plants and animals are often adorned with potentially conspicuous colours (e.g. red, yellow, orange, blue, purple). These include the dazzling colours of fruits and flowers, the brilliant warning colours of frogs, snakes, and invertebrates, and the spectacular sexually selected colours of insects, fish, birds, and lizards. Such signals are often thought to utilize pre-existing sensitivities in the receiver's visual systems. This raises the question: what was the initial function of conspicuous colouration and colour vision? Here, we review the origins of colour vision, fruit, flowers, and aposematic and sexually selected colouration. We find that aposematic colouration is widely distributed across animals but relatively young, evolving only in the last ~150 million years (Myr). Sexually selected colouration in animals appears confined to arthropods and chordates, and is also relatively young (generally <100 Myr). Colourful flowers likely evolved ~200 million years ago (Mya), whereas colourful fruits/seeds likely evolved ~300 Mya. Colour vision (sensu lato) appears to be substantially older, and likely originated ~400-500 Mya in both arthropods and chordates. Thus, colour vision may have evolved long before extant lineages with fruit, flowers, aposematism, and sexual colour signals. We also find that there appears to have been an explosion of colour within the last ~100 Myr, including >200 origins of aposematic colouration across nine animal phyla and >100 origins of sexually selected colouration among arthropods and chordates.

Conservation of OFD1 Protein Motifs: Implications for Discovery of Novel Interactors and the OFD1 Function

OFD1 is a protein involved in many cellular processes, including cilia biogenesis, mitotic spindle assembly, translation, autophagy and the repair of double-strand DNA breaks. Despite many potential interactors identified in high-throughput studies, only a few have been directly confirmed with their binding sites identified. We performed an analysis of the evolutionary conservation of the OFD1 sequence in three clades: 80 Tetrapoda, 144 Vertebrata or 26 Animalia species, and identified 59 protein-binding motifs localized in the OFD1 regions conserved in various clades. Our results indicate that OFD1 contains 14 potential post-translational modification (PTM) sites targeted by at least eight protein kinases, seven motifs bound by proteins recognizing phosphorylated aa residues and a binding site for phosphatase 2A. Moreover, OFD1 harbors both a motif that enables its phosphorylation by mitogen-activated protein kinases (MAPKs) and a specific docking site for these proteins. Generally, our results suggest that OFD1 forms a scaffold for interaction with many proteins and is tightly regulated by PTMs and ligands. Future research on OFD1 should focus on the regulation of OFD1 function and localization.

Sex Chromosome Turnovers and Stability in Snakes

For a long time, snakes were presented as a textbook example of a group with gradual differentiation of homologous ZZ/ZW sex chromosomes. However, recent advances revealed that the ZZ/ZW sex chromosomes characterize only caenophidian snakes and certain species of boas and pythons have nonhomologous XX/XY sex chromosomes. We used genome coverage analysis in four non-caenophidian species to identify their sex chromosomes, and we examined the homology of sex chromosomes across phylogenetically informative snake lineages. We identified sex chromosomes for the first time in 13 species of non-caenophidian snakes, providing much deeper insights into the evolutionary history of snake sex chromosomes. The evolution of sex chromosomes in snakes is more complex than previously thought. Snakes may have had ancestral XX/XY sex chromosomes, which are still present in a blind snake and some boas, and there were several transitions to derived XX/XY sex chromosomes with different gene content and two or even three transitions to ZZ/ZW sex chromosomes. However, we discuss more alternative scenarios. In any case, we document that (1) some genomic regions were likely repeatedly co-opted as sex chromosomes in phylogenetically distant lineages, even with opposite types of heterogamety; (2) snake lineages differ greatly in the rate of differentiation of sex chromosomes; (3) snakes likely originally possessed sex chromosomes prone to turnovers. The sex chromosomes became evolutionarily highly stable once their differentiation progressed in the megadiverse caenophidian snakes. Snakes thus provide an ideal system for studying the evolutionary factors that drive unequal rates of differentiation, turnovers and stability of sex chromosomes.

First Look at the Venoms of Two Sinomicrurus Snakes: Differences in Yield, Proteomic Profiles, and Immunorecognition by Commercial Antivenoms

Chinese coral snakes (Sinomicrurus) are highly neglected regarding their venom profiles and harm to humans, which impedes our ability to deeply understand their biological properties and explore their medicinal potential. In this study, we performed a comparative analysis to reveal the venom profiles of two Chinese coral snakes in terms of their venom yields, proteomic profiles, and immunorecognition by commercial antivenoms. The results showed that Sinomicrurus kelloggi expels more venom (lyophilized venom mass) than Sinomicrurus maccelellandi but possesses a similar solid venom content. These interspecific differences in venom yield were influenced by the snout-vent length. The venoms of these two species varied in their electrophoretic profiles, as well as in the presence or absence and relative abundance of protein families. They exhibited a 3-FTx-predominant phenotype, where the S. maccelellandi venom was dominated by 3-FTx (32.43%), SVMP (23.63%), PLA2 (19.88%), and SVSP (12.61%), while the S. kelloggi venom was dominated by 3-FTx (65.81%), LAAO (11.35%), and AMP (10.09%). While both the commercial Naja atra and Bungarus multicinctus antivenoms could immunorecognize these two Chinese coral snake venoms, the N. atra antivenom possessed a higher neutralization capability than the B. multicinctus antivenom for both species of coral snakes. Our findings show significant interspecific variations in the venom profiles of these Sinomicrurus snakes for the first time. We suggest screening or preparing specific antivenoms with high efficiency for the clinical treatment of envenomation caused by these snakes.

Heads and Tails: Comparative Osteology of Nearctic Dipsadid Snakes

Although numerous studies have addressed some aspects of the cranial osteology of Nearctic dipsadid species, only the species within the genera Heterodon and Carphophis have a formal published description of their skull. Similarly, vertebral data on such species are extremely scarce, and most of the available literature is focused on fossils. Such group has a complex phylogenetic history, being recovered as monophyletic or nonmonophyletic depending on the approach. In this paper, we provide detailed and comparative descriptions of the osteology of dipsadid species distributed in the Nearctic region based on 69 specimens of dry material and high-resolution computed tomography (CT) scans. Additionally, we explore the morphological variation of the skull and cervical vertebrae within the context of distinct phylogenetic hypotheses previously proposed. Only two suprageneric groups previously proposed shared exclusive morphological traits: (Carphophis amoenus + Contia tenuis), proposed by three studies, and (Diadophis punctatus (Ca. amoenus + Co. tenuis)), proposed by one study. Large and detailed studies on the skull, mandible, and vertebrae represent an important step toward the understanding of the evolution of species, especially when they also show intraspecific variation.

A review of age estimation methods in non-avian reptiles by growth marks in hard tissues

Age and growth-related data are basic biological parameters, essential in population ecology, evolution, and conservation biology. There is a growing body of published information on reptile demography derived from sclerochronology, a technique based on counting the growth layers deposited in bones (skeletochronology) and other hard body structures. Since the data are not always easily available, we compiled the existing published data, described the current status of knowledge, synthetized the conclusions of disparate studies, and identified patterns of research and information gaps, prioritizing the needs for future research. Our database includes the results of 468 published studies covering 236 reptile species from 41 families. These represent less than 2% of the total number of known extant species. Turtles and crocodiles are proportionally better studied, while snakes are the least examined group. The distribution of the research does not reflect conservation needs; we found an important geographic bias, with an overrepresentation of Northern temperate species. Only 23% of the studies checked the assumption of periodicity of growth marks deposition, and the method was found to be reliable or adequate in 79% of the cases. Overall, the data obtained through sclerochronology can be considered robust, especially if validation methods are employed, since the general goal is to characterize population parameters, trends, and dynamics, rather than determining the exact age of any specimen in particular.

PHYLOGENY AND SYSTEMATICS OF CYATHOCOTYLID DIGENEANS (DIGENEA: DIPLOSTOMOIDEA) PARASITIZING SNAKES WITH DESCRIPTION OF THREE NEW SPECIES OF GOGATEA FROM AUSTRALIA AND VIETNAM

The Cyathocotylidae Mühling, 1896 is a small but broadly distributed family of digeneans parasitic in a wide range of vertebrate definitive hosts, from fish to mammals. Despite existing taxonomic questions, only a few studies have generated DNA sequence data from cyathocotylids, and only a single publication included sequence data from a cyathocotylid parasitic in snakes. Four genera are known to infect snakes: GogateaLutz, 1935, Szidatia Dubois, 1938, MesostephanoidesDubois, 1951, and SerpentostephanusSudarikov, 1961. Members of these genera were known from only Asia and Africa. In the present study, we describe 2 new species of Gogatea from snakes in Australia and 1 from Vietnam. The new species from Vietnam described herein is the first member of the genus that lacks a ventral sucker. We used partial sequences of the nuclear large ribosomal subunit (28S) and mitochondrial cytochrome c oxidase subunit I (COI) genes to explore phylogenetic relationships among cyathocotylids and species differentiation. In addition, this is the first report of a cyathocotylid from snakes in Australia, thus posing interesting questions regarding the dispersal and historical biogeography of these parasites. Cyathocotylid genera from snakes have a long, convoluted taxonomic history. The genera Gogatea, Mesostephanoides, and Szidatia were distinguished from each other based on very small morphological differences. Therefore, the validity of Szidatia and Mesostephanoides was often questioned in the literature. Based on the detailed morphological analysis of our freshly collected high-quality specimens and comparison with published information, we synonymize Mesostephanoides and Szidatia with Gogatea.

Three new species of Ophiotaenia La Rue, 1911 (Cestoda: Proteocephalidae) from dipsadine snakes (Squamata: Colubridae) in Ecuador

The parasite fauna of Neotropical reptiles is poorly known, and the number of parasites described in these hosts does not seem to correspond to the actual species diversity in this zoogeographical region. This also applies to tapeworms such as proteocephalids, which are rarely found in reptiles and are strictly specific to their reptilian hosts. In the present paper, three new species of Ophiotaenia La Rue, 1911 are described from three dipsadine snake species (Squamata: Colubridae) in Ecuador, namely O. jeanmarctouzeti sp. n. from the Neotropical blunt-headed treesnake Imantodes cenchoa (Linnaeus), O. barraganae sp. n. from the beautiful calico snake Oxyrhopus formosus (Wied-Neuwied) and O. velascoae sp. n. from the forest flame snake Oxyrhopus petolarius (Linnaeus). The new species are characterised by type 1 uterine development, the number and distribution of testes, the size of the scolex and other metric features. As no molecular data are available on the specimens collected more than 35 years ago, the phylogenetic relationships of the individual taxa are not known.

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).

Mandibular gland proteomics of the Mexican alligator lizard, Abronia graminea, and the red-lipped arboreal alligator lizard, Abronia lythrochila

A useful approach to deepen our knowledge about the origin and evolution of venom systems in Reptilia has been exploring the vast biodiversity of this clade of vertebrates in search of orally produced proteins with toxic actions, as well as their corresponding delivery systems. The occurrence of toxins in anguimorph lizards has been demonstrated experimentally or inferred from reports of the toxic effects of the oral secretions of taxa within the Varanidae and Helodermatidae families. In the present study, we have focused on two alligator lizards of the Anguidae family, the Mexican alligator lizard, Abronia graminea, and the red-lipped arboreal alligator lizard, A. lythrochila. In addition, the fine morphology of teeth of the latter species is described. The presence of a conserved set of proteins, including B-type natriuretic peptides, cysteine-rich secretory proteins, group III phospholipase A2, and kallikrein, in submandibular gland extracts was demonstrated for both Abronia species. These proteins belong to toxin families found in oral gland secretions of venomous reptile species. This finding, along with previous demonstration of toxin-producing taxa in both paleo- and neoanguimorpha clades, provides further support for the existence of a handful of conserved toxin families in oral secretions across the 100+ million years of Anguimorpha cladogenesis.

Patterns of girdle shape and their correlates in Australian limb-reduced skinks

The evolution of limb reduction in squamates is a classic example of convergence, but the skeletal morphological patterns associated with it are underexplored. To provide insights on the biomechanical and developmental consequences of transitions to limb reduction, we use geometric morphometrics to examine the morphology of pectoral and pelvic girdles in 90 species of limb-reduced skinks and their fully limbed relatives. Clavicle shapes converge towards an acute anterior bend when forelimbs are lost but hindlimbs are retained-a morphology typical of sand-swimmers. This may either indicate functional adaptations to locomotion in fine substrates, or a developmental consequence of complete limb loss. The shape of limb-bearing elements of both girdles (coracoid and pelvis) instead closely mirrors limb reduction, becoming more simplified as undulation replaces limbed locomotion. Integration between girdles decreases in taxa lacking elements of the forelimbs but not hindlimbs, indicating differential selection on each girdle in response to distinct locomotory strategies. However, this pattern becomes less clear when considering phylogenetic history, perhaps because it is limited to one specific clade (Lerista). We show how the functional demands of locomotion can induce changes at different levels of organismal organization, including both external and internal structures.

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

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.

Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae)

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

The skull of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae)

We provide a detailed and first description of the skull, hyoid apparatus, and trachea of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae). Cyclura is a radiation of iguanas restricted to islands of the Caribbean Sea. Species of Cyclura have high rates of endemism, and all species are severely threatened with extinction. Our anatomical description of this threatened iguana is based on high-resolution computed tomography scans of one adult, one putative adult or near adult, and one juvenile specimen, and includes three-dimensional segmented renderings and visualizations. We discuss some observations of intraspecific and ontogenetic variation, and provide a brief comparison with specimens of another species of Cyclura and published descriptions of other iguanas. Our study provides a cranial osteological framework for Cyclura and augments the body of knowledge on iguana anatomy generally. Finally, we posit that our description and future studies may facilitate identification of fossil Cyclura, which could help understand the paleobiogeography of the genus.

Compound osteoderms preserved in amber reveal the oldest known skink

Scincidae is one of the most species-rich and cosmopolitan clades of squamate reptiles. Abundant disarticulated fossil material has also been attributed to this group, however, no complete pre-Cenozoic crown-scincid specimens have been found. A specimen in Burmite (99 MYA) is the first fossil that can be unambiguously referred to this clade. Our analyses place it as nested within extant skinks, supported by the presence of compound osteoderms formed by articulated small ostedermites. The specimen has a combination of dorsal and ventral compound osteoderms and overlapping cycloid scales that is limited to skinks. We propose that this type of osteoderm evolved as a response to an increased overlap of scales, and to reduced stiffness of the dermal armour. Compound osteoderms could be a key innovation that facilitated diversification in this megadiverse family.

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.

The First Identification of Homomorphic XY Sex Chromosomes by Integrating Cytogenetic and Transcriptomic Approaches in Plestiodon elegans (Scincidae)

The sex chromosomes of skinks are usually poorly differentiated and hardly distinguished by cytogenetic methods. Therefore, identifying sex chromosomes in species lacking easily recognizable heteromorphic sex chromosomes is necessary to fully understand sex chromosome diversity. In this paper, we employed cytogenetics, sex quantification of genes, and transcriptomic approaches to characterize the sex chromosomes in Plestiodon elegans. Cytogenetic examination of metaphases revealed a diploid number of 2n = 26, consisting of 12 macrochromosomes and 14 microchromosomes, with no significant heteromorphic chromosome pairs, speculating that the sex chromosomes may be homomorphic or poorly differentiated. The results of the sex quantification of genes showed that Calumenin (calu), COPI coat complex subunit γ 2 (copg2), and Smoothened (smo) were at half the dose in males as in females, suggesting that they are on the X chromosome. Transcriptomic data analysis from the gonads yielded the excess expression male-specific genes (n = 16), in which five PCR molecular markers were developed. Restricting the observed heterozygosity to males suggests the presence of homomorphic sex chromosomes in P. elegans, XX/XY. This is the first breakthrough in the study of the sex chromosomes of Plestiodon.

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.

Environmental (and Random?) Sex Determination in Endangered and Invasive Phelsuma Geckos

INTRODUCTION

Sex is a fundamental characteristic of an individual. It is therefore puzzling why in some systems sex is precisely determined by a genotype, while in others it is influenced by the environment or even subtle, not to say random, factors. Some stochasticity in sex determination would be expected if environmental conditions did not have a large sex-specific effect on fitness. Although data are only available for a small fraction of species, geckos show exceptional variability in sex determination.

METHODS

We tested the effects of three incubation temperatures on sex ratio and adult body size in the invasive gecko Phelsuma laticauda and the vulnerable gecko Phelsuma nigristriata.

RESULTS

We document temperature-dependent sex determination (TSD) in both species. Only females hatched at a low temperature (24°C), whereas male production peaked at an intermediate temperature (28°C) and declined, at least in P. laticauda, again at the highest temperature (31°C). Interestingly, full siblings hatched from eggs glued together during the incubation at temperatures producing both sexes are often of the opposite sex. We found no significant effect of incubation temperature on adult body length.

CONCLUSIONS

Documentation of TSD in the day geckos has implications for conservation practice in environmental management of endangered species or eradication of invasive species. However, it appears that a very subtle (random?) factor may also be involved in their sex determination. In line with this, we found no significant effect of incubation temperature on adult body length.

How seasonality influences the thermal biology of lizards with different thermoregulatory strategies: a meta-analysis

Ectotherms that maintain thermal balance in the face of varying climates should be able to colonise a wide range of habitats. In lizards, thermoregulation usually appears as a variety of behaviours that buffer external influences over physiology. Basking species rely on solar radiation to raise body temperatures and usually show high thermoregulatory precision. By contrast, species that do not bask are often constrained by climatic conditions in their habitats, thus having lower thermoregulatory precision. While much focus has been given to the effects of mean habitat temperatures, relatively less is known about how seasonality affects the thermal biology of lizards on a macroecological scale. Considering the current climate crisis, assessing how lizards cope with temporal variations in environmental temperature is essential to understand better how these organisms will fare under climate change. Activity body temperatures (Tb ) represent the internal temperature of an animal measured in nature during its active period (i.e. realised thermal niche), and preferred body temperatures (Tpref ) are those selected by an animal in a laboratory thermal gradient that lacks thermoregulatory costs (i.e. fundamental thermal niche). Both traits form the bulk of thermal ecology research and are often studied in the context of seasonality. In this study, we used a meta-analysis to test how environmental temperature seasonality influences the seasonal variation in the Tb and Tpref of lizards that differ in thermoregulatory strategy (basking versus non-basking). Based on 333 effect sizes from 137 species, we found that Tb varied over a greater magnitude than Tpref across seasons. Variations in Tb were not influenced by environmental temperature seasonality; however, body size and thermoregulatory strategy mediated Tb responses. Specifically, larger species were subjected to greater seasonal variations in Tb , and basking species endured greater seasonal variations in Tb compared to non-basking species. On the other hand, the seasonal variation in Tpref increased with environmental temperature seasonality regardless of body size. Thermoregulatory strategy also influenced Tpref , suggesting that behaviour has an important role in mediating Tpref responses to seasonal variations in the thermal landscape. After controlling for phylogenetic effects, we showed that Tb and Tpref varied significantly across lizard families. Taken together, our results support the notion that the relationship between thermal biology responses and climatic parameters can be taxon and trait dependent. Our results also showcase the importance of considering ecological and behavioural aspects in macroecological studies. We further highlight current systematic, geographical, and knowledge gaps in thermal ecology research. Our work should benefit those who aim to understand more fully how seasonality shapes thermal biology in lizards, ultimately contributing to the goal of elucidating the evolution of temperature-sensitive traits in ectotherms.

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.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

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.

Description of a New Species of the Pareas hamptoni Complex from Yunnan, China, with Confirmation of P. hamptoni Sensu Stricto in China (Squamata, Pareidae)

We describe a new species of the genus Pareas, based on three specimens collected from Guanyinshan Provincial Nature Reserve in Yuanyang County, Honghe Prefecture, Yunnan Province, China. The new species is distinguished from its congeners by one preocular, one postocular or postocular fused with subocular, loreal not bordering the orbit, one row enlarged vertebral scales, five rows keeled mid-dorsal scales at the middle of the body, 189-192 ventral scales and 72-89 subcaudal scales. The dorsal surfaces of the head and body are yellowish red or yellowish brown, and the belly and ventral surfaces of the head and tail are pinkish yellow or yellow with more or less small black spots. Phylogenetic analyses of mitochondrial DNA recovered the new species being the sister taxon to P. hamptoni sensu stricto. The genetic divergences between the new species and P. hamptoni sensu stricto were 4.2% in the Cyt b sequences and 5.0% in the ND4 sequences. In addition, based on specimens collected from Honghe and Wenshan prefectures, we confirmed that P. hamptoni sensu stricto is distributed in China.

A Palaeogene stem crotaphytid (Aciprion formosum) and the phylogenetic affinities of early fossil pleurodontan iguanians

Pleurodonta is an ancient, diverse clade of iguanian lizard distributed primarily in the Western Hemisphere. Although the clade is a frequent subject of systematic research, phylogenetic resolution among the major pleurodontan clades is elusive. That uncertainty has complicated the interpretations of many fossil pleurodontans. I describe a fossil skull of a pleurodontan lizard from the Palaeogene of Wyoming that was previously allocated to the puzzling taxon Aciprion formosum, and provide an updated morphological matrix for iguanian lizards. Phylogenetic analyses using Bayesian inference demonstrate that the fossil skull is the oldest and first definitive stem member of Crotaphytidae (collared and leopard lizards), establishing the presence of that clade in North America during the Palaeogene. I also discuss new or revised hypotheses for the relationships of several early pleurodontans. In particular, I examine potential evidence for crown-Pleurodonta in the Cretaceous of Mongolia (Polrussia), stem Pleurodonta in the Cretaceous of North America (Magnuviator) and a stem anole in the Eocene of North America (Afairiguana). I suggest that the placement of the fossil crotaphytid is stable to the uncertain phylogeny of Pleurodonta, but recognize the dynamic nature of fossil diagnosis and the potential for updated systematic hypotheses for the other fossils analysed here.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles

BACKGROUND

Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles.

RESULTS

The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species.

CONCLUSIONS

Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

Built for success: Distribution, morphology, ecology and life history of the world's skinks

In animals, the success of particular lineages can be measured in terms of their number of species, the extent of their geographic range, the breadth of their habitats and ecological niches, and the diversity of their morphological and life-history traits. Here, we review the distribution, ecology, morphology and life history of skinks, a diverse lineage of terrestrial vertebrates. We compared key traits between the three subfamilies of skinks, and between skinks and non-scincid lizards. There are currently 1743 described species of skink, which represent 24% of global lizard diversity. Since 2010, 16% of lizard descriptions have been of skinks. The centres of skink diversity are in Australia, New Guinea, southeast Asia, Oceania, Madagascar and central Africa. Compared with non-scincid lizards, skinks have larger distributional ranges, but smaller body sizes. Sexual size dimorphism is rare in skinks. Almost a quarter (23%) of skinks exhibit limb reduction or loss, compared with just 3% of non-scincid lizards. Skinks are more likely to be viviparous (34% of species) compared with non-scincids (13%), and have higher clutch/litter sizes than non-scincids. Although skinks mature later than non-scincids, their longevity is similar to that exhibited by other lizard groups. Most skinks (88%) are active foragers, and they are more likely to be carnivorous than non-scincids. Skinks are more likely to be diurnal or cathemeral than other lizard groups, but they generally have lower field body temperatures compared with non-scincids. The success of skinks appears to be both a result of them hitting upon a winning body plan and ecology, and their capacity to regularly deviate from this body plan and adapt their ecology and life history (e.g. repeated limb reduction and loss, transitions to viviparity) to prevailing conditions.

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.

Frog phylogeny: A time-calibrated, species-level tree based on hundreds of loci and 5,242 species

Large-scale, time-calibrated phylogenies from supermatrix studies have become crucial for evolutionary and ecological studies in many groups of organisms. However, in frogs (anuran amphibians), there is a serious problem with existing supermatrix estimates. Specifically, these trees are based on a limited number of loci (15 or fewer), and the higher-level relationships estimated are discordant with recent phylogenomic estimates based on much larger numbers of loci. Here, we attempted to rectify this problem by generating an expanded supermatrix and combining this with data from phylogenomic studies. To assist in aligning ribosomal sequences for this supermatrix, we developed a new program (TaxonomyAlign) to help perform taxonomy-guided alignments. The new combined matrix contained 5,242 anuran species with data from 307 markers, but with 95% missing data overall. This dataset represented a 71% increase in species sampled relative to the previous largest supermatrix analysis of anurans (adding 2,175 species). Maximum-likelihood analyses generated a tree in which higher-level relationships (and estimated clade ages) were generally concordant with those from phylogenomic analyses but were more discordant with the previous largest supermatrix analysis. We found few obvious problems arising from the extensive missing data in most species. We also generated a set of 100 time-calibrated trees for use in comparative analyses. Overall, we provide an improved estimate of anuran phylogeny based on the largest number of combined taxa and markers to date. More broadly, we demonstrate the potential to combine phylogenomic and supermatrix analyses in other groups of organisms.

Squamate metabolic rates decrease in winter beyond the effect of temperature

The reptilian form of hibernation (brumation) is much less studied than its mammalian and insect equivalents. Hibernation and brumation share some basic features but may differ in others. Evidence for hypometabolism in brumating reptiles beyond the effect of temperature is sporadic and often ignored. We calculated the standard metabolic rates (SMR, oxygen uptake during inactivity), in winter and/or summer, of 156 individuals representing 59 species of Israeli squamates across all 17 local families. For 32 species, we measured the same individuals during both seasons. We measured gas exchange continuously in a dark metabolic chamber, under the average January high and low temperatures (20°C and 12°C), during daytime and nighttime. We examined how SMR changes with season, biome, body size, temperature and time of day, using phylogenetic mixed models. Metabolic rates increased at sunrise in the diurnal species, despite no light or other external cues, while in nocturnal species the metabolic rates did not increase. Cathemeral species shifted from a diurnal-like diel pattern in winter to a nocturnal-like pattern in summer. Regardless of season, Mediterranean species SMRs were 30% higher than similar-sized desert species. Summer SMR of all species together scaled with body size with an exponent of 0.84 but dropped to 0.71 during brumation. Individuals measured during both seasons decreased their SMR between summer and winter by a 47%, on average, at 20°C and by 70% at 12°C. Q10 was 1.75 times higher in winter than in summer, possibly indicating an active suppression of metabolic processes under cold temperatures. Our results challenge the commonly held perception that squamate physiology is mainly shaped by temperature, with little role for intrinsic metabolic regulation. The patterns we describe indicate that seasonal, diel and geographic factors can trigger remarkable shifts in metabolism across squamate species.

De Novo Whole Genome Assemblies for Two Southern African Dwarf Chameleons (Bradypodion, Chamaeleonidae)

A complete and high-quality reference genome has become a fundamental tool for the study of functional, comparative, and evolutionary genomics. However, efforts to produce high-quality genomes for African taxa are lagging given the limited access to sufficient resources and technologies. The southern African dwarf chameleons (Bradypodion) are a relatively young lineage, with a large body of evidence demonstrating the highly adaptive capacity of these lizards. Bradypodion are known for their habitat specialization, with evidence of convergent phenotypes across the phylogeny. However, the underlying genetic architecture of these phenotypes remains unknown for Bradypodion, and without adequate genomic resources, many evolutionary questions cannot be answered. We present de novo assembled whole genomes for Bradypodion pumilum and Bradypodion ventrale, using Pacific Biosciences long-read sequencing data. BUSCO analysis revealed that 96.36% of single copy orthologs were present in the B. pumilum genome and 94% in B. ventrale. Moreover, these genomes boast scaffold N50 of 389.6 and 374.9 Mb, respectively. Based on a whole genome alignment of both Bradypodion genomes, B. pumilum is highly syntenic with B. ventrale. Furthermore, Bradypodion is also syntenic with Anolis lizards, despite the divergence between these lineages estimated to be nearly 170 Ma. Coalescent analysis of the genomic data also suggests that historical changes in effective population size for these species correspond to notable shifts in the southern African environment. These high-quality Bradypodion genome assemblies will support future research on the evolutionary history, diversification, and genetic underpinnings of adaptation in Bradypodion.

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.

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.

The revised reference genome of the leopard gecko (Eublepharis macularius) provides insight into the considerations of genome phasing and assembly

Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000.

Chemical signaling glands are unlinked to species diversification in lizards

Sexual selection has long been thought to increase species diversification. Sexually selected traits, such as sexual signals that contribute to reproductive isolation, were thought to promote diversification. However, studies exploring links between sexually selected traits and species diversification have thus far primarily focused on visual or acoustic signals. Many animals often employ chemical signals (i.e., pheromones) for sexual communications, but large-scale analyses on the role of chemical communications in driving species diversification have been missing. Here, for the first time, we investigate whether traits associated with chemical communications-the presence of follicular epidermal glands-promote diversification across 6,672 lizard species. In most analyses, we found no strong association between the presence of follicular epidermal glands and species diversification rates, either across all lizard species or at lower phylogenetic scales. Previous studies suggest that follicular gland secretions act as species recognition signals that prevent hybridization during speciation in lizards. However, we show that geographic range overlap was no different in sibling species pairs with and without follicular epidermal glands. Together, these results imply that either follicular epidermal glands do not primarily function in sexual communications or sexually selected traits in general (here chemical communication) have a limited effect on species diversification. In our additional analysis accounting for sex-specific differences in glands, we again found no detectable effect of follicular epidermal glands on species diversification rates. Thus, our study challenges the general role of sexually selected traits in broad-scale species diversification patterns.

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.