Phylogenomics resolves the puzzling phylogeny of banded newts (genus Ommatotriton)

Resolving the order of speciation events that occurred in rapid succession is inherently hard and typically requires a phylogenomic approach. A case in point concerns the previously unresolved phylogeny of the three species of banded newt (genus Ommatotriton). We obtain c. 7k nuclear DNA markers using target enrichment by sequence capture and analyze the dataset using maximum likelihood inference of concatenated data with RAxML, summary multi-species coalescent analysis with ASTRAL and Bayesian species tree inference using a diffusion model with SNAPPER, and use TreeMix and PhyloNet to test for interspecific gene flow. All analyses recover three distinct species with no evidence of interspecific gene flow. All analyses retrieved the topology (O. nesterovi, (O. ophryticus, O. vittatus)), with high support. SNAPPER did show the tendency to get stuck in a local optimum, resulting in a different but still highly supported topology. Furthermore, we notice that fewer SNAPPER runs get stuck in a local optimum when we include an outgroup. Therefore, we recommend the exploration of multiple independent runs and the use of an outgroup with this approach. The banded newt radiation illustrates the use of genome-wide data to tackle formerly unresolved phylogenies.

Ancient structural variants control sex-specific flowering time morphs in walnuts and hickories

Balanced mating type polymorphisms offer a distinct window into the forces shaping sexual reproduction strategies. Multiple hermaphroditic genera in Juglandaceae, including walnuts (Juglans) and hickories (Carya), show a 1:1 genetic dimorphism for male versus female flowering order (heterodichogamy). We map two distinct Mendelian inheritance mechanisms to ancient (>37 million years old) genus-wide structural DNA polymorphisms. The dominant haplotype for female-first flowering in Juglans contains tandem repeats of the 3' untranslated region of a gene putatively involved in trehalose-6-phosphate metabolism and is associated with increased cis gene expression in developing male flowers, possibly mediated by small RNAs. The Carya locus contains ~20 syntenic genes and shows molecular signatures of sex chromosome-like evolution. Inheritance mechanisms for heterodichogamy are deeply conserved, yet may occasionally turn over, as in sex determination.

Neo-Sex Chromosome Evolution in Treehoppers Despite Long-Term X Chromosome Conservation

Sex chromosomes follow distinct evolutionary trajectories compared to the rest of the genome. In many cases, sex chromosomes (X and Y or Z and W) significantly differentiate from one another resulting in heteromorphic sex chromosome systems. Such heteromorphic systems are thought to act as an evolutionary trap that prevents subsequent turnover of the sex chromosome system. For old, degenerated sex chromosome systems, chromosomal fusion with an autosome may be one way that sex chromosomes can "refresh" their sequence content. We investigated these dynamics using treehoppers (hemipteran insects of the family Membracidae), which ancestrally have XX/X0 sex chromosomes. We assembled the most complete reference assembly for treehoppers to date for Umbonia crassicornis and employed comparative genomic analyses of 12 additional treehopper species to analyze X chromosome variation across different evolutionary timescales. We find that the X chromosome is largely conserved, with one exception being an X-autosome fusion in Calloconophora caliginosa. We also compare the ancestral treehopper X with other X chromosomes in Auchenorrhyncha (the clade containing treehoppers, leafhoppers, spittlebugs, cicadas, and planthoppers), revealing X conservation across more than 300 million years. These findings shed light on chromosomal evolution dynamics in treehoppers and the role of chromosomal rearrangements in sex chromosome evolution.

Pivotal temperature is not for everyone: Evidence for temperature-dependent sex determination in three gecko species

The prevalence of environmental sex determination (ESD) in squamate reptiles is often overestimated in the literature. This is surprising because we have reliable data demonstrating ESD in only a few species. The documentation of ESD in three species of geckos presented here has significantly increased our knowledge, given that satisfactory evidence for ESD existed in only eight other gecko species. For the first time, we document the occurrence of ESD in the family Sphaerodactylidae. Our finding of unexpected variability in the shapes of reaction norms among geckos highlights that traditional descriptions using parameters such as pivotal temperature, that is, temperature producing a 50:50 sex ratio, are unsatisfactory. For example, the gecko Pachydactylus tigrinus lacks any pivotal temperature and its sex ratios are strongly female-biased across the entire range of viable temperatures. We argue for the effective capture of the relationship between temperature and sex ratio using specific nonlinear models rather than using classical simplistic descriptions and classifications of reaction norms.

Differentiated sex chromosomes, karyotype evolution, and spontaneous triploidy in carphodactylid geckos

