Female heterogamety in Madagascar chameleons (Squamata: Chamaeleonidae: Furcifer): differentiation of sex and neo-sex chromosomes

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

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

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.

Microchromosome fusions underpin convergent evolution of chameleon karyotypes

Evolutionary shifts in chromosome compositions (karyotypes) are major drivers of lineage and genomic diversification. Fusion of ancestral chromosomes is one hypothesized mechanism for the evolutionary reduction of the total chromosome number, a frequently implied karyotypic shift. Empirical tests of this hypothesis require model systems with variable karyotypes, known chromosome features, and a robust phylogeny. Here we used chameleons, diverse lizards with exceptionally variable karyotypes ($2n=20\text{-}62$), to test whether chromosomal fusions explain the repeated evolution of karyotypes with fewer chromosomes than ancestral karyotypes. Using a multidisciplinary approach including cytogenetic analyses and phylogenetic comparative methods, we found that a model of constant loss through time best explained chromosome evolution across the chameleon phylogeny. Next, we tested whether fusions of microchromosomes into macrochromosomes explained these evolutionary losses using generalized linear models. Multiple comparisons supported microchromosome fusions as the predominant agent of evolutionary loss. We further compared our results to various natural history traits and found no correlations. As such, we infer that the tendency of microchromosomes to fuse was a quality of the ancestral chameleon genome and that the genomic predisposition of ancestors is a more substantive predictor of chromosome change than the ecological, physiological, and biogeographical factors involved in their diversification.

Efficient and Highly Continuous Chromosome-Level Genome Assembly of the First Chameleon Genome

Most amniote genomes are diploid and moderate in size (∼1-6 Gbp) and contain a large proportion of repetitive sequences. The development of next-generation sequencing technology, especially the emergence of high-fidelity (HiFi) long-read data, has made it feasible to resolve high-quality genome assembly for nonmodel species efficiently. However, reference genomes for squamate reptiles have lagged behind other amniote lineages. Here, we de novo assembled the first genome from the Chameleonidae family, the panther chameleon (Furcifer pardalis). We obtained telomere-to-telomere contigs using only HiFi data, reaching a contig N50 of 158.72 Mbp. The final chromosome-level assembly is 1.61 Gbp in size, and 100% of primary contigs were placed to pseudochromosomes using Hi-C interaction data. We also found that sequencing depth > 30 folds can ensure both the integrity and accuracy of the genome, whereas insufficient depth led to false increase in genome size and proportion of duplicated genes. We provide a high-quality reference genome valuable for evolutionary and ecological studies in chameleons as well as provide comparative genomic resources for squamate reptiles.

Madagascar Leaf-Tail Geckos (Uroplatus spp.) Share Independently Evolved Differentiated ZZ/ZW Sex Chromosomes

Geckos are an excellent group to study the evolution of sex determination, as they possess a remarkable variability ranging from a complete absence of sex chromosomes to highly differentiated sex chromosomes. We explored sex determination in the Madagascar leaf-tail geckos of the genus Uroplatus. The cytogenetic analyses revealed highly heterochromatic W chromosomes in all three examined species (Uroplatus henkeli, U. alluaudi, U. sikorae). The comparative gene coverage analysis between sexes in U. henkeli uncovered an extensive Z-specific region, with a gene content shared with the chicken chromosomes 8, 20, 26 and 28. The genomic region homologous to chicken chromosome 28 has been independently co-opted for the role of sex chromosomes in several vertebrate lineages, including monitors, beaded lizards and monotremes, perhaps because it contains the amh gene, whose homologs were repeatedly recruited as a sex-determining locus. We demonstrate that all tested species of leaf-tail geckos share homologous sex chromosomes despite the differences in shape and size of their W chromosomes, which are not homologous to the sex chromosomes of other closely related genera. The rather old (at least 40 million years), highly differentiated sex chromosomes of Uroplatus geckos can serve as a great system to study the convergence of sex chromosomes evolved from the same genomic region.

ZZ/ZW Sex Chromosomes in the Madagascar Girdled Lizard, Zonosaurus madagascariensis (Squamata: Gerrhosauridae)

