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

Pivotal temperature is not for everyone: Evidence for temperature-dependent sex determination in three gecko species

The prevalence of environmental sex determination (ESD) in squamate reptiles is often overestimated in the literature. This is surprising because we have reliable data demonstrating ESD in only a few species. The documentation of ESD in three species of geckos presented here has significantly increased our knowledge, given that satisfactory evidence for ESD existed in only eight other gecko species. For the first time, we document the occurrence of ESD in the family Sphaerodactylidae. Our finding of unexpected variability in the shapes of reaction norms among geckos highlights that traditional descriptions using parameters such as pivotal temperature, that is, temperature producing a 50:50 sex ratio, are unsatisfactory. For example, the gecko Pachydactylus tigrinus lacks any pivotal temperature and its sex ratios are strongly female-biased across the entire range of viable temperatures. We argue for the effective capture of the relationship between temperature and sex ratio using specific nonlinear models rather than using classical simplistic descriptions and classifications of reaction norms.

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.

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.

A Brief Review of Meiotic Chromosomes in Early Spermatogenesis and Oogenesis and Mitotic Chromosomes in the Viviparous Lizard Zootoca vivipara (Squamata: Lacertidae) with Multiple Sex Chromosomes

This brief review is focused on the viviparous lizard Zootoca vivipara (Lichtenstein, 1823), of the family Lacertidae, which possesses female heterogamety and multiple sex chromosomes (male 2n = 36, Z₁Z₁Z₂Z₂/Z₁Z₂W, female 2n = 35, with variable W sex chromosome). Multiple sex chromosomes and their changes may influence meiosis and the female meiotic drive, and they may play a role in reproductive isolation. In two cryptic taxa of Z. vivipara with different W sex chromosomes, meiosis during early spermatogenesis and oogenesis proceeds normally, without any disturbances, with the formation of haploid spermatocytes, and in female meiosis with the formation of synaptonemal complexes (SCs) and the lampbrush chromosomes. In females, the SC number was constantly equal to 19 (according to the SC length, 16 SC autosomal bivalents plus three presumed SC sex chromosome elements). No variability in the chromosomes at the early stages of meiotic prophase I, and no significant disturbances in the chromosome segregation at the anaphase-telophase I stage, have been discovered, and haploid oocytes (n = 17) at the metaphase II stage have been revealed. There should be a factor/factors that maintain the multiple sex chromosomes, their equal transmission, and the course of meiosis in these cryptic forms of Z. vivipara.

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.

A novel neo-sex chromosome in Sylvietta brachyura (Macrosphenidae) adds to the extraordinary avian sex chromosome diversity among Sylvioidea songbirds

We report the discovery of a novel neo-sex chromosome in an African warbler, Sylvietta brachyura (northern crombec; Macrosphenidae). This species is part of the Sylvioidea superfamily, where four separate autosome-sex chromosome translocation events have previously been discovered via comparative genomics of 11 of the 22 families in this clade. Our discovery here resulted from analyses of genomic data of single species-representatives from three additional Sylvioidea families (Macrosphenidae, Pycnonotidae and Leiothrichidae). In all three species, we confirmed the translocation of a part of chromosome 4A to the sex chromosomes, which originated basally in Sylvioidea. In S. brachyura, we found that a part of chromosome 8 has been translocated to the sex chromosomes, forming a unique neo-sex chromosome in this lineage. Furthermore, the non-recombining part of 4A in S. brachyura is smaller than in other Sylvioidea species, which suggests that recombination continued along this region after the fusion event in the Sylvioidea ancestor. These findings reveal additional sex chromosome diversity among the Sylvioidea, where five separate translocation events are now confirmed.

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.

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.

Cytogenetically Elusive Sex Chromosomes in Scincoidean Lizards