Geckos exhibit derived karyotypes without a clear distinction between macrochromosomes and microchromosomes and intriguing diversity in sex determination mechanisms. We conducted cytogenetic analyses in six species from the genera Nephrurus, Phyllurus, and Saltuarius of the gecko family Carphodactylidae. We confirmed the presence of a female heterogametic system with markedly differentiated and heteromorphic sex chromosomes in all examined species, typically with the W chromosome notably larger than the Z chromosome. One species, Nephrurus cinctus, possesses unusual multiple Z1Z1Z2Z2/Z1Z2W sex chromosomes. The morphology of the sex chromosomes, along with repetitive DNA content, suggests that the differentiation or emergence of sex chromosomes occurred independently in the genus Phyllurus. Furthermore, our study unveils a case of spontaneous triploidy in a fully grown individual of Saltuarius cornutus (3n = 57) and explores its implications for reproduction in carphodactylid geckos. We revealed that most carphodactylids retain the putative ancestral gekkotan karyotype of 2n = 38, characterized by predominantly acrocentric chromosomes that gradually decrease in size. If present, biarmed chromosomes emerge through pericentric inversions, maintaining the chromosome (and centromere) numbers. However, Phyllurus platurus is a notable exception, with a karyotype of 2n = 22 chromosomes. Its eight pairs of biarmed chromosomes were probably formed by Robertsonian fusions of acrocentric chromosomes. The family underscores a remarkable instance of evolutionary stability in chromosome numbers, followed by a profound transformation through parallel interchromosomal rearrangements. Our study highlights the need to continue generating cytogenetic data in order to test long-standing ideas about reproductive biology and the evolution of genome and sex determination.

Environmental (and Random?) Sex Determination in Endangered and Invasive Phelsuma Geckos

INTRODUCTION

Sex is a fundamental characteristic of an individual. It is therefore puzzling why in some systems sex is precisely determined by a genotype, while in others it is influenced by the environment or even subtle, not to say random, factors. Some stochasticity in sex determination would be expected if environmental conditions did not have a large sex-specific effect on fitness. Although data are only available for a small fraction of species, geckos show exceptional variability in sex determination.

METHODS

We tested the effects of three incubation temperatures on sex ratio and adult body size in the invasive gecko Phelsuma laticauda and the vulnerable gecko Phelsuma nigristriata.

RESULTS

We document temperature-dependent sex determination (TSD) in both species. Only females hatched at a low temperature (24°C), whereas male production peaked at an intermediate temperature (28°C) and declined, at least in P. laticauda, again at the highest temperature (31°C). Interestingly, full siblings hatched from eggs glued together during the incubation at temperatures producing both sexes are often of the opposite sex. We found no significant effect of incubation temperature on adult body length.

CONCLUSIONS

Documentation of TSD in the day geckos has implications for conservation practice in environmental management of endangered species or eradication of invasive species. However, it appears that a very subtle (random?) factor may also be involved in their sex determination. In line with this, we found no significant effect of incubation temperature on adult body length.

Using a Handful of Transcriptomes to Detect Sex-Linked Markers and Develop Molecular Sexing Assays in a Species with Homomorphic Sex Chromosomes

To understand the biology of a species, it is often crucial to be able to differentiate males and females. However, many species lack easily identifiable sexually dimorphic traits. In those that possess sex chromosomes, molecular sexing offers a good alternative, and molecular sexing assays can be developed through the comparison of male and female genomic sequences. However, in many nonmodel species, sex chromosomes are poorly differentiated, and identifying sex-linked sequences and developing sexing assays can be challenging. In this study, we highlight a simple transcriptome-based procedure for the detection of sex-linked markers suitable for the development of sexing assays that circumvents limitations of more commonly used approaches. We apply it to the spotted snow skink Carinascincus ocellatus, a viviparous lizard with homomorphic XY chromosomes that has environmentally induced sex reversal. With transcriptomes from three males and three females alone, we identify thousands of putative Y-linked sequences. We confirm linkage through alignment of assembled transcripts to a distantly related lizard genome and readily design multiple single locus polymerase chain reaction primers to sex C. ocellatus and related species. Our approach also facilitates valuable comparisons of sex determining systems on a broad taxonomic scale.

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.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Heteromorphic ZZ/ZW sex chromosomes sharing gene content with mammalian XX/XY are conserved in Madagascan chameleons of the genus Furcifer

Chameleons are well-known lizards with unique morphology and physiology, but their sex determination has remained poorly studied. Madagascan chameleons of the genus Furcifer have cytogenetically distinct Z and W sex chromosomes and occasionally Z1Z1Z2Z2/Z1Z2W multiple neo-sex chromosomes. To identify the gene content of their sex chromosomes, we microdissected and sequenced the sex chromosomes of F. oustaleti (ZZ/ZW) and F. pardalis (Z1Z1Z2Z2/Z1Z2W). In addition, we sequenced the genomes of a male and a female of F. lateralis (ZZ/ZW) and F. pardalis and performed a comparative coverage analysis between the sexes. Despite the notable heteromorphy and distinctiveness in heterochromatin content, the Z and W sex chromosomes share approximately 90% of their gene content. This finding demonstrates poor correlation of the degree of differentiation of sex chromosomes at the cytogenetic and gene level. The test of homology based on the comparison of gene copy number variation revealed that female heterogamety with differentiated sex chromosomes remained stable in the genus Furcifer for at least 20 million years. These chameleons co-opted for the role of sex chromosomes the same genomic region as viviparous mammals, lacertids and geckos of the genus Paroedura, which makes these groups excellent model for studies of convergent and divergent evolution of sex chromosomes.

