CONSERVED SEX CHROMOSOMES ACROSS ADAPTIVELY RADIATEDANOLISLIZARDS

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.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

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.

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

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

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei – a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

A highly-complete chromosome-scale genome assembly of the brown anole, Anolis sagrei, provides insight into the evolution of sex chromosomes and is a crucial resource for this model lizard species.

Propagule size and sex ratio influence colonisation dynamics after introduction of a non-native lizard

The composition of founding populations plays an important role in colonisation dynamics and can influence population growth during early stages of biological invasion. Specifically, founding populations with small propagules (i.e., low number of founders) are vulnerable to the Allee effect and have reduced likelihood of establishment compared to those with large propagules. The founding sex ratio can also impact establishment via its influence on mating success and offspring production. Our goal was to test the effects of propagule size and sex ratio on offspring production and annual population growth following introductions of a non-native lizard species (Anolis sagrei). We manipulated propagule composition on nine small islands, then examined offspring production, population growth, and survival rate of founders and their descendants encompassing three generations. By the third reproductive season, per capita offspring production was higher on islands seeded with a relatively large propagule size, but population growth was not associated with propagule size. Propagule sex ratio did not affect offspring production, but populations with a female-biased propagule had positive growth, whereas those with a male-biased propagule had negative growth in the first year. Populations were not affected by propagule sex ratio in subsequent years, possibly due to rapid shifts towards balanced (or slightly female-biased) population sex ratios. Overall, we show that different components of population fitness have different responses to propagule size and sex ratio in ways that could affect early stages of biological invasion. Despite these effects, the short lifespan and high fecundity of A. sagrei likely helped small populations to overcome Allee effects and enabled all populations to successfully establish. Our rare experimental manipulation of propagule size and sex ratio can inform predictions of colonisation dynamics in response to different compositions of founding populations, which is critical in the context of population ecology and invasion dynamics.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

Background

High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.

Findings

We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.

Conclusions

These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Do male and female heterogamety really differ in expression regulation? Lack of global dosage balance in pygopodid geckos

Differentiation of sex chromosomes is thought to have evolved with cessation of recombination and subsequent loss of genes from the degenerated partner (Y and W) of sex chromosomes, which in turn leads to imbalance of gene dosage between sexes. Based on work with traditional model species, theory suggests that unequal gene copy numbers lead to the evolution of mechanisms to counter this imbalance. Dosage compensation, or at least achieving dosage balance in expression of sex-linked genes between sexes, has largely been documented in lineages with male heterogamety (XX/XY sex determination), while ZZ/ZW systems are assumed to be usually associated with the lack of chromosome-wide gene dose regulatory mechanisms. Here, we document that although the pygopodid geckos evolved male heterogamety with a degenerated Y chromosome 32-72 Ma, one species in particular, Burton's legless lizard (Lialis burtonis), does not possess dosage balance in the expression of genes in its X-specific region. We summarize studies on gene dose regulatory mechanisms in animals and conclude that there is in them no significant dichotomy between male and female heterogamety. We speculate that gene dose regulatory mechanisms are likely to be related to the general mechanisms of sex determination instead of type of heterogamety. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Are Geckos Special in Sex Determination? Independently Evolved Differentiated ZZ/ZW Sex Chromosomes in Carphodactylid Geckos

Amniotes possess astonishing variability in sex determination ranging from environmental sex determination (ESD) to genotypic sex determination (GSD) with highly differentiated sex chromosomes. Geckos are one of the few amniote groups with substantial variability in sex determination. What makes them special in this respect? We hypothesized that the extraordinary variability of sex determination in geckos can be explained by two alternatives: 1) unusual lability of sex determination, predicting that the current GSD systems were recently formed and are prone to turnovers; and 2) independent transitions from the ancestral ESD to later stable GSD, which assumes that geckos possessed ancestrally ESD, but once sex chromosomes emerged, they remain stable in the long term. Here, based on genomic data, we document that the differentiated ZZ/ZW sex chromosomes evolved within carphodactylid geckos independently from other gekkotan lineages and remained stable in the genera Nephrurus, Underwoodisaurus, and Saltuarius for at least 15 Myr and potentially up to 45 Myr. These results together with evidence for the stability of sex chromosomes in other gekkotan lineages support more our second hypothesis suggesting that geckos do not dramatically differ from the evolutionary transitions in sex determination observed in the majority of the amniote lineages.

