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

Comparative Cytogenetics of the Malagasy Ground Geckos of the Paroedura bastardi and Paroedura picta Species Groups

We present a comparative chromosome study of several taxa of the Malagasy ground geckos of the Paroedura bastardi and P. picta species groups. We employed a preliminary molecular analysis using a trait of the mitochondrial 16S rRNA gene (of about 570 bp) to assess the taxonomic status of the samples studied and a cytogenetic analysis with standard karyotyping (5% Giemsa solution), silver staining (Ag-NOR staining) and sequential C-banding (C-banding + Giemsa and + fluorochromes). Our results show that all the taxa studied of the P. bastardi group (P. ibityensis, P. rennerae and P. cf. guibeae) have a similar karyotype composed of 2n = 34 chromosomes, with two metacentric pairs (1 and 3) and all other pairs being acrocentric. Chromosome diversification in the P. bastardi group was mainly linked to the diversification of heteromorphic sex chromosome systems (ZZ/ZW) in P. ibityensis and P. rennerae, while no heteromorphic sex chromosome pair was found in P. cf. guibeae. The two taxa investigated of the P. picta species group (here named P. picta and P. cf. picta based on molecular data) showed the same chromosome number of 2n = 36, mostly acrocentric elements, but differed in the number of metacentric elements, probably as a result of an inversion at chromosome pair 2. We highlight that the genus Paroedura is characterized by the independent diversification of heterogametic sex chromosomes in different evolutionary lineages and, similarly to other phylogenetically related gecko genera, by a progressive formation of a biarmed element by means of tandem fusions and inversions of distinct pairs.

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.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.

Microchromosome fusions underpin convergent evolution of chameleon karyotypes

Evolutionary shifts in chromosome compositions (karyotypes) are major drivers of lineage and genomic diversification. Fusion of ancestral chromosomes is one hypothesized mechanism for the evolutionary reduction of the total chromosome number, a frequently implied karyotypic shift. Empirical tests of this hypothesis require model systems with variable karyotypes, known chromosome features, and a robust phylogeny. Here we used chameleons, diverse lizards with exceptionally variable karyotypes ($2n=20\text{-}62$), to test whether chromosomal fusions explain the repeated evolution of karyotypes with fewer chromosomes than ancestral karyotypes. Using a multidisciplinary approach including cytogenetic analyses and phylogenetic comparative methods, we found that a model of constant loss through time best explained chromosome evolution across the chameleon phylogeny. Next, we tested whether fusions of microchromosomes into macrochromosomes explained these evolutionary losses using generalized linear models. Multiple comparisons supported microchromosome fusions as the predominant agent of evolutionary loss. We further compared our results to various natural history traits and found no correlations. As such, we infer that the tendency of microchromosomes to fuse was a quality of the ancestral chameleon genome and that the genomic predisposition of ancestors is a more substantive predictor of chromosome change than the ecological, physiological, and biogeographical factors involved in their diversification.

A lizard is never late: Squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012 and 2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

Gene expression of male pathway genes sox9 and amh during early sex differentiation in a reptile departs from the classical amniote model

BACKGROUND

Sex determination is the process whereby the bipotential embryonic gonads become committed to differentiate into testes or ovaries. In genetic sex determination (GSD), the sex determining trigger is encoded by a gene on the sex chromosomes, which activates a network of downstream genes; in mammals these include SOX9, AMH and DMRT1 in the male pathway, and FOXL2 in the female pathway. Although mammalian and avian GSD systems have been well studied, few data are available for reptilian GSD systems.

RESULTS

We conducted an unbiased transcriptome-wide analysis of gonad development throughout differentiation in central bearded dragon (Pogona vitticeps) embryos with GSD. We found that sex differentiation of transcriptomic profiles occurs at a very early stage, before the gonad consolidates as a body distinct from the gonad-kidney complex. The male pathway genes dmrt1 and amh and the female pathway gene foxl2 play a key role in early sex differentiation in P. vitticeps, but the central player of the mammalian male trajectory, sox9, is not differentially expressed in P. vitticeps at the bipotential stage. The most striking difference from GSD systems of other amniotes is the high expression of the male pathway genes amh and sox9 in female gonads during development. We propose that a default male trajectory progresses if not repressed by a W-linked dominant gene that tips the balance of gene expression towards the female trajectory. Further, weighted gene expression correlation network analysis revealed novel candidates for male and female sex differentiation.