The genome of the invasive and broadly polyphagous Diaprepes root weevil, Diaprepes abbreviatus (Coleoptera), reveals an arsenal of putative polysaccharide-degrading enzymes

The Diaprepes root weevil (DRW), Diaprepes abbreviatus, is a broadly polyphagous invasive pest of agriculture in the southern United States and the Caribbean. Its genome was sequenced, assembled, and annotated to study genomic correlates of specialized plant-feeding and invasiveness and to facilitate the development of new methods for DRW control. The 1.69 Gb D. abbreviatus genome assembly was distributed across 653 contigs, with an N50 of 7.8 Mb and the largest contig of 62 Mb. Most of the genome was comprised of repetitive sequences, with 66.17% in transposable elements, 5.75% in macrosatellites, and 2.06% in microsatellites. Most expected orthologous genes were present and fully assembled, with 99.5% of BUSCO genes present and 1.5% duplicated. One hundred and nine contigs (27.19 Mb) were identified as putative fragments of the X and Y sex chromosomes, and homology assessment with other beetle X chromosomes indicated a possible sex chromosome turnover event. Genome annotation identified 18,412 genes, including 43 putative horizontally transferred (HT) loci. Notably, 258 genes were identified from gene families known to encode plant cell wall degrading enzymes and invertases, including carbohydrate esterases, polysaccharide lyases, and glycoside hydrolases (GH). GH genes were unusually numerous, with 239 putative genes representing 19 GH families. Interestingly, several other beetle species with large numbers of GH genes are (like D. abbreviatus) successful invasive pests of agriculture or forestry.

It's a trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex-determining gene in a lizard (Anguimorpha; Shinisauridae)

Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally to genetically influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades are not conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 with a novel ZW system on LG3. The linkage group, LG3, corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1 Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2, critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex-determining gene for Shinisaurus.

Histological Analysis of Gonadal Ridge Development and Sex Differentiation of Gonads in Three Gecko Species

Reptiles constitute a highly diverse group of vertebrates, with their evolutionary lineages having diverged relatively early. The types of sex determination exemplify the diversity of reptiles; however, there are limited data regarding the gonadal development in squamate reptiles. Geckos constitute a group that is increasingly used in research and that serves as a potential reptilian model organism. The aim of this study was to trace the changes in the structure of developing gonads in the embryos of three gecko species: the crested gecko, leopard gecko, and mourning gecko. These species represent different families of the Gekkota infraorder and exhibit different types of sex determination. Gonadal development was examined from the formation of the earliest gonadal ridges through the development of undifferentiated gonadal structures, sex differentiation of gonads, and the formation of testicular and ovarian structures. The study showed that the gonadal primordia of these three gecko species formed on the most dorsally located surface of the dorsal mesentery, and both the coelomic epithelium and the nephric mesenchyme contributed to their development. As in other reptile species, primordial germ cells settled in the gonadal ridges, and the undifferentiated gonad was composed of a cortex and a medulla. Ovarian differentiation started with the thickening of the gonadal cortex and proliferation of germ cells in this region. A characteristic feature of the developing gecko ovaries was the thickened crescent-shaped cortex on the medial and ventral surfaces of the ovaries. The ovarian medulla also grew and exhibited diverse tendencies to form cords. In the leopard gecko, advanced cord-like structures with lumens were observed in the ovaries, which were not seen in the crested gecko. Testicular differentiation was characterized by cortical thinning and the disappearance of germ cells in this region. In the medulla, the development of distinct cords with early lumen formation was noted. A characteristic feature of embryonic gonads was their growth in a horizontal plane. In this study, gonadal development was characterized by several features that are shared by geckos and other reptiles, along with features that are specific only to geckos.

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus)

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

It's a Trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex determining gene in a lizard (Anguimorpha; Shinisauridae)

Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally- to genetically-influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades aren't conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including: monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 chromosome with a novel ZW system on LG3. The linkage group LG3 corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2-critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex determining gene for Shinisaurus.

Chromosome evolution in Iberolacerta, a genus that deviates from the standard karyotype formula of Lacertidae

