Novel sex-determining genes in fish and sex chromosome evolution

Sex chromosome turnover and biodiversity in fishes

The impact of sex chromosomes and their turnover in speciation remains a subject of ongoing debate in the field of evolutionary biology. Fishes are the largest group of vertebrates, and they exhibit unparalleled sexual plasticity, as well as diverse sex-determining (SD) genes, sex chromosomes, and sex determination mechanisms. This diversity is hypothesized to be associated with the frequent turnover of sex chromosomes in fishes. Although it is evident that amh and amhr2 are repeatedly and independently recruited as SD genes, their relationship with the rapid turnover of sex chromosomes and the biodiversity of fishes remains unknown. We summarize the canonical models of sex chromosome turnover and highlight the vital roles of gene mutation and hybridization with empirical evidence. We revisit Haldane's rule and the large X-effect and propose the hypothesis that sex chromosomes accelerate speciation by multiplying genotypes via hybridization. By integrating recent findings on the turnover of SD genes, sex chromosomes, and sex determination systems in fish species, this review provides insights into the relationship between sex chromosome evolution and biodiversity in fishes.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

Comparative transcriptome analysis of ovaries and testes reveals sex-biased genes and pathways in zebrafish

Zebrafish (Danio rerio) is a widely recognized and extensively studied model organism in scientific research. The regulatory mechanism of gonadal development and differentiation of this species has aroused considerable attention. Nonetheless, the major sex-biased genes and pathways associated with gonadal development remain elusive. Therefore, to comprehend this intricate process, gonadal transcriptome sequencing was carried out to identify differentially expressed genes (DEGs) between the testes and ovaries of adult zebrafish. The preliminary assessment yielded a total of 23,529,272 and 23,521,368 clean reads from the cDNA libraries of ovaries and testes. Afterward, a comparative analysis of the transcriptome revealed 3,604 upregulated and 11,371 downregulated DEGs in the ovaries compared to the testes. Of these genes, 428 were exclusively expressed in females, while 3,516 were exclusively expressed in males. Additionally, further assessments were conducted to explore the functions associated with these DEGs in various biological processes. The data revealed their involvement in sex-biased pathways, such as progesterone-mediated oocyte maturation, oocyte meiosis, cytokine-cytokine receptor interaction, and cardiac muscle contraction. Finally, the expression levels of 14 sex-biased DEGs (cdc20, ccnb1, ypel3, chn1, bmp15, rspo1, tnfsf10, egfra, acta2, cox8a, gsdf, dmrt1, star, and cyp17a1) associated with the enriched pathways were subjected to further validation through qRT-PCR. The data acquired from these investigations offer valuable resources to support further exploration of the mechanisms governing sexual dimorphism and gonadal development in zebrafish.

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.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

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.

Estrogen Signaling Inhibits the Expression of anti-Müllerian hormone (amh) and gonadal-soma-derived factor (gsdf) during the Critical Time of Sexual Fate Determination in Zebrafish

The mechanism of fish gonadal sex differentiation is complex and regulated by multiple factors. It has been widely known that proper steroidogenesis in Leydig cells and sex-related genes in Sertoli cells play important roles in gonadal sex differentiation. In teleosts, the precise interaction of these signals during the sexual fate determination remains elusive, especially their effect on the bi-potential gonad during the critical stage of sexual fate determination. Recently, all-testis phenotypes have been observed in the cyp17a1-deficient zebrafish and common carp, as well as in cyp19a1a-deficient zebrafish. By mating cyp17a1-deficient fish with transgenic zebrafish Tg(piwil1:EGFP-nanos3UTR), germ cells in the gonads were labelled with enhanced green fluorescent protein (EGFP). We classified the cyp17a1-deficient zebrafish and their control siblings into primordial germ cell (PGC)-rich and -less groups according to the fluorescence area of the EGFP labelling. Intriguingly, the EGFP-labelled bi-potential gonads in cyp17a1+/+ fish from the PGC-rich group were significantly larger than those of the cyp17a1-/- fish at 23 days post-fertilization (dpf). Based on the transcriptome analysis, we observed that the cyp17a1-deficient fish of the PGC-rich group displayed a significantly upregulated expression of amh and gsdf compared to that of control fish. Likewise, the upregulated expressions of amh and gsdf were observed in cyp19a1a-deficient fish as examined at 23 dpf. This upregulation of amh and gsdf could be repressed by treatment with an exogenous supplement of estradiol. Moreover, tamoxifen, an effective antagonist of both estrogen receptor α and β (ERα and Erβ), upregulates the expression of amh and gsdf in wild-type (WT) fish. Using the cyp17a1- and cyp19a1a-deficient zebrafish, we provide evidence to show that the upregulated expression of amh and gsdf due to the compromised estrogen signaling probably determines their sexual fate towards testis differentiation. Collectively, our data suggest that estrogen signaling inhibits the expression of amh and gsdf during the critical time of sexual fate determination, which may broaden the scope of sex steroid hormones in regulating gonadal sex differentiation in fish.