CONCLUSION

Our data reveal that interpretation of putative mechanisms of GSD in reptiles cannot solely depend on lessons drawn from mammals.

A lizard is never late: squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012-2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

Chromosome-Level Genome Assembly of the Cape Cliff Lizard (Hemicordylus capensis)

Squamates represent a highly diverse and species-rich vertebrate group that is remarkably understudied from a genomic perspective. A scarcity of genomic data is particularly evident for scincomorph lizards, which encompass over 10% of all living squamates, and for which high-quality genomic resources are currently lacking. To address this knowledge gap, we present the first chromosome-level reference genome for this group, generated from a male Cape cliff lizard (Hemicordylus capensis), using highly accurate PacBio HiFi long-read sequencing data, long-range Omni-C chromosomal conformation capture data and transcriptomic data for annotation. The rHemCap1.1 genome assembly spans 2.29 Gb, with a scaffold N50 of 359.65 Mb, and includes 25,300 protein-coding genes, with a BUSCO completeness score of 95.5% (sauropsida_odb10). We have generated the most contiguous and complete chromosome-level squamate reference genome assembly publicly available to date. Furthermore, we used short-read resequencing of 35 males and females and applied a differential coverage approach to infer the sex-determination system of the species, which was previously unknown. Our results suggest this species has XX/XY sex chromosomes, representing the first evidence of sex determination in the family Cordylidae. This reference genome will help to establish this species as an evolutionary model for studying variation in body armor, a key trait in cordylids and other squamate groups. Lastly, this is the first squamate reference genome from a continental African species and, as such, represents a valuable resource not only for further evolutionary research in cordylids but also in closely related groups.

Cytogenetic Analysis of Seven Species of Gekkonid and Phyllodactylid Geckos

Geckos (Gekkota), the species-rich clade of reptiles with more than 2200 currently recognized species, demonstrate a remarkable variability in diploid chromosome numbers (2n = 16-48) and mode of sex determination. However, only a small fraction of gekkotan species have been studied with cytogenetic methods. Here, we applied both conventional (karyotype reconstruction and C-banding) and molecular (fluorescence in situ hybridization with probes for rDNA loci and telomeric repeats) cytogenetic analyses in seven species of geckos, namely Blaesodactylus boivini, Chondrodactylus laevigatus, Gekko badenii, Gekko cf. lionotum, Hemidactylus sahgali, Homopholis wahlbergii (Gekkonidae) and Ptyodactylus togoensis (Phyllodactylidae), in order to provide further insights into the evolution of karyotypes in geckos. Our analysis revealed the presence of interstitial telomeric repeats in four species, but we were not able to conclude if they are remnants of previous chromosome rearrangements or were formed by an accumulation of telomeric-like satellite motifs. Even though sex chromosomes were previously identified in several species from the genera Hemidactylus and Gekko by cytogenetic and/or genomic methods, they were not detected by us in any examined species. Our examined species either have poorly differentiated sex chromosomes or, possibly, environmental sex determination. Future studies should explore the effect of temperature and conduct genome-wide analyses in order to identify the mode of sex determination in these species.

Identification of Iguania Ancestral Syntenic Blocks and Putative Sex Chromosomes in the Veiled Chameleon (Chamaeleo calyptratus, Chamaeleonidae, Iguania)