This paper describes the preparation of flow-sorted chromosome paints from the Iberian Rock lizard Iberolacerta monticola, exemplifying their subsequent use in cross-species comparisons of chromosome painting. We carried out comparative analyses of chromosome evolution in the congeneric species I. galani and I. bonnali, as well as in two other species of Lacertini (Lacerta schreiberi and Timon lepidus) whose sex chromosomes were also studied through comparative genomic hybridization. Most species of Lacertini possess a diplod number of 2n = 38, with 36 acrocentric macrochromosomes and 2 microchromosomes. However, the nine species included in the genus Iberolacerta do not possess microchromosomes. Furthermore, very conspicuous differences from the standard Lacertini karyotype were observed in the three Pyrenean species of this genus, which included several biarmed metacentrics and a Z1Z2W multiple sex-chromosome system. With the possible exception of L. schreiberi, all the species of the family Lacertidae described to date appear to share homologous Z chromosomes, which date back to the last common ancestor of the whole group. We provide conclusive evidence that L. schreiberi should no longer be considered an exception to this rule, and demonstrate that the loss of microchromosomes in Iberolacerta was produced by their fusion to a middle-sized chromosome. Furthermore, we show that the multiple sex-chromosome system of the Pyrenean species of Iberolacerta originated from the fusion of the ancestral W chromosome with one of the shortest autosomes, and provide additional evidence of the fast evolution of DNA sequences linked to the W chromosome in Lacertini.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Robertsonian fusion triggers recombination suppression on sex chromosomes in Coleonyx geckos

The classical hypothesis proposes that the lack of recombination on sex chromosomes arises due to selection for linkage between a sex-determining locus and sexually antagonistic loci, primarily facilitated by inversions. However, cessation of recombination on sex chromosomes could be attributed also to neutral processes, connected with other chromosome rearrangements or can reflect sex-specific recombination patterns existing already before sex chromosome differentiation. Three Coleonyx gecko species share a complex X1X1X2X2/X1X2Y system of sex chromosomes evolved via a fusion of the Y chromosome with an autosome. We analyzed synaptonemal complexes and sequenced flow-sorted sex chromosomes to investigate the effect of chromosomal rearrangement on recombination and differentiation of these sex chromosomes. The gecko sex chromosomes evolved from syntenic regions that were also co-opted also for sex chromosomes in other reptiles. We showed that in male geckos, recombination is less prevalent in the proximal regions of chromosomes and is even further drastically reduced around the centromere of the neo-Y chromosome. We highlight that pre-existing recombination patterns and Robertsonian fusions can be responsible for the cessation of recombination on sex chromosomes and that such processes can be largely neutral.

The revised reference genome of the leopard gecko (Eublepharis macularius) provides insight into the considerations of genome phasing and assembly

Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000.

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.

Sex-Linked Loci on the W Chromosome in the Multi-Ocellated Racerunner (Eremias multiocellata) Confirm Genetic Sex-Determination Stability in Lacertid Lizards

The multi-ocellated racerunner, Eremias multiocellata, was considered to have temperature-dependent sex determination (TSD), as its sex ratio can be influenced at different temperatures. However, such an observation contrasts with recent findings that suggest TSD is less common than previously thought. Here, a genotyping-by-sequencing (GBS) approach was employed to identify sex-linked markers in the E. multiocellata, for which the mechanism choice of TSD or GSD is still controversial. We preliminarily identified 119 sex-linked markers based on sex-associated sex-specific sequences, 97% of which indicated female heterogamety. After eliminating the false positives, 38 sex-linked markers were recognized, all of which showed the ZW/ZZ system. Then, eight of the novel markers were verified by PCR amplification from 15 populations of E. multiocellata, which support the GSD in E. multiocellata without geographic variation. To test the conservation of sex chromosome in Eremias, the eight markers were further cross-tested by PCR amplification in 10 individuals of the Mongolian racerunner (Eremias argus), two of which exhibited cross-utility. The novel sex-linked markers could be mapped on the W chromosome of the sand lizard (Lacerta agilis). Our finding that the sex-linked markers are shared in closely related species, along with a conserved synteny of the W chromosome, further supports the homology and conservation of sex chromosomes in the lacertid lizards.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Divergent evolution of genetic sex determination mechanisms along environmental gradients

Sex determination (SD) is a crucial developmental process, but its molecular underpinnings are very diverse, both between and within species. SD mechanisms have traditionally been categorized as either genetic (GSD) or environmental (ESD), depending on the type of cue that triggers sexual differentiation. However, mixed systems, with both genetic and environmental components, are more prevalent than previously thought. Here, we show theoretically that environmental effects on expression levels of genes within SD regulatory mechanisms can easily trigger within-species evolutionary divergence of SD mechanisms. This may lead to the stable coexistence of multiple SD mechanisms and to spatial variation in the occurrence of different SD mechanisms along environmental gradients. We applied the model to the SD system of the housefly, a global species with world-wide latitudinal clines in the frequencies of different SD systems, and found that it correctly predicted these clines if specific genes in the housefly SD system were assumed to have temperature-dependent expression levels. We conclude that environmental sensitivity of gene regulatory networks may play an important role in diversification of SD mechanisms.