Scincoidea, the reptilian clade that includes girdled lizards, night lizards, plated lizards and skinks, are considered as a lineage with diversity in sex-determining systems. Recently, the hypothesis on the variability in sex determination in skinks and even more the absence of sex chromosomes in some of them has been rivalling. Homologous, evolutionary stable XX/XY sex chromosomes were documented to be widespread across skinks. However, sex determination in the other scincoidean families is highly understudied. ZZ/ZW sex chromosomes have been identified only in night lizards and a single species of plated lizards. It seems that although there are different sex chromosome systems among scincoidean lineages, they share one common trait: they are mostly poorly differentiated and often undetectable by cytogenetic methods. Here, we report one of the exceptions, demonstrating for the first time ZZ/ZW sex chromosomes in the plated lizard Zonosaurus madagascariensis. Its sex chromosomes are morphologically similar, but the W is clearly detectable by comparative genomic hybridization (CGH), suggesting that the Z and W chromosomes highly differ in sequence content. Our findings confirm the presence of female heterogamety in plated lizards and provides novel insights to expand our understanding of sex chromosome evolution in scincoidean lizards.

Identification of Iguania Ancestral Syntenic Blocks and Putative Sex Chromosomes in the Veiled Chameleon (Chamaeleo calyptratus, Chamaeleonidae, Iguania)

The veiled chameleon (Chamaeleo calyptratus) is a typical member of the family Chamaeleonidae and a promising object for comparative cytogenetics and genomics. The karyotype of C. calyptratus differs from the putative ancestral chameleon karyotype (2n = 36) due to a smaller chromosome number (2n = 24) resulting from multiple chromosome fusions. The homomorphic sex chromosomes of an XX/XY system were described recently using male-specific RADseq markers. However, the chromosomal pair carrying these markers was not identified. Here we obtained chromosome-specific DNA libraries of C. calyptratus by chromosome flow sorting that were assigned by FISH and sequenced. Sequence comparison with three squamate reptiles reference genomes revealed the ancestral syntenic regions in the C. calyptratus chromosomes. We demonstrated that reducing the chromosome number in the C. calyptratus karyotype occurred through two fusions between microchromosomes and four fusions between micro-and macrochromosomes. PCR-assisted mapping of a previously described Y-specific marker indicates that chromosome 5 may be the sex chromosome pair. One of the chromosome 5 conserved synteny blocks shares homology with the ancestral pleurodont X chromosome, assuming parallelism in the evolution of sex chromosomes from two basal Iguania clades (pleurodonts and acrodonts). The comparative chromosome map produced here can serve as the foundation for future genome assembly of chameleons and vertebrate-wide comparative genomic studies.

First Insights on the Karyotype Diversification of the Endemic Malagasy Leaf-Toed Geckos (Squamata: Gekkonidae: Uroplatus)

Simple Summary

The geckos of the genus Uroplatus include peculiar endemic species to Madagascar. Even though they have been the subject of several morphological and molecular studies, karyological analyses have been performed only on U. phantasticus, leaving the chromosomal diversity of the genus completely unexplored. In this study, we performed a preliminary molecular analysis and a comparative cytogenetic study providing the first karyotype description of eight species of Uroplatus and an assessment of their karyological variability. We found chromosome diversity in the species studied in terms of total chromosome number (2n = 34–38), localization of loci of Nucleolar Organizer Regions (NORs) (alternatively on the 2nd, 6th, 10th or 16th pair), heterochromatin composition and occurrence of heteromorphic sex chromosome pairs. Adding our newly generated data to those available from the literature, we show that in the genus Uroplatus, as well as in a larger group of phylogenetically related gecko genera, chromosome diversification mainly occurred toward a reduction in the chromosome number by means of chromosome fusions and translocation of NOR-bearing chromosomes. We also hypothesize that the diversification of sex chromosome systems occurred independently in different genera.

Abstract

