Sex chromosome evolution in snakes inferred from divergence patterns of two gametologous genes and chromosome distribution of sex chromosome-linked repetitive sequences

Sex Chromosome Turnovers and Stability in Snakes

For a long time, snakes were presented as a textbook example of a group with gradual differentiation of homologous ZZ/ZW sex chromosomes. However, recent advances revealed that the ZZ/ZW sex chromosomes characterize only caenophidian snakes and certain species of boas and pythons have nonhomologous XX/XY sex chromosomes. We used genome coverage analysis in four non-caenophidian species to identify their sex chromosomes, and we examined the homology of sex chromosomes across phylogenetically informative snake lineages. We identified sex chromosomes for the first time in 13 species of non-caenophidian snakes, providing much deeper insights into the evolutionary history of snake sex chromosomes. The evolution of sex chromosomes in snakes is more complex than previously thought. Snakes may have had ancestral XX/XY sex chromosomes, which are still present in a blind snake and some boas, and there were several transitions to derived XX/XY sex chromosomes with different gene content and two or even three transitions to ZZ/ZW sex chromosomes. However, we discuss more alternative scenarios. In any case, we document that (1) some genomic regions were likely repeatedly co-opted as sex chromosomes in phylogenetically distant lineages, even with opposite types of heterogamety; (2) snake lineages differ greatly in the rate of differentiation of sex chromosomes; (3) snakes likely originally possessed sex chromosomes prone to turnovers. The sex chromosomes became evolutionarily highly stable once their differentiation progressed in the megadiverse caenophidian snakes. Snakes thus provide an ideal system for studying the evolutionary factors that drive unequal rates of differentiation, turnovers and stability of sex chromosomes.

Sex chromosome cycle as a mechanism of stable sex determination

Recent advances in DNA sequencing technology have enabled the precise decoding of genomes in non-model organisms, providing a basis for unraveling the patterns and mechanisms of sex chromosome evolution. Studies of different species have yielded conflicting results regarding the traditional theory that sex chromosomes evolve from autosomes via the accumulation of deleterious mutations and degeneration of the Y (or W) chromosome. The concept of the 'sex chromosome cycle,' emerging from this context, posits that at any stage of the cycle (i.e., differentiation, degeneration, or loss), sex chromosome turnover can occur while maintaining stable sex determination. Thus, understanding the mechanisms that drive both the persistence and turnover of sex chromosomes at each stage of the cycle is crucial. In this review, we integrate recent findings on the mechanisms underlying maintenance and turnover, with a special focus on several organisms having unique sex chromosomes. Our review suggests that the diversity of sex chromosomes in the maintenance of stable sex determination is underappreciated and emphasizes the need for more research on the sex chromosome cycle.

Chromosomal mapping of repetitive DNA and retroelement sequences and its implications for the chromosomal evolution process in Ctenoluciidae (Characiformes)

Ctenoluciidae is a Neotropical freshwater fish family composed of two genera, Ctenolucius (C. beani and C. hujeta) and Boulengerella (B. cuvieri, B. lateristriga, B. lucius, B. maculata, and B. xyrekes), which present diploid number conservation of 36 chromosomes and a strong association of telomeric sequences with ribosomal DNAs. In the present study, we performed chromosomal mapping of microsatellites and transposable elements (TEs) in Boulengerella species and Ctenolucius hujeta. We aim to understand how those sequences are distributed in these organisms' genomes and their influence on the chromosomal evolution of the group. Our results indicate that repetitive sequences may had an active role in the karyotypic diversification of this family, especially in the formation of chromosomal hotspots that are traceable in the diversification processes of Ctenoluciidae karyotypes. We demonstrate that (GATA)n sequences also accumulate in the secondary constriction formed by the 18 S rDNA site, which shows consistent size heteromorphism between males and females in all Boulengerella species, suggesting an initial process of sex chromosome differentiation.

Following the Pathway of W Chromosome Differentiation in Triportheus (Teleostei: Characiformes)

In this work, we trace the dynamics of satellite DNAs (SatDNAs) accumulation and elimination along the pathway of W chromosome differentiation using the well-known Triportheus fish model. Triportheus stands out due to a conserved ZZ/ZW sex chromosome system present in all examined species. While the Z chromosome is conserved in all species, the W chromosome is invariably smaller and exhibits differences in size and morphology. The presumed ancestral W chromosome is comparable to that of T. auritus, and contains 19 different SatDNA families. Here, by examining five additional Triportheus species, we showed that the majority of these repetitive sequences were eliminated as speciation was taking place. The W chromosomes continued degeneration, while the Z chromosomes of some species began to accumulate some TauSatDNAs. Additional species-specific SatDNAs that made up the heterochromatic region of both Z and W chromosomes were most likely amplified in each species. Therefore, the W chromosomes of the various Triportheus species have undergone significant evolutionary changes in a short period of time (15-25 Myr) after their divergence.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

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.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

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

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