Massive expansion of sex-specific SNPs, transposon-related elements, and neocentromere formation shape the young W-chromosome from the mosquitofish Gambusia affinis

BACKGROUND

The Western mosquitofish, Gambusia affinis, is a model for sex chromosome organization and evolution of female heterogamety. We previously identified a G. affinis female-specific marker, orthologous to the aminomethyl transferase (amt) gene of the related platyfish (Xiphophorus maculatus). Here, we have analyzed the structure and differentiation of the G. affinis W-chromosome, using a cytogenomics and bioinformatics approach.

RESULTS

The long arm of the G. affinis W-chromosome (Wq) is highly enriched in dispersed repetitive sequences, but neither heterochromatic nor epigenetically silenced by hypermethylation. In line with this, Wq sequences are highly transcribed, including an active nucleolus organizing region (NOR). Female-specific SNPs and evolutionary young transposable elements were highly enriched and dispersed along the W-chromosome long arm, suggesting constrained recombination. Wq copy number expanded elements also include female-specific transcribed sequences from the amt locus with homology to TE. Collectively, the G. affinis W-chromosome is actively differentiating by sex-specific copy number expansion of transcribed TE-related elements, but not (yet) by extensive sequence divergence or gene decay.

CONCLUSIONS

The G. affinis W-chromosome exhibits characteristic genomic properties of an evolutionary young sex chromosome. Strikingly, the observed sex-specific changes in the genomic landscape are confined to the W long arm, which is separated from the rest of the W-chromosome by a neocentromere acquired during sex chromosome evolution and may thus have become functionally insulated. In contrast, W short arm sequences were apparently shielded from repeat-driven differentiation, retained Z-chromosome like genomic features, and may have preserved pseudo-autosomal properties.

Incomplete dosage balance and dosage compensation in the ZZ/ZW Gila monster (Heloderma suspectum) revealed by de novo genome assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z-chromosome scaffolds. These putative Z-chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

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.

The male-heterogametic sex determination system on chromosome 15 of Salix triandra and Salix arbutifolia reveals ancestral male heterogamety and subsequent turnover events in the genus Salix

Dioecious Salix evolved more than 45 million years ago, but have homomorphic sex chromosomes, suggesting that turnover event(s) prevented major differentiation. Sex chromosome turnover events have been inferred in the sister genus Populus. The genus Salix includes two main clades, Salix and Vetrix, with several previously studied Vetrix clade species having female-heterogametic (ZW) or male-heterogametic (XY) sex-determining systems (SDSs) on chromosome 15, while three Salix clade species have XY SDSs on chromosome 7. We here studied two basal taxa of the Vetrix clade, S. arbutifolia and S. triandra using S. purpurea as the reference genome. Analyses of whole genome resequencing data for genome-wide associations (GWAS) with the sexes and genetic differentiation between the sexes (FST values) showed that both species have male heterogamety with a sex-determining locus on chromosome 15, suggesting an early turnover event within the Vetrix clade, perhaps promoted by sexually antagonistic or (and) sex-ratio selection. Changepoint analysis based on FST values identified small sex-linked regions of ~3.33 Mb and ~2.80 Mb in S. arbutifolia and S. triandra, respectively. The SDS of S. arbutifolia was consistent with recent results that used its own genome as reference. Ancestral state reconstruction of SDS suggests that at least two turnover events occurred in Salix.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Madagascar Leaf-Tail Geckos (Uroplatus spp.) Share Independently Evolved Differentiated ZZ/ZW Sex Chromosomes

Geckos are an excellent group to study the evolution of sex determination, as they possess a remarkable variability ranging from a complete absence of sex chromosomes to highly differentiated sex chromosomes. We explored sex determination in the Madagascar leaf-tail geckos of the genus Uroplatus. The cytogenetic analyses revealed highly heterochromatic W chromosomes in all three examined species (Uroplatus henkeli, U. alluaudi, U. sikorae). The comparative gene coverage analysis between sexes in U. henkeli uncovered an extensive Z-specific region, with a gene content shared with the chicken chromosomes 8, 20, 26 and 28. The genomic region homologous to chicken chromosome 28 has been independently co-opted for the role of sex chromosomes in several vertebrate lineages, including monitors, beaded lizards and monotremes, perhaps because it contains the amh gene, whose homologs were repeatedly recruited as a sex-determining locus. We demonstrate that all tested species of leaf-tail geckos share homologous sex chromosomes despite the differences in shape and size of their W chromosomes, which are not homologous to the sex chromosomes of other closely related genera. The rather old (at least 40 million years), highly differentiated sex chromosomes of Uroplatus geckos can serve as a great system to study the convergence of sex chromosomes evolved from the same genomic region.

