The ZW sex microchromosomes of an Australian dragon lizard share no homology with those of other reptiles or birds

Evolution of sex determination in crustaceans

Sex determination (SD) involves mechanisms that determine whether an individual will develop into a male, female, or in rare cases, hermaphrodite. Crustaceans harbor extremely diverse SD systems, including hermaphroditism, environmental sex determination (ESD), genetic sex determination (GSD), and cytoplasmic sex determination (e.g., Wolbachia controlled SD systems). Such diversity lays the groundwork for researching the evolution of SD in crustaceans, i.e., transitions among different SD systems. However, most previous research has focused on understanding the mechanism of SD within a single lineage or species, overlooking the transition across different SD systems. To help bridge this gap, we summarize the understanding of SD in various clades of crustaceans, and discuss how different SD systems might evolve from one another. Furthermore, we review the genetic basis for transitions between different SD systems (i.e., Dmrt genes) and propose the microcrustacean Daphnia (clade Branchiopoda) as a model to study the transition from ESD to GSD.

Models and Molecular Markers of Spermatogonial Stem Cells in Vertebrates: To Find Models in Nonmammals

Spermatogonial stem cells (SSCs) are the germline stem cells that are essential for the maintenance of spermatogenesis in the testis. However, it has not been sufficiently understood in amphibians, reptiles, and fish because numerous studies have been focused mainly on mammals. The aim of this review is to discuss scientific ways to elucidate SSC models of nonmammals in the context of the evolution of testicular organization since rodent SSC models. To further understand the SSC models in nonmammals, we point out common markers of an SSC pool (undifferentiated spermatogonia) in various types of testes where the kinetics of the SSC pool appears. This review includes the knowledge of (1) common molecular markers of vertebrate type A spermatogonia including putative SSC markers, (2) localization of the markers on the spermatogonia that have been reported in previous studies, (3) highlighting the most common markers in vertebrates, and (4) suggesting ways of finding SSC models in nonmammals.

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.

Karyotype Characterisation of Two Australian Dragon Lizards (Squamata: Agamidae: Amphibolurinae) Reveals Subtle Chromosomal Rearrangements Between Related Species with Similar Karyotypes

Agamid lizards (Squamata: Agamidae) are karyotypically heterogeneous. Among the 101 species currently described from Australia, all are from the subfamily Amphibolurinae. This group is, with some exceptions, karyotypically conserved, and all species involving heterogametic sex show female heterogamety. Here, we describe the chromosomes of 2 additional Australian agamid lizards, Tympanocryptis lineata and Rankinia diemensis. These species are phylogenetically and cytogenetically sisters to the well-characterised Pogona vitticeps, but their sex chromosomes and other chromosomal characteristics are unknown. In this study, we applied advanced molecular cytogenetic techniques, such as fluorescence in situ hybridisation (FISH) and cross-species gene mapping, to characterise chromosomes and to identify sex chromosomes in these species. Our data suggest that both species have a conserved karyotype with P. vitticeps but with subtle rearrangements in the chromosomal landscapes. We could identify that T. lineata possesses a female heterogametic system (ZZ/ZW) with a pair of sex microchromosomes, while R. diemensis may have heterogametic sex chromosomes, but this requires further investigations. Our study shows the pattern of chromosomal rearrangements between closely related species, explaining the speciation within Australian agamid lizards of similar karyotypes.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction

Background

Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes.

Results

We generated 96.7 × 10⁹ 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 10⁹ paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females.

Conclusions

The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

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.

Cross-Species BAC Mapping Highlights Conservation of Chromosome Synteny across Dragon Lizards (Squamata: Agamidae)

Dragon lizards (Squamata: Agamidae) comprise about 520 species in six subfamilies distributed across Asia, Australasia and Africa. Only five species are known to have sex chromosomes. All of them possess ZZ/ZW sex chromosomes, which are microchromosomes in four species from the subfamily Amphibolurinae, but much larger in Phrynocephalus vlangalii from the subfamily Agaminae. In most previous studies of these sex chromosomes, the focus has been on Australian species from the subfamily Amphibolurinae, but only the sex chromosomes of the Australian central bearded dragon (Pogona vitticeps) are well-characterized cytogenetically. To determine the level of synteny of the sex chromosomes of P. vitticeps across agamid subfamilies, we performed cross-species two-colour FISH using two bacterial artificial chromosome (BAC) clones from the pseudo-autosomal regions of P. vitticeps. We mapped these two BACs across representative species from all six subfamilies as well as two species of chameleons, the sister group to agamids. We found that one of these BAC sequences is conserved in macrochromosomes and the other in microchromosomes across the agamid lineages. However, within the Amphibolurinae, there is evidence of multiple chromosomal rearrangements with one of the BACs mapping to the second-largest chromosome pair and to the microchromosomes in multiple species including the sex chromosomes of P. vitticeps. Intriguingly, no hybridization signal was observed in chameleons for either of these BACs, suggesting a likely agamid origin of these sequences. Our study shows lineage-specific evolution of sequences/syntenic blocks and successive rearrangements and reveals a complex history of sequences leading to their association with important biological processes such as the evolution of sex chromosomes and sex determination.

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.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