Large-scale comparative analysis of cytogenetic markers across Lepidoptera

Fluorescence in situ hybridization (FISH) allows identification of particular chromosomes and their rearrangements. Using FISH with signal enhancement via antibody amplification and enzymatically catalysed reporter deposition, we evaluated applicability of universal cytogenetic markers, namely 18S and 5S rDNA genes, U1 and U2 snRNA genes, and histone H3 genes, in the study of the karyotype evolution in moths and butterflies. Major rDNA underwent rather erratic evolution, which does not always reflect chromosomal changes. In contrast, the hybridization pattern of histone H3 genes was well conserved, reflecting the stable organisation of lepidopteran genomes. Unlike 5S rDNA and U1 and U2 snRNA genes which we failed to detect, except for 5S rDNA in a few representatives of early diverging lepidopteran lineages. To explain the negative FISH results, we used quantitative PCR and Southern hybridization to estimate the copy number and organization of the studied genes in selected species. The results suggested that their detection was hampered by long spacers between the genes and/or their scattered distribution. Our results question homology of 5S rDNA and U1 and U2 snRNA loci in comparative studies. We recommend the use of histone H3 in studies of karyotype evolution.

With or Without W? Molecular and Cytogenetic Markers are Not Sufficient for Identification of Environmentally-Induced Sex Reversal in the Bearded Dragon

Transitions from environmental sex determination (ESD) to genotypic sex determination (GSD) require an intermediate step of sex reversal, i.e., the production of individuals with a mismatch between the ancestral genotypic and the phenotypic sex. Among amniotes, the sole well-documented transition in this direction was shown in the laboratory in the central bearded dragon, Pogona vitticeps, where very high incubation temperatures led to the production of females with the male-typical (ZZ) genotype. These sex-reversed females then produced offspring whose sex depended on the incubation temperature. Sex-reversed animals identified by molecular and cytogenetic markers were also reported in the field, and their increasing incidence was speculated as a climate warming-driven transition in sex determination. We show that the molecular and cytogenetic markers normally sex-linked in P. vitticeps are also sex-linked in P. henrylawsoni and P. minor, which points to quite ancient sex chromosomes in this lineage. Nevertheless, we demonstrate, based on a crossing experiment with a male bearded dragon who possesses a mismatch between phenotypic sex and genotype, that the used cytogenetic and molecular markers might not be reliable for the identification of sex reversal. Sex reversal should not be considered as the only mechanism causing a mismatch between genetic sex-linked markers and phenotypic sex, which can emerge also by other processes, here most likely by a rare recombination between regions of sex chromosomes which are normally sex-linked. We warn that sex-linked, even apparently for a long evolutionary time, and sex-specific molecular and cytogenetic markers are not a reliable tool for the identification of sex-reversed individuals in a population and that sex reversal has to be verified by other approaches, particularly by observation of the sex ratio of the progeny.

Differences in Homomorphic Sex Chromosomes Are Associated with Population Divergence in Sex Determination in Carinascincus ocellatus (Scincidae: Lygosominae)

Sex determination directs development as male or female in sexually reproducing organisms. Evolutionary transitions in sex determination have occurred frequently, suggesting simple mechanisms behind the transitions, yet their detail remains elusive. Here we explore the links between mechanisms of transitions in sex determination and sex chromosome evolution at both recent and deeper temporal scales (<1 Myr; ~79 Myr). We studied a rare example of a species with intraspecific variation in sex determination, Carinascincus ocellatus, and a relative, Liopholis whitii, using c-banding and mapping of repeat motifs and a custom Y chromosome probe set to identify the sex chromosomes. We identified both unique and conserved regions of the Y chromosome among C. ocellatus populations differing in sex determination. There was no evidence for homology of sex chromosomes between C. ocellatus and L. whitii, suggesting independent evolutionary origins. We discuss sex chromosome homology between members of the subfamily Lygosominae and propose links between sex chromosome evolution, sex determination transitions, and karyotype evolution.

Independent Evolution of Sex Chromosomes in Eublepharid Geckos, A Lineage with Environmental and Genotypic Sex Determination