We provide here the first karyotype description of eight Uroplatus species and a characterization of their chromosomal diversity. We performed a molecular taxonomic assessment of several Uroplatus samples using the mitochondrial 12S marker and a comparative cytogenetic analysis with standard karyotyping, silver staining (Ag-NOR) and sequential C-banding + Giemsa, +Chromomycin A3 (CMA₃), +4′,6-diamidino-2-phenylindole (DAPI). We found chromosomal variability in terms of chromosome number (2n = 34–38), heterochromatin composition and number and localization of loci or Nucleolar Organizer Regions (NORs) (alternatively on the 2nd, 6th, 10th or 16th pair). Chromosome morphology is almost constant, with karyotypes composed of acrocentric chromosomes, gradually decreasing in length. C-banding evidenced a general low content of heterochromatin, mostly localized on pericentromeric and telomeric regions. Centromeric bands varied among the species studied, resulting in CMA₃ positive and DAPI negative or positive to both fluorochromes. We also provide evidence of a first putative heteromorphic sex chromosome system in the genus. In fact, in U. alluaudi the 10th pair was highly heteromorphic, with a metacentric, largely heterochromatic W chromosome, which was much bigger than the Z. We propose an evolutionary scenario of chromosome reduction from 2n = 38 to 2n = 34, by means of translocations of microchromosomes on larger chromosomes (often involving the NOR-bearing microchromosomes). Adding our data to those available from the literature, we show that similar processes characterized the evolutionary radiation of a larger gecko clade. Finally, we hypothesize that sex chromosome diversification occurred independently in different genera.

Cytogenetic Analysis of the Members of the Snake Genera Cylindrophis, Eryx, Python, and Tropidophis

The recent discovery of two independently evolved XX/XY sex determination systems in the snake genera Python and Boa sparked a new drive to study the evolution of sex chromosomes in poorly studied lineages of snakes, where female heterogamety was previously assumed. Therefore, we examined seven species from the genera Eryx, Cylindrophis, Python, and Tropidophis by conventional and molecular cytogenetic methods. Despite the fact that these species have similar karyotypes in terms of chromosome number and morphology, we detected variability in the distribution of heterochromatin, telomeric repeats, and rDNA loci. Heterochromatic blocks were mainly detected in the centromeric regions in all species, although accumulations were detected in pericentromeric and telomeric regions in a few macrochromosomes in several of the studied species. All species show the expected topology of telomeric repeats at the edge of all chromosomes, with the exception of Eryx muelleri, where additional accumulations were detected in the centromeres of three pairs of macrochromosomes. The rDNA loci accumulate in one pair of microchromosomes in all Eryx species and in Cylindrophis ruffus, in one macrochromosome pair in Tropidophis melanurus and in two pairs of microchromosomes in Python regius. Sex-specific differences were not detected, suggesting that these species likely have homomorphic, poorly differentiated sex chromosomes.

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.

Sex chromosome evolution among amniotes: is the origin of sex chromosomes non-random?

Sex chromosomes are a great example of a convergent evolution at the genomic level, having evolved dozens of times just within amniotes. An intriguing question is whether this repeated evolution was random, or whether some ancestral syntenic blocks have significantly higher chance to be co-opted for the role of sex chromosomes owing to their gene content related to gonad development. Here, we summarize current knowledge on the evolutionary history of sex determination and sex chromosomes in amniotes and evaluate the hypothesis of non-random emergence of sex chromosomes. The current data on the origin of sex chromosomes in amniotes suggest that their evolution is indeed non-random. However, this non-random pattern is not very strong, and many syntenic blocks representing putatively independently evolved sex chromosomes are unique. Still, repeatedly co-opted chromosomes are an excellent model system, as independent co-option of the same genomic region for the role of sex chromosome offers a great opportunity for testing evolutionary scenarios on the sex chromosome evolution under the explicit control for the genomic background and gene identity. Future studies should use these systems more to explore the convergent/divergent evolution of sex chromosomes. 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)'.

Lizards as Model Organisms of Sex Chromosome Evolution: What We Really Know from a Systematic Distribution of Available Data?

Lizards represent unique model organisms in the study of sex determination and sex chromosome evolution. Among tetrapods, they are characterized by an unparalleled diversity of sex determination systems, including temperature-dependent sex determination (TSD) and genetic sex determination (GSD) under either male or female heterogamety. Sex chromosome systems are also extremely variable in lizards. They include simple (XY and ZW) and multiple (X1X2Y and Z1Z2W) sex chromosome systems and encompass all the different hypothesized stages of diversification of heterogametic chromosomes, from homomorphic to heteromorphic and completely heterochromatic sex chromosomes. The co-occurrence of TSD, GSD and different sex chromosome systems also characterizes different lizard taxa, which represent ideal models to study the emergence and the evolutionary drivers of sex reversal and sex chromosome turnover. In this review, we present a synthesis of general genome and karyotype features of non-snakes squamates and discuss the main theories and evidences on the evolution and diversification of their different sex determination and sex chromosome systems. We here provide a systematic assessment of the available data on lizard sex chromosome systems and an overview of the main cytogenetic and molecular methods used for their identification, using a qualitative and quantitative approach.

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.

