The Diversity and Evolution of Sex Chromosomes in Frogs

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Evolution and regulation of animal sex chromosomes

Animal sex chromosomes typically carry the upstream sex-determining gene that triggers testis or ovary development and, in some species, are regulated by global dosage compensation in response to functional decay of the Y chromosome. Despite the importance of these pathways, they exhibit striking differences across species, raising fundamental questions regarding the mechanisms underlying their evolutionary turnover. Recent studies of non-model organisms, including insects, reptiles and teleosts, have yielded a broad view of the diversity of sex chromosomes that challenges established theories. Moreover, continued studies in model organisms with recently developed technologies have characterized the dynamics of sex determination and dosage compensation in three-dimensional nuclear space and at single-cell resolution. Here, we synthesize recent insights into sex chromosomes from a variety of species to review their evolutionary dynamics with respect to the canonical model, as well as their diverse mechanisms of regulation.

Rapid Sex Chromosome Turnover in African Clawed Frogs (Xenopus) and the Origins of New Sex Chromosomes

Sex chromosomes of some closely related species are not homologous, and sex chromosome turnover is often attributed to mechanisms that involve linkage to or recombination arrest around sex-determining loci. We examined sex chromosome turnover and recombination landscapes in African clawed frogs (genus Xenopus) with reduced representation genome sequences from 929 individuals from 19 species. We recovered extensive variation in sex chromosomes, including at least eight nonhomologous sex-associated regions-five newly reported here, with most maintaining female heterogamety, but two independent origins of Y chromosomes. Seven of these regions are found in allopolyploid species in the subgenus Xenopus, and all of these reside in one of their two subgenomes, which highlights functional asymmetry between subgenomes. In three species with chromosome-scale genome assemblies (Xenopus borealis, Xenopus laevis, and Xenopus tropicalis), sex-specific recombination landscapes have similar patterns of sex differences in rates and locations of recombination. Across these Xenopus species, sex-associated regions are significantly nearer chromosome ends than expected by chance, even though this is where the ancestral recombination rate is highest in both sexes before the regions became sex associated. As well, expansions of sex-associated recombination arrest occurred multiple times. New information on sex linkage along with among-species variation in female specificity of the sex-determining gene dm-w argues against a "jumping gene" model, where dm-w moves around the genome. The diversity of sex chromosomes in Xenopus raises questions about the roles of natural and sexual selection, polyploidy, the recombination landscape, and neutral processes in driving sex chromosome turnover in animal groups with mostly heterogametic females.

Yes, polygenic sex determination is a thing!

The process of sexual development in animals is modulated by a variety of mechanisms. Some species respond to environmental cues, while, in others, sex determination is thought to be controlled by a single 'master regulator' gene. However, many animals respond to a combination of environmental cues (e.g., temperature) and genetic factors (e.g., sex chromosomes). Even among species in which genetic factors predominate, there is a continuum between monofactorial and polygenic systems. The perception that polygenic systems are rare may result from experiments that lack the statistical power to detect multiple loci. Intellectual biases against the existence of polygenic sex determination (PSD) may further arise from misconceptions about the regulation of developmental processes and a misreading of theoretical results on the stability of polygenic systems of sex determination.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Genome Survey of Male Rana dybowskii to Further Understand the Sex Determination Mechanism

Rana dybowskii is one of the important aquaculture species in Northeast China. The fallopian tubes of female R. dybowskii are used to prepare oviductus ranae (an important traditional Chinese medicine). Therefore, R. dybowskii females have higher economical value than males. An increasing female R. dybowskii population can increase the benefits from R. dybowskii culture. However, the genome of amphibians is complex, making it difficult to investigate their sex determination mechanism. In this study, we analyzed the genome of male R. dybowskii using next-generation sequencing technology. A total of 200,046,452,400 bp of clean data were obtained, and the K-mer analysis indicated that the depth was 50×. The genome size of R. dybowskii was approximately 3585.05 M, with a heterozygosity rate, repeat sequence ratio, and genome GC content of 1.15%, 68.96%, and approximately 43.0%, respectively. In total, 270,785 contigs and 498 scaffolds were generated. The size of the contigs and scaffolds was 3,748,543,415 and 3,765,862,278 bp, respectively, with the N50 length of 31,988 and 336,385,783. The longest contig and scaffold were of the size 137,967,485 and 1,808,367,828 bp, respectively. The number of contigs and scaffolds > 10K nt was 99,620 and 451, respectively. Through annotation, 40,913 genes were obtained, including 156,609 CDS (i.e., 3.83 CDS per gene). Sequence alignment was performed with the assembled scaffolding genome in this study. Two and one fragment had high homology with two male-specific DNA molecular markers of R. dybowskii discovered previously (namely, MSM-222 and MSM-261, respectively). In addition, the Dmrt1 gene of R. dybowskii was obtained with a length of 18,893 bp by comparison and splicing. The forward primers amplifying MSM-222 and MSM-261 were located at 322-343 and 14,501-14,526 bp of Dmrt1, respectively. However, sequence alignment revealed that MSM-222 and MSM-261 were not located on Dmrt1, and only some homologous parts were observed. This indicated that in addition to Dmrt1, other important genes may play a crucial role in the sex determination mechanism of R. dybowskii. Our study provided a foundation for the subsequent high-quality genome construction and provided important genomic resources for future studies on R. dybowskii.