Geckos demonstrate a remarkable variability in sex determination systems, but our limited knowledge prohibits accurate conclusions on the evolution of sex determination in this group. Eyelid geckos (Eublepharidae) are of particular interest, as they encompass species with both environmental and genotypic sex determination. We identified for the first time the X-specific gene content in the Yucatán banded gecko, Coleonyx elegans, possessing X1X1X2X2/X1X2Y multiple sex chromosomes by comparative genome coverage analysis between sexes. The X-specific gene content of Coleonyx elegans was revealed to be partially homologous to genomic regions linked to the chicken autosomes 1, 6 and 11. A qPCR-based test was applied to validate a subset of X-specific genes by comparing the difference in gene copy numbers between sexes, and to explore the homology of sex chromosomes across eleven eublepharid, two phyllodactylid and one sphaerodactylid species. Homologous sex chromosomes are shared between Coleonyx elegans and Coleonyx mitratus, two species diverged approximately 34 million years ago, but not with other tested species. As far as we know, the X-specific gene content of Coleonyx elegans / Coleonyx mitratus was never involved in the sex chromosomes of other gecko lineages, indicating that the sex chromosomes in this clade of eublepharid geckos evolved independently.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

Corytophanids Replaced the Pleurodont XY System with a New Pair of XY Chromosomes

Almost all lizard families in the pleurodont clade share the same XY system. This system was meticulously studied in Anolis carolinensis, where it shows a highly degenerated Y chromosome and a male-specific X chromosome dosage compensation mechanism. Corytophanids (casque-headed lizards) have been proposed as the only family in the pleurodont clade to lack the XY system. In this study, we worked with extensive genomic and transcriptomic data from Basiliscus vittatus, a member of the Corytophanidae family that inhabits the tropical rainforests of Mexico. We confirmed that B. vittatus underwent a sex chromosome system turnover, which consisted in the loss of the pleurodont XY system and the gain of a new pair of XY chromosomes that are orthologous to chicken chromosome 17. We estimated the origin of the sex chromosome system to have occurred ∼63 Ma in the ancestor of corytophanids. Moreover, we identified 12 XY gametologues with particular attributes, such as functions related to the membrane and intracellular trafficking, very low expression levels, blood specificity, and incomplete dosage compensation in males.

First Report of Sex Chromosomes in Night Lizards (Scincoidea: Xantusiidae)

Squamate reptiles (lizards, snakes, and amphibians) are an outstanding group for studying sex chromosome evolution—they are old, speciose, geographically widespread, and exhibit myriad sex-determining modes. Yet, the vast majority of squamate species lack heteromorphic sex chromosomes. Cataloging the sex chromosome systems of species lacking easily identifiable, heteromorphic sex chromosomes, therefore, is essential before we are to fully understand the evolution of vertebrate sex chromosomes. Here, we use restriction site-associated DNA sequencing (RADseq) to classify the sex chromosome system of the granite night lizard, Xantusia henshawi. RADseq is an effective alternative to traditional cytogenetic methods for determining a species’ sex chromosome system (i.e., XX/XY or ZZ/ZW), particularly in taxa with non-differentiated sex chromosomes. Although many xantusiid lineages have been karyotyped, none possess heteromorphic sex chromosomes. We identified a ZZ/ZW sex chromosome system in X. henshawi—the first such data for this family. Furthermore, we report that the X. henshawi sex chromosome contains fragments of genes found on Gallus gallus chromosomes 7, 12, and 18 (which are homologous to Anolis carolinensis chromosome 2), the first vertebrate sex chromosomes to utilize this linkage group.

ZZ/ZW Sex Determination with Multiple Neo-Sex Chromosomes is Common in Madagascan Chameleons of the Genus Furcifer (Reptilia: Chamaeleonidae)

Chameleons are well-known, highly distinctive lizards characterized by unique morphological and physiological traits, but their karyotypes and sex determination system have remained poorly studied. We studied karyotypes in six species of Madagascan chameleons of the genus Furcifer by classical (conventional stain, C-banding) and molecular (comparative genomic hybridization, in situ hybridization with rDNA, microsatellite, and telomeric sequences) cytogenetic approaches. In contrast to most sauropsid lineages, the chameleons of the genus Furcifer show chromosomal variability even among closely related species, with diploid chromosome numbers varying from 2n = 22 to 2n = 28. We identified female heterogamety with cytogenetically distinct Z and W sex chromosomes in all studied species. Notably, multiple neo-sex chromosomes in the form Z₁Z₁Z₂Z₂/Z₁Z₂W were uncovered in four species of the genus (F. bifidus, F. verrucosus, F. willsii, and previously studied F. pardalis). Phylogenetic distribution and morphology of sex chromosomes suggest that multiple sex chromosomes, which are generally very rare among vertebrates with female heterogamety, possibly evolved several times within the genus Furcifer. Although acrodontan lizards (chameleons and dragon lizards) demonstrate otherwise notable variability in sex determination, it seems that female heterogamety with differentiated sex chromosomes remained stable in the chameleons of the genus Furcifer for about 30 million years.

