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

Molecular phylogenic trees of CTNNB1 gene. This figure shows neighbor-joining trees of CTNNB1 gene with the long alignment for 20 tetrapod species and the short alignment for 26 squamate species. (PDF 289 kb

Karyotypic analysis and FISH mapping of microsatellite motifs reveal highly differentiated XX/XY sex chromosomes in the pink-tailed worm-lizard (Aprasia parapulchella, Pygopodidae, Squamata)

BACKGROUND: The infraorder Gekkota is intriguing because it contains multiple chromosomal and environmental sex determination systems that vary even among closely related taxa. Here, we compare male and females karyotypes of the pink-tailed worm-lizard (Aprasia parapulchella), a small legless lizard belonging to the endemic Australian family Pygopodidae. RESULTS: We applied comparative genomic hybridization to reveal an XX/XY sex chromosome system in which the Y chromosome is highly differentiated from the X in both gross morphology and DNA sequence. In addition, FISH mapping has revealed that two microsatellite repeat motifs, (AGAT)n and (AC)n, have been amplified multiple times on the Y chromosome. CONCLUSION: XY karyotypes are found in other pygopodids (Delma inornata and Lialis burtonis), suggesting that the common ancestor of Pygopodidae also had XY sex chromosomes. However, the morphology and size of the Y chromosomes are different among the three species, suggesting that the processes underlying the evolution of sex chromosomes in the Pygopodidae involved chromosome rearrangements and accumulation and amplification of repeats

Highly Differentiated ZW Sex Microchromosomes in the Australian Varanus Species Evolved through Rapid Amplification of Repetitive Sequences

Transitions between sex determination systems have occurred in many lineages of squamates and it follows that novel sex chromosomes will also have arisen multiple times. The formation of sex chromosomes may be reinforced by inhibition of recombination and the accumulation of repetitive DNA sequences. The karyotypes of monitor lizards are known to be highly conserved yet the sex chromosomes in this family have not been fully investigated. Here, we compare male and female karyotypes of three Australian monitor lizards, Varanus acanthurus, V. gouldii and V. rosenbergi, from two different clades. V. acanthurus belongs to the acanthurus clade and the other two belong to the gouldii clade. We applied C-banding and comparative genomic hybridization to reveal that these species have ZZ/ZW sex micro-chromosomes in which the W chromosome is highly differentiated from the Z chromosome. In combination with previous reports, all six Varanus species in which sex chromosomes have been identified have ZZ/ZW sex chromosomes, spanning several clades on the varanid phylogeny, making it likely that the ZZ/ZW sex chromosome is ancestral for this family. However, repetitive sequences of these ZW chromosome pairs differed among species. In particular, an (AAT)n microsatellite repeat motif mapped by fluorescence in situ hybridization on part of W chromosome in V. acanthurus only, whereas a (CGG)n motif mapped onto the W chromosomes of V. gouldii and V. rosenbergi. Furthermore, the W chromosome probe for V. acanthurus produced hybridization signals only on the centromeric regions of W chromosomes of the other two species. These results suggest that the W chromosome sequences were not conserved between gouldii and acanthurus clades and that these repetitive sequences have been amplified rapidly and independently on the W chromosome of the two clades after their divergence

Diversity of reptile sex chromosome evolution revealed by cytogenetic and linked-read sequencing

Reptile sex determination is attracting much attention because the great diversity of sex-determination and dosage compensation mechanisms permits us to approach fundamental questions about mechanisms of sex chromosome turnover. Recent studies have made significant progress in better understanding diversity and conservation of reptile sex chromosomes, with however no reptile master sex determination genes identified. Here we describe an integrated genomics and cytogenetics pipeline, combining probes generated from the microdissected sex chromosomes with transcriptome and genome sequencing to explore the sex chromosome diversity in non-model Australian reptiles. We tested our pipeline on a turtle, two species of geckos, and a monitor lizard. Genes identified on sex chromosomes were compared to the chicken genome to identify homologous regions among the four species. We identified candidate sex determining genes within these regions, including conserved vertebrate sex-determining genes pdgfa, pdgfra amh and wt1, and demonstrated their testis or ovary-specific expression. All four species showed gene-by-gene rather than chromosome-wide dosage compensation. Our results imply that reptile sex chromosomes originated by independent acquisition of sex-determining genes on different autosomes, as well as translocations between different ancestral macro- and microchromosomes. We discuss the evolutionary drivers of the slow differentiation and turnover of reptile sex chromosomes

Sequence and gene content of a large fragment of a lizard sex chromosome and evaluation of candidate sex differentiating gene R-spondin 1

Background: Scant genomic information from non-avian reptile sex chromosomes is available, and for only a few lizards, several snakes and one turtle species, and it represents only a small fraction of the total sex chromosome sequences in these species. Results: We report a 352 kb of contiguous sequence from the sex chromosome of a squamate reptile, Pogona vitticeps, with a ZZ/ZW sex microchromosome system. This contig contains five protein coding genes (oprd1, rcc1, znf91, znf131, znf180), and major families of repetitive sequences with a high number of copies of LTR and non-LTR retrotransposons, including the CR1 and Bov-B LINEs. The two genes, oprd1 and rcc1 are part of a homologous syntenic block, which is conserved among amniotes. While oprd1 and rcc1 have no known function in sex determination or differentiation in amniotes, this homologous syntenic block in mammals and chicken also contains R-spondin 1 (rspo1), the ovarian differentiating gene in mammals. In order to explore the probability that rspo1 is sex determining in dragon lizards, genomic BAC and cDNA clones were mapped using fluorescence in situ hybridisation. Their location on an autosomal microchromosome pair, not on the ZW sex microchromosomes, eliminates rspo1 as a candidate sex determining gene in P. vitticeps. Conclusion: Our study has characterized the largest contiguous stretch of physically mapped sex chromosome sequence (352 kb) from a ZZ/ZW lizard species. Although this region represents only a small fraction of the sex chromosomes of P. vitticeps, it has revealed several features typically associated with sex chromosomes including the accumulation of large blocks of repetitive sequences

Chromosomal localization of the 18S-28S and 5S rRNA genes and (TTAGGG)n sequences of butterfly lizards (Leiolepis belliana belliana and Leiolepis boehmei, Agamidae, Squamata)

Chromosomal mapping of the butterfly lizards Leiolepis belliana belliana and L. boehmei was done using the 18S-28S and 5S rRNA genes and telomeric (TTAGGG)n sequences. The karyotype of L. b. belliana was 2n = 36, whereas that of L. boehmei was 2n = 34. The 18S-28S rRNA genes were located at the secondary constriction of the long arm of chromosome 1, while the 5S rRNA genes were found in the pericentromeric region of chromosome 6 in both species. Hybridization signals for the (TTAGGG)n sequence were observed at the telomeric ends of all chromosomes, as well as interstitially at the same position as the 18S-28S rRNA genes in L. boehmei. This finding suggests that in L. boehmei telomere-to-telomere fusion probably occurred between chromosome 1 and a microchromosome where the 18S-28S rRNA genes were located or, alternatively, at the secondary constriction of chromosome 1. The absence of telomeric sequence signals in chromosome 1 of L. b. belliana suggested that its chromosomes may have only a few copies of the (TTAGGG)n sequence or that there may have been a gradual loss of the repeat sequences during chromosomal evolution

Characterization of Squamate Olfactory Receptor Genes and Their Transcripts by the High-Throughput Sequencing Approach

The olfactory receptor (OR) genes represent the largest multigene family in the genome of terrestrial vertebrates. Here, the high-throughput next-generation sequencing (NGS) approach was applied to characterization of OR gene repertoires in the green anole lizard Anolis carolinensis and the Japanese four-lined ratsnake Elaphe quadrivirgata. Tagged polymerase chain reaction (PCR) products amplified from either genomic DNA or cDNA of the two species were used for parallel pyrosequencing, assembling, and screening for errors in PCR and pyrosequencing. Starting from the lizard genomic DNA, we accurately identified 56 of 136 OR genes that were identified from its draft genome sequence. These recovered genes were broadly distributed in the phylogenetic tree of vertebrate OR genes without severe biases toward particular OR families. Ninety-six OR genes were identified from the ratsnake genomic DNA, implying that the snake has more OR gene loci than the anole lizard in response to an increased need for the acuity of olfaction. This view is supported by the estimated number of OR genes in the Burmese python's draft genome (∼280), although squamates may generally have fewer OR genes than terrestrial mammals and amphibians. The OR gene repertoire of the python seems unique in that many class I OR genes are retained. The NGS approach also allowed us to identify candidates of highly expressed and silent OR gene copies in the lizard's olfactory epithelium. The approach will facilitate efficient and parallel characterization of considerable unbiased proportions of multigene family members and their transcripts from nonmodel organisms

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition