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.

Evolutionary paths to new phenotypes

Ecological model systems inform on innovative traits in plants and animals.

Incomplete dosage balance and 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.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Coordinating tiny limbs and long bodies: Geometric mechanics of lizard terrestrial swimming

Although typically possessing four limbs and short bodies, lizards have evolved diverse morphologies, including elongate trunks with tiny limbs. Such forms are hypothesized to aid locomotion in cluttered/fossorial environments but propulsion mechanisms (e.g., the use of body and/or limbs to interact with substrates) and potential body/limb coordination remain unstudied. Here, we use biological experiments, a geometric theory of locomotion, and robophysical models to investigate body-limb coordination in diverse lizards. Locomotor field studies in short-limbed, elongate lizards (Brachymeles and Lerista) and laboratory studies of fully limbed lizards (Uma scoparia and Sceloporus olivaceus) and a snake (Chionactis occipitalis) reveal that body-wave dynamics can be described by a combination of standing and traveling waves; the ratio of the amplitudes of these components is inversely related to the degree of limb reduction and body elongation. The geometric theory (which replaces laborious calculation with diagrams) helps explain our observations, predicting that the advantage of traveling-wave body undulations (compared with a standing wave) emerges when the dominant thrust-generation mechanism arises from the body rather than the limbs and reveals that such soil-dwelling lizards propel via "terrestrial swimming" like sand-swimming lizards and snakes. We test our hypothesis by inducing the use of traveling waves in stereotyped lizards via modulating the ground-penetration resistance. Study of a limbed/undulatory robophysical model demonstrates that a traveling wave is beneficial when propulsion is generated by body-environment interaction. Our models could be valuable in understanding functional constraints on the evolutionary processes of elongation and limb reduction as well as advancing robot designs.

Reproductive Biology, Sperm storage, and Sexual Maturity of Thamnodynastes strigatus (Serpentes: Dipsadidae)

Life history strategies determine and influence many aspects of species fitness. In this study, we describe the reproductive biology - reproductive cycle, sperm storage, and sexual maturity - of Thamnodynastes strigatus in South Brazil. We analyzed 49 individuals (25 males and 24 females) from herpetological collections. The reproductive cycle of males and females was described considering the morpho-anatomical and histological changes in the testes, ductus deferens, and kidney, as well in the ovary and oviduct. The age at the onset of sexual maturity was determined by skeletochronology of the caudal vertebra. The reproductive cycle is seasonal semi-synchronous and most individuals have a reproductive peak in spring and summer. The seasonal biennial reproductive cycle and viviparity are two phylogenetically conserved characters in Tachymenini snakes. Thamnodynastes Strigatus females store sperm in the utero-vaginal junction furrows during autumn. There were no differences between the ages of sexual maturity of males (4-11y) and females (4-12y). Females reach sexual maturity at larger body sizes, and this may confer an adaptive advantage due to a higher fecundity potential. Herein, we confirmed the previously described seasonal biennial reproductive cycle of T. strigatus through histological analysis.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Viviparous Reptile Regarded to Have Temperature-Dependent Sex Determination Has Old XY Chromosomes

The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12–15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79–116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

Evolutionary origins of viviparity consistent with palaeoclimate and lineage diversification

It is of fundamental importance for the field of evolutionary biology to understand when and why major evolutionary transitions occur. Live-bearing young (viviparity) is a major evolutionary change and has evolved from egg-laying (oviparity) independently in many vertebrate lineages and most abundantly in lizards and snakes. Although contemporary viviparous squamate species generally occupy cold climatic regions across the globe, it is not known whether viviparity evolved as a response to cold climate in the first place. Here, we used available published time-calibrated squamate phylogenies and parity data on 3,498 taxa. We compared the accumulation of transitions from oviparity to viviparity relative to background diversification and a simulated binary trait. Extracting the date of each transition in the phylogenies and informed by 65 my of global palaeoclimatic data, we tested the nonexclusive hypotheses that viviparity evolved under the following: (a) cold, (b) long-term stable climatic conditions and (c) with background diversification rate. We show that stable and long-lasting cold climatic conditions are correlated with transitions to viviparity across squamates. This correlation of parity mode and palaeoclimate is mirrored by background diversification in squamates, and simulations of a binary trait also showed a similar association with palaeoclimate, meaning that trait evolution cannot be separated from squamate lineage diversification. We suggest that parity mode transitions depend on environmental and intrinsic effects and that background diversification rate may be a factor in trait diversification more generally.

