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.

Diversity and Evolution of Highly Repetitive DNA Sequences Constituting Chromosome Site-Specific Heterochromatin in Two Gerbillinae Species.

Constitutive heterochromatin, consisting of repetitive sequences, diverges very rapidly; therefore, its nucleotide sequences and chromosomal distributions are often largely different, even between closely related species. The chromosome C-banding patterns of two Gerbillinae species, Meriones unguiculatus and Gerbillus perpallidus, vary greatly, even though they belong to the same subfamily. To understand the evolution of C-positive heterochromatin in these species, we isolated highly repetitive sequences, determined their nucleotide sequences, and characterized them using chromosomal and filter hybridization. We obtained a centromeric repeat (MUN-HaeIII) and a chromosome 13-specific repeat (MUN-EcoRI) from M. unguiculatus. We also isolated a centromeric/pericentromeric repeat (GPE-MBD) and an interspersed-type repeat that was predominantly amplified in the X and Y chromosomes (GPE-EcoRI) from G. perpallidus. GPE-MBD was found to contain a 17-bp motif that is essential for binding to the centromere-associated protein CENP-B. This indicates that it may play a role in the formation of a specified structure and/or function of centromeres. The nucleotide sequences of the three sequence families, except GPE-EcoRI, were conserved only in Gerbillinae. GPE-EcoRI was derived from the long interspersed nuclear elements 1 retrotransposon and showed sequence homology throughout Muridae and Cricetidae species, indicating that the repeat sequence occurred at least in the common ancestor of Muridae and Cricetidae. Due to a lack of assembly data of highly repetitive sequences constituting heterochromatin in whole-genome sequences of vertebrate species published to date, the knowledge obtained in this study provides useful information for a deep understanding of the evolution of repetitive sequences in not only rodents but also in mammals.

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.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation.

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

The PRDM14-CtBP1/2-PRC2 complex regulates transcriptional repression during the transition from primed to naïve pluripotency.

The pluripotency-associated transcriptional network is regulated by a core circuitry of transcription factors. The PR domain-containing protein PRDM14 maintains pluripotency by activating and repressing transcription in a target gene-dependent manner. However, the mechanisms underlying dichotomic switching of PRDM14-mediated transcriptional control remain elusive. Here, we identified C-terminal binding protein 1 and 2 (CtBP1 and CtBP2; generically referred to as CtBP1/2) as components of the PRDM14-mediated repressive complex. CtBP1/2 binding to PRDM14 depends on CBFA2T2, a core component of the PRDM14 complex. The loss of Ctbp1/2 impaired the PRDM14-mediated transcriptional repression required for pluripotency maintenance and transition from primed to naïve pluripotency. Furthermore, CtBP1/2 interacted with the PRC2 complexes, and the loss of Ctbp1/2 impaired Polycomb repressive complex 2 (PRC2) and H3K27me3 enrichment at target genes after Prdm14 induction. These results provide evidence that the target gene-dependent transcriptional activity of PRDM14 is regulated by partner switching to ensure the transition from primed to naïve pluripotency.This article has an associated First Person interview with the first author of the paper.

ZW Sex Chromosomes in Australian Dragon Lizards (Agamidae) Originated from a Combination of Duplication and Translocation in the Nucleolar Organising Region.