hsd17b1 is a key gene for ovarian differentiation of the Siberian sturgeon

This is the first work using gonads from undifferentiated, genetically-sexed Siberian sturgeon describing expression changes in genes related to steroid synthesis and female and male sex differentiation. One factor identified as relevant for ovarian differentiation was the gene coding for the enzyme Hsd17b1, which converts estrone into estradiol-17β. hsd17b1 was highly activated in female gonads at 2.5 months of age, around the onset of sex differentiation, preceding activation of two other genes involved in estrogen production (cyp19a1 and foxl2). hsd17b1 was also strongly repressed in males. Two known foxl2 paralogs are found in Siberian sturgeon-foxl2 and foxl2l-but only foxl2 appeared to be associated with ovarian differentiation. With regard to the male pathway, neither 11-oxygenated androgens nor classic male genes (amh, dmrt1, sox9, and dhh) were found to be involved in male sex differentiation, leaving open the question of which genes participate in early male gonad development in this ancient fish. Taken together, these results indicate an estrogen-dependence of female sex differentiation and 11-oxygenated androgen-independence of male sex differentiation.

Sex-Associated SNP Confirmation of Sex-Reversed Male Farmed Japanese Flounder Paralichthys olivaceus

Female Japanese flounder Paralichthys olivaceus grow more rapidly than the male. The goal of all-female commercial production requires an efficient method of genetic sex identification. We conducted genome-wide association analysis of female and male farmed Japanese flounder (n = 24 per phenotypic sex) and found all regions of chromosome 24 to be significantly associated with phenotypic sex, suggesting it as the sex chromosome. Genetic sex was identified based on single nucleotide polymorphisms (SNP) on chromosome 24 (n = 3568) using multidimensional scaling analysis, and individuals were clearly separated according to sex by the first dimension. The 61 SNPs most highly associated with sex were selected, and an amplicon-based SNP panel was developed. This was used to determine genetic sex of 39 females and 40 males. Eleven phenotypic males were assigned as female with XX genotype, suggesting sex reversal. Genetic sex was also assessed based on the indel of the amh gene promoter, which is the major candidate sex gene of Japanese flounder. We found four SNPs perfectly associated with genotypic sex in the sex-associated SNP panel, one of which was located in exon 2 of the amh gene. Along with the indel of the amh gene promoter, the sex-associated SNP panel will be of value in identifying genetic sex of farmed Japanese flounder. Molecular sexing will facilitate all-female production by breeding sex-reversed males.

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.

No evidence for sex chromosomes in natural populations of the cichlid fish Astatotilapia burtoni

Sex determination (SD) is not conserved among teleost fishes and can even differ between populations of the same species. Across the outstandingly species-rich fish family Cichlidae, more and more SD systems are being discovered. Still, the picture of SD evolution in this group is far from being complete. Lake Tanganyika and its affluent rivers are home to Astatotilapia burtoni, which belongs to the extremely successful East African cichlid lineage Haplochromini. Previously, in different families of an A. burtoni laboratory strain, an XYW system and an XY system have been described. The latter was also found in a second laboratory strain. In a laboratory-reared family descending from a population of the species' southern distribution, a second XY system was discovered. Yet, an analysis of sex chromosomes for the whole species distribution is missing. Here, we examined the genomes of 11 natural populations of A. burtoni, encompassing a wide range of its distribution, for sex-linked regions. We did not detect signs of differentiated sex chromosomes and also not the previously described sex chromosomal systems present in laboratory lines, suggesting different SD systems in the same species under natural and (long-term) artificial conditions. We suggest that SD in A. burtoni is more labile than previously assumed and consists of a combination of non-genetic, polygenic, or poorly differentiated sex chromosomes.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

Comparative gonad transcriptome analysis in cobia (Rachycentron canadum)