The veiled chameleon (Chamaeleo calyptratus) is a typical member of the family Chamaeleonidae and a promising object for comparative cytogenetics and genomics. The karyotype of C. calyptratus differs from the putative ancestral chameleon karyotype (2n = 36) due to a smaller chromosome number (2n = 24) resulting from multiple chromosome fusions. The homomorphic sex chromosomes of an XX/XY system were described recently using male-specific RADseq markers. However, the chromosomal pair carrying these markers was not identified. Here we obtained chromosome-specific DNA libraries of C. calyptratus by chromosome flow sorting that were assigned by FISH and sequenced. Sequence comparison with three squamate reptiles reference genomes revealed the ancestral syntenic regions in the C. calyptratus chromosomes. We demonstrated that reducing the chromosome number in the C. calyptratus karyotype occurred through two fusions between microchromosomes and four fusions between micro-and macrochromosomes. PCR-assisted mapping of a previously described Y-specific marker indicates that chromosome 5 may be the sex chromosome pair. One of the chromosome 5 conserved synteny blocks shares homology with the ancestral pleurodont X chromosome, assuming parallelism in the evolution of sex chromosomes from two basal Iguania clades (pleurodonts and acrodonts). The comparative chromosome map produced here can serve as the foundation for future genome assembly of chameleons and vertebrate-wide comparative genomic studies.

Sex reversal explains some, but not all, climate-mediated sex ratio variation within a viviparous reptile

Evolutionary transitions in sex-determining systems have occurred frequently yet understanding how they occur remains a major challenge. In reptiles, transitions from genetic to temperature-dependent sex determination can occur if the gene products that determine sex evolve thermal sensitivity, resulting in sex-reversed individuals. However, evidence of sex reversal is limited to oviparous reptiles. Here we used thermal experiments to test whether sex reversal is responsible for differences in sex determination in a viviparous reptile, Carinascincus ocellatus, a species with XY sex chromosomes and population-specific sex ratio response to temperature. We show that sex reversal is occurring and that its frequency is related to temperature. Sex reversal was unidirectional (phenotypic males with XX genotype) and observed in both high- and low-elevation populations. We propose that XX-biased genotypic sex ratios could produce either male- or female-biased phenotypic sex ratios as observed in low-elevation C. ocellatus under variable rates of XX sex reversal. We discuss reasons why sex reversal may not influence sex ratios at high elevation. Our results suggest that the mechanism responsible for evolutionary transitions from genotypic to temperature-dependent sex determination is more complex than can be explained by a single process such as sex reversal.

Findzx: an automated pipeline for detecting and visualising sex chromosomes using whole-genome sequencing data

BACKGROUND

Sex chromosomes have evolved numerous times, as revealed by recent genomic studies. However, large gaps in our knowledge of sex chromosome diversity across the tree of life remain. Filling these gaps, through the study of novel species, is crucial for improved understanding of why and how sex chromosomes evolve. Characterization of sex chromosomes in already well-studied organisms is also important to avoid misinterpretations of population genomic patterns caused by undetected sex chromosome variation.

RESULTS

Here we present findZX, an automated Snakemake-based computational pipeline for detecting and visualizing sex chromosomes through differences in genome coverage and heterozygosity between any number of males and females. A main feature of the pipeline is the option to perform a genome coordinate liftover to a reference genome of another species. This allows users to inspect sex-linked regions over larger contiguous chromosome regions, while also providing important between-species synteny information. To demonstrate its effectiveness, we applied findZX to publicly available genomic data from species belonging to widely different taxonomic groups (mammals, birds, reptiles, and fish), with sex chromosome systems of different ages, sizes, and levels of differentiation. We also demonstrate that the liftover method is robust over large phylogenetic distances (> 80 million years of evolution).

CONCLUSIONS

With findZX we provide an easy-to-use and highly effective tool for identification of sex chromosomes. The pipeline is compatible with both Linux and MacOS systems, and scalable to suit different computational platforms.

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.

Dynamics of epigenetic modifiers and environmentally sensitive proteins in a reptile with temperature induced sex reversal