The lizards of the species-rich clade Scincoidea including cordylids, gerrhosaurids, skinks, and xantusiids, show an almost cosmopolitan geographical distribution and a remarkable ecological and morphological divergence. However, previous studies revealed limited variability in cytogenetic traits. The sex determination mode was revealed only in a handful of gerrhosaurid, skink, and xantusiid species, which demonstrated either ZZ/ZW or XX/XY sex chromosomes. In this study, we explored the karyotypes of six species of skinks, two species of cordylids, and one gerrhosaurid. We applied conventional and molecular cytogenetic methods, including C-banding, fluorescence in situ hybridization with probes specific for telomeric motifs and rDNA loci, and comparative genomic hybridization. The diploid chromosome numbers are rather conserved among these species, but the chromosome morphology, the presence of interstitial telomeric sequences, and the topology of rDNA loci vary significantly. Notably, XX/XY sex chromosomes were identified only in Tiliqua scincoides, where, in contrast to the X chromosome, the Y chromosome lacks accumulations of rDNA loci. We confirm that within the lizards of the scincoidean clade, sex chromosomes remained in a generally poor stage of differentiation.

Chromosome-Centric View of Genome Organization and Evolution

Genetic material in all cellular organisms is packed into chromosomes, which represent essential units of inheritance, recombination, and evolution [...].

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.

Revisiting the Karyotypes of Alligators and Caimans (Crocodylia, Alligatoridae) after a Half-Century Delay: Bridging the Gap in the Chromosomal Evolution of Reptiles

Although crocodilians have attracted enormous attention in other research fields, from the cytogenetic point of view, this group remains understudied. Here, we analyzed the karyotypes of eight species formally described from the Alligatoridae family using differential staining, fluorescence in situ hybridization with rDNA and repetitive motifs as a probe, whole chromosome painting (WCP), and comparative genome hybridization. All Caimaninae species have a diploid chromosome number (2n) 42 and karyotypes dominated by acrocentric chromosomes, in contrast to both species of Alligatorinae, which have 2n = 32 and karyotypes that are predominantly metacentric, suggesting fusion/fission rearrangements. Our WCP results supported this scenario by revealing the homeology of the largest metacentric pair present in both Alligator spp. with two smaller pairs of acrocentrics in Caimaninae species. The clusters of 18S rDNA were found on one chromosome pair in all species, except for Paleosuchus spp., which possessed three chromosome pairs bearing these sites. Similarly, comparative genomic hybridization demonstrated an advanced stage of sequence divergence among the caiman genomes, with Paleosuchus standing out as the most divergent. Thus, although Alligatoridae exhibited rather low species diversity and some level of karyotype stasis, their genomic content indicates that they are not as conserved as previously thought. These new data deepen the discussion of cytotaxonomy in this family.

Anhydroicaritin Inhibits EMT in Breast Cancer by Enhancing GPX1 Expression: A Research Based on Sequencing Technologies and Bioinformatics Analysis

Background: Breast cancer (BC) is the leading cause of cancer-related deaths among women worldwide. The application of advanced technology has promoted accurate diagnosis and treatment of cancer. Anhydroicaritin (AHI) is a flavonoid with therapeutic potential in BC treatment. The current study aimed to determine AHI's mechanism in BC treatment via RNA sequencing, comprehensive bioinformatics analysis, and experimental verification. Methods: Network pharmacology and MTT (3-(4,5)-dimethylthiazolyl-3,5- diphenyltetrazolium bromide) experiments were conducted to first confirm AHI's anti-BC effect. RNA sequencing was performed to identify the genes affected by AHI. Differential expression analysis, survival analysis, gene set enrichment analysis, and immune infiltration analysis were performed via bioinformatics analysis. Western blot analysis, reverse transcription-polymerase chain reaction (RT-PCR) experiment, molecular docking, and drug affinity responsive target stability (DARTS) experiments were also performed to confirm AHI's direct effect on glutathione peroxidase 1 (GPX1) expression. Confocal immunofluorescence analysis was conducted to verify AHI's effect on the occurrence and development of epithelial-mesenchymal transition (EMT). Finally, BC nude mouse xenografts were established, and AHI's molecular mechanism on BC was explored. Results: Network pharmacology results demonstrated that AHI's therapeutic targets on BC were related to the proliferation, invasion, and metastasis of BC cells. AHI significantly inhibited the proliferation of 4T1 and MDA-MB-231 BC cells in the MTT experiments. RNA sequencing results showed that AHI upregulated the GPX1 expression in the 4T1 and MDA-MB-231 BC cells. Next, bioinformatics analysis revealed that GPX1 is less expressed in BC than in normal breast tissues. Patients with high GPX1 expression levels tended to have prolonged overall survival and disease-free survival than patients with low GPX1 expression levels in BC. Western blot and RT-PCR experiments revealed that AHI increased the protein and mRNA levels of GPX1. Molecular docking and DARTS experiments confirmed the direct binding combination between AHI and GPX1. After the evaluation of the EMT scores of 1,078 patients with BC, we found a potential anti-BC role of GPX1 possibly via suppression of the malignant EMT. The confocal immunofluorescence analysis showed that AHI increased E-cadherin expression levels and reduced vimentin expression levels in BC cells. Animal experiments showed that AHI significantly inhibited tumor growth. AHI also inhibited EMT by enhancing GPX1 and caspase3 cleavage, hence inhibiting EMT markers (i.e., N-cadherin and vimentin) and Ki-67. Conclusion: GPX1 plays a critical role in BC, which may be a biomarker for the prognosis. In addition, AHI suppressed EMT by increasing GPX1 expression, which may serve as a potential therapy for BC treatment.