ZZ/ZW Sex Chromosomes in the Madagascar Girdled Lizard, Zonosaurus madagascariensis (Squamata: Gerrhosauridae)

Scincoidea, the reptilian clade that includes girdled lizards, night lizards, plated lizards and skinks, are considered as a lineage with diversity in sex-determining systems. Recently, the hypothesis on the variability in sex determination in skinks and even more the absence of sex chromosomes in some of them has been rivalling. Homologous, evolutionary stable XX/XY sex chromosomes were documented to be widespread across skinks. However, sex determination in the other scincoidean families is highly understudied. ZZ/ZW sex chromosomes have been identified only in night lizards and a single species of plated lizards. It seems that although there are different sex chromosome systems among scincoidean lineages, they share one common trait: they are mostly poorly differentiated and often undetectable by cytogenetic methods. Here, we report one of the exceptions, demonstrating for the first time ZZ/ZW sex chromosomes in the plated lizard Zonosaurus madagascariensis. Its sex chromosomes are morphologically similar, but the W is clearly detectable by comparative genomic hybridization (CGH), suggesting that the Z and W chromosomes highly differ in sequence content. Our findings confirm the presence of female heterogamety in plated lizards and provides novel insights to expand our understanding of sex chromosome evolution in scincoidean lizards.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

New insights into Xenopus sex chromosome genomics from the Marsabit clawed frog X. borealis

In many groups, sex chromosomes change frequently but the drivers of their rapid evolution are varied and often poorly characterized. With an aim of further understanding sex chromosome turnover, we investigated the polymorphic sex chromosomes of the Marsabit clawed frog, Xenopus borealis, using genomic data and a new chromosome-scale genome assembly. We confirmed previous findings that 54.1 Mb of chromosome 8L is sex-linked in animals from east Kenya and a laboratory strain, but most (or all) of this region is not sex-linked in natural populations from west Kenya. Previous work suggests possible degeneration of the Z chromosomes in the east population because many sex-linked transcripts of this female heterogametic population have female-biased expression, and we therefore expected this chromosome to not be present in the west population. In contrast, our simulations support a model where most or all of the sex-linked portion of the Z chromosome from the east acquired autosomal segregation in the west, and where much genetic variation specific to the large sex-linked portion of the W chromosome from the east is not present in the west. These recent changes are consistent with the hot-potato model, wherein sex chromosome turnover is favoured by natural selection if it purges a (minimally) degenerate sex-specific sex chromosome, but counterintuitively suggest natural selection failed to purge a Z chromosome that has signs of more advanced and possibly more ancient regulatory degeneration. These findings highlight complex evolutionary dynamics of young, rapidly evolving Xenopus sex chromosomes and set the stage for mechanistic work aimed at pinpointing additional sex-determining genes in this group.

Sex Determination and Ovarian Development in Reptiles and Amphibians: From Genetic Pathways to Environmental Influences

BACKGROUND

Reptiles and amphibians provide untapped potential for discovering how a diversity of genetic pathways and environmental conditions are incorporated into developmental processes that can lead to similar functional outcomes. These groups display a multitude of reproductive strategies, and whereas many attributes are conserved within groups and even across vertebrates, several aspects of sexual development show considerable variation.

SUMMARY

In this review, we focus our attention on the development of the reptilian and amphibian ovary. First, we review and describe the events leading to ovarian development, including sex determination and ovarian maturation, through a comparative lens. We then describe how these events are influenced by environmental factors, focusing on temperature and exposure to anthropogenic chemicals. Lastly, we identify critical knowledge gaps and future research directions that will be crucial to moving forward in our understanding of ovarian development and the influences of the environment in reptiles and amphibians.

KEY MESSAGES

Reptiles and amphibians provide excellent models for understanding the diversity of sex determination strategies and reproductive development. However, a greater understanding of the basic biology of these systems is necessary for deciphering the adaptive and potentially disruptive implications of embryo-by-environment interactions in a rapidly changing world.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops

In sharp contrast to birds and mammals, in numerous cold-blooded vertebrates, sex chromosomes have been described as homomorphic. This sex chromosome homomorphy has been suggested to result from the high turnovers often observed across deeply diverged clades. However, little is known about the tempo and mode of sex chromosome evolution among the most closely related species. Here, we examined the evolution of sex chromosome among nine species of the torrent frog genus Amolops. We analyzed male and female GBS and RAD-seq from 182 individuals and performed PCR verification for 176 individuals. We identified signatures of sex chromosomes involving two pairs of chromosomes. We found that sex-chromosome homomorphy results from both turnover and X–Y recombination in the Amolops species, which simultaneously exhibits heterogeneous evolution on homologous and non-homologous sex chromosomes. A low turnover rate of non-homologous sex chromosomes exists in these torrent frogs. The ongoing X–Y recombination in homologous sex chromosomes will act as an indispensable force in preventing sex chromosomes from differentiating.