The mechanisms by which sex is determined, and how a sexual phenotype is stably maintained during adulthood, has been the focus of vigorous scientific inquiry. Resources common to the biomedical field (automated staining and imaging platforms) were leveraged to provide the first immunofluorescent data for a reptile species with temperature induced sex reversal. Two four-plex immunofluorescent panels were explored across three sex classes (sex reversed ZZf females, normal ZWf females, and normal ZZm males). One panel was stained for chromatin remodelling genes JARID2 and KDM6B, and methylation marks H3K27me3, and H3K4me3 (Jumonji Panel). The other CaRe panel stained for environmental response genes CIRBP and RelA, and H3K27me3 and H3K4me3. Our study characterised tissue specific expression and cellular localisation patterns of these proteins and histone marks, providing new insights to the molecular characteristics of adult gonads in a dragon lizard Pogona vitticeps. The confirmation that mammalian antibodies cross react in P. vitticeps paves the way for experiments that can take advantage of this new immunohistochemical resource to gain a new understanding of the role of these proteins during embryonic development, and most importantly for P. vitticeps, the molecular underpinnings of sex reversal.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

Whole-chromosome fusions in the karyotype evolution of Sceloporus (Iguania, Reptilia) are more frequent in sex chromosomes than autosomes

Whole-chromosome fusions play a major role in the karyotypic evolution of reptiles. It has been suggested that certain chromosomes tend to fuse with sex chromosomes more frequently than others. However, the comparative genomic synteny data are too scarce to draw strong conclusions. We obtained and sequenced chromosome-specific DNA pools of Sceloporus malachiticus, an iguanian species which has experienced many chromosome fusions. We found that four of seven lineage-specific fusions involved sex chromosomes, and that certain syntenic blocks which constitute the sex chromosomes, such as the homologues of the Anolis carolinensis chromosomes 11 and 16, are repeatedly involved in sex chromosome formation in different squamate species. To test the hypothesis that the karyotypic shift could be associated with changes in recombination patterns, we performed a synaptonemal complex analysis in this species and in Sceloporus variabilis (2n = 34). It revealed that the sex chromosomes in S. malachiticus had two distal pseudoautosomal regions and a medial differentiated region. We found that multiple fusions little affected the recombination rate in S. malachiticus. Our data confirm more frequent involvement of certain chromosomes in sex chromosome formation, but do not reveal a connection between the gonosome-autosome fusions and the evolution of recombination rate. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

With or Without W? Molecular and Cytogenetic Markers are Not Sufficient for Identification of Environmentally-Induced Sex Reversal in the Bearded Dragon

Transitions from environmental sex determination (ESD) to genotypic sex determination (GSD) require an intermediate step of sex reversal, i.e., the production of individuals with a mismatch between the ancestral genotypic and the phenotypic sex. Among amniotes, the sole well-documented transition in this direction was shown in the laboratory in the central bearded dragon, Pogona vitticeps, where very high incubation temperatures led to the production of females with the male-typical (ZZ) genotype. These sex-reversed females then produced offspring whose sex depended on the incubation temperature. Sex-reversed animals identified by molecular and cytogenetic markers were also reported in the field, and their increasing incidence was speculated as a climate warming-driven transition in sex determination. We show that the molecular and cytogenetic markers normally sex-linked in P. vitticeps are also sex-linked in P. henrylawsoni and P. minor, which points to quite ancient sex chromosomes in this lineage. Nevertheless, we demonstrate, based on a crossing experiment with a male bearded dragon who possesses a mismatch between phenotypic sex and genotype, that the used cytogenetic and molecular markers might not be reliable for the identification of sex reversal. Sex reversal should not be considered as the only mechanism causing a mismatch between genetic sex-linked markers and phenotypic sex, which can emerge also by other processes, here most likely by a rare recombination between regions of sex chromosomes which are normally sex-linked. We warn that sex-linked, even apparently for a long evolutionary time, and sex-specific molecular and cytogenetic markers are not a reliable tool for the identification of sex-reversed individuals in a population and that sex reversal has to be verified by other approaches, particularly by observation of the sex ratio of the progeny.