Weird Animals, Sex, and Genome Evolution

Making my career in Australia exposed me to the tyranny of distance, but it gave me opportunities to study our unique native fauna. Distantly related animal species present genetic variation that we can use to explore the most fundamental biological structures and processes. I have compared chromosomes and genomes of kangaroos and platypus, tiger snakes and emus, devils (Tasmanian) and dragons (lizards). I particularly love the challenges posed by sex chromosomes, which, apart from determining sex, provide stunning examples of epigenetic control and break all the evolutionary rules that we currently understand. Here I describe some of those amazing animals and the insights on genome structure, function, and evolution they have afforded us. I also describe my sometimes-random walk in science and the factors and people who influenced my direction. Being a woman in science is still not easy, and I hope others will find encouragement and empathy in my story.

Expression of a hindlimb-determining factor Pitx1 in the forelimb of the lizard Pogona vitticeps during morphogenesis

With over 9000 species, squamates, which include lizards and snakes, are the largest group of reptiles and second-largest order of vertebrates, spanning a vast array of appendicular skeletal morphology. As such, they provide a promising system for examining developmental and molecular processes underlying limb morphology. Using the central bearded dragon (Pogona vitticeps) as the primary study model, we examined limb morphometry throughout embryonic development and characterized the expression of three known developmental genes (GHR, Pitx1 and Shh) from early embryonic stage through to hatchling stage via reverse transcription quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC). In this study, all genes were found to be transcribed in both the forelimbs and hindlimbs of P. vitticeps. While the highest level of GHR expression occurred at the hatchling stage, Pitx1 and Shh expression was greatest earlier during embryogenesis, which coincides with the onset of the differentiation between forelimb and hindlimb length. We compared our finding of Pitx1 expression—a hindlimb-determining gene—in the forelimbs of P. vitticeps to that in a closely related Australian agamid lizard, Ctenophorus pictus, where we found Pitx1 expression to be more highly expressed in the hindlimb compared with the forelimb during early and late morphogenesis—a result consistent with that found across other tetrapods. Expression of Pitx1 in forelimbs has only rarely been documented, including via in situ hybridization in a chicken and a frog. Our findings from both RT-qPCR and IHC indicate that further research across a wider range of tetrapods is needed to more fully understand evolutionary variation in molecular processes underlying limb morphology.

Sex determination mode does not affect body or genital development of the central bearded dragon (Pogona vitticeps)

Background

The development of male- or female-specific phenotypes in squamates is typically controlled by either temperature-dependent sex determination (TSD) or chromosome-based genetic sex determination (GSD). However, while sex determination is a major switch in individual phenotypic development, it is unknownhow evolutionary transitions between GSD and TSD might impact on the evolution of squamate phenotypes, particularly the fast-evolving and diverse genitalia. Here, we take the unique opportunity of studying the impact of both sex determination mechanisms on the embryological development of the central bearded dragon (Pogona vitticeps). This is possible because of the transitional sex determination system of this species, in which genetically male individuals reverse sex at high incubation temperatures. This can trigger the evolutionary transition of GSD to TSD in a single generation, making P. vitticeps an ideal model organism for comparing the effects of both sex determination processes in the same species.

Results

We conducted four incubation experiments on 265 P. vitticeps eggs, covering two temperature regimes ("normal" at 28 °C and "sex reversing" at 36 °C) and the two maternal sexual genotypes (concordant ZW females or sex-reversed ZZ females). From this, we provide the first detailed staging system for the species, with a focus on genital and limb development. This was augmented by a new sex chromosome identification methodology for P. vitticeps that is non-destructive to the embryo. We found a strong correlation between embryo age and embryo stage. Aside from faster growth in 36 °C treatments, body and external genital development was entirely unperturbed by temperature, sex reversal or maternal sexual genotype. Unexpectedly, all females developed hemipenes (the genital phenotype of adult male P. vitticeps), which regress close to hatching.