Existence of Bov-B LINE Retrotransposons in Snake Lineages Reveals Recent Multiple Horizontal Gene Transfers with Copy Number Variation

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can “jump” and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10⁻⁶% to 5.63 × 10⁻²% per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution.

Genome-wide SNP analysis of Siamese cobra (Naja kaouthia) reveals the molecular basis of transitions between Z and W sex chromosomes and supports the presence of an ancestral super-sex chromosome in amniotes

Elucidation of the process of sex chromosome differentiation is necessary to understand the dynamics of evolutionary mechanisms in organisms. The W sex chromosome of the Siamese cobra (Naja kaouthia) contains a large number of repeats and shares amniote sex chromosomal linkages. Diversity Arrays Technology provides an effective approach to identify sex-specific loci that are epoch-making, to understand the dynamics of molecular transitions between the Z and W sex chromosomes in a snake lineage. From a total of 543 sex-specific loci, 90 showed partial homology with sex chromosomes of several amniotes and 89 loci were homologous to transposable elements. Two loci were confirmed as W-specific nucleotides after PCR amplification. These loci might result from a sex chromosome differentiation process and involve putative sex-determination regions in the Siamese cobra. Sex-specific loci shared linkage homologies among amniote sex chromosomes, supporting an ancestral super-sex chromosome.

Partial Amniote Sex Chromosomal Linkage Homologies Shared on Snake W Sex Chromosomes Support the Ancestral Super-Sex Chromosome Evolution in Amniotes