Background: Cobia (Rachycentron canadum) is a species of fish with high commercial potential particularly due to fast growth rates. The evidence of sexual size dimorphism favoring females indicate potential benefits in having a monosex culture. However, the involvement of genetic factors responsible for sexual development and gonadal maintenance that produces phenotypic sex in cobia is largely unknown. Methods: In the present study, we performed transcriptome sequencing of cobia to identify sex-biased significantly differentially expressed genes (DEGs) in testes and ovaries. The reliability of the gonad transcriptome data was validated by qPCR analysis of eight selected significantly differential expressed sex-related candidate genes. Results: This comparative gonad transcriptomic analysis revealed that 7,120 and 4,628 DEGs are up-regulated in testes or ovaries, respectively. Further functional annotation analyses identified 76 important candidate genes involved in sex determination cascades or sex differentiation, including 42 known testis-biased DEGs (dmrt1, amh and sox9 etc.), and 34 known ovary-biased DEGs (foxl2, sox3 and cyp19a etc.). Moreover, eleven significantly enriched pathways functionally related to sex determination and sex differentiation were identified, including Wnt signaling pathway, oocyte meiosis, the TGF-beta signaling pathway and MAPK signaling pathway. Conclusion: This work represents the first comparative gonad transcriptome study in cobia. The putative sex-associated DEGs and pathways provide an important molecular basis for further investigation of cobia's sex determination, gonadal development as well as potential control breeding of monosex female populations for a possible aquaculture setting.

Genomic imprinting-like monoallelic paternal expression determines sex of channel catfish

The X and Y chromosomes of channel catfish have the same gene contents. Here, we report allelic hypermethylation of the X chromosome within the sex determination region (SDR). Accordingly, the X-borne hydin-1 gene was silenced, whereas the Y-borne hydin-1 gene was expressed, making monoallelic expression of hydin-1 responsible for sex determination, much like genomic imprinting. Treatment with a methylation inhibitor, 5-aza-dC, erased the epigenetic marks within the SDR and caused sex reversal of genetic females into phenotypic males. After the treatment, hydin-1 and six other genes related to cell cycle control and proliferative growth were up-regulated, while three genes related to female sex differentiation were down-regulated in genetic females, providing additional support for epigenetic sex determination in catfish. This mechanism of sex determination provides insights into the plasticity of genetic sex determination in lower vertebrates and its connection with temperature sex determination where DNA methylation is broadly involved.

Trioecy is maintained as a time-stable mating system in the pink sea urchin Toxopneustes roseus from the Mexican Pacific

Trioecy is a sexual system that consists of the co-occurrence of females, males and hermaphrodites in a population and is common in plants; however, in animals it is uncommon and poorly understood. In echinoderms, trioecy had never been recorded until now. Frequencies of females, males, and hermaphrodites were evaluated and gametogenic development was histologically characterized in a population of Toxopneustes roseus inhabiting the Mexican Pacific. Trioecy in this population is functional and temporally stable, since the three sexes coexisted in each sampling month. The hermaphrodites presented similar gametogenic development as the females and males and participated during the spawning season, contributing to the population's reproductive process. Trioecy is considered an evolutionarily transitory state, and it is extremely difficult to explain its presence in a species. We hypothesize that continuous ocean warming represents a threat to the survival of this population of T. roseus, since its early developmental stages, which represent a population bottleneck, are more vulnerable to high temperatures than other sea urchins inhabiting the area, while its population density is significantly lower. These conditions generate a strongly stressed environment, which is the determining factor that maintains the stability of trioecy in the species in which it has been studied.

Repeated turnovers keep sex chromosomes young in willows

Background

Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively.

Results

Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the “deleterious mutation load” and “sexually antagonistic selection” theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration.

Conclusions

Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

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.

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.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish

SignificanceAlthough turnover of sex chromosomes is very common in many vertebrate lineages, the transition process is still elusive. We studied the sex-determining region (SDR) of 12 congeneric fish species of Takifugu fish that compose an important model for the study of genomics and sex determination. We found that while nine species retained their ancestral SDR, three species had acquired derived SDRs. Although the derived SDRs resided in three different chromosomes, they harbored a shared supergene flanked by two putative transposable elements. The results highlight the underestimated role of a mobile supergene in turnover of sex chromosomes in vertebrates.

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.

Generation of a chromosome‐level genome assembly for Pacific halibut ( Hippoglossus stenolepis ) and characterization of its sex‐determining genomic region