Cytogenetic Evidence for Sex Chromosomes and Karyotype Evolution in Anguimorphan Lizards

Anguimorphan lizards are a morphologically variable group of squamate reptiles with a wide geographical distribution. In spite of their importance, they have been cytogenetically understudied. Here, we present the results of the cytogenetic examination of 23 species from five anguimorphan families (Anguidae, Helodermatidae, Shinisauridae, Varanidae and Xenosauridae). We applied both conventional (Giemsa staining and C-banding) and molecular cytogenetic methods (fluorescence in situ hybridization with probes for the telomeric motifs and rDNA loci, comparative genome hybridization), intending to describe the karyotypes of previously unstudied species, to uncover the sex determination mode, and to reveal the distribution of variability in cytogenetic characteristics among anguimorphan lizards. We documented that karyotypes are generally quite variable across anguimorphan lineages, with anguids being the most varying. However, the derived chromosome number of 2n = 40 exhibits a notable long-term evolutionary stasis in monitors. Differentiated ZZ/ZW sex chromosomes were documented in monitors and helodermatids, as well as in the anguids Abronia lythrochila, and preliminary also in Celestus warreni and Gerrhonotus liocephalus. Several other anguimorphan species have likely poorly differentiated sex chromosomes, which cannot be detected by the applied cytogenetic methods, although the presence of environmental sex determination cannot be excluded. In addition, we uncovered a rare case of spontaneous triploidy in a fully grown Varanus primordius.

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.

Poorly differentiated XX/XY sex chromosomes are widely shared across skink radiation

Differentiated sex chromosomes are believed to be evolutionarily stable, while poorly differentiated sex chromosomes are considered to be prone to turnovers. With around 1700 currently known species forming ca 15% of reptile species diversity, skinks (family Scincidae) are a very diverse group of squamates known for their large ecological and morphological variability. Skinks generally have poorly differentiated and cytogenetically indistinguishable sex chromosomes, and their sex determination was suggested to be highly variable. Here, we determined X-linked genes in the common sandfish (Scincus scincus) and demonstrate that skinks have shared the same homologous XX/XY sex chromosomes across their wide phylogenetic spectrum for at least 85 million years, approaching the age of the highly differentiated ZZ/ZW sex chromosomes of birds and advanced snakes. Skinks thus demonstrate that even poorly differentiated sex chromosomes can be evolutionarily stable. The conservation of sex chromosomes across skinks allows us to introduce the first molecular sexing method widely applicable in this group.

Genome Complexity Reduction High-Throughput Genome Sequencing of Green Iguana (Iguana iguana) Reveal a Paradigm Shift in Understanding Sex-Chromosomal Linkages on Homomorphic X and Y Sex Chromosomes

The majority of lizards classified in the superfamily Iguanoidea have an XX/XY sex-determination system in which sex-chromosomal linkage shows homology with chicken (Gallus gallus) chromosome 15 (GGA15). However, the genomics of sex chromosomes remain largely unexplored owing to the presence of homomorphic sex chromosomes in majority of the species. Recent advances in high-throughput genome complexity reduction sequencing provide an effective approach to the identification of sex-specific loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA), and a better understanding of sex chromosome dynamics in Iguanoidea. In this study, we applied Diversity Arrays Technology (DArTseq^TM) in 29 phenotypic sex assignments (14 males and 15 females) of green iguana (Iguana iguana). We confirmed a male heterogametic (XX/XY) sex determination mode in this species, identifying 29 perfectly sex-linked SNP/PA loci and 164 moderately sex-linked SNP/PA loci, providing evidence probably indicative of XY recombination. Three loci from among the perfectly sex-linked SNP/PA loci showed partial homology with several amniote sex chromosomal linkages. The results support the hypothesis of an ancestral super-sex chromosome with overlaps of partial sex-chromosomal linkages. However, only one locus among the moderately sex-linked loci showed homology with GGA15, which suggests that the specific region homologous to GGA15 was located outside the non-recombination region but in close proximity to this region of the sex chromosome in green iguana. Therefore, the location of GGA15 might be further from the putative sex-determination locus in green iguana. This is a paradigm shift in understanding linkages on homomorphic X and Y sex chromosomes. The DArTseq platform provides an easy-to-use strategy for future research on the evolution of sex chromosomes in Iguanoidea, particularly for non-model species with homomorphic or highly cryptic sex chromosomes.