Chromosomal Analysis in Crotophaga ani (Aves, Cuculiformes) Reveals Extensive Genomic Reorganization and an Unusual Z-Autosome Robertsonian Translocation

Although cytogenetics studies in cuckoos (Aves, Cuculiformes) have demonstrated an interesting karyotype variation, such as variations in the chromosome morphology and diploid number, their chromosome organization and evolution, and relation with other birds are poorly understood. Hence, we combined conventional and molecular cytogenetic approaches to investigate chromosome homologies between chicken and the smooth-billed ani (Crotophaga ani). Our results demonstrate extensive chromosome reorganization in C. ani, with interchromosomal rearrangements involving macro and microchromosomes. Intrachromosomal rearrangements were observed in some macrochromosomes, including the Z chromosome. The most evolutionary notable finding was a Robertsonian translocation between the microchromosome 17 and the Z chromosome, a rare event in birds. Additionally, the simple short repeats (SSRs) tested here were preferentially accumulated in the microchromosomes and in the Z and W chromosomes, showing no relationship with the constitutive heterochromatin regions, except in the W chromosome. Taken together, our results suggest that the avian sex chromosome is more complex than previously postulated and revealed the role of microchromosomes in the avian sex chromosome evolution, especially cuckoos.

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.

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.

Origin, clonal diversity, and evolution of the parthenogenetic lizard Darevskia unisexualis

BACKGROUND

The hybridization of female D. raddei and male D. valentini gave rise to the parthenogenetic Caucasian rock lizard Darevskia unisexualis. A previously identified genetic polymorphism in the species consisted of one common and two allozyme clones. Analysis of microsatellites and single nucleotide polymorphisms (SNPs) from the three species yields estimates of clonal diversity and tests the hypothesis of a single origin for D. unisexualis.

RESULTS

Genotyping and sequencing of four microsatellite-containing loci for 109 specimens of D. unisexualis, 17 D. valentini, and 45 D. raddei nairensis identified 12 presumptive clones, including one widespread and 11 rare clones. Most individuals in some localities had a rare clone. Clone-specific alleles in D. unisexualis were compared with those of the parental species. The results inferred a single hybridization event. Post-formation mutations best explain the less common clones.

CONCLUSIONS

Interspecific analyses identify alleles inherited by D. unisexualis from its bisexual ancestors. SNP analyses fail to reject the hypothesis of a single interspecific origin of D. unisexualis, followed by microsatellite mutations in this initial clone. Microsatellites detect higher clonal diversity in D. unisexualis compared to allozymes and identify the likely origins of clones. Our approach may be applicable to other unisexual species whose origins involve interspecific hybridization.

First Report of Sex Chromosomes in Night Lizards (Scincoidea: Xantusiidae)

Squamate reptiles (lizards, snakes, and amphibians) are an outstanding group for studying sex chromosome evolution—they are old, speciose, geographically widespread, and exhibit myriad sex-determining modes. Yet, the vast majority of squamate species lack heteromorphic sex chromosomes. Cataloging the sex chromosome systems of species lacking easily identifiable, heteromorphic sex chromosomes, therefore, is essential before we are to fully understand the evolution of vertebrate sex chromosomes. Here, we use restriction site-associated DNA sequencing (RADseq) to classify the sex chromosome system of the granite night lizard, Xantusia henshawi. RADseq is an effective alternative to traditional cytogenetic methods for determining a species’ sex chromosome system (i.e., XX/XY or ZZ/ZW), particularly in taxa with non-differentiated sex chromosomes. Although many xantusiid lineages have been karyotyped, none possess heteromorphic sex chromosomes. We identified a ZZ/ZW sex chromosome system in X. henshawi—the first such data for this family. Furthermore, we report that the X. henshawi sex chromosome contains fragments of genes found on Gallus gallus chromosomes 7, 12, and 18 (which are homologous to Anolis carolinensis chromosome 2), the first vertebrate sex chromosomes to utilize this linkage group.

Reproductive phenotype predicts adult bite-force performance in sex-reversed dragons (Pogona vitticeps)