The Pacific halibut (Hippoglossus stenolepis) is a key species in the North Pacific Ocean and Bering Sea ecosystems, where it also supports important fisheries. However, the lack of genomic resources limits our understanding of evolutionary, environmental and anthropogenic forces affecting key life history characteristics of Pacific halibut and prevents the application of genomic tools in fisheries management and conservation efforts. In the present study, we report on the first generation of a high‐quality chromosome‐level assembly of the Pacific halibut genome, with an estimated size of 602 Mb, 24 chromosome‐length scaffolds that contain 99.8% of the assembly and a N₅₀ scaffold length of 27.3 Mb. In the first application of this important resource, we conducted genome‐wide analyses of sex‐specific genetic variation by pool sequencing and characterized a potential sex‐determining region in chromosome 9 with a high density of female‐specific SNPs. Within this region, we identified the bmpr1ba gene as a potential candidate for master sex‐determining (MSD) gene. bmpr1ba is a member of the TGF‐β family that in teleosts has provided the largest number of MSD genes, including a paralogue of this gene in Atlantic herring. The genome assembly constitutes an essential resource for future studies on Pacific halibut population structure and dynamics, evolutionary history and responses to environmental and anthropogenic influences. Furthermore, the genomic location of the sex‐determining region in Pacific halibut has been identified and a putative candidate MSD gene has been proposed, providing further support for the rapid evolution of sex‐determining mechanisms in teleost fish.

Identification of sex chromosome and sex-determining gene of southern catfish (Silurus meridionalis) based on XX, XY and YY genome sequencing

Teleosts are important models to study sex chromosomes and sex-determining (SD) genes because they present a variety of sex determination systems. Here, we used Nanopore and Hi-C technologies to generate a high-contiguity chromosome-level genome assembly of a YY southern catfish (Silurus meridionalis). The assembly is 750.0 Mb long, with contig N50 of 15.96 Mb and scaffold N50 of 27.22 Mb. We also sequenced and assembled an XY male genome with a size of 727.2 Mb and contig N50 of 13.69 Mb. We identified a candidate SD gene through comparisons to our previous assembly of an XX individual. By resequencing male and female pools, we characterized a 2.38 Mb sex-determining region (SDR) on Chr24. Analysis of read coverage and comparison of the X and Y chromosome sequences showed a Y specific insertion (approx. 500 kb) in the SDR which contained a male-specific duplicate of amhr2 (named amhr2y). amhr2y and amhr2 shared high-nucleotide identity (81.0%) in the coding region but extremely low identity in the promotor and intron regions. The exclusive expression in the male gonadal primordium and loss-of-function inducing male to female sex reversal confirmed the role of amhr2y in male sex determination. Our study provides a new example of amhr2 as the SD gene in fish and sheds light on the convergent evolution of the duplication of AMH/AMHR2 pathway members underlying the evolution of sex determination in different fish lineages.

Sexual development dysgenesis in interspecific hybrids of Medaka fish

Fish are amongst vertebrates the group with the highest diversity of known sex-determining genes. Particularly, the genus Oryzias is a suitable taxon to understand how different sex determination genetic networks evolved in closely related species. Two closely related species, O. latipes and O. curvinotus, do not only share the same XX/XY sex chromosome system, but also the same male sex-determining gene, dmrt1bY. We performed whole mRNA transcriptomes and morphology analyses of the gonads of hybrids resulting from reciprocal crosses between O. latipes and O. curvinotus. XY male hybrids, presenting meiotic arrest and no production of sperm were sterile, and about 30% of the XY hybrids underwent male-to-female sex reversal. Both XX and XY hybrid females exhibited reduced fertility and developed ovotestis while aging. Transcriptome data showed that male-related genes are upregulated in the XX and XY female hybrids. The transcriptomes of both types of female and of the male gonads are characterized by upregulation of meiosis and germ cell differentiation genes. Differences in the parental species in the downstream pathways of sexual development could explain sex reversal, sterility, and the development of intersex gonads in the hybrids. We hypothesize that male-to-female sex reversal may be connected to a different development time between species at which dmrt1bY expression starts. Our results provide molecular clues for the proximate mechanisms of hybrid incompatibility and Haldane’s rule.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Evaluation of four methods to identify the homozygotic sex chromosome in small populations

BACKGROUND

Whole genomes are commonly assembled into a collection of scaffolds and often lack annotations of autosomes, sex chromosomes, and organelle genomes (i.e., mitochondrial and chloroplast). As these chromosome types differ in effective population size and can have highly disparate evolutionary histories, it is imperative to take this information into account when analysing genomic variation. Here we assessed the accuracy of four methods for identifying the homogametic sex chromosome in a small population using two whole genome sequences (WGS) and 133 RAD sequences of white-tailed eagles (Haliaeetus albicilla): i) difference in read depth per scaffold in a male and a female, ii) heterozygosity per scaffold in a male and a female, iii) mapping to the reference genome of a related species (chicken) with annotated sex chromosomes, and iv) analysis of SNP-loadings from a principal components analysis (PCA), based on the low-depth RADseq data.

RESULTS