Interstitial Telomeric Repeats Are Rare in Turtles

Telomeres are nucleoprotein complexes protecting chromosome ends in most eukaryotic organisms. In addition to chromosome ends, telomeric-like motifs can be accumulated in centromeric, pericentromeric and intermediate (i.e., between centromeres and telomeres) positions as so-called interstitial telomeric repeats (ITRs). We mapped the distribution of (TTAGGG)_n repeats in the karyotypes of 30 species from nine families of turtles using fluorescence in situ hybridization. All examined species showed the expected terminal topology of telomeric motifs at the edges of chromosomes. We detected ITRs in only five species from three families. Combining our and literature data, we inferred seven independent origins of ITRs among turtles. ITRs occurred in turtles in centromeric positions, often in several chromosomal pairs, in a given species. Their distribution does not correspond directly to interchromosomal rearrangements. Our findings support that centromeres and non-recombining parts of sex chromosomes are very dynamic genomic regions, even in turtles, a group generally thought to be slowly evolving. However, in contrast to squamate reptiles (lizards and snakes), where ITRs were found in more than half of the examined species, and birds, the presence of ITRs is generally rare in turtles, which agrees with the expected low rates of chromosomal rearrangements and rather slow karyotype evolution in this group.

Deciphering Ancestral Sex Chromosome Turnovers Based on Analysis of Male Mutation Bias

The age of sex chromosomes is commonly obtained by comparing the substitution rates of XY gametologs. Coupled with phylogenetic reconstructions, one can refine the origin of a sex chromosome system relative to specific speciation events. However, these approaches are insufficient to determine the presence and duration of ancestral sex chromosome systems that were lost in some species. In this study, we worked with genomic and transcriptomic data from mammals and squamates and analyzed the effect of male mutation bias on X-linked sequences in these groups. We searched for signatures indicating whether monotremes shared the same sex chromosomes with placental mammals or whether pleurodonts and acrodonts had a common ancestral sex chromosome system. Our analyses indicate that platypus did not share the XY chromosomes with placental mammals, in agreement with previous work. In contrast, analyses of agamids showed that this lineage maintained the pleurodont XY chromosomes for several million years. We performed multiple simulations using different strengths of male mutation bias to confirm the results. Overall, our work shows that variations in substitution rates due to male mutation bias could be applied to uncover signatures of ancestral sex chromosome systems.

Sex is determined by XX/XY sex chromosomes in Australasian side-necked turtles (Testudines: Chelidae)

Turtles demonstrate variability in sex determination and, hence, constitute an excellent model for the evolution of sex chromosomes. Notably, the sex determination of the freshwater turtles from the family Chelidae, a species-rich group with wide geographical distribution in the southern hemisphere, is still poorly explored. Here we documented the presence of an XX/XY sex determination system in seven species of the Australasian chelid genera Chelodina, Emydura, and Elseya by conventional (karyogram reconstruction, C-banding) and molecular cytogenetic methods (comparative genome hybridization, in situ hybridization with probes specific for GATA microsatellite motif, the rDNA loci, and the telomeric repeats). The sex chromosomes are microchromosomes in all examined species of the genus Chelodina. In contrast, the sex chromosomes are the 4^th largest pair of macrochromosomes in the genera Emydura and Elseya. Their X chromosomes are submetacentric, while their Y chromosomes are metacentric. The chelid Y chromosomes contain a substantial male-specific genomic region with an accumulation of the GATA microsatellite motif, and occasionally, of the rDNA loci and telomeric repeats. Despite morphological differences between sex chromosomes, we conclude that male heterogamety was likely already present in the common ancestor of Chelodina, Emydura and Elseya in the Mesozoic period.