Fang evolution in venomous snakes: Adaptation of 3D tooth shape to the biomechanical properties of their prey

Venomous snakes are among the world's most specialized predators. During feeding, they use fangs to penetrate the body tissues of their prey, but the success of this penetration depends on the shape of these highly specialized teeth. Here, we examined the evolution of fang shape in a wide range of snakes using 3D geometric morphometrics (3DGM) and cross-sectional tooth sharpness measurements. We investigated the relationship of these variables with six diet categories based on the prey's biomechanical properties, and tested for evolutionary convergence using two methods. Our results show that slender elongate fangs with sharp tips are used by snakes that target soft-skinned prey (e.g., mammals), whereas fangs become more robust and blunter as the target's skin becomes scaly (e.g., fish and reptiles) and eventually hard-shelled (e.g., crustaceans), both with and without correction for evolutionary allometry. Convergence in fang shape is present, indicating that fangs of snakes with the same diet are more similar than those of closely related species with different diets. Establishing the relationship between fang morphology and diet helps to explain how snakes became adapted to different lifestyles, while also providing a proxy to infer diet in lesser known species or extinct snakes from the fossil record.

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

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

Evolutionary Dynamics and Consequences of Parthenogenesis in Vertebrates

Parthenogenesis is asexual reproduction without any required participation from males and, as such, is a null model for sexual reproduction. In a comparative context, we can expand our understanding of the evolution and ecology of sex by investigating the consequences of parthenogenesis. In this review, we examine the theoretical predictions of and empirical results on the evolution of asexual reproduction in vertebrates, focusing on recent studies addressing the origins and geographic spread of parthenogenetic lineages and the genomic consequences of an asexual life history. With advances in computational methods and genome technologies, researchers are poised to make rapid and significant progress in studying the origin and evolution of parthenogenesis in vertebrates, thus providing an important perspective on understanding biodiversity patterns of both asexual and sexual populations.

Chromosome-Level Assembly of the Common Lizard (Zootoca vivipara) Genome

Squamate reptiles exhibit high variation in their phenotypic traits and geographical distributions and are therefore fascinating taxa for evolutionary and ecological research. However, genomic resources are very limited for this group of species, consequently inhibiting research efforts. To address this gap, we assembled a high-quality genome of the common lizard, Zootoca vivipara (Lacertidae), using a combination of high coverage Illumina (shotgun and mate-pair) and PacBio sequencing data, coupled with RNAseq data and genetic linkage map generation. The 1.46-Gb genome assembly has a scaffold N50 of 11.52 Mb with N50 contig size of 220.4 kb and only 2.96% gaps. A BUSCO analysis indicates that 97.7% of the single-copy Tetrapoda orthologs were recovered in the assembly. In total, 19,829 gene models were annotated to the genome using a combination of ab initio and homology-based methods. To improve the chromosome-level assembly, we generated a high-density linkage map from wild-caught families and developed a novel analytical pipeline to accommodate multiple paternity and unknown father genotypes. We successfully anchored and oriented almost 90% of the genome on 19 linkage groups. This annotated and oriented chromosome-level reference genome represents a valuable resource to facilitate evolutionary studies in squamate reptiles.

Using 3D-digital photogrammetry to examine scaling of the body axis in burrowing skinks