Satellite DNAs and the evolution of the multiple X1X2Y sex chromosomes in the wolf fish Hoplias malabaricus (Teleostei; Characiformes)

Multiple sex chromosomes usually arise from chromosomal rearrangements which involve ancestral sex chromosomes. There is a fundamental condition to be met for their long-term fixation: the meiosis must function, leading to the stability of the emerged system, mainly concerning the segregation of the sex multivalent. Here, we sought to analyze the degree of differentiation and meiotic pairing properties in the selected fish multiple sex chromosome system present in the wolf-fish Hoplias malabaricus (HMA). This species complex encompasses seven known karyotype forms (karyomorphs) where the karyomorph C (HMA-C) exhibits a nascent XY sex chromosomes from which the multiple X1X2Y system evolved in karyomorph HMA-D via a Y-autosome fusion. We combined genomic and cytogenetic approaches to analyze the satellite DNA (satDNA) content in the genome of HMA-D karyomorph and to investigate its potential contribution to X1X2Y sex chromosome differentiation. We revealed 56 satDNA monomers of which the majority was AT-rich and with repeat units longer than 100 bp. Seven out of 18 satDNA families chosen for chromosomal mapping by fluorescence in situ hybridization (FISH) formed detectable accumulation in at least one of the three sex chromosomes (X1, X2 and neo-Y). Nine satDNA monomers showed only two hybridization signals limited to HMA-D autosomes, and the two remaining ones provided no visible FISH signals. Out of seven satDNAs located on the HMA-D sex chromosomes, five mapped also to XY chromosomes of HMA-C. We showed that after the autosome-Y fusion event, the neo-Y chromosome has not substantially accumulated or eliminated satDNA sequences except for minor changes in the centromere-proximal region. Finally, based on the obtained FISHpatterns, we speculate on the possible contribution of satDNA to sex trivalent pairing and segregation.

The importance of reproductive isolation in driving diversification and speciation within Peruvian mimetic poison frogs (Dendrobatidae)

To explain how populations with distinct warning signals coexist in close parapatry, we experimentally assessed intrinsic mechanisms acting as reproductive barriers within three poison-frog species from the Peruvian Amazon belonging to a Müllerian mimetic ring (Ranitomeya variabilis, Ranitomeya imitator and Ranitomeya fantastica). We tested the role of prezygotic and postzygotic isolation barriers between phenotypically different ecotypes of each species, using no-choice mating experiments and offspring survival analysis. Our results show that prezygotic mating preference did not occur except for one specific ecotype of R. imitator, and that all three species were able to produce viable inter-population F1 hybrids. However, while R. variabilis and R. imitator hybrids were able to produce viable F2 generations, we found that for R. fantastica, every F1 hybrid males were sterile while females remained fertile. This unexpected result, echoing with Haldane's rule of speciation, validated phylogenetic studies which tentatively diagnose these populations of R. fantastica as two different species. Our work suggests that postzygotic genetic barriers likely participate in the extraordinary phenotypic diversity observed within Müllerian mimetic Ranitomeya populations, by maintaining species boundaries.

Comparative molecular and conventional cytogenetic analyses of three species of Rhinella (Anura; Bufonidae)