Highly rapid and diverse sex chromosome evolution in the Odorrana frog species complex

Sex chromosomes in poikilothermal vertebrates are characterized by rapid and diverse evolution at the species or population level. Our previous study revealed that the Taiwanese frog Odorrana swinhoana (2n = 26) has a unique system of multiple sex chromosomes created by three sequential translocations among chromosomes 1, 3, and 7. To reveal the evolutionary history of sex chromosomes in the Odorrana species complex, we first identified the original, homomorphic sex chromosomes, prior to the occurrence of translocations, in the ancestral-type population of O. swinhoana. Then, we extended the investigation to a closely related Japanese species, Odorrana utsunomiyaorum, which is distributed on two small islands. We used a high-throughput nuclear genomic approach to analyze single-nucleotide polymorphisms and identify the sex-linked markers. Those isolated from the O. swinhoana ancestral-type population were found to be aligned to chromosome 1 and showed male heterogamety. In contrast, almost all the sex-linked markers isolated from O. utsunomiyaorum were heterozygous in females and homozygous in males and were aligned to chromosome 9. Morphologically, we confirmed chromosome 9 to be heteromorphic in females, showing a ZZ-ZW sex determination system, in which the W chromosomes were heterochromatinized in a stripe pattern along the chromosome axis. These results indicated that after divergence of the two species, the ancestral homomorphic sex chromosome 1 underwent highly rapid and diverse evolution, i.e., sequential translocations with two autosomes in O. swinhoana, and turnover to chromosome 9 in O. utsunomiyaorum, with a transition from XY to ZW heterogamety and change to heteromorphy.

Genetic determination and JARID2 over-expression in a thermal incubation experiment in Casque-Headed Lizard

Non-avian reptiles, unlike mammals and birds, have undergone numerous sex determination changes. Casque-Headed Lizards have replaced the ancestral XY system shared across pleurodonts with a new pair of XY chromosomes. However, the evolutionary forces that triggered this transition have remained unclear. An interesting hypothesis suggests that species with intermediate states, with sex chromosomes but also thermal-induced sex reversal at specific incubation temperatures, could be more susceptible to sex determination turnovers. We contrasted genotypic data (presence/absence of the Y chromosome) against the histology of gonads of embryos from stages 35–37 incubated at various temperatures, including typical male-producing (26°C) and female-producing (32°C) temperatures. Our work apparently reports for the first time the histology of gonads, including morphological changes, from stages 35–37 of development in the family Corytophanidae. We also observed that all embryos developed hemipenes, suggesting sex-linked developmental heterochrony. We observed perfect concordance between genotype and phenotype at all temperatures. However, analysis of transcriptomic data from embryos incubated at 26°C and 32°C identified transcript variants of the chromatin modifiers JARID2 and KDM6B that have been linked to temperature-dependent sex determination in other reptiles. Our work tested the validity of a mixed sex determination system in the family Corytophanidae. We found that XY chromosomes are dominant; however, our work supports the hypothesis of a conserved transcriptional response to incubation temperatures across non-avian reptiles that could be a reminiscence of an ancestral sex determination system.

Cytogenetic Analysis of the Members of the Snake Genera Cylindrophis, Eryx, Python, and Tropidophis

The recent discovery of two independently evolved XX/XY sex determination systems in the snake genera Python and Boa sparked a new drive to study the evolution of sex chromosomes in poorly studied lineages of snakes, where female heterogamety was previously assumed. Therefore, we examined seven species from the genera Eryx, Cylindrophis, Python, and Tropidophis by conventional and molecular cytogenetic methods. Despite the fact that these species have similar karyotypes in terms of chromosome number and morphology, we detected variability in the distribution of heterochromatin, telomeric repeats, and rDNA loci. Heterochromatic blocks were mainly detected in the centromeric regions in all species, although accumulations were detected in pericentromeric and telomeric regions in a few macrochromosomes in several of the studied species. All species show the expected topology of telomeric repeats at the edge of all chromosomes, with the exception of Eryx muelleri, where additional accumulations were detected in the centromeres of three pairs of macrochromosomes. The rDNA loci accumulate in one pair of microchromosomes in all Eryx species and in Cylindrophis ruffus, in one macrochromosome pair in Tropidophis melanurus and in two pairs of microchromosomes in Python regius. Sex-specific differences were not detected, suggesting that these species likely have homomorphic, poorly differentiated sex chromosomes.

Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

Identification of a cis-sex chromosome transition in banded geckos (Coleonyx, Eublepharidae, Gekkota)