Squamate reptile chromosome 2 (SR2) is thought to be an important remnant of an ancestral amniote super-sex chromosome, but a recent study showed that the Siamese cobra W sex chromosome is also a part of this larger ancestral chromosome. To confirm the existence of an ancestral amniote super-sex chromosome and understand the mechanisms of amniote sex chromosome evolution, chromosome maps of two snake species [Russell's viper: Daboia russelii (DRU) and the common tiger snake: Notechis scutatus (NSC)] were constructed using bacterial artificial chromosomes (BACs) derived from chicken and zebra finch libraries containing amniote sex chromosomal linkages. Sixteen BACs were mapped on the W sex chromosome of DRU and/or NSC, suggesting that these BACs contained a common genomic region shared with the W sex chromosome of these snakes. Two of the sixteen BACs were co-localized to DRU2 and NSC2, corresponding to SR2. Prediction of genomic content from all BACs mapped on snake W sex chromosomes revealed a large proportion of long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) retrotransposons. These results led us to predict that amplification of LINE and SINE may have occurred on snake W chromosomes during evolution. Genome compartmentalization, such as transposon amplification, might be the key factor influencing chromosome structure and differentiation. Multiple sequence alignments of all BACs mapped on snake W sex chromosomes did not reveal common sequences. Our findings indicate that the SR2 and snake W sex chromosomes may have been part of a larger ancestral amniote super-sex chromosome, and support the view of sex chromosome evolution as a colorful myriad of situations and trajectories in which many diverse processes are in action.

Landscape of snake' sex chromosomes evolution spanning 85 MYR reveals ancestry of sequences despite distinct evolutionary trajectories

Most of snakes exhibit a ZZ/ZW sex chromosome system, with different stages of degeneration. However, undifferentiated sex chromosomes and unique Y sex-linked markers, suggest that an XY system has also evolved in ancestral lineages. Comparative cytogenetic mappings revealed that several genes share ancestry among X, Y and Z chromosomes, implying that XY and ZW may have undergone transitions during serpent's evolution. In this study, we performed a comparative cytogenetic analysis to identify homologies of sex chromosomes across ancestral (Henophidia) and more recent (Caenophidia) snakes. Our analysis suggests that, despite ~ 85 myr of independent evolution, henophidians and caenophidians retained conserved synteny over much of their genomes. However, our findings allowed us to discover that ancestral and recent lineages of snakes do not share the same sex chromosome and followed distinct pathways for sex chromosomes evolution.

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.

Evolutionary Insights of the ZW Sex Chromosomesin Snakes: A New Chapter Added by the AmazonianPuffing Snakes of the Genus Spilotes

Amazonian puffing snakes (Spilotes; Colubridae) are snakes widely distributed in the Neotropical region. However, chromosomal data are scarce in this group and, when available, are only limited to karyotype description using conventional staining. In this paper, we focused on the process of karyotype evolution and trends for sex chromosomes in two Amazonian Puffer Snakes (S. pulllatus and S. sulphureus). We performed an extensive karyotype characterization using conventional and molecular cytogenetic approaches. The karyotype of S. sulphureus (presented here for the first time) exhibits a 2n = 36, similar to that previously described in S. pullatus. Both species have highly differentiated ZZ/ZW sex chromosomes, where the W chromosome is highly heterochromatic in S. pullatus but euchromatic in S. sulphureus. Both W chromosomes are homologous between these species as revealed by cross-species comparative genomic hybridization, even with heterogeneous distributions of several repetitive sequences across their genomes, including on the Z and on the W chromosomes. Our study provides evidence that W chromosomes in these two species have shared ancestry.

Chromosome map of the Siamese cobra: did partial synteny of sex chromosomes in the amniote represent "a hypothetical ancestral super-sex chromosome" or random distribution?

BACKGROUND

Unlike the chromosome constitution of most snakes (2n=36), the cobra karyotype shows a diploid chromosome number of 38 with a highly heterochromatic W chromosome and a large morphologically different chromosome 2. To investigate the process of sex chromosome differentiation and evolution between cobras, most snakes, and other amniotes, we constructed a chromosome map of the Siamese cobra (Naja kaouthia) with 43 bacterial artificial chromosomes (BACs) derived from the chicken and zebra finch libraries using the fluorescence in situ hybridization (FISH) technique, and compared it with those of the chicken, the zebra finch, and other amniotes.

RESULTS

We produced a detailed chromosome map of the Siamese cobra genome, focusing on chromosome 2 and sex chromosomes. Synteny of the Siamese cobra chromosome 2 (NKA2) and NKAZ were highly conserved among snakes and other squamate reptiles, except for intrachromosomal rearrangements occurring in NKA2. Interestingly, twelve BACs that had partial homology with sex chromosomes of several amniotes were mapped on the heterochromatic NKAW as hybridization signals such as repeat sequences. Sequence analysis showed that most of these BACs contained high proportions of transposable elements. In addition, hybridization signals of telomeric repeat (TTAGGG)n and six microsatellite repeat motifs ((AAGG)8, (AGAT)8, (AAAC)8, (ACAG)8, (AATC)8, and (AAAAT)6) were observed on NKAW, and most of these were also found on other amniote sex chromosomes.

CONCLUSIONS

The frequent amplification of repeats might involve heterochromatinization and promote sex chromosome differentiation in the Siamese cobra W sex chromosome. Repeat sequences are also shared among amniote sex chromosomes, which supports the hypothesis of an ancestral super-sex chromosome with overlaps of partial syntenies. Alternatively, amplification of microsatellite repeat motifs could have occurred independently in each lineage, representing convergent sex chromosomal differentiation among amniote sex chromosomes.

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 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.

Global Dosage Compensation Is Ubiquitous in Lepidoptera, but Counteracted by the Masculinization of the Z Chromosome

While chromosome-wide dosage compensation of the X chromosome has been found in many species, studies in ZW clades have indicated that compensation of the Z is more localized and/or incomplete. In the ZW Lepidoptera, some species show complete compensation of the Z chromosome, while others lack full equalization, but what drives these inconsistencies is unclear. Here, we compare patterns of male and female gene expression on the Z chromosome of two closely related butterfly species, Papilio xuthus and Papilio machaon, and in multiple tissues of two moths species, Plodia interpunctella and Bombyx mori, which were previously found to differ in the extent to which they equalize Z-linked gene expression between the sexes. We find that, while some species and tissues seem to have incomplete dosage compensation, this is in fact due to the accumulation of male-biased genes and the depletion of female-biased genes on the Z chromosome. Once this is accounted for, the Z chromosome is fully compensated in all four species, through the up-regulation of Z expression in females and in some cases additional down-regulation in males. We further find that both sex-biased genes and Z-linked genes have increased rates of expression divergence in this clade, and that this can lead to fast shifts in patterns of gene expression even between closely related species. Taken together, these results show that the uneven distribution of sex-biased genes on sex chromosomes can confound conclusions about dosage compensation and that Z chromosome-wide dosage compensation is not only possible but ubiquitous among Lepidoptera.

A novel method of caenophidian snake sex identification using molecular markers based on two gametologous genes

Sex identification provides important information for ecological and evolutionary studies, as well as benefiting snake conservation management. Traditional methods such as cloacal probing or cloacal popping are counterproductive for sex identification concerning very small species, resulting in difficulties in the management of their breeding programs. In this study, the nucleotide sequences of gametologous genes (CTNNB1 and WAC genes) were used for the development of molecular sexing markers in caenophidian snakes. Two candidate markers were developed with the two primer sets, and successfully amplified by a single band on the agarose gel in male (ZZ) and two bands, differing in fragment sizes, in female (ZW) of 16 caenophidian snakes for CTNNB1 and 12 caenophidian snakes for WAC. Another candidate marker was developed with the primer set to amplify the specific sequence for CTNNB1W homolog, and the PCR products were successfully obtained in a female‐specific 250‐bp DNA bands. The three candidate PCR sexing markers provide a simple sex identification method based on the amplification of gametologous genes, and they can be used to facilitate effective caenophidian snake conservation and management programs.