Discovery of Putative XX/XY Male Heterogamety in Emydura subglobosa Turtles Exposes a Novel Trajectory of Sex Chromosome Evolution in Emydura

The discovery of sex chromosome systems in non-model organisms has elicited growing recognition that sex chromosomes evolved via diverse paths that are not fully elucidated. Lineages with labile sex determination, such as turtles, hold critical cues, yet data are skewed toward hide-neck turtles (suborder Cryptodira) and scant for side-neck turtles (suborder Pleurodira). Here, we used classic and molecular cytogenetics to investigate Emydura subglobosa (ESU), an unstudied side-neck turtle with genotypic sex determination from the family Chelidae, where extensive morphological divergence exists among XX/XY systems. Our data represent the first cytogenetic description for ESU. Similarities were found between ESU and E. macquarii (EMA), such as identical chromosome number (2n = 50), a single and dimorphic nucleolus organizer region (NOR) localized in a microchromosome pair (ESU14) of both sexes (detected via FISH of 18S rDNA). Only the larger NOR is active (detected by silver staining). As in EMA, comparative genome hybridization revealed putative macro XX/XY chromosomes in ESU (the 4th largest pair). Our comparative analyses and revaluation of previous data strongly support the hypothesis that Emydura's XX/XY system evolved via fusion of an ancestral micro-Y (retained by Chelodina longicollis) onto a macro-autosome. This evolutionary trajectory differs from the purported independent evolution of XX/XY from separate ancestral autosomes in Chelodina and Emydura that was previously reported. Our data permit dating this Y-autosome fusion to at least the split of Emydura around 45 Mya and add critical information about the evolution of the remarkable diversity of sex-determining mechanisms in turtles, reptiles, and vertebrates.

Conserved sex chromosomes and karyotype evolution in monitor lizards (Varanidae)

Despite their long history with the basal split dating back to the Eocene, all species of monitor lizards (family Varanidae) studied so far share the same chromosome number of 2n = 40. However, there are differences in the morphology of the macrochromosome pairs 5-8. Further, sex determination, which revealed ZZ/ZW sex microchromosomes, was studied only in a few varanid species and only with techniques that did not test their homology. The aim of this study was to (i) test if cryptic interchromosomal rearrangements of larger chromosomal blocks occurred during the karyotype evolution of this group, (ii) contribute to the reconstruction of the varanid ancestral karyotype, and (iii) test homology of sex chromosomes among varanids. We investigated these issues by hybridizing flow sorted chromosome paints from Varanus komodoensis to metaphases of nine species of monitor lizards. The results show that differences in the morphology of the chromosome pairs 5-8 can be attributed to intrachromosomal rearrangements, which led to transitions between acrocentric and metacentric chromosomes in both directions. We also documented the first case of spontaneous triploidy among varanids in Varanus albigularis. The triploid individual was fully grown, which demonstrates that polyploidization is compatible with life in this lineage. We found that the W chromosome differs between species in size and heterochromatin content. The varanid Z chromosome is clearly conserved in all the analyzed species. Varanids, in addition to iguanas, caenophidian snakes, and lacertid lizards, are another squamate group with highly conserved sex chromosomes over a long evolutionary time.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Distribution and amplification of interstitial telomeric sequences (ITSs) in Australian dragon lizards support frequent chromosome fusions in Iguania