The genus Rhinella corresponds to a group of anurans characterized by numerous taxonomic and systemic challenges, leading to their organization into species complexes. Cytogenetic data for this genus thus far are limited to the diploid number and chromosome morphology, which remain highly conserved among the species. In this study, we analyse the karyotypes of three species of the genus Rhinella (Rhinella granulosa, Rhinella margaritifera, and Rhinella marina) using both classical (conventional staining and C-banding) and molecular (FISH-fluorescence in situ hybridization with 18S rDNA, telomeric sequences, and microsatellite probes) cytogenetic approaches. The aim of this study is to provide data that can reveal variations in the distribution of repetitive sequences that can contribute to understanding karyotypic diversification in these species. The results revealed a conserved karyotype across the species, with 2n = 22 and FN = 44, with metacentric and submetacentric chromosomes. C-banding revealed heterochromatic blocks in the pericentromeric region for all species, with a proximal block on the long arms of pairs 3 and 6 in R. marina and on the short arms of pairs 4 and 6 in R. margaritifera. Additionally, 18S rDNA probes hybridized to pair 5 in R. granulosa, to pair 7 in R. marina, and to pair 10 in R. margaritifera. Telomeric sequence probes displayed signals exclusively in the distal region of the chromosomes, while microsatellite DNA probes showed species-specific patterns. These findings indicate that despite a conserved karyotypical macrostructure, chromosomal differences exist among the species due to the accumulation of repetitive sequences. This variation may be attributed to chromosome rearrangements or differential accumulation of these sequences, highlighting the dynamic role of repetitive sequences in the chromosomal evolution of Rhinella species. Ultimately, this study emphasizes the importance of the role of repetitive DNAs in chromosomal rearrangements to elucidate the evolutionary mechanisms leading to independent diversification in the distinct phylogenetic groups of Rhinella.

Altering traits and fates of wild populations with Mendelian DNA sequence modifying Allele Sails

Population-scale genome modification can alter the composition or fate of wild populations. Synthetic gene drives provide one set of tools, but their use is complicated by scientific, regulatory, and social issues associated with transgene persistence and flow. Here we propose an alternative approach. An Allele Sail consists of a genome editor (the Wind) that introduces DNA sequence edits, and is inherited in a Mendelian fashion. Meanwhile, the edits (the Sail) experience an arithmetic, Super-Mendelian increase in frequency. We model this system and identify contexts in which a single, low frequency release of an editor brings edits to a very high frequency. We also identify conditions in which manipulation of sex determination can bring about population suppression. In regulatory frameworks that distinguish between transgenics (GMO) and their edited non-transgenic progeny (non-GMO) Allele Sails may prove useful since the spread and persistence of the GM component can be limited.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

Homology in Sex Determination in Two Distant Spiny Frogs, Nanorana quadranus and Quasipaa yei

Sex determination is remarkably diverse, with frequent transitions between sex chromosomes, in amphibians. Under these transitions, some chromosomes are more likely to be recurrently co-opted as sex chromosomes, as they are often observed across deeply divergent taxa. However, little is known about the pattern of sex chromosome evolution among closely related groups. Here, we examined sex chromosome and sex determination in two spiny frogs, Nanorana quadranus and Quasipaa yei. We conducted an analysis of genotyping-by-sequencing (GBS) data from a total of 34 individuals to identify sex-specific makers, with the results verified by PCR. The results suggest that chromosome 1 is a homologous sex chromosome with an XY pattern in both species. This chromosome has been evolutionarily conserved across these closely related groups within a period of time. The DMRT1 gene is proposed to be implicated in homology across two distantly related spiny frog species as a putative candidate sex-determining gene. Harboring the DMRT1 gene, chromosome 1 would have been independently co-opted for sex determination in deeply divergent groups of anurans.

Understanding the genetic sex-determining mechanism in Hyla eximia treefrog inferred from H-Y antigen

Genetic sex-determining mechanisms have been extensively elucidated in mammals; however, the sex chromosomes, sex-determining genes, and gene regulatory networks involved in sex differentiation remain poorly understood in amphibians. In this study, we investigated the sex-determining mechanism in the Hyla eximia treefrog based on karyotypic analysis and identification of H-Y antigen, a sex-linked peptide that is present in the gonads of the heterogametic sex (XY or ZW) in all vertebrates. Results show a diploid chromosome number 2n = 24 with homomorphic sex chromosomes. The heterogametic sex, ZW-female, were hypothesized based on H-Y antigen mRNA expression in female gonads (24,ZZ/24,ZW). The treefrog H-Y peptide exhibited a high percentage of identity with other vertebrate sequences uploaded to GenBank database. To obtain gene expression profiles, we also obtained the coding sequence of the housekeeping Actb gene. High H-Y antigen expression levels were further confirmed in ovaries using real-time polymerase chain reaction (RT-PCR) during non-breeding season, we noted a decrease in the expression of the H-Y antigen during breeding season. This study provides evidence that sex hormones might suppress H-Y antigen expression in the gonads of heterogametic females 24,ZW during the breeding season. These findings suggest that H-Y gene expression is a well-suited model for studying heterogametic sex by comparing the male and female gonads.