Sex-related differences in morphology and behavior are well documented, but the relative contributions of genes and environment to these traits are less well understood. Species that undergo sex reversal, such as the central bearded dragon (Pogona vitticeps), offer an opportunity to better understand sexually dimorphic traits because sexual phenotypes can exist on different chromosomal backgrounds. Reproductively female dragons with a discordant sex chromosome complement (sex reversed), at least as juveniles, exhibit traits in common with males (e.g., longer tails and greater boldness). However, the impact of sex reversal on sexually dimorphic traits in adult dragons is unknown. Here, we investigate the effect of sex reversal on bite-force performance, which may be important in resource acquisition (e.g., mates and/or food). We measured body size, head size, and bite force of the three sexual phenotypes in a colony of captive animals. Among adults, we found that males (ZZm) bite more forcefully than either chromosomally concordant females (ZWf) or sex-reversed females (ZZf), and this difference is associated with having relatively larger head dimensions. Therefore, adult sex-reversed females, despite apparently exhibiting male traits as juveniles, do not develop the larger head and enhanced bite force of adult male bearded dragons. This pattern is further illustrated in the full sample by a lack of positive allometry of bite force in sex-reversed females that is observed in males. The results reveal a close association between reproductive phenotype and bite force performance, regardless of sex chromosome complement.

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.

Escaping the evolutionary trap? Sex chromosome turnover in basilisks and related lizards (Corytophanidae: Squamata)

Most pleurodont lizard families (anoles, iguanas and their relatives), with the exception of the basilisks and casquehead lizards (family Corytophanidae), share homologous XX/XY sex chromosomes, syntenic with chicken chromosome 15. Here, we used a suite of methods (i.e. RADseq, RNAseq and qPCR) to identify corytophanid sex chromosomes for the first time. We reveal that all examined corytophanid species have partially degenerated XX/XY sex chromosomes, syntenic with chicken chromosome 17. Transcriptomic analyses showed that the expression of X-linked genes in the corytophanid, Basiliscus vittatus, is not balanced between the sexes, which is rather exceptional under male heterogamety, and unlike the dosage-balanced sex chromosomes in other well-studied XX/XY systems, including the green anole, Anolis carolinensis. Corytophanid sex chromosomes may represent a rare example of a turnover away from stable, differentiated sex chromosomes. However, because of poor phylogenetic resolution among pleurodont families, we cannot reject the alternative hypothesis that corytophanid sex chromosomes evolved independently from an unknown ancestral system.

Chromosomics: Bridging the Gap between Genomes and Chromosomes

The recent advances in DNA sequencing technology are enabling a rapid increase in the number of genomes being sequenced. However, many fundamental questions in genome biology remain unanswered, because sequence data alone is unable to provide insight into how the genome is organised into chromosomes, the position and interaction of those chromosomes in the cell, and how chromosomes and their interactions with each other change in response to environmental stimuli or over time. The intimate relationship between DNA sequence and chromosome structure and function highlights the need to integrate genomic and cytogenetic data to more comprehensively understand the role genome architecture plays in genome plasticity. We propose adoption of the term 'chromosomics' as an approach encompassing genome sequencing, cytogenetics and cell biology, and present examples of where chromosomics has already led to novel discoveries, such as the sex-determining gene in eutherian mammals. More importantly, we look to the future and the questions that could be answered as we enter into the chromosomics revolution, such as the role of chromosome rearrangements in speciation and the role more rapidly evolving regions of the genome, like centromeres, play in genome plasticity. However, for chromosomics to reach its full potential, we need to address several challenges, particularly the training of a new generation of cytogeneticists, and the commitment to a closer union among the research areas of genomics, cytogenetics, cell biology and bioinformatics. Overcoming these challenges will lead to ground-breaking discoveries in understanding genome evolution and function.

Distribution and amplification of interstitial telomeric sequences (ITSs) in Australian dragon lizards support frequent chromosome fusions in Iguania