Telomeric sequences are generally located at the ends of chromosomes; however, they can also be found in non-terminal chromosomal regions when they are known as interstitial telomeric sequences (ITSs). Distribution of ITSs across closely related and divergent species elucidates karyotype evolution and speciation as ITSs provide evolutionary evidence for chromosome fusion. In this study, we performed physical mapping of telomeric repeats by fluorescence in situ hybridisation (FISH) in seven Australian dragon lizards thought to represent derived karyotypes of squamate reptiles and a gecko lizard with considerably different karyotypic feature. Telomeric repeats were present at both ends of all chromosomes in all species, while varying numbers of ITSs were also found on microchromosomes and in pericentromeric or centromeric regions on macrochromosomes in five lizard species examined. This suggests that chromosomal rearrangements from ancestral squamate reptiles to Iguania occurred mainly by fusion between ancestral types of acrocentric chromosomes and/or between microchromosomes, leading to appearance of bi-armed macrochromosomes, and in the reduction of microchromosome numbers. These results support the previously proposed hypothesis of karyotype evolution in squamate reptiles. In addition, we observed the presence of telomeric sequences in the similar regions to heterochromatin of the W microchromosome in Pogona barbata and Doporiphora nobbi, while sex chromosomes for the two species contained part of the nucleolar organiser regions (NORs). This likely implies that these ITSs are a part of the satellite DNA and not relics of chromosome fusions. Amplification of telomeric repeats may have involved heterochromatinisation of sex-specific W chromosomes and play a role in the organisation of the nucleolus.

Evolution of Dosage Compensation in Anolis carolinensis, a Reptile with XX/XY Chromosomal Sex Determination

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

Evolutionary Dynamics of the W Chromosome in Caenophidian Snakes

The caenophidian (assigned also as “advanced”) snakes are traditionally viewed as a group of reptiles with a limited karyotypic variation and stable ZZ/ZW sex chromosomes. The W chromosomes of the caenophidian snakes are heterochromatic, and pioneering studies demonstrated that they are rich in repetitive elements. However, a comparative study of the evolutionary dynamics of the repetitive content of the W chromosome across the whole lineage is missing. Using molecular-cytogenetic techniques, we explored the distribution of four repetitive motifs (microsatellites GATA, GACA, AG and telomeric-like sequences), which are frequently accumulated in differentiated sex chromosomes in vertebrates, in the genomes of 13 species of the caenophidian snakes covering a wide phylogenetic spectrum of the lineage. The results demonstrate a striking variability in the morphology and the repetitive content of the W chromosomes even between closely-related species, which is in contrast to the homology and long-term stability of the gene content of the caenophidian Z chromosome. We uncovered that the tested microsatellite motifs are accumulated on the degenerated, heterochromatic W chromosomes in all tested species of the caenophidian snakes with the exception of the Javan file snake representing a basal clade. On the other hand, the presence of the accumulation of the telomeric-like sequences on the caenophidian W chromosome is evolutionary much less stable. Moreover, we demonstrated that large accumulations of telomeric-like motifs on the W chromosome contribute to sexual differences in the number of copies of the telomeric and telomeric-like repeats estimated by quantitative PCR, which might be confusing and incorrectly interpreted as sexual differences in telomere length.

Evolution of Karyotypes in Chameleons

The reconstruction of the evolutionary dynamics of karyotypes and sex determining systems in squamate reptiles is precluded by the lack of data in many groups including most chameleons (Squamata: Acrodonta: Chamaeleonidae). We performed cytogenetic analysis in 16 species of chameleons from 8 genera covering the phylogenetic diversity of the family and also phylogenetic reconstruction of karyotype evolution in this group. In comparison to other squamates, chameleons demonstrate rather variable karyotypes, differing in chromosome number, morphology and presence of interstitial telomeric signal (ITS). On the other hand, the location of rDNA is quite conserved among chameleon species. Phylogenetic analysis combining our new results and previously published data tentatively suggests that the ancestral chromosome number for chameleons is 2n = 36, which is the same as assumed for other lineages of the clade Iguania, i.e., agamids and iguanas. In general, we observed a tendency for the reduction of chromosome number during the evolution of chameleons, however, in Rieppeleon brevicaudatus, we uncovered a chromosome number of 2n = 62, very unusual among squamates, originating from a number of chromosome splits. Despite the presence of the highly differentiated ZZ/ZW sex chromosomes in the genus Furcifer, we did not detect any unequivocal sexual differences in the karyotypes of any other studied species of chameleons tested using differential staining and comparative genomic hybridization, suggesting that sex chromosomes in most chameleons are only poorly differentiated.

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno's original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster. Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways.

Localization and distribution of gonadal proteins in the oviparous lizard Sceloporus aeneus (Squamata: Phrynosomatidae)