Y-chromosome Degeneration due to Speciation and Founder Effect

The Y chromosome in the XY sex-determination system is often shorter than its X counterpart, a condition attributed to degeneration after Y recombination ceases. Contrary to the traditional view of continuous, gradual degeneration, our study reveals stabilization within large mating populations. In these populations, we demonstrate that both mutant and active alleles on the Y chromosome can reach equilibrium through a mutation-selection balance. However, the emergence of a new species, particularly through the founder effect, can disrupt this equilibrium. Specifically, if the male founders of a new species carry only a mutant allele for a particular Y-linked gene, this allele becomes fixed, leading to the loss of the corresponding active gene on the Y chromosome. Our findings suggest that the rate of Y-chromosome degeneration may be linked to the frequency of speciation events associated with single-male founder events.

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.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Evolution of a ZW sex chromosome system in willows

Transitions in the heterogamety of sex chromosomes (e.g., XY to ZW or vice versa) fundamentally alter the genetic basis of sex determination, however the details of these changes have been studied in only a few cases. In an XY to ZW transition, the X is likely to give rise to the W because they both carry feminizing genes and the X is expected to harbour less genetic load than the Y. Here, using a new reference genome for Salix exigua, we trace the X, Y, Z, and W sex determination regions during the homologous transition from an XY system to a ZW system in willow (Salix). We show that both the W and the Z arose from the Y chromosome. We find that the new Z chromosome shares multiple homologous putative masculinizing factors with the ancestral Y, whereas the new W lost these masculinizing factors and gained feminizing factors. The origination of both the W and Z from the Y was permitted by an unexpectedly low genetic load on the Y and this indicates that the origins of sex chromosomes during homologous transitions may be more flexible than previously considered.

The first linkage map for Australo-Papuan Treefrogs (family: Pelodryadidae) reveals the sex-determination system of the Green-eyed Treefrog (Litoria serrata)

Amphibians represent a useful taxon to study the evolution of sex determination because of their highly variable sex-determination systems. However, the sex-determination system for many amphibian families remains unknown, in part because of a lack of genomic resources. Here, using an F1 family of Green-eyed Treefrogs (Litoria serrata), we produce the first genetic linkage map for any Australo-Papuan Treefrogs (family: Pelodryadidae). The resulting linkage map contains 8662 SNPs across 13 linkage groups. Using an independent set of sexed adults, we identify a small region in linkage group 6 matching an XY sex-determination system. These results suggest Litoria serrata possesses a male heterogametic system, with a candidate sex-determination locus on linkage group 6. Furthermore, this linkage map represents the first genomic resource for Australo-Papuan Treefrogs, an ecologically diverse family of over 220 species.

Homeology of sex chromosomes in Amazonian Harttia armored catfishes supports the X-fission hypothesis for the X1X2Y sex chromosome system origin

The Neotropical monophyletic catfish genus Harttia represents an excellent model to study karyotype and sex chromosome evolution in teleosts. Its species split into three phylogenetic clades distributed along the Brazilian territory and they differ widely in karyotype traits, including the presence of standard or multiple sex chromosome systems in some members. Here, we investigate the chromosomal rearrangements and associated synteny blocks involved in the origin of a multiple X1X2Y sex chromosome system present in three out of six sampled Amazonian-clade species. Using 5S and 18S ribosomal DNA fluorescence in situ hybridization and whole chromosome painting with probes corresponding to X1 and X2 chromosomes of X1X2Y system from H. punctata, we confirm previous assumptions that X1X2Y sex chromosome systems of H. punctata, H. duriventris and H. villasboas represent the same linkage groups which also form the putative XY sex chromosomes of H. rondoni. The shared homeology between X1X2Y sex chromosomes suggests they might have originated once in the common ancestor of these closely related species. A joint arrangement of mapped H. punctata X1 and X2 sex chromosomes in early diverging species of different Harttia clades suggests that the X1X2Y sex chromosome system may have formed through an X chromosome fission rather than previously proposed Y-autosome fusion.

