New insights into sex chromosome evolution in anole lizards (Reptilia, Dactyloidae)

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.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

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

Homeology of sex chromosomes in Amazonian Harttia armored catfishes supports the X-fission hypothesis for the X1X2Y sex chromosome system origin

The Neotropical monophyletic catfish genus Harttia represents an excellent model to study karyotype and sex chromosome evolution in teleosts. Its species split into three phylogenetic clades distributed along the Brazilian territory and they differ widely in karyotype traits, including the presence of standard or multiple sex chromosome systems in some members. Here, we investigate the chromosomal rearrangements and associated synteny blocks involved in the origin of a multiple X1X2Y sex chromosome system present in three out of six sampled Amazonian-clade species. Using 5S and 18S ribosomal DNA fluorescence in situ hybridization and whole chromosome painting with probes corresponding to X1 and X2 chromosomes of X1X2Y system from H. punctata, we confirm previous assumptions that X1X2Y sex chromosome systems of H. punctata, H. duriventris and H. villasboas represent the same linkage groups which also form the putative XY sex chromosomes of H. rondoni. The shared homeology between X1X2Y sex chromosomes suggests they might have originated once in the common ancestor of these closely related species. A joint arrangement of mapped H. punctata X1 and X2 sex chromosomes in early diverging species of different Harttia clades suggests that the X1X2Y sex chromosome system may have formed through an X chromosome fission rather than previously proposed Y-autosome fusion.

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.

Chromosome Evolution of the Liolaemus monticola (Liolaemidae) Complex: Chromosomal and Molecular Aspects

Chromosomal rearrangements can directly influence population differentiation and speciation. The Liolaemus monticola complex in Chile is a unique model consisting of several chromosome races arranged in a latitudinal sequence of increasing karyotype complexity from south to north. Here, we compared chromosomal and mitochondrial cytochrome b data from 15 localities across the northern geographic distribution of L. monticola. We expanded the distribution of the previously described Multiple Fissions race (re-described as MF2), in the Coastal range between the Aconcagua River and the Petorca River, and described a new Multiple Fissions 1 (MF1) race in the Andean range. Both races present centric fissions in pairs 1 and 2, as well as a pericentric inversion in one fission product of pair 2 that changes the NOR position. Additionally, we detected a new chromosomal race north of the Petorca River, the Northern Modified 2 (NM2) race, which is polymorphic for novel centric fissions in pairs 3 and 4. Our results increase the number of chromosomal races in L. monticola to seven, suggesting a complex evolutionary history of chromosomal rearrangements, population isolation by barriers, and hybridization. These results show the relevant role of chromosome mutations in evolution, especially for highly speciose groups such as Liolaemus lizards.

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.

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.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish 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.

Whole-chromosome fusions in the karyotype evolution of Sceloporus (Iguania, Reptilia) are more frequent in sex chromosomes than autosomes

Whole-chromosome fusions play a major role in the karyotypic evolution of reptiles. It has been suggested that certain chromosomes tend to fuse with sex chromosomes more frequently than others. However, the comparative genomic synteny data are too scarce to draw strong conclusions. We obtained and sequenced chromosome-specific DNA pools of Sceloporus malachiticus, an iguanian species which has experienced many chromosome fusions. We found that four of seven lineage-specific fusions involved sex chromosomes, and that certain syntenic blocks which constitute the sex chromosomes, such as the homologues of the Anolis carolinensis chromosomes 11 and 16, are repeatedly involved in sex chromosome formation in different squamate species. To test the hypothesis that the karyotypic shift could be associated with changes in recombination patterns, we performed a synaptonemal complex analysis in this species and in Sceloporus variabilis (2n = 34). It revealed that the sex chromosomes in S. malachiticus had two distal pseudoautosomal regions and a medial differentiated region. We found that multiple fusions little affected the recombination rate in S. malachiticus. Our data confirm more frequent involvement of certain chromosomes in sex chromosome formation, but do not reveal a connection between the gonosome-autosome fusions and the evolution of recombination rate. 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)'.

Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae)

Heteromorphic W and Y sex chromosomes often experience gene loss and heterochromatinization, which is frequently viewed as their "degeneration". However, the evolutionary trajectories of the heterochromosomes are in fact more complex since they may not only lose but also acquire new sequences. Previously, we found that the heterochromatic W chromosome of a lizard Eremias velox (Lacertidae) is decondensed and thus transcriptionally active during the lampbrush stage. To determine possible sources of this transcription, we sequenced DNA from a microdissected W chromosome sample and a total female DNA sample and analyzed the results of reference-based and de novo assembly. We found a new repetitive sequence, consisting of fragments of an autosomal protein-coding gene ATF7IP2, several SINE elements, and sequences of unknown origin. This repetitive element is distributed across the whole length of the W chromosome, except the centromeric region. Since it retained only 3 out of 10 original ATF7IP2 exons, it remains unclear whether it is able to produce a protein product. Subsequent studies are required to test the presence of this element in other species of Lacertidae and possible functionality. Our results provide further evidence for the view of W and Y chromosomes as not just "degraded" copies of Z and X chromosomes but independent genomic segments in which novel genetic elements may arise.