Among vertebrates, several specific proteins are involved in the function and development of gonads. Several genes such as SOX9, FOXL2, DDX4, IFITM3, and DPPA3, are active during embryonic differentiation and maintain their expression in adult tissues, playing important roles in the function and development of the line cell, where these are produced. Among reptiles, molecular mechanisms for sex differentiation have been analyzed in turtles, crocodiles, and some lizards, while in adult stages such studies are scarce. The aim of this study was to locate and analyze the distribution of important gonadal proteins in adult and embryonic ovaries and testes of the oviparous lizard Sceloporus aeneus (Squamata: Phrynosomatidae). Adult specimens and embryos of the lizard S. aeneus were collected in Milpa Alta, a suburb located Southwest of Mexico City. Expression of gonadal proteins was analyzed using immunofluorescent staining and confocal microscopy. Our results showed that SOX9 is located in Sertoli cells of embryonic and adult testes. FOXL2 is expressed in follicular cells of adult ovaries. DDX4 and IFITM3 are located in germ line cells as well as in follicular cells of adult ovaries. DPPA3 was observed in somatic and germ line cells of adult and embryonic gonads. Our observations show that important molecules of vertebrate ovaries and testes are conserved in S. aeneus and it is suggested that these may have a similar role during gonadal development and function.

Stable Cretaceous sex chromosomes enable molecular sexing in softshell turtles (Testudines: Trionychidae)

Turtles demonstrate variability in sex determination ranging from environmental sex determination (ESD) to highly differentiated sex chromosomes. However, the evolutionary dynamics of sex determining systems in this group is not well known. Differentiated ZZ/ZW sex chromosomes were identified in two species of the softshell turtles (Trionychidae) from the subfamily Trionychinae and Z-specific genes were identified in a single species. We tested Z-specificity of a subset of these genes by quantitative PCR comparing copy gene numbers in male and female genomes in 10 species covering the phylogenetic diversity of trionychids. We demonstrated that differentiated ZZ/ZW sex chromosomes are conserved across the whole family and that they were already present in the common ancestor of the extant trionychids. As the sister lineage, Carettochelys insculpta, possess ESD, we can date the origin of the sex chromosomes in trionychids between 200 Mya (split of Trionychidae and Carettochelyidae) and 120 Mya (basal splitting of the recent trionychids). The results support the evolutionary stability of differentiated sex chromosomes in some lineages of ectothermic vertebrates. Moreover, our approach determining sex-linkage of protein coding genes can be used as a reliable technique of molecular sexing across trionychids useful for effective breeding strategy in conservation projects of endangered species.

Evolutionary stability of sex chromosomes in snakes

Amniote vertebrates possess various mechanisms of sex determination, but their variability is not equally distributed. The large evolutionary stability of sex chromosomes in viviparous mammals and birds was believed to be connected with their endothermy. However, some ectotherm lineages seem to be comparably conserved in sex determination, but previously there was a lack of molecular evidence to confirm this. Here, we document a stability of sex chromosomes in advanced snakes based on the testing of Z-specificity of genes using quantitative PCR (qPCR) across 37 snake species (our qPCR technique is suitable for molecular sexing in potentially all advanced snakes). We discovered that at least part of sex chromosomes is homologous across all families of caenophidian snakes (Acrochordidae, Xenodermatidae, Pareatidae, Viperidae, Homalopsidae, Colubridae, Elapidae and Lamprophiidae). The emergence of differentiated sex chromosomes can be dated back to about 60 Ma and preceded the extensive diversification of advanced snakes, the group with more than 3000 species. The Z-specific genes of caenophidian snakes are (pseudo)autosomal in the members of the snake families Pythonidae, Xenopeltidae, Boidae, Erycidae and Sanziniidae, as well as in outgroups with differentiated sex chromosomes such as monitor lizards, iguanas and chameleons. Along with iguanas, advanced snakes are therefore another example of ectothermic amniotes with a long-term stability of sex chromosomes comparable with endotherms.

Mixed-Up Sex Chromosomes: Identification of Sex Chromosomes in the X1X1X2X2/X1X2Y System of the Legless Lizards of the Genus Lialis (Squamata: Gekkota: Pygopodidae)