HerpSexDet: the herpetological database of sex determination and sex reversal

Wildlife exhibits various sex-determination systems where sex chromosomes and environmental temperatures may both contribute to individual sexual development. The causes and consequences of this variability are important questions for evolutionary ecology, especially in light of ongoing environmental change. Amphibians and reptiles are emerging as a key group for studying these questions, with new data accumulating acceleratingly. We collected empirical data from earlier databases, reviews and primary literature to create the most up-to-date database on herpetological sex determination. We named our database HerpSexDet, which currently features data on genetic and temperature-dependent sex determination as well as reports on sex reversal for a total of 192 amphibian and 697 reptile species. This dataset, which we will regularly update in the future, facilitates interspecific comparative studies on the evolution of sex determination and its consequences for species-specific traits such as life history and conservation status, and may also help guiding future research by identifying species or higher taxa that are potentially most enlightening for the study of environmentally driven sex reversal.

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log₁₀ (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log₁₀ (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log₁₀ (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log₁₀ (Foxl2) and Log₁₀ (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest

Satellite DNAs (satDNAs) are one of the most abundant elements in genomes. Characterized as tandemly organized sequences that can be amplified into multiple copies, mainly in heterochromatic regions. The frog P. boiei (2n = 22, ZZ♂/ZW♀) is found in the Brazilian Atlantic forest and has an atypical pattern of heterochromatin distribution when compared to other anuran amphibians, with large pericentromeric blocks on all chromosomes. In addition, females of Proceratophrys boiei have a metacentric sex chromosome W showing heterochromatin in all chromosomal extension. In this work, we performed high-throughput genomic, bioinformatic, and cytogenetic analyses to characterize the satellite DNA content (satellitome) in P. boiei, mainly due to high amount of C-positive heterochromatin and the highly heterochromatic W sex chromosome. After all the analyses, it is remarkable that the satellitome of P. boiei is composed of a high number of satDNA families (226), making P. boiei the frog species with the highest number of satellites described so far. Consistent with the observation of large centromeric C-positive heterochromatin blocks, the genome of P. boiei is enriched with high copy number of repetitive DNAs, with total satDNA abundance comprising 16.87% of the genome. We successfully mapped via Fluorescence in situ hybridization the two most abundant repeats in the genome, PboSat01-176 and PboSat02-192, highlighting the presence of certain satDNAs sequences in strategic chromosomal regions (e.g., centromere and pericentromeric region), which leads to their participation in crucial processes for genomic organization and maintenance. Our study reveals a great diversity of satellite repeats that are driving genomic organization in this frog species. The characterization and approaches regarding satDNAs in this species of frog allowed the confirmation of some insights from satellite biology and a possible relationship with the evolution of sex chromosomes, especially in anuran amphibians, including P. boiei, for which data were not available.

Sex-linked markers in an Australian frog Platyplectrum ornatum (Limnodynastidae) with a small genome and homomorphic sex chromosomes

Amphibians have highly diverse sex-determining modes leading to a notable interest in vertebrate sex determination and sex chromosome evolution. The identification of sex-determining systems in amphibians, however, is often difficult as a vast majority consist of homomorphic sex chromosomes making them hard to distinguish. In this study, we used Diversity Array Technology sequencing (DArTseq) to identify the sex-determining system in the ornate burrowing frog from Australia, Platyplectrum ornatum. We applied DArTseq to 44 individuals, 19 males and 25 females, collected from two locations to develop sex-linked markers. Unexpectedly, these 44 individuals were classified into two distinct population clusters based on our SNP analyses, 36 individuals in cluster 1, and 8 individuals in cluster 2. We then performed sex-linkage analyses separately in each cluster. We identified 35 sex-linked markers from cluster 1, which were all associated with maleness. Therefore, P. ornatum cluster 1 is utilising a male heterogametic (XX/XY) sex-determining system. On the other hand, we identified 210 sex-linked markers from cluster 2, of which 89 were male specific, i.e., identifying XX/XY sex determining system and 111 were female specific, i.e., identifying ZZ/ZW sex determining system, suggesting existence of either male or female heterogametic sex determining system in cluster 2. We also performed cytogenetic analyses in 1 male and 1 female from cluster 1; however, we did not detect any visible differentiation between the X and Y sex chromosomes. We also mapped sex-linked markers from the two clusters against the P. ornatum genome and our comparative analysis indicated that the sex chromosomes in both clusters shared homologies to chromosome 10 (autosome) of Rana temporaria and ZWY sex chromosome of Xenopus tropicalis. Our preliminary data suggest that it is plausible that the cluster 2 has a potential to be either male or female heterogamety in sex determination, requiring further investigation.

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.

Whole-Genome Inter-Sex Variation in Russian Sturgeon (Acipenser gueldenstaedtii)

The Russian sturgeon (Acipenser gueldenstaedtii, AG) is an endangered fish species increasingly raised on fish farms for black caviar. Understanding the process of sex determination in AG is, therefore, of scientific and commercial importance. AG lacks sexual dimorphism until sexual maturation and has a predominantly octoploid genome without a definite sex chromosome. A conserved short female-specific genomic sequence was recently described, leading to the development of a genetic sex marker. However, no biological function has been reported for this sequence. Thus, the mechanism of sex determination and the overall inter-sex genomic variation in AG are still unknown. To comprehensively analyze the inter-sex genomic variation and assess the overall inter-species variation between AG and A. ruthenus (AR, sterlet), a related tetraploid sturgeon species, we performed whole-genome sequencing on DNA from 10 fish-farm-raised adult AG (5 males and 5 females). We produced a partially assembled, ~2390 MBp draft genome for AG. We validated in AG the female-specific region previously described in AR. We identified ~2.8 million loci (SNP/indels) varying between the species, but only ~7400 sex-associated loci in AG. We mapped the sex-associated AG loci to the AR genome and identified 15 peaks of sex-associated variation (10 kb segments with 30 or more sex-associated variants), 1 of which matched the previously reported sex-variable region. Finally, we identified 14 known and predicted genes in proximity to these peaks. Our analysis suggests that one or more of these genes may have functional roles in sex determination and/or sexual differentiation in sturgeons. Further functional studies are required to elucidate these roles.

Ecology and the evolution of sex chromosomes

Sex chromosomes are common features of animal genomes, often carrying a sex determination gene responsible for initiating the development of sexually dimorphic traits. The specific chromosome that serves as the sex chromosome differs across taxa as a result of fusions between sex chromosomes and autosomes, along with sex chromosome turnover-autosomes becoming sex chromosomes and sex chromosomes 'reverting' back to autosomes. In addition, the types of genes on sex chromosomes frequently differ from the autosomes, and genes on sex chromosomes often evolve faster than autosomal genes. Sex-specific selection pressures, such as sexual antagonism and sexual selection, are hypothesized to be responsible for sex chromosome turnovers, the unique gene content of sex chromosomes and the accelerated evolutionary rates of genes on sex chromosomes. Sex-specific selection has pronounced effects on sex chromosomes because their sex-biased inheritance can tilt the balance of selection in favour of one sex. Despite the general consensus that sex-specific selection affects sex chromosome evolution, most population genetic models are agnostic as to the specific sources of these sex-specific selection pressures, and many of the details about the effects of sex-specific selection remain unresolved. Here, I review the evidence that ecological factors, including variable selection across heterogeneous environments and conflicts between sexual and natural selection, can be important determinants of sex-specific selection pressures that shape sex chromosome evolution. I also explain how studying the ecology of sex chromosome evolution can help us understand important and unresolved aspects of both sex chromosome evolution and sex-specific selection.

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.

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.

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.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Sex-related differences in aging rate are associated with sex chromosome system in amphibians

Sex‐related differences in mortality are widespread in the animal kingdom. Although studies have shown that sex determination systems might drive lifespan evolution, sex chromosome influence on aging rates have not been investigated so far, likely due to an apparent lack of demographic data from clades including both XY (with heterogametic males) and ZW (heterogametic females) systems. Taking advantage of a unique collection of capture–recapture datasets in amphibians, a vertebrate group where XY and ZW systems have repeatedly evolved over the past 200 million years, we examined whether sex heterogamy can predict sex differences in aging rates and lifespans. We showed that the strength and direction of sex differences in aging rates (and not lifespan) differ between XY and ZW systems. Sex‐specific variation in aging rates was moderate within each system, but aging rates tended to be consistently higher in the heterogametic sex. This led to small but detectable effects of sex chromosome system on sex differences in aging rates in our models. Although preliminary, our results suggest that exposed recessive deleterious mutations on the X/Z chromosome (the “unguarded X/Z effect”) or repeat‐rich Y/W chromosome (the “toxic Y/W effect”) could accelerate aging in the heterogametic sex in some vertebrate clades.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.