Detection of cryptic diversity in lizards (Squamata) from two Biosphere Reserves in Mesoamerica

A combined approach based on karyology and DNA taxonomy allowed us to characterize the taxonomic peculiarities in 10 Mesoamerican lizard species, belonging to six genera and five families, inhabiting two Biosphere Reserve in Chiapas, Mexico: La Sepultura Biosphere Reserve, and Montes Azules Biosphere. The karyotypes of four species, Phyllodactylus sp. 3 (P. tuberculosus species group) (2n = 38), Holcosus festivus (Lichtenstein et von Martens, 1856) (2n = 50), Anolis lemurinus Cope, 1861 (2n = 40), and A. uniformis Cope, 1885 (2n = 29-30) are described for the first time, the last one showing a particular X₁X₁X₂X₂/X₁X₂Y condition. In Aspidoscelis deppii (Wiegmann, 1834) (2n = 50) and Anolis capito Peters, 1863 (2n = 42), we found a different karyotype from the ones previously reported for these species. Moreover, in A. capito, the cytogenetic observation is concurrent with a considerable genetic divergence (9%) at the studied mtDNA marker (MT-ND2), which is indicative of a putative new cryptic species. The skink Scincella cherriei (Cope, 1893), showed high values of genetic divergence (5.2% at 16S gene) between the specimens from Montes Azules and those from Costa Rica and Nicaragua, comparable to the values typical of sister species in skinks. A lower level of genetic divergence, compatible with an intraspecific phylogeographic structure, has been identified in Lepidophyma flavimaculatum Duméril, 1851. These new data identify taxa that urgently require more in-depth taxonomic studies especially in these areas where habitat alteration is proceeding at an alarming rate.

Chromosome Painting Does Not Support a Sex Chromosome Turnover in Lacerta agilis Linnaeus, 1758

Reptiles show a remarkable diversity of sex determination mechanisms and sex chromosome systems, derived from different autosomal pairs. The origin of the ZW sex chromosomes of Lacerta agilis, a widespread Eurasian lizard species, is a matter of discussion: is it a small macrochromosome from the 11-18 group common to all lacertids, or does this species have a unique ZW pair derived from the large chromosome 5? Using independent molecular cytogenetic methods, we investigated the karyotype of L. agilis exigua from Siberia, Russia, to identify the sex chromosomes. FISH with a flow-sorted chromosome painting probe derived from L. strigata and specific to chromosomes 13, 14, and Z confirmed that the Z chromosome of L. agilis is a small macrochromosome, the same as in L. strigata. FISH with the telomeric probe showed an extensive accumulation of the telomere-like repeat in the W chromosome in agreement with previous studies, excluding the possibility that the lineages of L. agilis studied in different works could have different sex chromosome systems due to a putative intra-species polymorphism. Our results reinforce the idea of the stability of the sex chromosomes and lack of evidence for sex-chromosome turnovers in known species of Lacertidae.

Repeated sex chromosome evolution in vertebrates supported by expanded avian sex chromosomes

Sex chromosomes have evolved from the same autosomes multiple times across vertebrates, suggesting that selection for recombination suppression has acted repeatedly and independently on certain genetic backgrounds. Here, we perform comparative genomics of a bird clade (larks and their sister lineage; Alaudidae and Panuridae) where multiple autosome-sex chromosome fusions appear to have formed expanded sex chromosomes. We detected the largest known avian sex chromosome (195.3 Mbp) and show that it originates from fusions between parts of four avian chromosomes: Z, 3, 4A and 5. Within these four chromosomes, we found evidence of five evolutionary strata where recombination had been suppressed at different time points, and show that stratum age explained the divergence rate of Z-W gametologs. Next, we analysed chromosome content and found that chromosome 3 was significantly enriched for genes with predicted sex-related functions. Finally, we demonstrate extensive homology to sex chromosomes in other vertebrate lineages: chromosomes Z, 3, 4A and 5 have independently evolved into sex chromosomes in fish (Z), turtles (Z, 5), lizards (Z, 4A), mammals (Z, 4A) and frogs (Z, 3, 4A, 5). Our results provide insights into and support for repeated evolution of sex chromosomes in vertebrates.

Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes)