Conclusions

The tight correlation between embryo age and embryo stage allows the precise targeting of specific developmental periods in the emerging field of molecular research on P. vitticeps. The stability of genital development in all treatments suggests that the two sex-determining mechanisms have little impact on genital evolution, despite their known role in triggering genital development. Hemipenis retention in developing female P. vitticeps, together with frequent occurrences of hemipenis-like structures during development in other squamate species, raises the possibility of a bias towards hemipenis formation in the ancestral developmental programme for squamate genitalia.

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno's original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster. Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways.

Highly conserved Z and molecularly diverged W chromosomes in the fish genus Triportheus (Characiformes, Triportheidae)

The main objectives of this study were to test: (1) whether the W-chromosome differentiation matches to species' evolutionary divergence (phylogenetic concordance) and (2) whether sex chromosomes share a common ancestor within a congeneric group. The monophyletic genus Triportheus (Characiformes, Triportheidae) was the model group for this study. All species in this genus so far analyzed have ZW sex chromosome system, where the Z is always the largest chromosome of the karyotype, whereas the W chromosome is highly variable ranging from almost homomorphic to highly heteromorphic. We applied conventional and molecular cytogenetic approaches including C-banding, ribosomal DNA mapping, comparative genomic hybridization (CGH) and cross-species whole chromosome painting (WCP) to test our questions. We developed Z- and W-chromosome paints from T. auritus for cross-species WCP and performed CGH in a representative species (T. signatus) to decipher level of homologies and rates of differentiation of W chromosomes. Our study revealed that the ZW sex chromosome system had a common origin, showing highly conserved Z chromosomes and remarkably divergent W chromosomes. Notably, the W chromosomes have evolved to different shapes and sequence contents within ~15-25 Myr of divergence time. Such differentiation highlights a dynamic process of W-chromosome evolution within congeneric species of Triportheus.

Cytogenetic Insights into the Evolution of Chromosomes and Sex Determination Reveal Striking Homology of Turtle Sex Chromosomes to Amphibian Autosomes

Turtle karyotypes are highly conserved compared to other vertebrates; yet, variation in diploid number (2n = 26-68) reflects profound genomic reorganization, which correlates with evolutionary turnovers in sex determination. We evaluate the published literature and newly collected comparative cytogenetic data (G- and C-banding, 18S-NOR, and telomere-FISH mapping) from 13 species spanning 2n = 28-68 to revisit turtle genome evolution and sex determination. Interstitial telomeric sites were detected in multiple lineages that underwent diploid number and sex determination turnovers, suggesting chromosomal rearrangements. C-banding revealed potential interspecific variation in centromere composition and interstitial heterochromatin at secondary constrictions. 18S-NORs were detected in secondary constrictions in a single chromosomal pair per species, refuting previous reports of multiple NORs in turtles. 18S-NORs are linked to ZW chromosomes in Apalone and Pelodiscus and to X (not Y) in Staurotypus. Notably, comparative genomics across amniotes revealed that the sex chromosomes of several turtles, as well as mammals and some lizards, are homologous to components of Xenopus tropicalis XTR1 (carrying Dmrt1). Other turtle sex chromosomes are homologous to XTR4 (carrying Wt1). Interestingly, all known turtle sex chromosomes, except in Trionychidae, evolved via inversions around Dmrt1 or Wt1. Thus, XTR1 appears to represent an amniote proto-sex chromosome (perhaps linked ancestrally to XTR4) that gave rise to turtle and other amniote sex chromosomes.

Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome

Background

Squamates (lizards and snakes) are a speciose lineage of reptiles displaying considerable karyotypic diversity, particularly among lizards. Understanding the evolution of this diversity requires comparison of genome organisation between species. Although the genomes of several squamate species have now been sequenced, only the green anole lizard has any sequence anchored to chromosomes. There is only limited gene mapping data available for five other squamates. This makes it difficult to reconstruct the events that have led to extant squamate karyotypic diversity. The purpose of this study was to anchor the recently sequenced central bearded dragon (Pogona vitticeps) genome to chromosomes to trace the evolution of squamate chromosomes. Assigning sequence to sex chromosomes was of particular interest for identifying candidate sex determining genes.

Results