Telomeric sequences are generally located at the ends of chromosomes; however, they can also be found in non-terminal chromosomal regions when they are known as interstitial telomeric sequences (ITSs). Distribution of ITSs across closely related and divergent species elucidates karyotype evolution and speciation as ITSs provide evolutionary evidence for chromosome fusion. In this study, we performed physical mapping of telomeric repeats by fluorescence in situ hybridisation (FISH) in seven Australian dragon lizards thought to represent derived karyotypes of squamate reptiles and a gecko lizard with considerably different karyotypic feature. Telomeric repeats were present at both ends of all chromosomes in all species, while varying numbers of ITSs were also found on microchromosomes and in pericentromeric or centromeric regions on macrochromosomes in five lizard species examined. This suggests that chromosomal rearrangements from ancestral squamate reptiles to Iguania occurred mainly by fusion between ancestral types of acrocentric chromosomes and/or between microchromosomes, leading to appearance of bi-armed macrochromosomes, and in the reduction of microchromosome numbers. These results support the previously proposed hypothesis of karyotype evolution in squamate reptiles. In addition, we observed the presence of telomeric sequences in the similar regions to heterochromatin of the W microchromosome in Pogona barbata and Doporiphora nobbi, while sex chromosomes for the two species contained part of the nucleolar organiser regions (NORs). This likely implies that these ITSs are a part of the satellite DNA and not relics of chromosome fusions. Amplification of telomeric repeats may have involved heterochromatinisation of sex-specific W chromosomes and play a role in the organisation of the nucleolus.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Multiple interspecific hybridization and microsatellite mutations provide clonal diversity in the parthenogenetic rock lizard Darevskia armeniaca

Background

The parthenogenetic Caucasian rock lizard Darevskia armeniaca, like most other parthenogenetic vertebrate species, originated through interspecific hybridization between the closely related sexual Darevskia mixta and Darevskia valentini. Darevskia armeniaca was shown to consist of one widespread allozyme clone and a few rare ones, but notwithstanding the origin of clonal diversity remains unclear. We conduct genomic analysis of D. armeniaca and its parental sexual species using microsatellite and SNP markers to identify the origin of parthenogenetic clonal lineages.

Results

Four microsatellite-containing loci were genotyped for 111 specimens of D. armeniaca, 17 D. valentini, and four D. mixta. For these species, a total of 47 alleles were isolated and sequenced. Analysis of the data revealed 13 genotypes or presumptive clones in parthenogenetic D. armeniaca, including one widespread clone, two apparently geographically restricted clones, and ten rare clones. Comparisons of genotype-specific markers in D. armeniaca with those of its parental species revealed three founder-events including a common and two rare clones. All other clones appeared to have originated via post-formation microsatellite mutations in the course of evolutionary history of D. armeniaca.

Conclusion

Our new approach to microsatellite genotyping reveals allele-specific microsatellite and SNP markers for each locus studied. Interspecies comparison of these markers identifies alleles inherited by parthenospecies from parental species, and provides new information on origin and evolution of clonal diversity in D. armeniaca. SNP analyses reveal at least three interspecific origins of D. armeniaca, and microsatellite mutations in these initial clones give rise to new clones. Thus, we first establish multiple origins of D. armeniaca. Our study identifies the most effective molecular markers for elucidating the origins of clonal diversity in other unisexual species that arose via interspecific hybridization.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5359-5) contains supplementary material, which is available to authorized users.

Developmental asynchrony and antagonism of sex determination pathways in a lizard with temperature-induced sex reversal