Three-dimensional (3D) modeling techniques have been increasingly utilized across disciplines for the visualization and analysis of complex structures. We employ 3D-digital photogrammetry for understanding the scaling of the body axis of 12 species of scincid lizards in the genus Brachymeles. These skinks represent a diverse radiation which shows tremendous variation in body size and degree of axial elongation. Because of the complex nature of the body axis, 3D-methods are important for understanding how the body axis evolves. 3D-digital photogrammetry presents a flexible, inexpensive, and portable system for the reconstruction of biological forms. As body size increased among species, the cross-sectional area and circumference of the head and other portions of the body axis increased isometrically, which indicates that species of differing sizes possess proportionally similar head and body shapes. These results suggest that there are no substantial head and body shape changes with body size among the sampled species, but further comparative studies with larger sample sizes and functional studies of size and morphology effects on burrowing or above-ground locomotion are needed.

What makes a fang? Phylogenetic and ecological controls on tooth evolution in rear-fanged snakes

BACKGROUND

Fangs are a putative key innovation that revolutionized prey capture and feeding in snakes, and - along with their associated venom phenotypes - have made snakes perhaps the most medically-significant vertebrate animals. Three snake clades are known for their forward-positioned fangs, and these clades (Elapidae, Viperidae, and Atractaspidinae) contain the majority of snakes that are traditionally considered venomous. However, many other snakes are "rear-fanged": they possess potentially venom-delivering teeth situated at the rear end of the upper jaw. Quantification of fang phenotypes - and especially those of rear-fanged species - has proved challenging or impossible owing to the small size and relative rarity of many such snakes. Consequently, it has been difficult to understand the evolutionary history of both venom and prey-capture strategies across extant snakes. We quantified variation in the dentition of 145 colubriform ("advanced") snake species using microCT scanning and compared dental characters with ecological data on species' diet and prey capture method(s) to understand broader patterns in snake fang evolution.

RESULTS

Dental traits such as maxilla length, tooth number, and fang size show strong phylogenetic signal across Colubriformes. We find extreme heterogeneity and evolutionary lability in the rear-fanged phenotype in colubrid (colubrine, dipsadine, and natricine lineages) and lamprophiid snakes, in contrast to relative uniformity in the front fanged phenotypes of other groups (vipers and, to a lesser extent, elapids). Fang size and position are correlated with venom-use in vipers, elapids, and colubrid snakes, with the latter group shifting fangs anteriorly by shortening the entire maxillary bone. We find that maxilla length and tooth number may also be correlated with the evolution of dietary specialization. Finally, an ancestral state reconstruction suggests that fang loss is a recurring phenomenon in colubrid snakes, likely accompanied by shifts in diet and prey capture mode.

CONCLUSIONS

Our study provides a framework for quantifying the complex morphologies associated with venom use in snakes. Our results suggest that fang phenotypes, and particularly the rear-fanged phenotype, in snakes are both diverse and labile, facilitating a wide range of ecological strategies and contributing to spectacular radiations of these organisms in tropical and subtropical biomes worldwide.

Reversals in complex traits uncovered as reticulation events: Lessons from the evolution of parity-mode, chromosome morphology, and maternal resource transfer

Complex traits include, among many others, the evolution of eyes, wings, body forms, reproductive modes, human intelligence, social behavior, diseases, and chromosome morphology. Dollo's law states that the evolution of complex traits is irreversible. However, potential exceptions have been proposed. Here, we investigated whether reticulation, a simple and elegant means by which complex characters may be reacquired, could account for suggested reversals in the evolution of complex characters using two datasets with sufficient genetic coverage and a total of five potential reversals. Our analyses uncovered a potential reversal in the evolution of parity mode and a potential reversal in the evolution of placentotrophy of fish (Cyprinodontiformes) as reticulation events. Moreover, in a reptile that exhibits a potential reversal from viviparity to oviparity (Zootoca vivipara), reticulation provided the most parsimonious explanation for sex chromosome evolution. Therefore, three of the five studied potential reversals were unraveled as reticulation events. This constitutes the first evidence that accounting for reticulation can fundamentally influence the interpretation of the evolution of complex traits, that testing for reticulation is crucial for obtaining robust phylogenies, and that complex ancestral characters may be reacquired through hybridization with a lineage that still exhibits the trait. Hybridization, rather than reappearance of ancestral traits by means of small evolutionary steps, may thus account for suggested exceptions to Dollo's law. Consequently, ruling out reticulation is required to claim the evolutionary reversal of complex characters and potential exceptions to Dollo's rule.