Oplegnathus fasciatus and O. punctatus (Teleostei: Centrarchiformes: Oplegnathidae), are commercially important rocky reef fishes, endemic to East Asia. Both species present an X1X2Y sex chromosome system. Here, we investigated the evolutionary forces behind the origin and differentiation of these sex chromosomes, with the aim to elucidate whether they had a single or convergent origin. To achieve this, conventional and molecular cytogenetic protocols, involving the mapping of repetitive DNA markers, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) were applied. Both species presented similar 2n, karyotype structure and hybridization patterns of repetitive DNA classes. 5S rDNA loci, besides being placed on the autosomal pair 22, resided in the terminal region of the long arms of both X1 chromosomes in females, and on the X1 and Y chromosomes in males. Furthermore, WCP experiments with a probe derived from the Y chromosome of O. fasciatus (OFAS-Y) entirely painted the X1 and X2 chromosomes in females and the X1, X2, and Y chromosomes in males of both species. CGH failed to reveal any sign of sequence differentiation on the Y chromosome in both species, thereby suggesting the shared early stage of neo-Y chromosome differentiation. Altogether, the present findings confirmed the origin of the X1X2Y sex chromosomes via Y-autosome centric fusion and strongly suggested their common origin.

Squamate Chromosome Size and GC Content Assessed by Flow Karyotyping

Chromosome homologies in reptiles have been investigated extensively by gene mapping and chromosome painting. Relative chromosome size can be estimated roughly from conventional karyotypes, but chromosome GC content cannot be evaluated by any of these approaches. However, GC content can be obtained by whole-genome sequencing, although complete data are available only for a limited number of reptilian species. Chromosomes can be characterized by size and GC content in bivariate flow karyotypes, in which the distribution of peaks represents the differences. We have analysed flow karyotypes from 9 representative squamate species and show chromosome profiles for each species based on the relationship between size and GC content. Our results reveal that the GC content of macrochromosomes is invariable in the 9 species. A higher GC content was found in microchromosomes, similar to profiles previously determined in crocodile, turtle, and chicken. The findings suggest that karyotype evolution in reptiles is characterized by unique features of chromosome GC content.

First report on B chromosome content in a reptilian species: the case of Anolis carolinensis

Supernumerary elements of the genome are often called B chromosomes. They usually consist of various autosomal sequences and, because of low selective pressure, are mostly pseudogenized and contain many repeats. There are numerous reports on B chromosomes in mammals, fish, invertebrates, plants, and fungi, but only a few of them have been studied using sequencing techniques. However, reptilian supernumerary chromosomes have been detected only cytogenetically and never sequenced or analyzed at the molecular level. One model squamate species with available genome sequence is Anolis carolinensis. The scope of the present article is to describe the genetic content of A. carolinensis supernumerary chromosomes. In this article, we confirm the presence of B chromosomes in this species by reverse painting and synaptonemal complex analysis. We applied low-pass high-throughput sequencing to analyze flow-sorted B chromosomes. Anole B chromosomes exhibit similar traits to other supernumerary chromosomes from different taxons: they contain two genes related to cell division control (INCENP and SPIRE2), are enriched in specific repeats, and show a high degree of pseudogenization. Therefore, the present study confirms that reptilian B chromosomes resemble supernumerary chromosomes of other taxons.

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.

Immunocytological analysis of meiotic recombination in two anole lizards (Squamata, Dactyloidae)

Although the evolutionary importance of meiotic recombination is not disputed, the significance of interspecies differences in the recombination rates and recombination landscapes remains under-appreciated. Recombination rates and distribution of chiasmata have been examined cytologically in many mammalian species, whereas data on other vertebrates are scarce. Immunolocalization of the protein of the synaptonemal complex (SYCP3), centromere proteins and the mismatch-repair protein MLH1 was used, which is associated with the most common type of recombination nodules, to analyze the pattern of meiotic recombination in the male of two species of iguanian lizards, Anolis carolinensis Voigt, 1832 and Deiroptyx coelestinus (Cope, 1862). These species are separated by a relatively long evolutionary history although they retain the ancestral iguanian karyotype. In both species similar and extremely uneven distributions of MLH1 foci along the macrochromosome bivalents were detected: approximately 90% of crossovers were located at the distal 20% of the chromosome arm length. Almost total suppression of recombination in the intermediate and proximal regions of the chromosome arms contradicts the hypothesis that "homogenous recombination" is responsible for the low variation in GC content across the anole genome. It also leads to strong linkage disequilibrium between the genes located in these regions, which may benefit conservation of co-adaptive gene arrays responsible for the ecological adaptations of the anoles.