Sex chromosomes in some reptiles share synteny with distantly related amniotes in regions orthologous to squamate chromosome 2. The latter finding suggests that chromosome 2 was formerly part of a larger ancestral (amniote) super-sex chromosome and raises questions about how sex chromosomes are formed and modified in reptiles. Australian dragon lizards (Agamidae) are emerging as an excellent model for studying these processes. In particular, they exhibit both genotypic (GSD) and temperature-dependent (TSD) sex determination, show evidence of transitions between the two modes and have evolved non-homologous ZW sex microchromosomes even within the same evolutionary lineage. They therefore represent an excellent group to probe further the idea of a shared ancestral super-sex chromosome and to investigate mechanisms for transition between different sex chromosome forms. Here, we compare sex chromosome homology among eight dragon lizard species from five genera to identify key cytological differences and the mechanisms that may be driving sex chromosome evolution in this group. We performed fluorescence in situ hybridisation to physically map bacterial artificial chromosome (BAC) clones from the bearded dragon, Pogona vitticeps' ZW sex chromosomes and a nucleolar organising region (NOR) probe in males and females of eight Agamid species exhibiting either GSD or TSD. We show that the sex chromosome derived BAC clone hybridises near the telomere of chromosome 2q in all eight species examined. This clone also hybridises to the sex microchromosomes of three species (P vitticeps, P. barbata and Diporiphora nobbi) and a pair of microchromosomes in three others (Ctenophorus pictus, Amphibolurus norrisi and Amphibolurus muricatus). No other chromosomes are marked by the probe in two species from the closely related genus Physignathus. A probe bearing nucleolar organising region (NOR) sequences maps close to the telomere of chromosome 2q in all eight species, and to the ZW pair in P. vitticeps and P. barbata, the W microchromosome in D. nobbi, and several microchromosomes in P. cocincinus. Our findings provide evidence of sequence homology between chromosome 2 and the sex chromosomes of multiple agamids. These data support the hypothesis that there was an ancestral sex chromosome in amniotes that gave rise to squamate chromosome 2 and raises the prospect that some particular property of this chromosome has favoured its role as a sex chromosome in amniotes. It is likely that the amplification of repetitive sequences associated with this region has driven the high level of heterochromatinisation of the sex-specific chromosomes in three species of agamid. Our data suggest a possible mechanism for chromosome rearrangement, including inversion and duplication near the telomeric regions of the ancestral chromosome 2 and subsequent translocation to the ZW sex microchromosomes in three agamid species. It is plausible that these chromosome rearrangements involving sex chromosomes also drove speciation in this group.

Karyotype Analysis of Four Blind Snake Species (Reptilia: Squamata: Scolecophidia) and Karyotypic Changes in Serpentes.

The suborder Serpentes is divided into 2 infraorders, Scolecophidia and Alethinophidia, which diverged at an early stage of snake diversification. In this study, we examined karyotypes of 4 scolecophidian species (Letheobia simonii, Xerotyphlops vermicularis, Indotyphlops braminus, and Myriopholis macrorhyncha) and performed FISH with 18S-28S rDNA as well as microchromosomal and Z chromosome-linked genes of Elaphe quadrivirgata (Alethinophidia) to investigate the karyotype evolution in the scolecophidian lineage. Diploid chromosome numbers of X. vermicularis and L. simonii were 30 (16 macrochromosomes and 14 microchromosomes) and 32 (16 macrochromosomes and 16 microchromosomes), respectively. The karyotype of a female M. macrorhyncha consisted of 15 macrochromosomes and 19 microchromosomes, including a heterochromatic microchromosome, indicating the presence of a heteromorphic chromosome pair. E. quadrivirgata Z-linked genes mapped to chromosome 4 of M. macrorhyncha, not to the heteromorphic pair. Therefore, M. macrorhyncha may have differentiated ZW sex chromosomes which are not homologous to those of E. quadrivirgata. One of the E. quadrivirgata microchromosomal genes mapped to the terminal region of chromosome 4q in X. vermicularis, suggesting that fusions between microchromosomes and macrochromosomes occurred in this species. rDNA was localized in different macrochromosomal pairs in the 2 diploid scolecophidian snakes examined here, whereas the gene location in a microchromosomal pair was conserved in 5 alethinophidian species examined. These results might imply the occurrence of chromosome fusions in the scolecophidian lineages. In I. braminus, a unique parthenogenetic snake with a triploid karyotype (21 macrochromosomes and 21 microchromosomes), morphological heteromorphisms were identified in chromosomes 1 and 7. Such heteromorphisms in 2 chromosomes were also observed in individuals from distant locations in the broad distribution range of this species, suggesting that the heteromorphisms were fixed in the genome at an early stage of its speciation.

Co-option of the PRDM14-CBFA2T complex from motor neurons to pluripotent cells during vertebrate evolution.