Geckos in general show extensive variability in sex determining systems, but only male heterogamety has been demonstrated in the members of their legless family Pygopodidae. In the pioneering study published more than 45 years ago, multiple sex chromosomes of the type X1X1X2X2/X1X2Y were described in Burton's legless lizard (Lialisburtonis) based on conventional cytogenetic techniques. We conducted cytogenetic analyses including comparative genomic hybridization and fluorescence in situ hybridization (FISH) with selected cytogenetic markers in this species and the previously cytogenetically unstudied Papua snake lizard (Lialis jicari) to better understand the nature of these sex chromosomes and their differentiation. Both species possess male heterogamety with an X1X1X2X2/X1X2Y sex chromosome system; however, the Y and one of the X chromosomes are not small chromosomes as previously reported in L. burtonis, but the largest macrochromosomal pair in the karyotype. The Y chromosomes in both species have large heterochromatic blocks with extensive accumulations of GATA and AC microsatellite motifs. FISH with telomeric probe revealed an exclusively terminal position of telomeric sequences in L. jicari (2n = 42 chromosomes in females), but extensive interstitial signals, potentially remnants of chromosomal fusions, in L.burtonis (2n = 34 in females). Our study shows that even largely differentiated and heteromorphic sex chromosomes might be misidentified by conventional cytogenetic analyses and that the application of more sensitive cytogenetic techniques for the identification of sex chromosomes is beneficial even in the classical examples of multiple sex chromosomes.

Cytogenetic Insights into the Evolution of Chromosomes and Sex Determination Reveal Striking Homology of Turtle Sex Chromosomes to Amphibian Autosomes

Turtle karyotypes are highly conserved compared to other vertebrates; yet, variation in diploid number (2n = 26-68) reflects profound genomic reorganization, which correlates with evolutionary turnovers in sex determination. We evaluate the published literature and newly collected comparative cytogenetic data (G- and C-banding, 18S-NOR, and telomere-FISH mapping) from 13 species spanning 2n = 28-68 to revisit turtle genome evolution and sex determination. Interstitial telomeric sites were detected in multiple lineages that underwent diploid number and sex determination turnovers, suggesting chromosomal rearrangements. C-banding revealed potential interspecific variation in centromere composition and interstitial heterochromatin at secondary constrictions. 18S-NORs were detected in secondary constrictions in a single chromosomal pair per species, refuting previous reports of multiple NORs in turtles. 18S-NORs are linked to ZW chromosomes in Apalone and Pelodiscus and to X (not Y) in Staurotypus. Notably, comparative genomics across amniotes revealed that the sex chromosomes of several turtles, as well as mammals and some lizards, are homologous to components of Xenopus tropicalis XTR1 (carrying Dmrt1). Other turtle sex chromosomes are homologous to XTR4 (carrying Wt1). Interestingly, all known turtle sex chromosomes, except in Trionychidae, evolved via inversions around Dmrt1 or Wt1. Thus, XTR1 appears to represent an amniote proto-sex chromosome (perhaps linked ancestrally to XTR4) that gave rise to turtle and other amniote sex chromosomes.

How to make a sex chromosome

Sex chromosomes can evolve once recombination is halted between a homologous pair of chromosomes. Owing to detailed studies using key model systems, we have a nuanced understanding and a rich review literature of what happens to sex chromosomes once recombination is arrested. However, three broad questions remain unanswered. First, why do sex chromosomes stop recombining in the first place? Second, how is recombination halted? Finally, why does the spread of recombination suppression, and therefore the rate of sex chromosome divergence, vary so substantially across clades? In this review, we consider each of these three questions in turn to address fundamental questions in the field, summarize our current understanding, and highlight important areas for future work.

Sex chromosome evolution begins when recombination between a homologous pair of chromosomes is halted. Here, Wright et al. review our current understanding of the causes and mechanisms of recombination suppression between incipient sex chromosomes and suggest future directions for the field.

Sex determination, longevity, and the birth and death of reptilian species