The best performing approach was the reference mapping (method iii), which identified 98.12% of the expected homogametic sex chromosome (Z). Read depth per scaffold (method i) identified 86.41% of the homogametic sex chromosome with few false positives. SNP-loading scores (method iv) identified 78.6% of the Z-chromosome and had a false positive discovery rate of more than 10%. Heterozygosity per scaffold (method ii) did not provide clear results due to a lack of diversity in both the Z and autosomal chromosomes, and potential interference from the heterogametic sex chromosome (W). The evaluation of these methods also revealed 10 Mb of putative PAR and gametologous regions.

CONCLUSION

Identification of the homogametic sex chromosome in a small population is best accomplished by reference mapping or examining differences in read depth between sexes.

Heterochiasmy and the establishment of gsdf as a novel sex determining gene in Atlantic halibut

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly from a male and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male-only and female-only meiotic recombination regions (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4.5 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. Our findings indicate that heterochiasmy may facilitate the evolution of genetic sex determination systems relying on linkage of sexually antagonistic loci to a sex-determining factor.

Even closely related fish species can have different sex chromosomes, but this turn-over of sex determination systems is poorly understood. Here, we used large-scale genome sequencing to determine the DNA sequence of the Atlantic halibut chromosomes and compared sequencing data from males and females to identify the sex chromosomes. We show that males have much higher gene activity of the gene gonadal somatic cell derived factor (gsdf), which is located on the sex chromosomes and has a role in testicular development. The genome contains many mobile DNA sequences, transposable elements (TEs), one placed in front of gsdf, enhancing its activity. This made gsdf the sex determining factor, thereby creating a new Y-chromosome. We further describe how all Atlantic halibut chromosomes behave similar to sex chromosomes in that most regions only recombine in one sex. This phenomenon may contribute to the rapid turn-over of genetic sex determination systems in fish. Our results highlight the molecular events creating a new Y-chromosome and show that the new Atlantic halibut Y was formed less than 4.5 million years ago. Future studies in Atlantic halibut and closely related species can shed light on mechanisms contributing to sex chromosome evolution in fish.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish

Background

Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens.

Results

We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication.

Conclusions

Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01205-y.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Reference genome of Lumpfish Cyclopterus lumpus Linnaeus provides evidence of male heterogametic sex determination through the AMH pathway

Teleosts exhibit extensive diversity of sex determination (SD) systems and mechanisms, providing the opportunity to study the evolution of sex determination and sex chromosomes. Here we sequenced the genome of the Common Lumpfish (Cyclopterus lumpus Linnaeus), a species of increasing importance to aquaculture, and identified the SD region and master SD locus using a 70K SNP array and tissue-specific expression data. The chromosome-level assembly identified 25 diploid chromosomes with a total size of 572.89 Mb, a scaffold N50 of 23.86 Mb, and genome annotation predicted 21,480 protein-coding genes. Genome wide association analysis located a highly sex-associated region on chromosome 13, suggesting that anti-Müllerian hormone (AMH) is the putative SD factor. Linkage disequilibrium and heterozygosity across chromosome 13 support a proto-XX/XY system, with an absence of widespread chromosome divergence between sexes. We identified three copies of AMH in the Lumpfish primary and alternate haplotype assemblies localized in the SD region. Comparison to sequences from other teleosts suggested a monophyletic relationship and conservation within the Cottioidei. One AMH copy showed similarity to AMH/AMHY in a related species and was also the only copy with expression in testis tissue, suggesting this copy may be the functional copy of AMH in Lumpfish. The two other copies arranged in tandem inverted duplication were highly similar, suggesting a recent duplication event. This study provides a resource for the study of early sex chromosome evolution and novel genomic resources that benefits Lumpfish conservation management and aquaculture.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

Control of testes mass by androgen receptor paralogs in a cichlid

Steroid hormones play numerous important and diverse roles in the differentiation and development of vertebrate primary and secondary reproductive characteristics. However, the exact role of androgen receptors-which bind circulating androgens-in this regulatory pathway is unclear. Teleost fishes further complicate this question by having two paralogs of the androgen receptor (ARα and ARβ) resulting from a duplication of their ancestral genome. We investigated the functional role of these two ARs on adult testes mass, by eliminating receptor function of one or both paralogs using CRISPR/Cas9 genome edited Astatotilapia burtoni, an African cichlid fish. Fish with two or more functional AR alleles were more likely to be male compared to fish with one or fewer, suggesting that the two paralogs may play redundant roles in the A. burtoni sex determination system. We replicated previous work showing that fish lacking functional ARβ possess testes smaller than wild-type fish, while fish lacking ARα possess testes larger than wild-type fish. However, we found novel evidence supporting a complex relationship between the two AR paralogs in the regulation of testes mass. For instance, the effects of ARα mutation on testes mass are eliminated in homozygous ARβ mutants but the reverse is not true. These results suggest a dynamic relationship between these two AR paralogs where ARβ functions may be permissive to ARα functions in the control of testes mass. This mechanism may contribute to the robust physiological plasticity displayed by A. burtoni and other social teleost fishes.

Sex Determination and Differentiation in Teleost: Roles of Genetics, Environment, and Brain

The fish reproductive system is a complex biological system. Nonetheless, reproductive organ development is conserved, which starts with sex determination and then sex differentiation. The sex of a teleost is determined and differentiated from bipotential primordium by genetics, environmental factors, or both. These two processes are species-specific. There are several prominent genes and environmental factors involved during sex determination and differentiation. At the cellular level, most of the sex-determining genes suppress the female pathway. For environmental factors, there are temperature, density, hypoxia, pH, and social interaction. Once the sexual fate is determined, sex differentiation takes over the gonadal developmental process. Environmental factors involve activation and suppression of various male and female pathways depending on the sexual fate. Alongside these factors, the role of the brain during sex determination and differentiation remains elusive. Nonetheless, GnRH III knockout has promoted a male sex-biased population, which shows brain involvement during sex determination. During sex differentiation, LH and FSH might not affect the gonadal differentiation, but are required for regulating sex differentiation. This review discusses the role of prominent genes, environmental factors, and the brain in sex determination and differentiation across a few teleost species.

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

A 180 Myr-old female-specific genome region in sturgeon reveals the oldest known vertebrate sex determining system with undifferentiated sex chromosomes

Several hypotheses explain the prevalence of undifferentiated sex chromosomes in poikilothermic vertebrates. Turnovers change the master sex determination gene, the sex chromosome or the sex determination system (e.g. XY to WZ). Jumping master genes stay main triggers but translocate to other chromosomes. Occasional recombination (e.g. in sex-reversed females) prevents sex chromosome degeneration. Recent research has uncovered conserved heteromorphic or even homomorphic sex chromosomes in several clades of non-avian and non-mammalian vertebrates. Sex determination in sturgeons (Acipenseridae) has been a long-standing basic biological question, linked to economical demands by the caviar-producing aquaculture. Here, we report the discovery of a sex-specific sequence from sterlet (Acipenser ruthenus). Using chromosome-scale assemblies and pool-sequencing, we first identified an approximately 16 kb female-specific region. We developed a PCR-genotyping test, yielding female-specific products in six species, spanning the entire phylogeny with the most divergent extant lineages (A. sturio, A. oxyrinchus versus A. ruthenus, Huso huso), stemming from an ancient tetraploidization. Similar results were obtained in two octoploid species (A. gueldenstaedtii, A. baerii). Conservation of a female-specific sequence for a long period, representing 180 Myr of sturgeon evolution, and across at least one polyploidization event, raises many interesting biological questions. We discuss a conserved undifferentiated sex chromosome system with a ZZ/ZW-mode of sex determination and potential alternatives.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Epigenetic mechanisms in sex determination and in the evolutionary transitions between sexual systems

The hypothesis that epigenetic mechanisms of gene expression regulation have two main roles in vertebrate sex is presented. First, and within a given generation, by contributing to the acquisition and maintenance of (i) the male or female function once during the lifetime in individuals of gonochoristic species; and (ii) the male and female function in the same individual, either at the same time in simultaneous hermaphrodites, or first as one sex and then as the other in sequential hermaphrodites. Second, if environmental conditions change, epigenetic mechanisms may have also a role across generations, by providing the necessary phenotypic plasticity to facilitate the transition: (i) from one sexual system to another, or (ii) from one sex-determining mechanism to another. Furthermore, if the environmental change lasts enough time, epimutations could facilitate assimilation into genetic changes that stabilize the new sexual system or sex-determining mechanism. Examples supporting these assertions are presented, caveats or difficulties and knowledge gaps identified, and possible ways to test this hypothesis suggested.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Female-specific genomic regions and molecular sex identification of the clearhead icefish (Protosalanx hyalocranius)

Background

The clearhead icefish, Protosalanx hyalocranius, is an economically important fishery species in China. Since 1980s, P. hyalocranius was widely introduced into lakes and reservoirs of northern China for aquaculture. However, the lack of a rapid and cost-effective sex identification method based on sex specific genetic markers has hindered study on sex determination mechanisms and breeding applications.

Results

Female-specific genomic regions were discovered by comparing whole genome re-sequencing data of both males and females. Two female-specific genomic regions larger than 50 bp were identified, and one (598 bp) contained a putative FOXI gene, which was paralogous to another FOXI gene with sex-associated SNPs. The two FOXI sequences displayed significant length difference with nine deletions of total length of 230 bp. This deletion-type structural variation could be easily and efficiently detected by traditional PCR and agarose gel electrophoresis with one 569 bp band for males and two bands (569 and 339 bp) for females, which were validated in 50 females and 40 males with known phenotypic sexes.

Conclusions

The results provided structural genomic evidence for the ZZ/ZW sex determination system in P. hyalocranius discovered in our previous study with association analysis of SNPs. Moreover, the female-specific markers and rapid and cost-effective PCR-based genetic sex identification method should have applications in further studies of sex determination mechanism for this species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07830-9.

Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY. By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

A First Insight into the Gonad Transcriptome of Hong Kong Catfish (Clarias fuscus)

Hong Kong catfish (Clarias fuscus) exhibit sexual dimorphism, particularly in body size. Due to the fast growth rate of males, the sexual size dimorphism of Hong Kong catfish has become an economically important trait. However, limited knowledge is known about the molecular mechanisms of sex determination and sex differentiation in this species. In this study, a first de novo transcriptome sequencing analysis of testes and ovaries was performed to identify sex-biased genes in Hong Kong catfish. The results showed that a total of 290,291 circular consensus sequences (CCSs) were obtained, from which 248,408 full-length non-chimeric (FLNC) reads were generated. After non-redundant analysis, a total of 37,305 unigenes were predicted, in which 34,342 unigenes were annotated with multiple public databases. Comparative transcriptomic analysis identified 5750 testis-biased differentially expressed genes (DEGs) and 6991 ovary-biased DEGs. The enrichment analysis showed that DEGs were classified into 783 Gene Ontology (GO) terms and 16 Kyoto Encyclopedia of Gene and Genome (KEGG) pathways. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, androgen secretion, gonadal development and steroid biosynthesis pathways. In addition, the expression levels of 23 unigenes were confirmed to validate the transcriptomic data by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first investigation into the transcriptome of Hong Kong catfish testes and ovaries. This study provides an important molecular basis for the sex determination and sex control breeding of Hong Kong catfish.

The Diversity and Evolution of Sex Chromosomes in Frogs

Frogs are ideal organisms for studying sex chromosome evolution because of their diversity in sex chromosome differentiation and sex-determination systems. We review 222 anuran frogs, spanning ~220 Myr of divergence, with characterized sex chromosomes, and discuss their evolution, phylogenetic distribution and transitions between homomorphic and heteromorphic states, as well as between sex-determination systems. Most (~75%) anurans have homomorphic sex chromosomes, with XY systems being three times more common than ZW systems. Most remaining anurans (~25%) have heteromorphic sex chromosomes, with XY and ZW systems almost equally represented. There are Y-autosome fusions in 11 species, and no W-/Z-/X-autosome fusions are known. The phylogeny represents at least 19 transitions between sex-determination systems and at least 16 cases of independent evolution of heteromorphic sex chromosomes from homomorphy, the likely ancestral state. Five lineages mostly have heteromorphic sex chromosomes, which might have evolved due to demographic and sexual selection attributes of those lineages. Males do not recombine over most of their genome, regardless of which is the heterogametic sex. Nevertheless, telomere-restricted recombination between ZW chromosomes has evolved at least once. More comparative genomic studies are needed to understand the evolutionary trajectories of sex chromosomes among frog lineages, especially in the ZW systems.

Genome wide identification, phylogeny, and synteny analysis of sox gene family in common carp (Cyprinus carpio)

•

27 SOX (high-mobility group HMG-box) genes were identified in the C. carp genome.

•

SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative protein series from 307 to 509 amino acids.

•

Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %.

•

Chromosomal location and synteny analysis display all SOX gene are located on different chromosomes.

Common carp (Cyprinus carpio) is a commercial fish species valuable for nutritious components and plays a vital role in human healthy nutrition. The SOX (SRY-related genes systematically characterized by a high-mobility group HMG-box) encoded important gene regulatory proteins, a family of transcription factors found in a broad range of animal taxa and extensively known for its contribution in multiple developmental processes including contribution in sex determination across phyla. In our current study, we initially accomplished a genome-wide analysis to report the SOX gene family in common carp fish based on available genomic sequences of zebrafish retrieved from gene repository databases, we focused on the global identification of the Sox gene family in Common carp among wide range of vertebrates and teleosts based on bioinformatics tools and techniques and explore the evolutionary relationships. In our results, a total of 27 SOX (high-mobility group HMG-box) domain genes were identified in the C. carp genome. The full length sequences of SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative proteins series from 307 to 509 amino acids and all gene having exon number expect SOX9 and SOX13. All the SOX proteins contained at least one conserved DNA-binding HMG-box domain and two (SOX7 and SOX18) were found C terminal. The Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %, average in biological regulation 45.188 %, biosynthetic process (19.992 %), regulation of cellular process 39.68, 45.508 % organic substance metabolic process, multicellular organismal process 23.23 %,developmental process 21.74 %, system development 16.59 %, gene expression 16.05 % and 14.337 % of RNA metabolic process. Chromosomal location and syntanic analysis show all SOX gene are located on different chromosomes and apparently does not fallow the unique pattern. The maximum linkage of chromosome is (2) on Unplaced Scaffold region. Finally, our results provide important genomic suggestion for upcoming studies of biochemical, physiological, and phylogenetic understanding on SOX genes among teleost.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

Males can adjust offspring sex ratio in an adaptive fashion through different mechanisms

Sex allocation research has primarily focused on offspring sex-ratio adjustment by mothers. Yet, fathers also benefit from producing more of the sex with greater fitness returns. Here, we review the state-of-the art in the study of male-driven sex allocation and, counter to the current paradigm, we propose that males can adaptively influence offspring sex ratio through a wide variety of mechanisms. This includes differential production and motility of X- versus Y-bearing sperms in mammals, variation in seminal fluid composition in haplo-diploid invertebrates, and epigenetic mechanisms in some fish and lizards exhibiting environmental sex determination. Conflicts of interest between mothers and fathers over offspring sex ratios can emerge, although many more studies are needed in this area. While many studies of sex allocation have focused on adaptive explanations with little attention to mechanisms, and vice versa, the integration of these two topics is essential for understanding male-driven sex allocation.

The rise and fall of the ancient northern pike master sex-determining gene

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

Epigenetic regulation of gonadal and brain aromatase expression in a cichlid fish with environmental sex determination

A fit animal must develop testes or ovaries, with brain and physiology to match. In species with alternative male morphs this coordination of development across tissues operates within sexes as well as between. For Pelvicachromis pulcher, an African cichlid in which early pH exposure influences both sex and alternative male morph, we sequence both copies of aromatase (cyp19a1), a key gene for sex determination. We analyze gene expression and epigenetic state, comparing gonad and brain tissue from females, alternative male morphs, and fry. Relative to brain, we find elevated expression of the A-copy in the ovaries but not testes. Methylation analysis suggests strong epigenetic regulation, with one region specifying sex and another specifying tissue. We find elevated brain expression of the B-copy with no sex or male morph differences. B-copy methylation follows that of the A-copy rather than corresponding to B-copy expression. In 30-day old fry, we see elevated B-copy expression in the head, but we do not see the expected elevated A-copy expression in the trunk that would reflect ovarian development. Interestingly, the A-copy epialleles that distinguish ovaries from testes are among the most explanatory patterns for variation among fry, suggesting epigenetic marking of sex prior to differentiation and thus laying the groundwork for mechanistic studies of epigenetic regulation of sex and morph differentiation.

Exposure to 4-nonylphenol induces a shift in the gene expression of gsdf and testis-ova formation and sex reversal in Japanese medaka (Oryzias latipes)

The branched isomer mixture 4-nonylphenol (4-NP) has been used worldwide as a surfactant, and can have endocrine-disrupting effects on aquatic organisms. For instance, 4-NP induces the formation of testis-ova (i.e., testicular and ovarian tissue in the same gonad) or male to female sex reversal of various teleost fishes. Recently, our group revealed that altered gsdf gene expression is associated with disruption of gonadal differentiation in Japanese medaka (Oryzias latipes) embryos exposed to methyltestosterone or bisphenol A, suggesting that gsdf might be useful as a biomarker for predicting the impact of endocrine-disrupting chemicals (EDCs) on gonadal differentiation. Here, we used 4-NP to examine further whether gsdf expression at the embryo stage is useful for predicting EDC impact on gonadal sex differentiation. When fertilized medaka eggs were exposed to 32 or 100 μg/L 4-NP, testis-ova in genetic males and sex reversal from genetic male to phenotypic female were observed. At stage 38 (just before hatching), 4-NP exposure at 1-100 μg/L did not affect gsdf expression in XX embryos compared with the nontreated control; however, in XY embryos, the gsdf expression in the 100 μg/L-exposed group was significantly lower than that in the controls. The 4-NP concentration at which gsdf expression was suppressed was equal to that at which testis-ova and sex reversal were induced. These results indicate that expression of the gsdf gene at the embryonic stage in medaka is a useful biomarker for predicting the impact of EDCs on sexual differentiation.