By using two different approaches to map conserved blocks of genes, we were able to anchor approximately 42 % of the dragon genome sequence to chromosomes. We constructed detailed comparative maps between dragon, anole and chicken genomes, and where possible, made broader comparisons across Squamata using cytogenetic mapping information for five other species. We show that squamate macrochromosomes are relatively well conserved between species, supporting findings from previous molecular cytogenetic studies. Macrochromosome diversity between members of the Toxicofera clade has been generated by intrachromosomal, and a small number of interchromosomal, rearrangements. We reconstructed the ancestral squamate macrochromosomes by drawing upon comparative cytogenetic mapping data from seven squamate species and propose the events leading to the arrangements observed in representative species. In addition, we assigned over 8 Mbp of sequence containing 219 genes to the Z chromosome, providing a list of genes to begin testing as candidate sex determining genes.

Conclusions

Anchoring of the dragon genome has provided substantial insight into the evolution of squamate genomes, enabling us to reconstruct ancestral macrochromosome arrangements at key positions in the squamate phylogeny, demonstrating that fusions between macrochromosomes or fusions of macrochromosomes and microchromosomes, have played an important role during the evolution of squamate genomes. Assigning sequence to the sex chromosomes has identified NR5A1 as a promising candidate sex determining gene in the dragon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2774-3) contains supplementary material, which is available to authorized users.

How Australian mammals contributed to our understanding of sex determination and sex chromosomes

Marsupials and monotremes can be thought of as independent experiments in mammalian evolution. The discovery of the human male-determining gene, SRY, how it works, how it evolved and defined our sex chromosomes, well illustrates the value of comparing distantly related animals and the folly of relying on humans and mice for an understanding of the most fundamental aspects of mammalian biology. The 25th anniversary of the discovery of SRY seems a good time to review the contributions of Australian mammals to these discoveries. The discovery of the mammalian sex determining gene, SRY, was a milestone in the history of human genetics. SRY opened up investigations into the pathway by which the genital ridge (bipotential gonad) becomes a testis. Studies of Australian mammals were important in the story of the discovery of SRY, not only in refuting the qualifications of the first candidate sex-determining gene, but also in confirming the ubiquity of SRY and raising questions as to how it works. Studies in marsupials also led to understanding of how SRY evolved from a gene on an autosome with functions in the brain and germ cells, and to identifying the ancestors of other genes on the human Y. The discovery that platypus have sex chromosomes homologous, not to the human XY, but to the bird ZW, dated the origin of the therian SRY and the XY chromosomes it defined. This led to important new models of how our sex chromosomes function, how they evolved, and what might befall this gene and the Y chromosome it defines.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0085-2) contains supplementary material, which is available to authorized users.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

Nascent ZW sex chromosomes in Thecadactylus rapicauda (Reptilia, Squamata, Phyllodactylidae)

The chromosomes of the turnip-tailed gecko Thecadactylus rapicauda from the Falcón State in northern Venezuela were examined by means of conventional staining, a variety of banding techniques and in situ hybridization with an 18S + 28S rDNA probe. In female specimens, C-banding analyses detected a cryptic W sex chromosome-associated interstitial heterochromatic segment which is absent in the Z sex chromosome. These ZW sex chromosomes are considered to be in a nascent stage of morphological differentiation and are absent in T. rapicauda collected in Guatemala. The amount, location and fluorochrome affinities of constitutive heterochromatin, the position of the nucleolus organizer region, and the genome sizes of female and male individuals were determined. The previously published cytogenetic data on T. rapicauda are discussed.

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.

Sex chromosomes and karyotype of the (nearly) mythical creature, the Gila monster, Heloderma suspectum (Squamata: Helodermatidae)

A wide variety of sex determination systems exist among squamate reptiles. They can therefore serve as an important model for studies of evolutionary transitions among particular sex determination systems. However, we still have only a limited knowledge of sex determination in certain important lineages of squamates. In this respect, one of the most understudied groups is the family Helodermatidae (Anguimorpha) encompassing the only two venomous species of lizards which are potentially lethal to human beings. We uncovered homomorphic ZZ/ZW sex chromosomes in the Gila monster (Heloderma suspectum) with a highly heterochromatic W chromosome. The sex chromosomes are morphologically similar to the ZZ/ZW sex chromosomes of monitor lizards (Varanidae). If the sex chromosomes of helodermatids and varanids are homologous, female heterogamety may be ancestral for the whole Anguimorpha group. Moreover, we found that the karyotype of the Gila monster consists of 2n = 36 chromosomes (14 larger metacentric chromosomes and 22 acrocentric microchromosomes). 2n = 36 is the widely distributed chromosomal number among squamates. In his pioneering works representing the only previous cytogenetic examination of the family Helodermatidae, Matthey reported the karyotype as 2n = 38 and suggested a different chromosomal morphology for this species. We believe that this was probably erroneously. We also discovered a strong accumulation of telomeric sequences on several pairs of microchromosomes in the Gila monster, which is a trait documented relatively rarely in vertebrates. These new data fill an important gap in our understanding of the sex determination and karyotype evolution of squamates.