Gene regulatory networks underlying cellular pluripotency are controlled by a core circuitry of transcription factors in mammals, including POU5F1. However, the evolutionary origin and transformation of pluripotency-related transcriptional networks have not been elucidated in deuterostomes. PR domain-containing protein 14 (PRDM14) is specifically expressed in pluripotent cells and germ cells, and is required for establishing embryonic stem cells (ESCs) and primordial germ cells in mice. Here, we compared the functions and expression patterns of PRDM14 orthologues within deuterostomes. Amphioxus PRDM14 and zebrafish PRDM14, but not sea urchin PRDM14, compensated for mouse PRDM14 function in maintaining mouse ESC pluripotency. Interestingly, sea urchin PRDM14 together with sea urchin CBFA2T, an essential partner of PRDM14 in mouse ESCs, complemented the self-renewal defect in mouse Prdm14 KO ESCs. Contrary to the Prdm14 expression pattern in mouse embryos, Prdm14 was expressed in motor neurons of amphioxus embryos, as observed in zebrafish embryos. Thus, Prdm14 expression in motor neurons was conserved in non-tetrapod deuterostomes and the co-option of the PRDM14-CBFA2T complex from motor neurons into pluripotent cells may have maintained the transcriptional network for pluripotency during vertebrate evolution.This article has an associated 'The people behind the papers' interview.

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.

Identification of genetic linkage group 1-linked sequences in Japanese eel (Anguilla japonica) by single chromosome sorting and sequencing.

Japanese eel (Anguilla japonica) constitutes one of the most important food fish in Japan; accordingly, genome sequencing and linkage mapping have been conducted for the purpose of artificial cultivation. In the next stage, integration of genomic sequences within linkage groups (LG) is required to construct higher-resolution genetic markers for quantitative trait loci mapping and selective breeding of beneficial traits in farming. In order to identify LG1-linked scaffolds from the draft genome assembly (323,776 scaffolds) reported previously, we attempted to isolate chromosomes corresponding to LG1 by flow sorting and subsequent analyses. Initially, single chromosomes were randomly collected by chromosome sorting and subjected to whole-genome amplification (WGA). A total of 60 WGA samples were screened by PCR with primers for a known LG1-linked scaffold, and five positive WGA samples were sequenced by next-generation sequencing (NGS). Following reference mapping analysis of the NGS reads, four of the five WGA samples were found to be enriched by LG1-linked sequences. These samples were cytogenetically assigned to chromosome 5 by fluorescence in situ hybridization. Using blastn searches with 82,081 contigs constructed from the NGS reads of the four WGA samples as queries, 2323 scaffolds were identified as putative LG1-linked scaffolds from the draft genome assembly. The total length of the putative LG1-linked scaffolds was 99.0 Mb, comparable to the estimated DNA amounts of chromosome 5 (91.1 Mb). These results suggest that the methodology developed herein is applicable to isolate specific chromosome DNAs and integrate unanchored scaffold sequences onto a particular LG and chromosome even in teleost fishes, in which isolation of specific chromosomes by flow sorting is generally difficult owing to similar morphologies, sizes, and GC-contents among chromosomes in the genome. The putative LG1-linked scaffolds of Japanese eel contain a total of 6833 short tandem repeats which will be available for higher-resolution linkage mapping.

Lack of satellite DNA species-specific homogenization and relationship to chromosomal rearrangements in monitor lizards (Varanidae, Squamata).

BACKGROUND: Satellite DNAs (stDNAs) are highly repeated sequences that constitute large portions of any genome. The evolutionary dynamics of stDNA (e.g. copy number, nucleotide sequence, location) can, therefore, provide an insight into genome organization and evolution. We investigated the evolutionary origin of VSAREP stDNA in 17 monitor lizards (seven Asian, five Australian, and five African) at molecular and cytogenetic level. RESULTS: Results revealed that VSAREP is conserved in the genome of Asian and Australian varanids, but not in African varanids, suggesting that these sequences are either differentiated or lost in the African varanids. Phylogenetic and arrangement network analyses revealed the existence of at least four VSAREP subfamilies. The similarity of each sequence unit within the same VSAREP subfamily from different species was higher than those of other VSAREP subfamilies belonging to the same species. Additionally, all VSAREP subfamilies isolated from the three Australian species (Varanus rosenbergi, V. gouldii, and V. acanthurus) were co-localized near the centromeric or pericentromeric regions of the macrochromosomes, except for chromosomes 3 and 4 in each Australian varanid. However, their chromosomal arrangements were different among species. CONCLUSIONS: The VSAREP stDNA family lack homogenized species-specific nucleotide positions in varanid lineage. Most VSAREP sequences were shared among varanids within the four VSAREP subfamilies. This suggests that nucleotide substitutions in each varanid species accumulated more slowly than homogenization rates in each VSAREP subfamily, resulting in non-species-specific evolution of stDNA profiles. Moreover, changes in location of VSAREP stDNA in each Australian varanid suggests a correlation with chromosomal rearrangements, leading to karyotypic differences among these species.

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.

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