Distinct telomere differences within a reproductively bimodal common lizard population

Different strategies of reproductive mode, either oviparity (egg-laying) or viviparity (live-bearing), will be associated with a range of other life-history differences that are expected to affect patterns of ageing and longevity. It is usually difficult to compare the effects of alternative reproductive modes because of evolutionary and ecological divergence. However, the very rare exemplars of reproductive bimodality, in which different modes exist within a single species, offer an opportunity for robust and controlled comparisons.One trait of interest that could be associated with life history, ageing and longevity is the length of the telomeres, which form protective caps at the chromosome ends and are generally considered a good indicator of cellular health. The shortening of these telomeres has been linked to stressful conditions; therefore, it is possible that differing reproductive costs will influence patterns of telomere loss. This is important because a number of studies have linked a shorter telomere length to reduced survival.Here, we have studied maternal and offspring telomere dynamics in the common lizard (Zootoca vivipara). Our study has focused on a population where oviparous and viviparous individuals co-occur in the same habitat and occasionally interbreed to form admixed individuals.While viviparity confers many advantages for offspring, it might also incur substantial costs for the mother, for example require more energy. Therefore, we predicted that viviparous mothers would have relatively shorter telomeres than oviparous mothers, with admixed mothers having intermediate telomere lengths. There is thought to be a heritable component to telomere length; therefore, we also hypothesized that offspring would follow the same pattern as the mothers.Contrary to our predictions, the viviparous mothers and offspring had the longest telomeres, and the oviparous mothers and offspring had the shortest telomeres. The differing telomere lengths may have evolved as an effect of the life-history divergence between the reproductive modes, for example due to the increased growth rate that viviparous individuals may undergo to reach a similar size at reproduction. A free http://onlinelibrary.wiley.com/doi/10.1111/1365-2435.13408/suppinfo can be found within the Supporting Information of this article.

Differential reproductive investment in co-occurring oviparous and viviparous common lizards (Zootoca vivipara) and implications for life-history trade-offs with viviparity

Live-bearing reproduction (viviparity) has evolved from egg-laying (oviparity) independently many times and most abundantly in squamate reptiles. Studying life-history trade-offs between the two reproductive modes is an inherently difficult task, as most transitions to viviparity are evolutionarily old and/or are confounded by environmental effects. The common lizard (Zootoca vivipara) is one of very few known reproductively bimodal species, in which some populations are oviparous and others viviparous. Oviparous and viviparous populations can occur in sympatry in the same environment, making this a unique system for investigating alternative life-history trade-offs between oviparous and viviparous reproduction. We find that viviparous females exhibit larger body size, smaller clutch sizes, a larger reproductive investment, and a higher hatching success rate than oviparous females. We find that offspring size and weight from viviparous females was lower compared to offspring from oviparous females, which may reflect space constraints during pregnancy. We suggest that viviparity in common lizards is associated with increased reproductive burden for viviparous females and speculate that this promoted the evolution of larger body size to create more physical space for developing embryos. In the context of life-history trade-offs in the evolution of viviparity, we suggest that the extent of correlation between reproductive traits, or differences between reproductive modes, may also depend on the time since the transition occurred.

Alternative reproductive adaptations predict asymmetric responses to climate change in lizards

