Resolving the phylogeny of lizards and snakes (Squamata) with extensive sampling of genes and species

Squamate reptiles (lizards and snakes) are one of the most diverse groups of terrestrial vertebrates. Recent molecular analyses have suggested a very different squamate phylogeny relative to morphological hypotheses, but many aspects remain uncertain from molecular data. Here, we analyse higher-level squamate phylogeny with a molecular dataset of unprecedented size, including 161 squamate species for up to 44 nuclear genes each (33 717 base pairs), using both concatenated and species-tree methods for the first time. Our results strongly resolve most squamate relationships and reveal some surprising results. In contrast to most other recent studies, we find that dibamids and gekkotans are together the sister group to all other squamates. Remarkably, we find that the distinctive scolecophidians (blind snakes) are paraphyletic with respect to other snakes, suggesting that snakes were primitively burrowers and subsequently re-invaded surface habitats. Finally, we find that some clades remain poorly supported, despite our extensive data. Our analyses show that weakly supported clades are associated with relatively short branches for which individual genes often show conflicting relationships. These latter results have important implications for all studies that attempt to resolve phylogenies with large-scale phylogenomic datasets.     

Effect of taxon sampling on recovering the phylogeny of squamate reptiles based on complete mitochondrial genome and nuclear gene sequence data

The complete nucleotide sequences of the mitochondrial (mt) genomes of three species of squamate lizards: Blanus cinereus (Amphisbaenidae), Anguis fragilis (Anguidae), and Tarentola mauritanica (Geckkonidae) were determined anew. The deduced amino acid sequences of all 13 mt protein-coding genes were combined into a single data set and phylogenetic relationships among main squamate lineages were analyzed under maximum likelihood (ML) and Bayesian Inference (BI). Within Squamata, the monophyly of Iguanidae, Anguimorpha, Amphisbaenia, Gekkota, Serpentes, and Acrodonta received high statistical support with both methods. It is particularly striking that this is the first molecular analysis that recovers the monophyly of Scincomorpha (including Scincidae, Xantusiidae, Cordylidae, and Lacertidae), although it is only supported in the Bayesian analysis, and it is sensitive to changes in the outgroup (see below). Phylogenetic relationships among the main squamate lineages could not be resolved with ML but received strong support with BI (above 95%). The newly reconstructed phylogeny of squamates does not support the Iguania-Scleroglossa split. Acrodonta and Serpentes form a clade, which is the sister group of the remaining squamate lineages. Within these, Gekkota were the first branching out, followed by Amphisbaenia, and a clade including Anguimorpha as sister group of Scincomorpha + Iguanidae. The recovered topology differed substantially from previously reported hypotheses on squamate relationships, and the relative effect of using different outgroups, genes, and taxon samplings were explored. The sister group relationship of Serpentes + Acrodonta, and their relative basal position within Squamata could be due to a long-branch attraction artifact. Phylogenetic relationships among Scincomorpha, Amphisbaenia, and Anguimorpha were found to be rather unresolved. Future improving of squamate phylogenetic relationships would rely on finding snake and acrodont species with slower mt evolutionary rates, ensuring thorough taxon coverage of squamate diversity, and incorporating more nuclear genes with appropriate evolutionary rates.     

A mitogenomic study on the phylogenetic position of snakes

Phylogenetic relationships of squamates (lizards, amphisbaenians and snakes) have received considerable attention, although no consensus has been reached concerning some basal divergences. This paper focuses on the Serpentes (snakes), whose phylogenetic position within the Squamata remains uncertain despite a number of morphological and molecular studies. Some mitogenomic studies have suggested a sister-group relationship between snakes and varanid lizards, while other studies have identified snakes and lizards as sister groups. However, recent studies using nuclear data have presented a different scenario, with snakes being more closely related to anguimorph and iguanian lizards. In this mitogenomic study we have examined the above hypotheses with the inclusion of amphisbaenians, one gekkotan and one acrodont lizard, taxa not represented in previous mitogenomic studies. To this end we have also extended the representation of snakes by sequencing five additional snake genomes: two scolecophidians ( Ramphotyphlops australis and Typhlops mirus ) two henophidians ( Eunectes notaeus and Boa constrictor ) and one caenophidian ( Elaphe guttata ). The phylogenetic analysis recovered snakes and amphisbaenians as sister groups, thereby differing from previous hypotheses. In addition to a discussion on previous morphological and molecular studies in light of the results presented here, the current study also provides some details regarding features of the new snake mitochondrial genomes described.     