Tracing the evolution of amniote chromosomes

A great deal of diversity in chromosome number and arrangement is observed across the amniote phylogeny. Understanding how this diversity is generated is important for determining the role of chromosomal rearrangements in generating phenotypic variation and speciation. Gaining this understanding is achieved by reconstructing the ancestral genome arrangement based on comparisons of genome organization of extant species. Ancestral karyotypes for several amniote lineages have been reconstructed, mainly from cross-species chromosome painting data. The availability of anchored whole genome sequences for amniote species has increased the evolutionary depth and confidence of ancestral reconstructions from those made solely from chromosome painting data. Nonetheless, there are still several key lineages where the appropriate data required for ancestral reconstructions is lacking. This review highlights the progress that has been made towards understanding the chromosomal changes that have occurred during amniote evolution and the reconstruction of ancestral karyotypes.

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.

Molecular cytogenetic map of the central bearded dragon, Pogona vitticeps (Squamata: Agamidae)

Reptiles, as the sister group to birds and mammals, are particularly valuable for comparative genomic studies among amniotes. The Australian central bearded dragon (Pogona vitticeps) is being developed as a reptilian model for such comparisons, with whole-genome sequencing near completion. The karyotype consists of 6 pairs of macrochromosomes and 10 pairs microchromosomes (2n = 32), including a female heterogametic ZW sex microchromosome pair. Here, we present a molecular cytogenetic map for P. vitticeps comprising 87 anchor bacterial artificial chromosome clones that together span each macro- and microchromosome. It is the first comprehensive cytogenetic map for any non-avian reptile. We identified an active nucleolus organizer region (NOR) on the sub-telomeric region of 2q by mapping 18S rDNA and Ag-NOR staining. We identified interstitial telomeric sequences in two microchromosome pairs and the W chromosome, indicating that microchromosome fusion has been a mechanism of karyotypic evolution in Australian agamids within the last 21 to 19 million years. Orthology searches against the chicken genome revealed an intrachromosomal rearrangement of P. vitticeps 1q, identified regions orthologous to chicken Z on P. vitticeps 2q, snake Z on P. vitticeps 6q and the autosomal microchromosome pair in P. vitticeps orthologous to turtle Pelodiscus sinensis ZW and lizard Anolis carolinensis XY. This cytogenetic map will be a valuable reference tool for future gene mapping studies and will provide the framework for the work currently underway to physically anchor genome sequences to chromosomes for this model Australian squamate.

Transitions between sex-determining systems in reptiles and amphibians

Important technological advances in genomics are driving a new understanding of the evolution of sex determination in vertebrates. In particular, comparative chromosome mapping in reptiles has shown an intriguing distribution of homology in sex chromosomes across reptile groups. When this new understanding is combined with the widespread distribution of genetic and temperature-dependent sex-determination mechanisms among reptiles, it is apparent that transitions between modes have occurred many times, as they have for amphibians (particularly between male and female heterogamety). It is also likely that thermosensitivity in sex determination is a key factor in those transitions in reptiles, and possibly in amphibians too. New models of sex determination involving temperature thresholds are providing the framework for the investigation of transitions and making possible key predictions about the homologies and sex-determination patterns expected among taxa in these groups. Molecular cytogenetics and other genomic approaches are essential to providing the fundamental material necessary to make advances in this field.

Strong conservation of the bird Z chromosome in reptilian genomes is revealed by comparative painting despite 275 million years divergence

The divergence of lineages leading to extant squamate reptiles (lizards, snakes, and amphisbaenians) and birds occurred about 275 million years ago. Birds, unlike squamates, have karyotypes that are typified by the presence of a number of very small chromosomes. Hence, a number of chromosome rearrangements might be expected between bird and squamate genomes. We used chromosome-specific DNA from flow-sorted chicken (Gallus gallus) Z sex chromosomes as a probe in cross-species hybridization to metaphase spreads of 28 species from 17 families representing most main squamate lineages and single species of crocodiles and turtles. In all but one case, the Z chromosome was conserved intact despite very ancient divergence of sauropsid lineages. Furthermore, the probe painted an autosomal region in seven species from our sample with characterized sex chromosomes, and this provides evidence against an ancestral avian-like system of sex determination in Squamata. The avian Z chromosome synteny is, therefore, conserved albeit it is not a sex chromosome in these squamate species.