Anthropogenic climate change ranks among the major global-scale threats to modern biodiversity. Extinction risks are known to increase via the interactions between rapid climatic alterations and environmentally-sensitive species traits that fail to adapt to those changes. Accumulating evidence reveals the influence of ecophysiological, ecological and phenological factors as drivers underlying demographic collapses that lead to population extinctions. However, the extent to which life-history traits influence population responses to climate change remains largely unexplored. The emerging ‘cul-de-sac hypothesis’ predicts that reptilian viviparity (‘live-bearing’ reproduction), a ‘key innovation’ facilitating historical invasions of cold climates, increases extinction risks under progressively warming climates compared to oviparous reproduction – as warming advances polewards/mountainwards, historically cold-climates shrink, leading viviparous species to face demographic collapses. We present the first large-scale test of this prediction based on multiple lizard radiations and on future projections of climate-based ecological niche models. Viviparous species were found to experience stronger elevational range shifts (and potentially increased extinctions) in coming decades, compared to oviparous lizards. Therefore, our analyses support the hypothesis’s fundamental prediction that elevational shifts are more severe in viviparous species, and highlight the role that life-history adaptations play in the responses of biodiversity to ongoing climate change.

Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Phylogenomic analyses of more than 4000 nuclear loci resolve the origin of snakes among lizard families

Squamate reptiles (lizards and snakes) are the most diverse group of terrestrial vertebrates, with more than 10 000 species. Despite considerable effort to resolve relationships among major squamates clades, some branches have remained difficult. Among the most vexing has been the placement of snakes among lizard families, with most studies yielding only weak support for the position of snakes. Furthermore, the placement of iguanian lizards has remained controversial. Here we used targeted sequence capture to obtain data from 4178 nuclear loci from ultraconserved elements from 32 squamate taxa (and five outgroups) including representatives of all major squamate groups. Using both concatenated and species-tree methods, we recover strong support for a sister relationship between iguanian and anguimorph lizards, with snakes strongly supported as the sister group of these two clades. These analyses strongly resolve the difficult placement of snakes within squamates and show overwhelming support for the contentious position of iguanians. More generally, we provide a strongly supported hypothesis of higher-level relationships in the most species-rich tetrapod clade using coalescent-based species-tree methods and approximately 100 times more loci than previous estimates.

Appendages and gene regulatory networks: Lessons from the limbless

Among squamate reptiles, dozens of lineages have independently evolved complete or partial limb reduction. This remarkable convergence of limbless and limb-reduced phenotypes provides multiple natural replicates of different ages to explore the evolution and development of the vertebrate limb and the gene regulatory network that controls its formation. The most successful and best known of the limb-reduced squamates are snakes, which evolved a limb-reduced body form more than 100 million years ago. Recent studies have revealed the unexpected finding that many ancient limb enhancers are conserved in the genomes of snakes. Analyses in limbed animals show that many of these limb enhancers are also active during development of the phallus, suggesting that these enhancers may have been retained in snakes due their importance in regulating transcription in the external genitalia. This hypothesis is substantiated by functional tests of snake enhancers, which demonstrate that snake enhancer elements have lost limb function while retaining genital enhancer function. The large degree of overlap in the gene regulatory networks deployed during limb and phallus development may act to constrain the divergence of shared gene network components and the evolution of appendage morphology. Future studies will reveal whether limb regulatory elements have undergone similar functional changes in other lineages of limb-reduced squamates.

Helminth parasites of Phyllodactylidae and Gekkonidae lizards in a Caatinga ecological station, northeastern Brazil

Abstract We investigated the parasites of five lizard species belonging to Phyllodactylidae (Phyllopezus pollicaris and Gymnodactylus geckoides) and Gekkonidae (Hemidactylus agrius, Lygodactylus klugei and Hemidactylus brasilianus) families in a semiarid region of Brazil. Six nematode species were identified: Parapharyngodon alvarengai and Spauligodon oxkutzcabiensis (Pharyngodonidae), Physaloptera lutzi (Physalopteridae), Skrjabinelazia intermedia (Seuratidae), Trichospirura sp. (Rhabdochonidae) and Piratuba sp. (Onchocercidae), and a cestode species, Oochoristica sp. (Linstowiidae). The most prevalent species were Spauligodon oxkutzcabiensis, which infected P. pollicaris (75%), and Parapharyngodon alvarengai, which infected G. geckoides (29%). South American lizards were identified as being new hosts for the Trichospirura genus (a usual parasite of mammals), and there were 16 new occurrences of parasite species in the five lizard species studied herein.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