Vertebrate sex differentiation follows a conserved suite of developmental events: the bipotential gonads differentiate and shortly thereafter sex specific traits become dimorphic. However, this may not apply to squamates, a diverse vertebrate lineage comprising of many species with thermosensitive sexual development. Of the three species with data on the relative timing of gonad differentiation and genital dimorphism, the females of two (Niveoscincus ocellatus and Barisia imbricata) exhibit a phase of temporary pseudohermaphroditism or TPH (gonads have differentiated well before genital dimorphism). We report a third example of TPH in Pogona vitticeps, an agamid with temperature-induced male to female sex reversal. These findings suggest that for female squamates, genital and gonad development may not be closely synchronised, so that TPH may be common. We further observed a high frequency of ovotestes, a usually rare gonadal phenotype characterised by a mix of male and female structures, exclusively associated with temperature-induced sex reversal. We propose that ovotestes are evidence of a period of antagonism between male and female sex-determining pathways during sex reversal. Female sexual development in squamates is considerably more complex than has been appreciated, providing numerous avenues for future exploration of the genetic and hormonal cues that govern sexual development.

Adaptation and conservation insights from the koala genome

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

The origin of multiple clones in the parthenogenetic lizard species Darevskia rostombekowi

The all-female Caucasian rock lizard Darevskia rostombekowi and other unisexual species of this genus reproduce normally via true parthenogenesis. Typically, diploid parthenogenetic reptiles exhibit some amount of clonal diversity. However, allozyme data from D. rostombekowi have suggested that this species consists of a single clone. Herein, we test this hypothesis by evaluating variation at three variable microsatellite loci for 42 specimens of D. rostombekowi from four populations in Armenia. Analyses based on single nucleotide polymorphisms of each locus reveal five genotypes or presumptive clones in this species. All individuals are heterozygous at the loci. The major clone occurs in 24 individuals and involves three populations. Four rare clones involve one or several individuals from one or two populations. Most variation owes to parent-specific single nucleotide polymorphisms, which occur as heterozygotes. This result fails to reject the hypothesis of a single hybridization founder event that resulted in the initial formation of one major clone. The other clones appear to have originated via post-formation microsatellite mutations of the major clone.

Recent advances and perspectives in the design and development of polymyxins

Covering: 1947-early 2017, particularly from 2005-early 2017The rise of bacterial pathogens with acquired resistance to almost all available antibiotics is becoming a serious public health issue. Polymyxins, antibiotics that were mostly abandoned a few decades ago because of toxicity concerns, are ultimately considered as a last-line therapy to treat infections caused by multi-drug resistant Gram-negative bacteria. This review surveys the progress in understanding polymyxin structure, and their chemistry, mechanisms of antibacterial activity and nephrotoxicity, biomarkers, synergy and combination with other antimicrobial agents and antibiofilm properties. An update of recent efforts in the design and development of a new generation of polymyxin drugs is also discussed. A novel approach considering the modification of the scaffold of polymyxins to integrate metabolism and detoxification issues into the drug design process is a promising new line to potentially prevent accumulation in the kidneys and reduce nephrotoxicity.

Immunocytological analysis of meiotic recombination in two anole lizards (Squamata, Dactyloidae)

Although the evolutionary importance of meiotic recombination is not disputed, the significance of interspecies differences in the recombination rates and recombination landscapes remains under-appreciated. Recombination rates and distribution of chiasmata have been examined cytologically in many mammalian species, whereas data on other vertebrates are scarce. Immunolocalization of the protein of the synaptonemal complex (SYCP3), centromere proteins and the mismatch-repair protein MLH1 was used, which is associated with the most common type of recombination nodules, to analyze the pattern of meiotic recombination in the male of two species of iguanian lizards, Anolis carolinensis Voigt, 1832 and Deiroptyx coelestinus (Cope, 1862). These species are separated by a relatively long evolutionary history although they retain the ancestral iguanian karyotype. In both species similar and extremely uneven distributions of MLH1 foci along the macrochromosome bivalents were detected: approximately 90% of crossovers were located at the distal 20% of the chromosome arm length. Almost total suppression of recombination in the intermediate and proximal regions of the chromosome arms contradicts the hypothesis that "homogenous recombination" is responsible for the low variation in GC content across the anole genome. It also leads to strong linkage disequilibrium between the genes located in these regions, which may benefit conservation of co-adaptive gene arrays responsible for the ecological adaptations of the anoles.