Vertebrate sex-determining mechanisms (SDMs) are triggered by the genotype (GSD), by temperature (TSD), or occasionally, by both. The causes and consequences of SDM diversity remain enigmatic. Theory predicts SDM effects on species diversification, and life-span effects on SDM evolutionary turnover. Yet, evidence is conflicting in clades with labile SDMs, such as reptiles. Here, we investigate whether SDM is associated with diversification in turtles and lizards, and whether alterative factors, such as lifespan's effect on transition rates, could explain the relative prevalence of SDMs in turtles and lizards (including and excluding snakes). We assembled a comprehensive dataset of SDM states for squamates and turtles and leveraged large phylogenies for these two groups. We found no evidence that SDMs affect turtle, squamate, or lizard diversification. However, SDM transition rates differ between groups. In lizards TSD-to-GSD surpass GSD-to-TSD transitions, explaining the predominance of GSD lizards in nature. SDM transitions are fewer in turtles and the rates are similar to each other (TSD-to-GSD equals GSD-to-TSD), which, coupled with TSD ancestry, could explain TSD's predominance in turtles. These contrasting patterns can be explained by differences in life history. Namely, our data support the notion that in general, shorter lizard lifespan renders TSD detrimental favoring GSD evolution in squamates, whereas turtle longevity permits TSD retention. Thus, based on the macro-evolutionary evidence we uncovered, we hypothesize that turtles and lizards followed different evolutionary trajectories with respect to SDM, likely mediated by differences in lifespan. Combined, our findings revealed a complex evolutionary interplay between SDMs and life histories that warrants further research that should make use of expanded datasets on unexamined taxa to enable more conclusive analyses.

Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome

Background

Squamates (lizards and snakes) are a speciose lineage of reptiles displaying considerable karyotypic diversity, particularly among lizards. Understanding the evolution of this diversity requires comparison of genome organisation between species. Although the genomes of several squamate species have now been sequenced, only the green anole lizard has any sequence anchored to chromosomes. There is only limited gene mapping data available for five other squamates. This makes it difficult to reconstruct the events that have led to extant squamate karyotypic diversity. The purpose of this study was to anchor the recently sequenced central bearded dragon (Pogona vitticeps) genome to chromosomes to trace the evolution of squamate chromosomes. Assigning sequence to sex chromosomes was of particular interest for identifying candidate sex determining genes.

Results

By using two different approaches to map conserved blocks of genes, we were able to anchor approximately 42 % of the dragon genome sequence to chromosomes. We constructed detailed comparative maps between dragon, anole and chicken genomes, and where possible, made broader comparisons across Squamata using cytogenetic mapping information for five other species. We show that squamate macrochromosomes are relatively well conserved between species, supporting findings from previous molecular cytogenetic studies. Macrochromosome diversity between members of the Toxicofera clade has been generated by intrachromosomal, and a small number of interchromosomal, rearrangements. We reconstructed the ancestral squamate macrochromosomes by drawing upon comparative cytogenetic mapping data from seven squamate species and propose the events leading to the arrangements observed in representative species. In addition, we assigned over 8 Mbp of sequence containing 219 genes to the Z chromosome, providing a list of genes to begin testing as candidate sex determining genes.

Conclusions

Anchoring of the dragon genome has provided substantial insight into the evolution of squamate genomes, enabling us to reconstruct ancestral macrochromosome arrangements at key positions in the squamate phylogeny, demonstrating that fusions between macrochromosomes or fusions of macrochromosomes and microchromosomes, have played an important role during the evolution of squamate genomes. Assigning sequence to the sex chromosomes has identified NR5A1 as a promising candidate sex determining gene in the dragon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2774-3) contains supplementary material, which is available to authorized users.

Mammalian X homolog acts as sex chromosome in lacertid lizards

Among amniotes, squamate reptiles are especially variable in their mechanisms of sex determination; however, based largely on cytogenetic data, some lineages possess highly evolutionary stable sex chromosomes. The still very limited knowledge of the genetic content of squamate sex chromosomes precludes a reliable reconstruction of the evolutionary history of sex determination in this group and consequently in all amniotes. Female heterogamety with a degenerated W chromosome typifies the lizards of the family Lacertidae, the widely distributed Old World clade including several hundreds of species. From the liver transcriptome of the lacertid Takydromus sexlineatus female, we selected candidates for Z-specific genes as the loci lacking single-nucleotide polymorphisms. We validated the candidate genes through the comparison of the copy numbers in the female and male genomes of T. sexlineatus and another lacertid species, Lacerta agilis, by quantitative PCR that also proved to be a reliable technique for the molecular sexing of the studied species. We suggest that this novel approach is effective for the detection of Z-specific and X-specific genes in lineages with degenerated W, respectively Y chromosomes. The analyzed gene content of the Z chromosome revealed that lacertid sex chromosomes are not homologous with those of other reptiles including birds, but instead the genes have orthologs in the X-conserved region shared by viviparous mammals. It is possible that this part of the vertebrate genome was independently co-opted for the function of sex chromosomes in viviparous mammals and lacertids because of its content of genes involved in gonad differentiation.