Molecular Cytogenetic Characterization of the Sicilian Endemic Pond Turtle Emys trinacris and the Yellow-Bellied Slider Trachemys scripta scripta (Testudines, Emydidae)

Turtles, a speciose group consisting of more than 300 species, demonstrate karyotypes with diploid chromosome numbers ranging from 2n = 26 to 2n = 68. However, cytogenetic analyses have been conducted only to 1/3rd of the turtle species, often limited to conventional staining methods. In order to expand our knowledge of the karyotype evolution in turtles, we examined the topology of the (TTAGGG)_n telomeric repeats and the rDNA loci by fluorescence in situ hybridization (FISH) on the karyotypes of two emydids: the Sicilian pond turtle, Emys trinacris, and the yellow-bellied slider, Trachemys scripta scripta (family Emydidae). Furthermore, AT-rich and GC-rich chromosome regions were detected by DAPI and CMA₃ stains, respectively. The cytogenetic analysis revealed that telomeric sequences are restricted to the terminal ends of all chromosomes and the rDNA loci are localized in one pair of microchromosomes in both species. The karyotype of the Sicilian endemic E. trinacris with diploid number 2n = 50, consisting of 13 pairs of macrochromosomes and 12 pairs of microchromosomes, is presented here for first time. Our comparative examination revealed similar cytogenetic features in Emys trinacris and the closely related E. orbicularis, as well as to other previously studied emydid species, demonstrating a low rate of karyotype evolution, as chromosomal rearrangements are rather infrequent in this group of turtles.

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.

Evolutionary Variability of W-Linked Repetitive Content in Lacertid Lizards

Lacertid lizards are a widely radiated group of squamate reptiles with long-term stable ZZ/ZW sex chromosomes. Despite their family-wide homology of Z-specific gene content, previous cytogenetic studies revealed significant variability in the size, morphology, and heterochromatin distribution of their W chromosome. However, there is little evidence about the accumulation and distribution of repetitive content on lacertid chromosomes, especially on their W chromosome. In order to expand our knowledge of the evolution of sex chromosome repetitive content, we examined the topology of telomeric and microsatellite motifs that tend to often accumulate on the sex chromosomes of reptiles in the karyotypes of 15 species of lacertids by fluorescence in situ hybridization (FISH). The topology of the above-mentioned motifs was compared to the pattern of heterochromatin distribution, as revealed by C-banding. Our results show that the topologies of the examined motifs on the W chromosome do not seem to follow a strong phylogenetic signal, indicating independent and species-specific accumulations. In addition, the degeneration of the W chromosome can also affect the Z chromosome and potentially also other parts of the genome. Our study provides solid evidence that the repetitive content of the degenerated sex chromosomes is one of the most evolutionary dynamic parts of the genome.