Karyological characterization of the common chameleon (Chamaeleo chamaeleon) provides insights on the evolution and diversification of sex chromosomes in Chamaeleonidae

Chameleons display high karyological diversity in chromosome number (from 2n = 20 to 62), morphology, heterochromatin distribution and location of specific chromosomal markers, making them unique study models in evolutionary cytogenetics. However, most available cytogenetic data are limited to the description of the chromosome number and morphology. Concerning sex chromosomes, our knowledge is limited to ZZ/ZW and Z1Z1Z2Z2/Z1Z2W systems in the genus Furcifer and the isolation of sex-linked, male-specific, sequences in Chamaeleo calyptratus, but the putative XY chromosomes have still to be identified in Chamaeleo and the conservation of male heterogamety in the genus needs confirmation from other species. In this study we performed a molecular and a cytogenetic analysis on C. chamaeleon, using standard, banding methods and molecular cytogenetics to provide a throughout karyological characterization of the species and to identify and locate the putative XY chromosomes. We confirm that the chromosome formula of the species is 2n = 24, with 12 metacentric macrochromosomes, 12 microchromosomes and NORs on the second chromosome pair. Heterochromatin was detected as weak C-bands on centromeric regions, differently from what was previously reported for C. calyptratus. Fluorescence in situ hybridization (FISH) showed the occurrence of interspersed telomeric signals on most macrochromosomes, suggesting that ancient chromosome fusions may have led to a reduction of the chromosome number. Using a combination of molecular and FISH analyses, we proved that male specific Restriction site-Associated DNA sequences (RADseq) isolated in C. calyptratus are conserved in C. chamaeleon and located the putative XY chromosomes on the second chromosome pair. We also identified different transposable elements in the focal taxa, which are highly interspersed on most chromosome pairs.

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.

Cytogenetic Analysis Did Not Reveal Differentiated Sex Chromosomes in Ten Species of Boas and Pythons (Reptilia: Serpentes)

Homologous and differentiated ZZ/ZW sex chromosomes (or derived multiple neo-sex chromosomes) were often described in caenophidian snakes, but sex chromosomes were unknown until recently in non-caenophidian snakes. Previous studies revealed that two species of boas (Boa imperator, B. constrictor) and one species of python (Python bivittatus) independently evolved XX/XY sex chromosomes. In addition, heteromorphic ZZ/ZW sex chromosomes were recently revealed in the Madagascar boa (Acrantophis sp. cf. dumerili) and putatively also in the blind snake Myriopholis macrorhyncha. Since the evolution of sex chromosomes in non-caenophidian snakes seems to be more complex than previously thought, we examined ten species of pythons and boas representing the families Boidae, Calabariidae, Candoiidae, Charinidae, Pythonidae, and Sanziniidae by conventional and molecular cytogenetic methods, aiming to reveal their sex chromosomes. Our results show that all examined species do not possess sex-specific differences in their genomes detectable by the applied cytogenetic methods, indicating the presence of poorly differentiated sex chromosomes or even the absence of sex chromosomes. Interestingly, fluorescence in situ hybridization with telomeric repeats revealed extensive distribution of interstitial telomeric repeats in eight species, which are likely a consequence of intra-chromosomal rearrangements.

Little evidence for switches to environmental sex determination and turnover of sex chromosomes in lacertid lizards

Amniotes possess variability in sex determination, from environmental sex determination (ESD), where no sex chromosomes are present, to genotypic sex determination (GSD) with highly differentiated sex chromosomes. Some evolutionary scenarios postulate high stability of differentiated sex chromosomes and rare transitions from GSD to ESD. However, sex chromosome turnovers and two independent transitions from highly differentiated ZZ/ZW sex chromosomes to ESD were previously reported in the lacertid lizards. Here, we examined the homology of sex chromosomes in the wide phylogenetic spectrum of lacertids and their outgroups by comparing gene copy numbers between sexes in genes previously found to be Z-specific in some lacertids. Our current sampling covers 45 species from 26 genera including lineages supposed to possess a derived sex determining systems. We found that all tested lacertids share homologous differentiated ZZ/ZW sex chromosomes, which were present already in their common ancestor living around 85 million years ago. These differentiated sex chromosomes are not present in amphisbaenians and teiid lizards, the close relatives of lacertids. Our study demonstrates how inaccuracies in data can influence the outcome of phylogenetic reconstructions of evolution of sex determination, in this case they overestimated the number of shifts from GSD to ESD and the rate in turnovers of sex chromosomes.