Sex-determination systems are highly variable amongst vertebrate groups, and the prevalence of genomic data has greatly expanded our knowledge of how diverse some groups truly are. Gecko lizards are known to possess a variety of sex-determination systems, and each new study increases our knowledge of this diversity. Here, we used RADseq to identify male-specific markers in the banded gecko Coleonyx brevis, indicating this species has a XX/XY sex-determination system. Furthermore, we show that these sex-linked regions are not homologous to the XX/XY sex chromosomes of two related Coleonyx species, C. elegans and C. mitratus, suggesting that a cis-sex chromosome turnover-a change in sex chromosomes without a concomitant change in heterogamety-has occurred within the genus. These findings demonstrate the utility of genome-scale data to uncover novel sex chromosomes and further highlight the diversity of gecko sex chromosomes.

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Sex-specific markers undetected in green sunfish Lepomis cyanellus using restriction-site associated DNA sequencing

We used restriction-site associated DNA sequencing for SNP discovery and genotyping of known-sex green sunfish Lepomis cyanellus DNA samples to search for sex-diagnostic single nucleotide polymorphisms (SNPs) and restriction-site associated sequences present in one sex and absent in the other. The bioinformatic analyses discovered candidate SNPs and sex-specific restriction-site associated sequences that fit patterns of male or female heterogametic sex determination systems. However, when primers were developed and tested, no candidates reliably identified phenotypic sex. The top performing SNP candidate (ZW_218) correlated with phenotypic sex 63.0% of the time and the presence-absence loci universally amplified in both sexes. We recommend further investigations that interrogate a larger fraction of the L. cyanellus genome. Additionally, studies on the effect of temperature and rearing density on sex determination, as well as breeding of sex-reversed individuals, could provide more insights into the sex determination system of L. cyanellus.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Switches, stability and reversals in the evolutionary history of sexual systems in fish

Sexual systems are highly diverse and have profound consequences for population dynamics and resilience. Yet, little is known about how they evolved. Using phylogenetic Bayesian modelling and a sample of 4614 species, we show that gonochorism is the likely ancestral condition in teleost fish. While all hermaphroditic forms revert quickly to gonochorism, protogyny and simultaneous hermaphroditism are evolutionarily more stable than protandry. In line with theoretical expectations, simultaneous hermaphroditism does not evolve directly from gonochorism but can evolve slowly from sequential hermaphroditism, particularly protandry. We find support for the predictions from life history theory that protogynous, but not protandrous, species live longer than gonochoristic species and invest the least in male gonad mass. The distribution of teleosts' sexual systems on the tree of life does not seem to reflect just adaptive predictions, suggesting that adaptations alone may not fully explain why some sexual forms evolve in some taxa but not others (Williams' paradox). We propose that future studies should incorporate mating systems, spawning behaviours, and the diversity of sex determining mechanisms. Some of the latter might constrain the evolution of hermaphroditism, while the non-duality of the embryological origin of teleost gonads might explain why protogyny predominates over protandry in teleosts.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Chromosome-level genome assembly of Asian yellow pond turtle (Mauremys mutica) with temperature-dependent sex determination system

Knowledge of sex determination has important implications in physiology, ecology and genetics, but the evolutionary mechanisms of sex determination systems in turtles have not been fully elucidated, due to a lack of reference genomes. Here, we generate a high-quality genome assembly of Asian yellow pond turtle (Mauremys mutica) using continuous long-read (PacBio platform), Illumina, and high-throughput chromatin conformation capture (Hi-C) technologies. The M. mutica haplotype has a genome size of 2.23 Gb with a contig N50 of 8.53 Mb and scaffold N50 of 141.98 Mb. 99.98% sequences of the total assembly are anchored to 26 pseudochromosomes. Comparative genomics analysis indicated that the lizard-snake-tuatara clade diverged from the bird-crocodilian-turtle clade at approximately 267.0-312.3 Mya. Intriguingly, positive selected genes are mostly enriched in the calcium signaling pathway and neuroactive ligand-receptor interaction, which are involved in the process of temperature-dependent sex determination. These findings provide important evolutionary insights into temperature-dependent sex determination system.

Asymmetrical sex reversal: Does the type of heterogamety predict propensity for sex reversal?

Sex reversal, a mismatch between phenotypic and genetic sex, can be induced by chemical and thermal insults in ectotherms. Therefore, climate change and environmental pollution may increase sex-reversal frequency in wild populations, with wide-ranging implications for sex ratios, population dynamics, and the evolution of sex determination. We propose that reconsidering the half-century old theory "Witschi's rule" should facilitate understanding the differences between species in sex-reversal propensity and thereby predicting their vulnerability to anthropogenic environmental change. The idea is that sex reversal should be asymmetrical: more likely to occur in the homogametic sex, assuming that sex-reversed heterogametic individuals would produce new genotypes with reduced fitness. A review of the existing evidence shows that while sex reversal can be induced in both homogametic and heterogametic individuals, the latter seem to require stronger stimuli in several cases. We provide guidelines for future studies on sex reversal to facilitate data comparability and reliability.

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.