Diversification rates are more strongly related to microhabitat than climate in squamate reptiles (lizards and snakes)

Patterns of species richness among clades can be directly explained by the ages of clades or their rates of diversification. The factors that most strongly influence diversification rates remain highly uncertain, since most studies typically consider only a single predictor variable. Here, we explore the relative impacts of macroclimate (i.e., occurring in tropical vs. temperate regions) and microhabitat use (i.e., terrestrial, fossorial, arboreal, aquatic) on diversification rates of squamate reptile clades (lizards and snakes). We obtained data on microhabitat, macroclimatic distribution, and phylogeny for >4000 species. We estimated diversification rates of squamate clades (mostly families) from a time‐calibrated tree, and used phylogenetic methods to test relationships between diversification rates and microhabitat and macroclimate. Across 72 squamate clades, the best‐fitting model included microhabitat but not climatic distribution. Microhabitat explained ～37% of the variation in diversification rates among clades, with a generally positive impact of arboreal microhabitat use on diversification, and negative impacts of fossorial and aquatic microhabitat use. Overall, our results show that the impacts of microhabitat on diversification rates can be more important than those of climate, despite much greater emphasis on climate in previous studies.     

Hidden support from unpromising data sets strongly unites snakes with anguimorph ‘lizards’

A combined analysis of nuclear, mitochondrial and morphological data robustly resolves snakes as the sister taxon to anguimorph ‘lizards’. Analysed in isolation, nuclear DNA (nDNA) produces a trichotomy between snakes, iguanians and anguimorphs, mitochondrial DNA (mtDNA) is largely uninformative at deeper levels, and morphology tends to nest snakes deep within anguimorphs or with various legless squamate groups. When analysed simultaneously, the nuclear signal is sufficiently strong that mtDNA and morphology are constrained to choose between alternative resolutions of the iguanian–anguimorph–snake trichotomy (generated by the nDNA) – and both support the snake–anguimorph solution. Combined analyses of fast‐evolving or idiosyncratically evolving markers (mtNDA, morphology) with conservative traits (e.g. nuclear genes) might be the best way to resolve ancient, closely spaced divergences. Fast or idiosyncratic markers potentially provide the most information about short, ancient internodes, but can converge on spurious trees if analysed in isolation. However, if constrained to only choosing between plausible trees, such data can contribute unique and valuable phylogenetic signal that resolves such problematic divergences.     

Interactions between locomotion,feeding, and bodily elongation during the evolution of snakes

Information from lizard lineages that have evolved a highly elongate (snake‐like) body form may clarify the selective forces important in the early evolution of snakes. Lizards have evolved bodily elongation via two distinct routes: as an adaptation to burrowing underground or to rapid locomotion above ground. These two routes involve diametrically opposite modifications to the body plan. Burrowing lizards have elongate trunks, small heads, short tails, and relatively constant body widths, whereas surface‐active taxa typically have shorter trunks, wider heads, longer tails, and more variable body widths. Snakes resemble burrowing rather than surface‐active (or aquatic) lizards in these respects, suggesting that snakes evolved from burrowing lizards. The trunk elongation of burrowing lizards increases the volume of the alimentary tract, so that an ability to ingest large meals (albeit consisting of small individual prey items) was present in the earliest snakes. Subsequent shifts to ingestion of wide‐bodied prey came later, after selection dismantled other gape‐constraining morphological attributes, some of which may also have arisen as adaptations to burrowing through hard soil (e.g. relatively small heads, rigid skulls). Adaptations of snake skulls to facilitate ingestion of large prey have evolved to compensate for the reduction of relative head size accompanying bodily elongation; relative to predator body mass, maximum sizes of prey taken by snakes may not be much larger than those of many lizards. This adaptive scenario suggests novel functional links between traits, and a series of testable predictions about the relationships between squamate morphology, habitat, and trophic ecology. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 293–304.     

Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.     