Proteomic Profile of Mabuya sp. (Squamata: Scincidae) Ovary and Placenta During Gestation

Reptiles are one of the most diverse groups of vertebrates, providing an integrated system for comparative studies on metabolic, animal physiology, and developmental biology. However, the molecular data available are limited and only recently have started to call attention in the "omics" sciences. Mabuya sp. is a viviparous placentrotrophic skink with particular reproductive features, including microlecithal eggs, early luteolysis, prolonged gestation, and development of a highly specialized placenta. This placenta is responsible for respiratory exchange and the transference of all nutrients necessary for embryonic development. Our aim was to identify differentially expressed proteins in the ovary and placenta of Mabuya sp. during early, mid, and late gestation; their possible metabolic pathways; and biological processes. We carried out a comparative proteomic analysis during gestation in both tissues by sodium dodecyl sulfate polyacrylamide gel electrophoresis, two-dimensional gel electrophoresis, and matrix-assisted laser desorption/ionization. Differential protein expression in both tissues (Student's t-test P < 0.05) was related to several processes such as cell structure, cell movement, and energy. Proteins found in ovary are mainly associated with follicular development and its regulation. In the placenta, particularly during mid and late gestation, protein expression is involved in nutrient metabolism, transport, protein synthesis, and embryonic development. This work provides new insights about the proteins expressed and their physiological mechanisms in Mabuya sp. placenta and ovary during gestation.

The Primary Structure of β(I)-Chain of Hemoglobin from Snake Sindhi Krait (Bungarus sindanus sindanus)

The amino acid sequence of β(I)-globin chain from Sindhi Krait (Bungarus sindanus sindanus) was determined to study the molecular evolution among snakes. The hemoglobin was isolated from the red blood cells and was analyzed by ion-exchange chromatography (IEX). The crude globin was subjected to reversed phased-high performance liquid chromatography (RP-HPLC) using C4 column. The N-terminal sequences of intact globin chains and tryptic peptides were determined by Edman degradation in a pulsed liquid gas phase sequencer using an online Phenylthiohydantoin analyzer. Sindhi Krait is expected to express three hemoglobin components that are composed of β(II), β(I), α(D) and α(A)-globin chains, as apparent by IEX, RP-HPLC and N-terminal sequence analyses. Sequence alignment and phylogenetic analyses of β(I) globin chain from Sindhi Krait showed closest relationship with β(I) globin chain from Rattlesnake, Water snake and Indigo snake. Interestingly, comparison of primary sequence of β(I) globin chain of Sindhi Krait with human β chain revealed 63 % similarity along with the retention of all heme contact points. Variations among the two sequences were prominent at αβ contact points and in regions directly not important for function.

How Should Genes and Taxa be Sampled for Phylogenomic Analyses with Missing Data? An Empirical Study in Iguanian Lizards

Targeted sequence capture is becoming a widespread tool for generating large phylogenomic data sets to address difficult phylogenetic problems. However, this methodology often generates data sets in which increasing the number of taxa and loci increases amounts of missing data. Thus, a fundamental (but still unresolved) question is whether sampling should be designed to maximize sampling of taxa or genes, or to minimize the inclusion of missing data cells. Here, we explore this question for an ancient, rapid radiation of lizards, the pleurodont iguanians. Pleurodonts include many well-known clades (e.g., anoles, basilisks, iguanas, and spiny lizards) but relationships among families have proven difficult to resolve strongly and consistently using traditional sequencing approaches. We generated up to 4921 ultraconserved elements with sampling strategies including 16, 29, and 44 taxa, from 1179 to approximately 2.4 million characters per matrix and approximately 30% to 60% total missing data. We then compared mean branch support for interfamilial relationships under these 15 different sampling strategies for both concatenated (maximum likelihood) and species tree (NJst) approaches (after showing that mean branch support appears to be related to accuracy). We found that both approaches had the highest support when including loci with up to 50% missing taxa (matrices with ~40-55% missing data overall). Thus, our results show that simply excluding all missing data may be highly problematic as the primary guiding principle for the inclusion or exclusion of taxa and genes. The optimal strategy was somewhat different for each approach, a pattern that has not been shown previously. For concatenated analyses, branch support was maximized when including many taxa (44) but fewer characters (1.1 million). For species-tree analyses, branch support was maximized with minimal taxon sampling (16) but many loci (4789 of 4921). We also show that the choice of these sampling strategies can be critically important for phylogenomic analyses, since some strategies lead to demonstrably incorrect inferences (using the same method) that have strong statistical support. Our preferred estimate provides strong support for most interfamilial relationships in this important but phylogenetically challenging group.

Group living in squamate reptiles: a review of evidence for stable aggregations

How sociality evolves and is maintained remains a key question in evolutionary biology. Most studies to date have focused on insects, birds, and mammals but data from a wider range of taxonomic groups are essential to identify general patterns and processes. The extent of social behaviour among squamate reptiles is under-appreciated, yet they are a promising group for further studies. Living in aggregations is posited as an important step in the evolution of more complex sociality. We review data on aggregations among squamates and find evidence for some form of aggregations in 94 species across 22 families. Of these, 18 species across 7 families exhibited 'stable' aggregations that entail overlapping home ranges and stable membership in long-term (years) or seasonal aggregations. Phylogenetic analysis suggests that stable aggregations have evolved multiple times in squamates. We: (i) identify significant gaps in our understanding; (ii) outline key traits which should be the focus of future research; and (iii) outline the potential for utilising reproductive skew theory to provide insights into squamate sociality.

Distinct patterns of desynchronized limb regression in malagasy scincine lizards (squamata, scincidae)

Scincine lizards in Madagascar form an endemic clade of about 60 species exhibiting a variety of ecomorphological adaptations. Several subclades have adapted to burrowing and convergently regressed their limbs and eyes, resulting in a variety of partial and completely limbless morphologies among extant taxa. However, patterns of limb regression in these taxa have not been studied in detail. Here we fill this gap in knowledge by providing a phylogenetic analysis of DNA sequences of three mitochondrial and four nuclear gene fragments in an extended sampling of Malagasy skinks, and microtomographic analyses of osteology of various burrowing taxa adapted to sand substrate. Based on our data we propose to (i) consider Sirenoscincus Sakata & Hikida, 2003, as junior synonym of Voeltzkowia Boettger, 1893; (ii) resurrect the genus name Grandidierina Mocquard, 1894, for four species previously included in Voeltzkowia; and (iii) consider Androngo Brygoo, 1982, as junior synonym of Pygomeles Grandidier, 1867. By supporting the clade consisting of the limbless Voeltzkowia mira and the forelimb-only taxa V. mobydick and V. yamagishii, our data indicate that full regression of limbs and eyes occurred in parallel twice in the genus Voeltzkowia (as hitherto defined) that we consider as a sand-swimming ecomorph: in the Voeltzkowia clade sensu stricto the regression first affected the hindlimbs and subsequently the forelimbs, whereas the Grandidierina clade first regressed the forelimbs and subsequently the hindlimbs following the pattern prevalent in squamates. Timetree reconstructions for the Malagasy Scincidae contain a substantial amount of uncertainty due to the absence of suitable primary fossil calibrations. However, our preliminary reconstructions suggest rapid limb regression in Malagasy scincids with an estimated maximal duration of 6 MYr for a complete regression in Paracontias, and 4 and 8 MYr respectively for complete regression of forelimbs in Grandidierina and hindlimbs in Voeltzkowia.

Integrated analyses resolve conflicts over squamate reptile phylogeny and reveal unexpected placements for fossil taxa

Squamate reptiles (lizards and snakes) are a pivotal group whose relationships have become increasingly controversial. Squamates include >9000 species, making them the second largest group of terrestrial vertebrates. They are important medicinally and as model systems for ecological and evolutionary research. However, studies of squamate biology are hindered by uncertainty over their relationships, and some consider squamate phylogeny unresolved, given recent conflicts between molecular and morphological results. To resolve these conflicts, we expand existing morphological and molecular datasets for squamates (691 morphological characters and 46 genes, for 161 living and 49 fossil taxa, including a new set of 81 morphological characters and adding two genes from published studies) and perform integrated analyses. Our results resolve higher-level relationships as indicated by molecular analyses, and reveal hidden morphological support for the molecular hypothesis (but not vice-versa). Furthermore, we find that integrating molecular, morphological, and paleontological data leads to surprising placements for two major fossil clades (Mosasauria and Polyglyphanodontia). These results further demonstrate the importance of combining fossil and molecular information, and the potential problems of estimating the placement of fossil taxa from morphological data alone. Thus, our results caution against estimating fossil relationships without considering relevant molecular data, and against placing fossils into molecular trees (e.g. for dating analyses) without considering the possible impact of molecular data on their placement.

Zebrin II / aldolase C expression in the cerebellum of the western diamondback rattlesnake (Crotalus atrox)

Aldolase C, also known as Zebrin II (ZII), is a glycolytic enzyme that is expressed in cerebellar Purkinje cells of the vertebrate cerebellum. In both mammals and birds, ZII is expressed heterogeneously, such that there are sagittal stripes of Purkinje cells with high ZII expression (ZII+), alternating with stripes of Purkinje cells with little or no expression (ZII-). The patterns of ZII+ and ZII- stripes in the cerebellum of birds and mammals are strikingly similar, suggesting that it may have first evolved in the stem reptiles. In this study, we examined the expression of ZII in the cerebellum of the western diamondback rattlesnake (Crotalus atrox). In contrast to birds and mammals, the cerebellum of the rattlesnake is much smaller and simpler, consisting of a small, unfoliated dome of cells. A pattern of alternating ZII+ and ZII- sagittal stripes cells was not observed: rather all Purkinje cells were ZII+. This suggests that ZII stripes have either been lost in snakes or that they evolved convergently in birds and mammals.

Locomotor benefits of being a slender and slick sand swimmer

Squamates classified as 'subarenaceous' possess the ability to move long distances within dry sand; body elongation among sand and soil burrowers has been hypothesized to enhance subsurface performance. Using X-ray imaging, we performed the first kinematic investigation of the subsurface locomotion of the long, slender shovel-nosed snake (Chionactis occipitalis) and compared its biomechanics with those of the shorter, limbed sandfish lizard (Scincus scincus). The sandfish was previously shown to maximize swimming speed and minimize the mechanical cost of transport during burial. Our measurements revealed that the snake also swims through sand by propagating traveling waves down the body, head to tail. Unlike the sandfish, the snake nearly followed its own tracks, thus swimming in an approximate tube of self-fluidized granular media. We measured deviations from tube movement by introducing a parameter, the local slip angle, βs, which measures the angle between the direction of movement of each segment and body orientation. The average βs was smaller for the snake than for the sandfish; granular resistive force theory (RFT) revealed that the curvature utilized by each animal optimized its performance. The snake benefits from its slender body shape (and increased vertebral number), which allows propagation of a higher number of optimal curvature body undulations. The snake's low skin friction also increases performance. The agreement between experiment and RFT combined with the relatively simple properties of the granular 'frictional fluid' make subarenaceous swimming an attractive system to study functional morphology and bauplan evolution.