BACKGROUND: The discovery of differentially organized sex chromosome systems suggests that heteromorphic sex chromosomes evolved from a pair of homologous chromosomes. Whereas karyotypes are highly conserved in alethinophidian snakes, the degeneration status of the W chromosomes varies among species. The Z and W chromosomes are morphologically homomorphic in henophidian species, whereas in snakes belonging to caenophidian families the W chromosomes are highly degenerated. Snakes therefore are excellent animal models in which to study sex chromosome evolution. Herein, we investigated the differentiation processes for snake sex chromosomes using both coding and repetitive sequences. We analyzed phylogenetic relationships of CTNNB1 and WAC genes, localized to the centromeric and telomeric regions, respectively, of the long arms on snake sex chromosomes, and chromosome distribution of sex chromosome-linked repetitive sequences in several henophidian and caenophidian species. RESULTS: Partial or full-length coding sequences of CTNNB1 and WAC were identified for Z homologs of henophidian species from Tropidophiidae, Boidae, Cylindrophiidae, Xenopeltidae, and Pythonidae, and for Z and W homologs of caenophidian species from Acrochordidae, Viperidae, Elapidae, and Colubridae. Female-specific sequences for the two genes were not found in the henophidian (boid and pythonid) species examined. Phylogenetic trees constructed using each gene showed that the Z and W homologs of the caenophidian species cluster separately. The repetitive sequence isolated from the W chromosome heterochromatin of the colubrid Elaphe quadrivirgata and a microsatellite motif (AGAT)8 were strongly hybridized with W chromosomes of the viperid and colubrid species examined. CONCLUSION: Our phylogenetic analyses suggest that the cessation of recombination between the Z and W homologs of CTNNB1 and WAC predated the diversification of the caenophidian families. As the repetitive sequences on the W chromosomes were shared among viperid and colubrid species, heterochromatinization of the proto-W chromosome appears to have occurred before the splitting of these two groups. These results collectively suggest that differentiation of the proto-Z and proto-W chromosomes extended to wide regions on the sex chromosomes in the common ancestor of caenophidian families during a relatively short period.

Comparison of the primary structures, cytotoxicities, and affinities to phospholipids of five kinds of cytotoxins from the venom of Indian cobra, Naja naja.

The molecular mechanism underlying the hemolytic and cytolytic processes of cobra cytotoxins (CTXs) is not yet fully elucidated. To examine this, we analyzed the amino acid sequences, hemolytic and cytotoxic activities, and affinities to phospholipids of the five major CTXs purified from the venom of Indian cobra, Naja naja. CTX2, CTX7, and CTX8 belonged to S-type, and CTX9 and CTX10 to P-type. Comparisons of CTX7 with CTX8 and CTX9 with CTX10 revealed similar primary structures and hemolytic and cytolytic activities. CTX2, whose primary structure was rather different from the others, showed several times weaker hemolytic and cytolytic biological activities than the others. The comparison of CTX2 with CTX7 suggested the importance of Lys30 in loop II for the strong hemolytic and cytolytic activities of S-type CTXs. Cloning of 12 CTX cDNAs from the Naja naja venom cDNA library revealed that 18 out of 23 substitutions found in CTX cDNAs were nonsynonymous. This clearly indicated the accelerated evolution of CTX genes. Multiple sequence alignment of 51 kinds of CTX cDNAs and calculations of nonsynonymous and synonymous substitutions indicated that the codons coding the three loops' regions, which may interact with the hydrophobic tails of phospholipids, have undergone an accelerated evolution. In contrast, the codons coding for amino acid residues considered to participate in the recognition and binding of the hydrophilic head groups of phospholipids, eight Cys residues, and those likely stabilizing ß core structure, were all conserved.

Amplification of microsatellite repeat motifs is associated with the evolutionary differentiation and heterochromatinization of sex chromosomes in Sauropsida.