Globin gene structure in a reptile supports the transpositional model for amniote α- and β-globin gene evolution

The haemoglobin protein, required for oxygen transportation in the body, is encoded by α- and β-globin genes that are arranged in clusters. The transpositional model for the evolution of distinct α-globin and β-globin clusters in amniotes is much simpler than the previously proposed whole genome duplication model. According to this model, all jawed vertebrates share one ancient region containing α- and β-globin genes and several flanking genes in the order MPG-C16orf35-(α-β)-GBY-LUC7L that has been conserved for more than 410 million years, whereas amniotes evolved a distinct β-globin cluster by insertion of a transposed β-globin gene from this ancient region into a cluster of olfactory receptors flanked by CCKBR and RRM1. It could not be determined whether this organisation is conserved in all amniotes because of the paucity of information from non-avian reptiles. To fill in this gap, we examined globin gene organisation in a squamate reptile, the Australian bearded dragon lizard, Pogona vitticeps (Agamidae). We report here that the α-globin cluster (HBK, HBA) is flanked by C16orf35 and GBY and is located on a pair of microchromosomes, whereas the β-globin cluster is flanked by RRM1 on the 3' end and is located on the long arm of chromosome 3. However, the CCKBR gene that flanks the β-globin cluster on the 5' end in other amniotes is located on the short arm of chromosome 5 in P. vitticeps, indicating that a chromosomal break between the β-globin cluster and CCKBR occurred at least in the agamid lineage. Our data from a reptile species provide further evidence to support the transpositional model for the evolution of β-globin gene cluster in amniotes.

Sex chromosome evolution in amniotes: applications for bacterial artificial chromosome libraries

Variability among sex chromosome pairs in amniotes denotes a dynamic history. Since amniotes diverged from a common ancestor, their sex chromosome pairs and, more broadly, sex-determining mechanisms have changed reversibly and frequently. These changes have been studied and characterized through the use of many tools and experimental approaches but perhaps most effectively through applications for bacterial artificial chromosome (BAC) libraries. Individual BAC clones carry 100–200 kb of sequence from one individual of a target species that can be isolated by screening, mapped onto karyotypes, and sequenced. With these techniques, researchers have identified differences and similarities in sex chromosome content and organization across amniotes and have addressed hypotheses regarding the frequency and direction of past changes. Here, we review studies of sex chromosome evolution in amniotes and the ways in which the field of research has been affected by the advent of BAC libraries.

A W-linked palindrome and gene conversion in New World sparrows and blackbirds

A hallmark feature of the male-specific region of the human Y chromosome is the presence of large and near-identical palindromes. These palindromes are maintained in a state of near identity via gene conversion between the arms of the palindrome, and both neutral and selection-based theories have been proposed to explain their enrichment on the human Y and X chromosomes. While those proposed theories would be applicable to sex chromosomes in other species, it has not been established whether near-identical palindromes are a common feature of sex chromosomes in a broader range of taxa, including other tetrapods. Here, we report the genomic sequencing and features of a 279-kb region of the non-recombining portion of the W chromosome spanning the CHD1W locus in a New World sparrow, the white-throated sparrow (Zonotrichia albicollis), and the corresponding region on the Z chromosome. As has been observed for other Y and W chromosomes, we detected a high repetitive element content (51%) and low gene content on the white-throated sparrow W chromosome. In addition, we identified a 22-kb near-identical (>99%) palindrome on the W chromosome that flanks the 5' end of the CHD1W gene. Signatures of gene conversion were readily detected between the arms of this palindrome, as was the presence of this palindrome in other New World sparrows and blackbirds. Near-identical palindromes are therefore present on the avian W chromosome and may persist due to the same forces proposed for the enrichment of these elements on the human sex chromosomes.

Are homologies in vertebrate sex determination due to shared ancestry or to limited options?

The same candidate genes and the same autosomes are repeatedly used as sex chromosomes in vertebrates. Are these systems identical by descent, or are some genes or chromosomes intrinsically better at triggering the first steps of sex determination?