Complete mitochondrial genomes of three lizard species and the systematic position of the Lacertidae (Squamata)*

Complete mitochondrial (mt) genomes were sequenced from representatives of three lacertid lizards: Podarcis siculus, Podarcis muralis and Phoenicolacerta kulzeri. In all three genomes the arrangement of the 22 tRNAs, the two rRNAs and the 13 protein‐coding genes conforms to the common vertebrate arrangement. The phylogenetic position of Lacertidae within the order Squamata was determined through sequence analyses based on large sections of complete mt genomes. The number of nucleotide sites used for tree construction was 9234 when outgroup taxa were included, and 10 499 when only Squamata were compared. The phylogenetic analyses confirmed the sister group relationship between Lacertidae and Amphisbaenia as previously proposed on the basis of molecular data. Additionally, Bayesian analysis revealed a well supported clade comprising (Gekkonidae (Lacertidae + Amphisbaenia)), which is not in accordance with the traditional morphological view and most of the previous molecular studies. It confirms, however, the close relationship between Gekkonidae and Amphisbaenia as revealed in a recent study based on complete mt genomes from a smaller number of taxa. Intra‐ and intergeneric sequence comparisons of six commonly used marker genes showed rather high levels of divergence within the Lacertidae. In the intrageneric comparison the control region proved to be considerably more conserved than the protein coding genes.     

A family of short retroposons (Squam1) from squamate reptiles (Reptilia: Squamata): Structure,evolution, and correlation with phylogeny

Sequences of the SINE family specific to squamate reptiles have been isolated from the genomes of lacertid lizards and sequenced. These retroposons, which we called Squam1, are 360–390 bp long and contain a region similar to the tRNA gene sequence at the 5’ end. This family has also been detected in representatives of other reptile families (varanids, iguanids (Anolis), gekkonids, and snakes), being absent from the genomes of crocodiles as well as amphibians, birds, and mammals. The primary structures of Squam1 copies have been comprehensively analyzed and compared with GenBank sequences. The genomes of most taxa contain two to three SINE subfamilies with specific diagnostic features in their primary structures. Individual similarity between the copies within each taxon is about 85%, with intrageneric similarity being only slightly higher. A comparison of consensus sequences between different lizard families has shown that Squam1 may be a convenient phylogenetic marker for this group of reptiles, having a number of both apomorphic and more or less pronounced synapomorphic features. By this criterion, snakes slightly differ from lizards but obviously belong to the same clade. However, they show no special affinity to varanids as the putative closest relatives of snakes, compared to other lizards.     

Phylogenetic Analysis of Reptilian Hemoglobins: Trees, Rates, and Divergences

Phylogenetic relationships among reptiles were examined using previously published and newly determined hemoglobin sequences. Trees reconstructed from these sequences using maximum-parsimony, neighbor-joining, and maximum-likelihood algorithms were compared with a phylogenetic tree of Amniota, which was assembled on the basis of published morphological data. All analyses differentiated α chains into α^A and α^D types, which are present in all reptiles except crocodiles, where only α^A chains are expressed. The occurrence of the α^D chain in squamates (lizards and snakes only in this study) appears to be a general characteristic of these species. Lizards and snakes also express two types of β chains (βI and βII), while only one type of β chain is present in birds and crocodiles. Reconstructed hemoglobin trees for both α and β sequences did not yield the monophyletic Archosauria (i.e., crocodilians + birds) and Lepidosauria (i.e., Sphenodon+ squamates) groups defined by the morphology tree. This discrepancy, as well as some other poorly resolved nodes, might be due to substantial heterogeneity in evolutionary rates among single hemoglobin lineages. Estimation of branch lengths based on uncorrected amino acid substitutions and on distances corrected for multiple substitutions (PAM distances) revealed that relative rates for squamate α^A and α^D chains and crocodilian β chains are at least twice as high as those of the rest of the chains considered. In contrast to these rate inequalities between reptilian orders, little variation was found within squamates, which allowed determination of absolute evolutionary rates for this subset of hemoglobins. Rate estimates for hemoglobins of lizards and snakes yielded 1.7 (α^A) and 3.3 (β) million years/PAM when calibrated with published divergence time vs. PAM distance correlates for several speciation events within snakes and for the squamate ↔ sphenodontid split. This suggests that hemoglobin chains of squamate reptiles evolved ∼3.5 (α^A) or ∼1.7 times (β) faster than their mammalian equivalents. These data also were used to obtain a first estimate of some intrasquamate divergence times. Received: 15 September 1997 / Accepted: 4 February 1998     

Patterns of interspecific variation in the heart rates of embryonic reptiles

New non-invasive technologies allow direct measurement of heart rates (and thus, developmental rates) of embryos. We applied these methods to a diverse array of oviparous reptiles (24 species of lizards, 18 snakes, 11 turtles, 1 crocodilian), to identify general influences on cardiac rates during embryogenesis. Heart rates increased with ambient temperature in all lineages, but (at the same temperature) were faster in lizards and turtles than in snakes and crocodilians. We analysed these data within a phylogenetic framework. Embryonic heart rates were faster in species with smaller adult sizes, smaller egg sizes, and shorter incubation periods. Phylogenetic changes in heart rates were negatively correlated with concurrent changes in adult body mass and residual incubation period among the lizards, snakes (especially within pythons) and crocodilians. The total number of embryonic heart beats between oviposition and hatching was lower in squamates than in turtles or the crocodilian. Within squamates, embryonic iguanians and gekkonids required more heartbeats to complete development than did embryos of the other squamate families that we tested. These differences plausibly reflect phylogenetic divergence in the proportion of embryogenesis completed before versus after laying.     

Molecular systematics of New World lampropeltinine snakes (Colubridae): implications for biogeography and evolution of food habits

We used mitochondrial gene sequences to infer phylogenetic relationships among North American snakes of the colubrid tribe Lampropeltini (Arizona, Bogertophis, Cemophora , New World Elaphe, Lampropellis, Pituophis, Rhinocheilus, Senticolis, Stilosoma) , and assessed the implications of our findings for the biogeography and evolution of food habits among these serpents. The maximum likelihood phylogeny identified Rhinocheilus as the sister taxon to all other lampropeltinines, and supported the monophyly of Lampropeltis (including Stilosoma) , New World Elaphe , and Pituophis , but not that of Bogertophis. This phylogeny also suggested a sister group relationship between Cemophora and Lampropeltis , and between New World Elaphe and Pituophis , and strongly supported that Sentkolis belongs within Lampropeltini, thus contradicting previous suggestions that Senticolis is not a lampropeltinine. Using a method for approximating ancestral areas of clades, we determined that western North America was most likely the ancestral area of lampropeltinines. Our survey of published studies, combined with unpublished data, indicated that lampropeltinines as a group feed mainly on mammals, less frequently on lizards, birds, and bird eggs, and only rarely on squamate eggs, snakes, anurans, and insects. Some individual species indeed emphasize mammals in their diets, but others most frequently eat lizards, squamate eggs, bird eggs, or snakes, whereas others take two prey types with similar frequency. Our reconstruction of the evolution of food habits among lampropeltinines suggests that a diet emphasizing lizards is ancestral, and therefore diets that mosdy consist of mammals, squamate and bird eggs, and snakes are derived within the clade. In at least some species, smaller individuals prey mostly on lizards and larger ones add mammals to their diets.     

Soft anatomy, diffuse homoplasy, and the relationships of lizards and snakes

Most previous phylogenetic analyses of squamates (‘lizards’ and snakes) employing large character sets have focused on osteology. Soft anatomical traits bearing on this problem have usually been considered in small subsets. Here, a comprehensive phylogenetic analysis of squamate soft anatomy is attempted. 126 informative characters are assessed for 23 squamate lineages, representing snakes, amphisbaenians, dibamids, and all the traditionally recognized ‘families’ of lizards. The traditionally recognized groupings Iguania, Scleroglossa, Gekkota, Scincomorpha, Anguimorpha and Varanoidea are corroborated in this analysis. More controversial taxa are resolved as follows. Xantusiids, amphisbaenians and dibamids cluster with gekkotans, and snakes are strongly allied with anguimorphs in general, and varanids in particular. Nearly all these clades are congruent with those found in a recent comprehensive osteological analysis; the strong support for snake‐varanid relationships found in both studies is particularly notable. This congruence is surprising given that previous studies of soft anatomy tended to give differing and often heterodox results. These previous results can be attributed to overrepresentation of misleading characters in small isolated data sets. Such misleading signals are minimized when data sets are combined. For instance, the snake‐varanid clade is contradicted by many characters, and analyses of particular organ systems therefore give differing results. However, characters that are incongruent with the snake‐varanid clade also disagree with each other (diffuse homoplasy), rather than forming coherent support for some particular alternative clade (concerted homoplasy). In a combined analysis these incongruent but diffuse characters cancel each other out to leave a very strong (and orthodox) phylogenetic signal. These results underscore the view that the raw amount of homoplasy — as revealed by consistency and retention indices — is not the only determinant of phylogenetic signal; the distribution of that homoplasy is also important. Thus, questioning a phylogenetic hypothesis (e.g. the snake‐varanid clade) by identifying numerous conflicting characters is insufficient — the structure of the conflicting characters should be assessed in a rigorous phylogenetic analysis.     

The root of the mammalian tree inferred from whole mitochondrial genomes

Morphological and molecular data are currently contradictory over the position of monotremes with respect to marsupial and placental mammals. As part of a re-evaluation of both forms of data we examine complete mitochondrial genomes in more detail. There is a particularly large discrepancy in the frequencies of thymine and cytosine (T-C) between mitochondrial genomes that appears to affect some deep divergences in the mammalian tree. We report that recoding nucleotides to RY-characters, and partitioning maximum-likelihood analyses among subsets of data reduces such biases, and improves the fit of models to the data, respectively. RY-coding also increases the signal on the internal branches relative to external, and thus increases the phylogenetic signal. In contrast to previous analyses of mitochondrial data, our analyses favor Theria (marsupials plus placentals) over Marsupionta (monotremes plus marsupials). However, a short therian stem lineage is inferred, which is at variance with the traditionally deep placement of monotremes on morphological data.     

Phylogenomics reveal a robust fungal tree of life

Our understanding of the tree of life (TOL) is still fragmentary. Until recently, molecular phylogeneticists have built trees based on ribosomal RNA sequences and selected protein sequences, which, however, usually suffered from lack of support for the deeper branches and inconsistencies probably due to limited subsampling of the entire genome. Now, phylogenetic hypotheses can be based on the analysis of full genomes. We used available complete genome data as well as the eukaryote orthologous group (KOG) proteins to reconstruct with confidence basal branches of the fungal TOL. Phylogenetic analysis of a core of 531 KOGs shared among 21 fungal genomes, three animal genomes and one plant genome showed a single tree with high support resulting from four different methods of phylogenetic reconstruction. The single tree that we inferred from our dataset showed excellent nodal support for each branch, suggesting that it reflects the true phylogenetic relationships of the species involved.     

A tiny lizard (Lepidosauria,Squamata) from the Lower Cretaceous of Spain

Abstract: The smallest living amniotes are all lizards, but the fossil history of this size trait in Squamata is difficult to follow because small skeletons have low preservation potential and are often hard to detect in the field. A new squamate taxon, Jucaraseps grandipes gen. et sp. nov., is here described on the basis of an articulated skeleton from the Early Cretaceous Spanish lagerstätten of Las Hoyas. It differs from other known Mesozoic lizards in combining very small body size with a short rostrum, low maxillary tooth count, a relatively slender and elongated body, and short limbs with large hind feet. Phylogenetic analysis using TNT places it on the stem of a clade encompassing scincomorphs, gekkotans, snakes, amphisbaenians and anguimorphs. Comparison with modern lizards suggests it was probably a cryptic surface or subsurface ground dweller but not a burrower.     

Coalescent-Based Genome Analyses Resolve the Early Branches of the Euarchontoglires