The sex chromosomes in Sauropsida (reptiles and birds) have evolved independently many times. They show astonishing diversity in morphology ranging from cryptic to highly differentiated sex chromosomes with male (XX/XY) and female heterogamety (ZZ/ZW). Comparing such diverse sex chromosome systems thus provides unparalleled opportunities to capture evolution of morphologically differentiated sex chromosomes in action. Here, we describe chromosomal mapping of 18 microsatellite repeat motifs in eight species of Sauropsida. More than two microsatellite repeat motifs were amplified on the sex-specific chromosome, W or Y, in five species (Bassiana duperreyi, Aprasia parapulchella, Notechis scutatus, Chelodina longicollis, and Gallus gallus) of which the sex-specific chromosomes were heteromorphic and heterochromatic. Motifs (AAGG)n and (ATCC)n were amplified on the W chromosome of Pogona vitticeps and the Y chromosome of Emydura macquarii, respectively. By contrast, no motifs were amplified on the W chromosome of Christinus marmoratus, which is not much differentiated from the Z chromosome. Taken together with previously published studies, our results suggest that the amplification of microsatellite repeats is tightly associated with the differentiation and heterochromatinization of sex-specific chromosomes in sauropsids as well as in other taxa. Although some motifs were common between the sex-specific chromosomes of multiple species, no correlation was observed between this commonality and the species phylogeny. Furthermore, comparative analysis of sex chromosome homology and chromosomal distribution of microsatellite repeats between two closely related chelid turtles, C. longicollis and E. macquarii, identified different ancestry and differentiation history. These suggest multiple evolutions of sex chromosomes in the Sauropsida.

No Interstitial Telomeres on Autosomes but Remarkable Amplification of Telomeric Repeats on the W Sex Chromosome in the Sand Lizard (Lacerta agilis).

Telomeres are repeat (TTAGGG) n sequences that form terminal ends of chromosomes and have several functions, such as protecting the coding DNA from erosion at mitosis. Due to chromosomal rearrangements through evolutionary history (e.g., inversions and fusions), telomeric sequences are also found between the centromere and the terminal ends (i.e., at interstitial telomeric sites, ITSs). ITS telomere sequences have been implicated in heritable disease caused by genomic instability of ITS polymorphic variants, both with respect to copy number and sequence. In the sand lizard (Lacerta agilis), we have shown that telomere length is predictive of lifetime fitness in females but not males. To assess whether this sex specific fitness effect could be traced to ITSs differences, we mapped (TTAGGG) n sequences using fluorescence in situ hybridization in fibroblast cells cultured from 4 specimens of known sex. No ITSs could be found on autosomes in either sex. However, females have heterogametic sex chromosomes in sand lizards (ZW, 2n = 38) and the female W chromosome showed degeneration and remarkable (TTAGGG) n amplification, which was absent in the Z chromosomes. This work warrants further research on sex chromosome content, in particular of the degenerate W chromosome, and links to female fitness in sand lizards.

High-coverage sequencing and annotated assembly of the genome of the Australian dragon lizard Pogona vitticeps.

BACKGROUND: The lizards of the family Agamidae are one of the most prominent elements of the Australian reptile fauna. Here, we present a genomic resource built on the basis of a wild-caught male ZZ central bearded dragon Pogona vitticeps. FINDINGS: The genomic sequence for P. vitticeps, generated on the Illumina HiSeq 2000 platform, comprised 317 Gbp (179X raw read depth) from 13 insert libraries ranging from 250 bp to 40 kbp. After filtering for low-quality and duplicated reads, 146 Gbp of data (83X) was available for assembly. Exceptionally high levels of heterozygosity (0.85 % of single nucleotide polymorphisms plus sequence insertions or deletions) complicated assembly; nevertheless, 96.4 % of reads mapped back to the assembled scaffolds, indicating that the assembly included most of the sequenced genome. Length of the assembly was 1.8 Gbp in 545,310 scaffolds (69,852 longer than 300 bp), the longest being 14.68 Mbp. N50 was 2.29 Mbp. Genes were annotated on the basis of de novo prediction, similarity to the green anole Anolis carolinensis, Gallus gallus and Homo sapiens proteins, and P. vitticeps transcriptome sequence assemblies, to yield 19,406 protein-coding genes in the assembly, 63 % of which had intact open reading frames. Our assembly captured 99 % (246 of 248) of core CEGMA genes, with 93 % (231) being complete. CONCLUSIONS: The quality of the P. vitticeps assembly is comparable or superior to that of other published squamate genomes, and the annotated P. vitticeps genome can be accessed through a genome browser available at https://genomics.canberra.edu.au.