The transcriptome of the veiled chameleon (Chamaeleo calyptratus): A resource for studying the evolution and development of vertebrates

PURPOSE

The veiled chameleon (Chamaeleo calyptratus) is an emerging model system for studying functional morphology and evolutionary developmental biology (evo-devo). Chameleons possess body plans that are highly adapted to an arboreal life style, featuring laterally compressed bodies, split hands/ft for grasping, a projectile tongue, turreted independently moving eyes, and a prehensile tail. Despite being one of the most phenotypically divergent clades of tetrapods, genomic resources for chameleons are severely lacking.

METHODS

To address this lack of resources, we used RNAseq to generate 288 million raw Illumina sequence reads from four adult tissues (male and female eyes and gonads) and whole embryos at three distinct developmental stages. We used these data to assemble a largely complete de novo transcriptome consisting of only 82 952 transcripts. In addition, a majority of assembled transcripts (67%) were successfully annotated.

RESULTS

We then demonstrated the utility of these data in the context of studying visual system evolution by examining the content of veiled chameleon opsin genes to show that chameleons possess all five ancestral tetrapod opsins.

CONCLUSION

We present this de novo, annotated, multi-tissue transcriptome assembly for the Veiled Chameleon, Chamaeleo calyptratus, as a resource to address a range of evolutionary and developmental questions. The associated raw reads and final annotated transcriptome assembly are freely available for use on NCBI and Figshare, respectively.

Turtles of the genera Geoemyda and Pangshura (Testudines: Geoemydidae) lack differentiated sex chromosomes: the end of a 40-year error cascade for Pangshura

For a long time, turtles of the family Geoemydidae have been considered exceptional because representatives of this family were thought to possess a wide variety of sex determination systems. In the present study, we cytogenetically studied Geoemyda spengleri and G. japonica and re-examined the putative presence of sex chromosomes in Pangshura smithii. Karyotypes were examined by assessing the occurrence of constitutive heterochromatin, by comparative genome hybridization and in situ hybridization with repetitive motifs, which are often accumulated on differentiated sex chromosomes in reptiles. We found similar karyotypes, similar distributions of constitutive heterochromatin and a similar topology of tested repetitive motifs for all three species. We did not detect differentiated sex chromosomes in any of the species. For P. smithii, a ZZ/ZW sex determination system, with differentiated sex chromosomes, was described more than 40 years ago, but this finding has never been re-examined and was cited in all reviews of sex determination in reptiles. Here, we show that the identification of sex chromosomes in the original report was based on the erroneous pairing of chromosomes in the karyogram, causing over decades an error cascade regarding the inferences derived from the putative existence of female heterogamety in geoemydid turtles.

Transitions in sex determination and sex chromosomes across vertebrate species

Despite the prevalence of sexual reproduction across eukaryotes, there is a remarkable diversity of sex-determination mechanisms. The underlying causes of this diversity remain unclear, and it is unknown whether there are convergent trends in the directionality of turnover in sex-determination mechanisms. We used the recently assembled Tree of Sex database to assess patterns in the evolution of sex-determination systems in the remarkably diverse vertebrate clades of teleost fish, squamate reptiles and amphibians. Contrary to theoretical predictions, we find no evidence that the evolution of separate sexes is irreversible, as transitions from separate sexes to hermaphroditism occur at higher rates than the reverse in fish. We also find that transitions from environmental sex determination to genetic sex determination occur at higher rates than the reverse in both squamates and fish, suggesting that genetic sex determination is more stable. However, our data are not consistent with the hypothesis that heteromorphic sex chromosomes are an "evolutionary trap." Rather, we find similar transition rates between homomorphic and heteromorphic sex chromosomes in both fish and amphibians, and to environmental sex determination from heteromorphic vs. homomorphic sex chromosome systems in fish. Finally, we find that transitions between male and female heterogamety occur at similar rates in amphibians and squamates, while transitions to male heterogamety occur at higher rates in fish. Together, these results provide the most comprehensive view to date of the evolution of vertebrate sex determination in a phylogenetic context, providing new insight into long-standing questions about the evolution of sexual reproduction.