Despite numerous large-scale phylogenomic studies, certain parts of the mammalian tree are extraordinarily difficult to resolve. We used the coding regions from 19 completely sequenced genomes to study the relationships within the super-clade Euarchontoglires (Primates, Rodentia, Lagomorpha, Dermoptera and Scandentia) because the placement of Scandentia within this clade is controversial. The difficulty in resolving this issue is due to the short time spans between the early divergences of Euarchontoglires, which may cause incongruent gene trees. The conflict in the data can be depicted by network analyses and the contentious relationships are best reconstructed by coalescent-based analyses. This method is expected to be superior to analyses of concatenated data in reconstructing a species tree from numerous gene trees. The total concatenated dataset used to study the relationships in this group comprises 5,875 protein-coding genes (9,799,170 nucleotides) from all orders except Dermoptera (flying lemurs). Reconstruction of the species tree from 1,006 gene trees using coalescent models placed Scandentia as sister group to the primates, which is in agreement with maximum likelihood analyses of concatenated nucleotide sequence data. Additionally, both analytical approaches favoured the Tarsier to be sister taxon to Anthropoidea, thus belonging to the Haplorrhine clade. When divergence times are short such as in radiations over periods of a few million years, even genome scale analyses struggle to resolve phylogenetic relationships. On these short branches processes such as incomplete lineage sorting and possibly hybridization occur and make it preferable to base phylogenomic analyses on coalescent methods.     

Mitogenomics of Vetigastropoda: insights into the evolution of pallial symmetry

The nucleotide sequences of the complete or nearly complete mitochondrial (mt) genomes of seven vetigastropods were determined: Angaria neglecta (Angarioidea), Phasianella solida (Phasianelloidea), Granata lyrata (Seguenzioidea), Tegula lividomaculata and Bolma rugosa (Trochoidea), Diodora graeca (Fissurelloidea) and Lepetodrilus schrolli (Lepetodriloidea). While the mt genomes of the superfamilies Angarioidea, Phasianelloidea, Seguenzioidea and Trochoidea conform generally to the ancestral gene order of Vetigastropoda and Gastropoda, those of the superfamilies Fissurelloidea and Lepetodriloidea have suffered important rearrangements. The gene order of the mtDNA of Chrysomallon squamiferum, a representative of Neomphalina, was also analysed since it has been proposed to be closely related to Vetigastropoda, and showed a distinct arrangement. The reconstructed phylogenies recovered Neomphalina as a distinct gastropod lineage that is the sister group (only with moderate bootstrap support) of a clade including Vetigastropoda and Neritimorpha + Caeno‐gastropoda while the relative position of Heterobranchia and Patellogastropoda in the gastropod tree could not be determined definitively due to their long branches. Within the monophyletic Vetigastropoda, the superfamily Fissurelloidea was recovered as the sister group of two lineages, one including Lepetodriloidea as the sister group of Seguenzioidea + Halitoidea, the other including Phasianelloidea, Angarioidea and Trochoidea without resolved relationships. The long branches of Fissurelloidea were found to introduce significant tree instability in phylogenetic reconstruction. The new phylogeny supports that the loss of the right pallial gill occurred multiple times in vetigastropod evolution as previously suggested and that Phasianelloidea, Angarioidea and Trochoidea radiated from a common asymmetric (single‐gilled) ancestor that lived in the middle Palaeozoic.     

Strong conservation of the bird Z chromosome in reptilian genomes is revealed by comparative painting despite 275 million years divergence

The divergence of lineages leading to extant squamate reptiles (lizards, snakes, and amphisbaenians) and birds occurred about 275 million years ago. Birds, unlike squamates, have karyotypes that are typified by the presence of a number of very small chromosomes. Hence, a number of chromosome rearrangements might be expected between bird and squamate genomes. We used chromosome-specific DNA from flow-sorted chicken (Gallus gallus) Z sex chromosomes as a probe in cross-species hybridization to metaphase spreads of 28 species from 17 families representing most main squamate lineages and single species of crocodiles and turtles. In all but one case, the Z chromosome was conserved intact despite very ancient divergence of sauropsid lineages. Furthermore, the probe painted an autosomal region in seven species from our sample with characterized sex chromosomes, and this provides evidence against an ancestral avian-like system of sex determination in Squamata. The avian Z chromosome synteny is, therefore, conserved albeit it is not a sex chromosome in these squamate species.