Sex reversal triggers the rapid transition from genetic to temperature-dependent sex.

Sex determination in animals is amazingly plastic. Vertebrates display contrasting strategies ranging from complete genetic control of sex (genotypic sex determination) to environmentally determined sex (for example, temperature-dependent sex determination). Phylogenetic analyses suggest frequent evolutionary transitions between genotypic and temperature-dependent sex determination in environmentally sensitive lineages, including reptiles. These transitions are thought to involve a genotypic system becoming sensitive to temperature, with sex determined by gene-environment interactions. Most mechanistic models of transitions invoke a role for sex reversal. Sex reversal has not yet been demonstrated in nature for any amniote, although it occurs in fish and rarely in amphibians. Here we make the first report of reptile sex reversal in the wild, in the Australian bearded dragon (Pogona vitticeps), and use sex-reversed animals to experimentally induce a rapid transition from genotypic to temperature-dependent sex determination. Controlled mating of normal males to sex-reversed females produces viable and fertile offspring whose phenotypic sex is determined solely by temperature (temperature-dependent sex determination). The W sex chromosome is eliminated from this lineage in the first generation. The instantaneous creation of a lineage of ZZ temperature-sensitive animals reveals a novel, climate-induced pathway for the rapid transition between genetic and temperature-dependent sex determination, and adds to concern about adaptation to rapid global climate change.

Molecular cloning and characterization of satellite DNA sequences from constitutive heterochromatin of the habu snake (Protobothrops flavoviridis, Viperidae) and the Burmese python (Python bivittatus, Pythonidae).

Highly repetitive DNA sequences of the centromeric heterochromatin provide valuable molecular cytogenetic markers for the investigation of genomic compartmentalization in the macrochromosomes and microchromosomes of sauropsids. Here, the relationship between centromeric heterochromatin and karyotype evolution was examined using cloned repetitive DNA sequences from two snake species, the habu snake (Protobothrops flavoviridis, Crotalinae, Viperidae) and Burmese python (Python bivittatus, Pythonidae). Three satellite DNA (stDNA) families were isolated from the heterochromatin of these snakes: 168-bp PFL-MspI from P. flavoviridis and 196-bp PBI-DdeI and 174-bp PBI-MspI from P. bivittatus. The PFL-MspI and PBI-DdeI sequences were localized to the centromeric regions of most chromosomes in the respective species, suggesting that the two sequences were the major components of the centromeric heterochromatin in these organisms. The PBI-MspI sequence was localized to the pericentromeric region of four chromosome pairs. The PFL-MspI and the PBI-DdeI sequences were conserved only in the genome of closely related species, Gloydius blomhoffii (Crotalinae) and Python molurus, respectively, although their locations on the chromosomes were slightly different. In contrast, the PBI-MspI sequence was also in the genomes of P. molurus and Boa constrictor (Boidae), and additionally localized to the centromeric regions of eight chromosome pairs in B. constrictor, suggesting that this sequence originated in the genome of a common ancestor of Pythonidae and Boidae, approximately 86 million years ago. The three stDNA sequences showed no genomic compartmentalization between the macrochromosomes and microchromosomes, suggesting that homogenization of the centromeric and/or pericentromeric stDNA sequences occurred in the macrochromosomes and microchromosomes of these snakes.

Non-homologous sex chromosomes in two geckos (Gekkonidae: Gekkota) with female heterogamety.

Evaluating homology between the sex chromosomes of different species is an important first step in deducing the origins and evolution of sex-determining mechanisms in a clade. Here, we describe the preparation of Z and W chromosome paints via chromosome microdissection from the Australian marbled gecko (Christinus marmoratus) and their subsequent use in evaluating sex chromosome homology with the ZW chromosomes of the Kwangsi gecko (Gekko hokouensis) from eastern Asia. We show that the ZW sex chromosomes of C. marmoratus and G. hokouensis are not homologous and represent independent origins of female heterogamety within the Gekkonidae. We also show that the C. marmoratus Z and W chromosomes are genetically similar to each other as revealed by C-banding, comparative genomic hybridization, and the reciprocal painting of Z and W chromosome probes. This implies that sex chromosomes in C. marmoratus are at an early stage of differentiation, suggesting a recent origin.