ZW, XY, and yet ZW: Sex chromosome evolution in snakes even more complicated

Snakes are historically important in the formulation of several central concepts on the evolution of sex chromosomes. For over 50 years, it was believed that all snakes shared the same ZZ/ZW sex chromosomes, which are homomorphic and poorly differentiated in "basal" snakes such as pythons and boas, while heteromorphic and well differentiated in "advanced" (caenophidian) snakes. Recent molecular studies revealed that differentiated sex chromosomes are indeed shared among all families of caenophidian snakes, but that boas and pythons evolved likely independently male heterogamety (XX/XY sex chromosomes). The historical report of heteromorphic ZZ/ZW sex chromosomes in a boid snake was previously regarded as ambiguous. In the current study, we document heteromorphic ZZ/ZW sex chromosomes in a boid snake. A comparative approach suggests that these heteromorphic sex chromosomes evolved very recently and that they are poorly differentiated at the sequence level. Interestingly, two snake lineages with confirmed male heterogamety possess homomorphic sex chromosomes, but heteromorphic sex chromosomes are present in both snake lineages with female heterogamety. We point out that this phenomenon is more common across squamates. The presence of female heterogamety in non-caenophidian snakes indicates that the evolution of sex chromosomes in this lineage is much more complex than previously thought, making snakes an even better model system for the evolution of sex chromosomes.

Evolutionary insight on localization of 18S, 28S rDNA genes on homologous chromosomes in Primates genomes

We explored the topology of 18S and 28S rDNA units by fluorescence in situ hybridization (FISH) in the karyotypes of thirteen species representatives from major groups of Primates and Tupaia minor (Günther, 1876) (Scandentia), in order to expand our knowledge of Primate genome reshuffling and to identify the possible dispersion mechanisms of rDNA sequences. We documented that rDNA probe signals were identified on one to six pairs of chromosomes, both acrocentric and metacentric ones. In addition, we examined the potential homology of chromosomes bearing rDNA genes across different species and in a wide phylogenetic perspective, based on the DAPI-inverted pattern and their synteny to human. Our analysis revealed an extensive variability in the topology of the rDNA signals across studied species. In some cases, closely related species show signals on homologous chromosomes, thus representing synapomorphies, while in other cases, signal was detected on distinct chromosomes, leading to species specific patterns. These results led us to support the hypothesis that different mechanisms are responsible for the distribution of the ribosomal DNA cluster in Primates.

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.

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.

Differentiation of Sex Chromosomes and Karyotype Characterisation in the Dragonsnake Xenodermus javanicus (Squamata: Xenodermatidae)

Highly differentiated heteromorphic ZZ/ZW sex chromosomes with a heterochromatic W are a basic principle among advanced snakes of the lineage Colubroidea, while other snake lineages generally lack these characteristics. For the first time, we cytogenetically examined the dragonsnake, Xenodermus javanicus, a member of the family Xenodermatidae, which is phylogenetically nested between snake lineages with and without differentiated sex chromosomes. Although most snakes have a karyotype with a stable chromosomal number of 2n = 36, the dragonsnake has an unusual, derived karyotype with 2n = 32 chromosomes. We found that heteromorphic ZZ/ZW sex chromosomes with a heterochromatic W are present in the dragonsnake, which suggests that the emergence of a highly differentiated W sex chromosome within snakes predates the split of Xenodermatidae and the clade including families Pareatidae, Viperidae, Homalopsidae, Lamprophiidae, Elapidae, and Colubridae. Although accumulations of interstitial telomeric sequences have not been previously reported in snakes, by using FISH with a telomeric probe we discovered them in 6 pairs of autosomes as well as in the W sex chromosome of the dragonsnake. Similarly to advanced snakes, the sex chromosomes of the dragonsnake have a significant accumulation of repeats containing a (GATA)n sequence. The results facilitate the dating of the differentiation of sex chromosomes within snakes back to the split between Xenodermatidae and other advanced snakes, i.e. around 40-75 mya.