A SNP in a Steroidogenic Enzyme Is Associated with Phenotypic Sex in Seriola 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.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Turnovers of Sex-Determining Mutation in the Golden Pompano and Related Species Provide Insights into Microevolution of Undifferentiated Sex Chromosome

The suppression of recombination is considered a hallmark of sex chromosome evolution. However, previous research has identified undifferentiated sex chromosomes and sex determination by single SNP in the greater amberjack (Seriola dumerili). We observed the same phenomena in the golden pompano (Trachinotus ovatus) of the same family Carangidae and discovered a different sex-determining SNP within the same gene Hsd17b1. We propose an evolutionary model elucidating the turnover of sex-determining mutations by highlighting the contrasting dynamics between purifying selection, responsible for maintaining W-linked Hsd17b1, and neutral evolution, which drives Z-linked Hsd17b1. Additionally, sporadic loss-of-function mutations in W-linked Hsd17b1 contribute to the conversion of W chromosomes into Z chromosomes. This model was directly supported by simulations, closely related species, and indirectly by zebrafish mutants. These findings shed new light on the early stages of sex chromosome evolution.

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.

Characterization of eight types of 17β-hydroxysteroid dehydrogenases from the Japanese sardine Sardinops melanostictus: The probable role of type 12a in ovarian estradiol synthesis

17β-hydroxysteroid dehydrogenases (Hsd17bs) play a critical role in sex steroid biosynthesis. Although multiple types of Hsd17b have been found in fish, there is limited research on their expression and function. Recently, we succeeded in identifying eight types of Hsd17b (types 3, 4, 7, 8, 10, 12a, 12b, and 14) by RNA sequencing in the Japanese sardine Sardinops melanostictus, a commercially important clupeoid fish; however, a homologous sequence of Hsd17b1, which catalyzes the key reaction of estradiol-17β (E2) synthesis, was absent. Here, we aimed to identify the Hsd17b type that plays a major role in E2 synthesis during ovarian development in Japanese sardine. The cDNAs encoding those eight types of Hsd17b were cloned and sequenced. The expressions of hsd17b3, hsd17b12a, and hsd17b12b were higher in ovary than in testis. In particular, hsd17b12a was predominantly expressed in the ovary. Expression of hsd17b3, hsd17b4, hsd17b12a, and hsd17b12b in the ovary increased during ovarian development. The enzymatic activities of Hsd17b3, Hsd17b12a, and Hsd17b12b were evaluated by expressing their recombinants in human embryonic kidney 293T cells. Hsd17b12a and Hsd17b12b catalyzed the conversion of androstenedione (AD) to testosterone (T) and estrone (E1) to E2. The results of in vitro bioassays using sardine ovaries indicated that E2 is synthesized from pregnenolone via AD and T, but not E1. These results suggest that Hsd17b12a plays a major role in E2 synthesis in sardine ovary by catalyzing the conversion of AD to T.

Comparative transcriptomic analysis of the different developmental stages of ovary in the cuttlefish Sepia pharaonis

The cuttlefish, Sepia pharaonis, is characterized by rapid growth and strong disease resistance, making it an important commercially farmed cephalopod species in the southeastern coastal regions of China. However, in the reproductive process of S. pharaonis, there are challenges such as a low output of eggs, poor quality, and low survival rates of newly hatched juveniles. Therefore, there is an urgent need to study the molecular mechanisms underlying ovarian development in this species. In this study, we conducted the first transcriptomic analysis of the ovary at four developmental stages: the undeveloped stage, developing stage, nearly-ripe stage, and ripe stage, and compared the transcriptomics among these four stages using Illumina sequencing technology. The total numbers of clean reads of the four stages ranged from 40,890,772 to 52,055,714 reads. A total of 136,829 DEGs were obtained, GC base ratios of raw data were between 38.44 and 44.59%, and the number of uniquely mapped reads spanned from 88.08 to 95.90%. The Pearson correlation coefficient demonstrated a strong correlation among different samples within the same group, PCA and Anosim analysis also revealed that the grouping of these four stages was feasible, and each stage could be distinguished from the others. GO enrichment analysis demonstrated that ovarian follicle growth, sex differentiation, and transforming growth factor beta receptor, played a foreshadowing role at the early ovarian development stage, and the terms of small molecule metabolic process, peptide metabolic process, and catalytic activity were prominent at the mature stage. Meanwhile, KEGG analysis showed that the early ovarian development of S. pharaonis was mainly associated with the cell cycle, DNA replication, and carbon metabolism, while the mid-late ovarian development was involved with the signal transduction, endocrine system, and reproduction pathway. RT-qPCR further confirmed the consistent expression patterns of genes such as 17β-HSD, GH, VGS, NFR, and NYR in the ovaries of S. pharaonis, exhibiting elevated levels of expression during the maturation stage. Conversely, ER and OM exhibited high expression levels during the early stages of ovarian development. These transcriptomic data provide insights into the molecular mechanisms of S. pharaonis ovarian development. The findings of this study will contribute to improving the reproduction and development of cuttlefish and enriching the bioinformatics knowledge of cephalopods.

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.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

The male and female genomes of golden pompano (Trachinotus ovatus) provide insights into the sex chromosome evolution and rapid growth

INTRODUCTION

Golden pompano (Trachinotus ovatus) is economically significant important for offshore cage aquaculture in China and Southeast Asian countries. Lack of high-quality genomic data and accurate gene annotations greatly restricts its genetic breeding progress.

OBJECTIVES

To decode the mechanisms of sex determination and rapid growth in golden pompano and facilitate the sex- and growth-aimed genetic breeding.

METHODS

Genome assemblies of male and female golden pompano were generated using Illumina, PacBio, BioNano, genetic maps and Hi-C sequencing data. Genomic comparisons, whole genome re-sequencing of 202 F1 individuals, QTL mapping and gonadal transcriptomes were used to analyze the sex determining region, sex chromosome evolution, SNP loci, and growth candidate genes. Zebrafish model was used to investigate the functions of growth candidate gene.

RESULTS

Female (644.45 Mb) and male (652.12 Mb) genomes of golden pompano were assembled and annotated at the chromosome level. Both genomes are highly conserved and no new or highly differentiated sex chromosomes occur. A 3.5 Mb sex determining region on LG15 was identified, where Hsd17b1, Micall2 and Lmx1a were putative candidates for sex determination. Three SNP loci significantly linked to growth were pinpointed, and a growth-linked gene gpsstr1 was identified by locus BSNP1369 (G→C, 17489695, Chr23). Loss of sstr1a (homologue of gpsstr1) in zebrafish caused growth retardation.

CONCLUSION

This study provides insights into sex chromosome evolution, sex determination and rapid growth of golden pompano.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis , Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene ( amhr2bY ), previously suggested to be the master sex determining (MSD) gene in P. flavescens . Phylogenetically related and structurally similar a mhr2 duplications ( amhr2b ) were found in P. schrenkii and Sander lucioperca , potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus , this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca . Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens . In P. fluviatilis , a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three ( c18h1orf198 , hsdl1 , tbc1d32 ) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.

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.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Gonadal transcriptomes reveal sex-biased expression genes associated with sex determination and differentiation in red-tail catfish (Hemibagrus wyckioides)

BACKGROUND

Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish.

RESULTS

In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways.

CONCLUSIONS

Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

Rxrs and their partner receptor genes inducing masculinization plausibly mediated by endocrine disruption in Paralichthys olivaceus

Retinoid X receptors (RXRs) can form homo- or heterodimers with orphan receptors involved in multiple intertwined signaling pathways. However, there is limited study on the formation of sex phenotypes and the regulation of steroidogenesis by RXRs in fish. Here, in Paralichthys olivaceus, we first indicated that PPARγ::RXRα was predictably a transcription factor for steroidogenesis genes, and Foxl2 and Dmrt1 were also transcription factors for rxrs and their partner receptor genes. When the flounder fry were exposed to LG100268 (LG, RXRs agonist, 50 mg/kg diet), the percentage of males increased from 50% to 71.4%. Before histological differentiation of the flounder ovary (3 cm TL) and testis (6 cm TL), the transcripts of rar β and rar γ (P < 0.05) were activated, and the steroidogenesis gene Hsd3b1 was down-regulated (P < 0.05). The ratios of testosterone (T)/17β-estradiol (E2) were all greatly increased (P < 0.05), and the ratio of 11-ketotestosterone (11-KT)/E2 was elevated at 3 cm TL. Moreover, LG was used to treat the cultured gonads in vitro (10 μM) and the fish with intraperitoneal injection in vivo (12 mg/kg body weight), respectively. LG was able to up-regulate rxr γ, rar γ, and ppar δ, and Hsd3b1 was significantly up-regulated (P < 0.05). The ratios of 11-KT/E2 in the culture medium and in the ovaries of the fish were decreased. Furthermore, the recombinant flounder Foxl2 protein was able to significantly down-regulate ppar γ (P < 0.05) and tr β (P < 0.01) in the ovaries in vitro, and the result of the Dmrt1 in the testes was opposite to that of the Foxl2, probably indicating a feedback loop between RXRs' partner receptors and Foxl2/Dmrt1. These findings introduce for the first time the mode of action of RXRs on the flounder steroidogenesis and provide important data to learn the potential function of RXRs in fish sex differentiation and the potential role of RXRs in aquatic animals in the presence of water pollutants.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

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.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Partial sex linkage and linkage disequilibrium on the guppy sex chromosome

The guppy Y chromosome has been considered a model system for the evolution of suppressed recombination between sex chromosomes, and it has been proposed that complete sex-linkage has evolved across about 3 Mb surrounding this fish's sex-determining locus, followed by recombination suppression across a further 7 Mb of the 23 Mb XY pair, forming younger "evolutionary strata". Sequences of the guppy genome show that Y is very similar to the X chromosome. Knowing which parts of the Y are completely nonrecombining, and whether there is indeed a large completely nonrecombining region, are important for understanding its evolution. Here, we describe analyses of PoolSeq data in samples from within multiple natural populations from Trinidad, yielding new results that support previous evidence for occasional recombination between the guppy Y and X. We detected recent demographic changes, notably that downstream populations have higher synonymous site diversity than upstream ones and other expected signals of bottlenecks. We detected evidence of associations between sequence variants and the sex-determining locus, rather than divergence under a complete lack of recombination. Although recombination is infrequent, it is frequent enough that associations with SNPs can suggest the region in which the sex-determining locus must be located. Diversity is elevated across a physically large region of the sex chromosome, conforming to predictions for a genome region with infrequent recombination that carries one or more sexually antagonistic polymorphisms. However, no consistently male-specific variants were found, supporting the suggestion that any completely sex-linked region may be very small.

Y-specific amh allele, amhy, is the master sex-determining gene in Japanese flounder Paralichthys olivaceus

Japanese flounder (Paralichthys olivaceus) is an important marine fish species of both fisheries and aquaculture in Northeast Asia. The commercial interest for all-female progenies due to several sex-related traits has prompted basic research on the mechanisms of sex determination in this species. By conducting a linkage analysis of the sex-determining locus, we initially identified 12 microsatellite markers linked to sex in 11 scaffolds, whose localization was restricted to a specific region of linkage group 9. Sequence analysis of this region identified 181 genes based on the UniProt database annotations. Among them, the amh gene was considered a potential candidate for sex determination because this gene is known to have taken over the role of sex determination in many teleosts. An in-depth sequence analysis of both the coding and non-coding regions of amh in XX and XY individuals detected nine SNPs linked with maleness. However, because these substitutions were synonymous, the upstream and downstream regions of amh were also investigated and a male-specific variant with deletions in the promoter region was detected. This truncated Y-specific amh variant was named amhy, and the amh shared by both sexes was named amhx. The association analysis using both females and males of the genotypic sex inferred by the presence/absence of amhy found complete association with phenotypic sex and genotype. Gene expression analysis in larvae derived from a single-pair progeny by quantitative real-time PCR detected amhy transcripts in the larval trunks between 20 and 100 days after hatching only in XY larvae. Localization of amhy by in situ hybridization was detected in presumptive Sertoli cells of XY gonads. Expression of amhx was almost undetectable in both XX and XY genotypes. Loss of Amh function by CRISPR-Cas9 induced male-to-female sex reversal, indicating that this gene was necessary for the masculinization of XY individuals. In conclusion, the complete linkage of amhy with males, its early expression in XY gonads before testicular differentiation, and the induction of sex reversal by loss-of-function mutation support the view that amhy is the sex-determining gene in this species.

Transcriptome Profiling and Expression Localization of Key Sex-Related Genes in a Socially-Controlled Hermaphroditic Clownfish, Amphiprion clarkii

Clownfish can be an excellent research model for investigating the socially-controlled sexual development of sequential hermaphrodite teleosts. However, the molecular cascades underlying the social cues that orchestrate the sexual development process remain poorly understood. Here, we performed a comparative transcriptomic analysis of gonads from females, males, and nonbreeders of Amphiprion clarkii, which constitute a complete social group, allowing us to investigate the molecular regulatory network under social control. Our analysis highlighted that the gonads of nonbreeders and males exhibited high similarities but were far from females, both in global transcriptomic profiles and histological characteristics, and identified numerous candidate genes involved in sexual development, some well-known and some novel. Significant upregulation of cyp19a1a, foxl2, nr5a1a, wnt4a, hsd3b7, and pgr in females provides strong evidence for the importance of steroidogenesis in ovarian development and maintenance, with cyp19a1a playing a central role. Amh and sox8 are two potential key factors that may regulate testicular tissue development in early and late stages, respectively, as they are expressed at higher levels in males than in females, but with slightly different expression timings. Unlike previous descriptions in other fishes, the unique expression pattern of dmrt1 in A. clarkii implied its potential function in both male and female gonads, and we speculated that it might play promoting roles in the early development of both testicular and ovarian tissues.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

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.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

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.

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

Some thoughts about the words we use for thinking about sex chromosome evolution

Sex chromosomes are familiar to most biologists since they first learned about genetics. However, research over the past 100 years has revealed that different organisms have evolved sex-determining systems independently. The differences in the ages of systems, and in how they evolved, both affect whether sex chromosomes have evolved. However, the diversity means that the terminology used tends to emphasize either the similarities or the differences, sometimes causing misunderstandings. In this article, I discuss some concepts where special care is needed with terminology. The following four terms regularly create problems: ‘sex chromosome’, ‘master sex-determining gene’, ‘evolutionary strata’ and ‘genetic degeneration’. There is no generally correct or wrong use of these words, but efforts are necessary to make clear how they are to be understood in specific situations. I briefly outline some widely accepted ideas about sex chromosomes, and then discuss these ‘problem terms’, highlighting some examples where careful use of the words helps bring to light current uncertainties and interesting questions for future work.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Morphological and Molecular Gonadal Sex Differentiation in the Wild Japanese eel Anguilla japonica

Most cultured Japanese eels (Anguilla japonica) show male sex differentiation; however, natural gonadal sex differentiation has not been evaluated. In this study, this process was characterized in wild eels. Differentiated ovaries and testes were observed after the eels grew to 320 and 300 mm in total length, respectively. The youngest ovary and testis appeared at 3 and 4 years old, respectively; however, undifferentiated gonads were found up to 7 years, suggesting that sex differentiation was triggered by growth rather than aging. gsdf, amh, foxl2b and foxl3b were highly expressed in the testes, whereas figla, sox3, foxn5, zar1, and zp3 were highly expressed in the ovaries. The expression of cyp19a1a and foxl2a did not differ significantly between the testis and ovary. In the ovaries, the cyp19a1a and foxl2a levels were highest in the early stages, suggesting that their function is limited to early ovarian differentiation. The foxn5, zar1 and zp3 levels tended to increase in the later stages, suggesting that they function after the initiation of ovarian differentiation. In undifferentiated gonads, dimorphic gene expression was not observed, suggesting that the molecular sex differentiation phase is short and difficult to detect. These findings provide the first demonstration of the whole course of natural gonadal sex differentiation in eels at molecular and morphological levels.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

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.

Prediction of the Sex-Associated Genomic Region in Tunas (Thunnus Fishes)

Fish species have a variety of sex determination systems. Tunas (genus Thunnus) have an XY genetic sex determination system. However, the Y chromosome or responsible locus has not yet been identified in males. In a previous study, a female genome of Pacific bluefin tuna (T. orientalis) was sequenced, and candidates for sex-associated DNA polymorphisms were identified by a genome-wide association study using resequencing data. In the present study, we sequenced a male genome of Pacific bluefin tuna by long-read and linked-read sequencing technologies and explored male-specific loci through a comparison with the female genome. As a result, we found a unique region carrying the male-specific haplotype, where a homolog of estrogen sulfotransferase gene was predicted to be encoded. The genome-wide mapping of previously resequenced data indicated that, among the functionally annotated genes, only this gene, named sult1st6y, was paternally inherited in the males of Pacific bluefin tuna. We reviewed the RNA-seq data of southern bluefin tuna (T. maccoyii) in the public database and found that sult1st6y of southern bluefin tuna was expressed in all male testes, but absent or suppressed in the female ovary. Since estrogen sulfotransferase is responsible for the inactivation of estrogens, it is reasonable to assume that the expression of sult1st6y in gonad cells may inhibit female development, thereby inducing the individuals to become males. Thus, our results raise a promising hypothesis that sult1st6y is the sex determination gene in Thunnus fishes or at least functions at a crucial point in the sex-differentiation cascade.

A single intronic single nucleotide polymorphism in splicing site of steroidogenic enzyme hsd17b1 is associated with phenotypic sex in oyster pompano, Trachinotus anak

Teleosts show varied master sex determining (MSD) genes and sex determination (SD) mechanisms, with frequent turnovers of sex chromosomes. Tracing the origins of MSD genes and turnovers of sex chromosomes in a taxonomic group is of particular interest in evolutionary biology. Oyster pompano (Trachinotus anak), a marine fish, belongs to the family Carangidae, in which 17b-hydroxysteroid dehydrogenase 1 (hsd17b1) has repeatedly evolved to an MSD gene. Whole-genome resequencing identified a single nucleotide polymorphism (SNP) at chromosome 24 to be strictly associated with phenotypic sex, with females being the heterozygous sex. This SNP is located in a splicing site at the first exon/intron boundary of hsd17b1. The Z-linked SNP results in malfunction of all spliced isoforms, whereas the W-linked isoforms were predicted to have open reading frames that are conserved among vertebrates, suggesting that hsd17b1 is a female-determining gene. The differential alternative splicing patterns of ZZ and ZW genotypes were consistently observed both in undifferentiated stages and differentiated gonads. We observed elevated recombination around the SD locus and no differentiation between Z and W chromosomes. The extreme diversity of mutational mechanisms that hsd17b1 evolves to an MSD gene highlights frequent in situ turnovers between sex chromosomes in the Carangidae.

The genome of New Zealand trevally (Carangidae: Pseudocaranx georgianus) uncovers a XY sex determination locus

BACKGROUND

The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers.

RESULTS

The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females).

CONCLUSIONS

The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

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 duplicated amh is the master sex-determining gene for Sebastes rockfish in the Northwest Pacific

Teleost fish are the most diverse group of vertebrates and provide opportunities to study the evolution of sex determination (SD) systems. Using genomic and functional analyses, we identified a male-specific duplication of anti-Müllerian hormone (amh) gene as the male master sex-determining (MSD) gene in Sebastes schlegelii. By resequencing 10 males and 10 females, we characterized a 5 kb-long fragment in HiC_Scaffold_12 as a male-specific region, which contained an amh gene (named amhy). We then demonstrated that amhy is a duplication of autosomal amh that was later translocated to the ancestral Y chromosome. amha and amhy shared high-nucleotide identity with the most significant difference being two insertions in intron 4 of amhy. Furthermore, amhy overexpression triggered female-to-male sex reversal in S. schlegelii, displaying its fundamental role in driving testis differentiation. We developed a PCR assay which successfully identified sexes in two species of northwest Pacific rockfish related to S. schlegelii. However, the PCR assay failed to distinguish the sexes in a separate clade of northeast Pacific rockfish. Our study provides new examples of amh as the MSD in fish and sheds light on the convergent evolution of amh duplication as the driving force of sex determination in different fish taxa.

Lessons from an unusual vertebrate sex-determining gene

So far, very few sex-determining genes have been identified in vertebrates and most of them, the so-called 'usual suspects', evolved from genes which fulfil essential functions during sexual development and are thus already tightly linked to the process that they now govern. The single exception to this 'usual suspects' rule in vertebrates so far is the conserved salmonid sex-determining gene, sdY (sexually dimorphic on the Y chromosome), that evolved from a gene known to be involved in regulation of the immune response. It is contained in a jumping sex locus that has been transposed or translocated into different ancestral autosomes during the evolution of salmonids. This special feature of sdY, i.e. being inserted in a 'jumping sex locus', could explain how salmonid sex chromosomes remain young and undifferentiated to escape degeneration. Recent knowledge on the mechanism of action of sdY demonstrates that it triggers its sex-determining action by deregulating oestrogen synthesis that is a conserved and crucial pathway for ovarian differentiation in vertebrates. This result suggests that sdY has evolved to cope with a pre-existing sex differentiation regulatory network. Therefore, 'limited options' for the emergence of new master sex-determining genes could be more constrained by their need to tightly interact with a conserved sex differentiation regulatory network rather than by being themselves 'usual suspects', already inside this sex regulatory network. 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)'.

Evolution of master sex determiners: TGF-β signalling pathways at regulatory crossroads

To date, more than 20 different vertebrate master sex-determining genes have been identified on different sex chromosomes of mammals, birds, frogs and fish. Interestingly, six of these genes are transcription factors (Dmrt1- or Sox3- related) and 13 others belong to the TGF-β signalling pathway (Amh, Amhr2, Bmpr1b, Gsdf and Gdf6). This pattern suggests that only a limited group of factors/signalling pathways are prone to become top regulators again and again. Although being clearly a subordinate member of the sex-regulatory network in mammals, the TGF-β signalling pathway made it to the top recurrently and independently. Facing this rolling wave of TGF-β signalling pathways, this review will decipher how the TGF-β signalling pathways cope with the canonical sex gene regulatory network and challenge the current evolutionary concepts accounting for the diversity of sex-determining mechanisms. 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)'.

A novel c.1759T>G variant in follicle-stimulating hormone-receptor gene is concordant with male determination in the flathead grey mullet (Mugil cephalus)

Various master key regulators (MKRs) that control a binary switch of sex determination (SD) have been found in fish; these provide an excellent model for the study of vertebrate genetic SD. The SD region in flathead grey mullet has been previously mapped to a 1 Mbp region harboring 27 genes, of which one is follicle-stimulating hormone receptor (fshr). Although this gene is involved in gonad differentiation and function, it has not been considered as an MKR of SD. We systematically investigated polymorphism in mullet fshr using DNA shotgun sequences, and compared them between males and females. Capable of encoding nonconservative amino acid substitutions, c.1732G>A and c.1759T>G exhibited association with sex on a population level (N = 83; P ≤ 6.7 × 10-19). Hence, 1732 A and 1759 G represent a male-specific haplotype of the gene, designated as "fshry." Additional flanking SNPs showed a weaker degree of association with sex, delimiting the SD critical region to 143 nucleotides on exon 14. Lack of homozygotes for fshry, and the resulting divergence from Hardy-Weinberg equilibrium (N = 170; P ≤ 3.9 × 10-5), were compatible with a male heterogametic model (XY/XX). Capable of replacing a phenylalanine with valine, c.1759T>G alters a conserved position across the sixth transmembrane domain of vertebrate FSHRs. Amino acid substitutions in this position in vertebrates are frequently associated with constant receptor activation and consequently with FSH/FSHR signaling alteration; thus, indicating a potential role of fshr as an MKR of SD.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BACKGROUND

Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.

RESULTS

The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.

CONCLUSIONS

The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Establishment of a new sex-determining allele driven by sexually antagonistic selection

The turnover of sex-determining loci has repeatedly occurred in a number of species, rather than having a diverged pair of sex chromosomes. We model the turnover process by considering a linked locus under sexually antagonistic selection. The entire process of a turnover may be divided into two phases, which are referred to as the stochastic and deterministic phases. The stochastic phase is when a new sex-determining allele just arises and is still rare and random genetic drift plays an important role. In the deterministic phase, the new allele further increases in frequency by positive selection. The theoretical results currently available are for the deterministic phase, which demonstrated that a turnover of a newly arisen sex-determining locus could benefit from selection at a linked locus under sexually antagonistic selection, by assuming that sexually antagonistic selection works in a form of balancing selection. In this work, we provide a comprehensive theoretical description of the entire process from the stochastic phase to the deterministic phase. In addition to balancing selection, we explore several other modes of selection on the linked locus. Our theory allows us make a quantitative argument on the rate of turnover and the effect of the mode of selection at the linked locus. We also performed simulations to explore the pattern of polymorphism around the new sex-determining locus. We find that the pattern of polymorphism is informative to infer how selection worked through the turnover process.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

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.

Activation of stress response axis as a key process in environment-induced sex plasticity in fish

The determination of sex is an important hallmark in the life cycle of organisms, in which the fate of gonads and then the individual sex are defined. In gonochoristic teleost fish, this process is characterized by a high plasticity, considering that in spite of genotypic sex many environmental factors can cause shifts from one to another molecular pathway, resulting in organisms with mismatching genotypic and phenotypic sexes. Interestingly, in most instances, both female-to-male or male-to-female sex-reversed individuals develop functional gonads with normal gametogenesis and respective progenies with full viability. The study of these mechanisms is being spread to other non-model species or to those inhabiting more extreme environmental conditions. Although water temperature is an important mechanism involved in sex determination, there are other environmental stressors affected by the climate change which are also implicated in stress response-induced masculinization in fish. In this regard, the brain has emerged as the transducer of the environment input that can influence the gonadal fate. Furthermore, the evaluation of other environmental stressors or their synergic effect on sex determination at conditions that simulate the natural environments is growing gradually. Within such scope, the concerns related to climate change impacts rely on the fact that many of biotic and abiotic parameters reported to affect sex ratios are expected to increase concomitantly as a result of increased greenhouse gas emissions and, particularly worrying, many of them are related to male bias in the populations, such as high temperature, hypoxia, and acidity. These environmental changes can also generate epigenetic changes in sex-related genes affecting their expression, with implications on sex differentiation not only of exposed individuals but also in following generations. The co-analysis of multi-stressors with potential inter- and transgenerational effects is essential to allow researchers to perform long-term predictions on climate change impacts in wild populations and for establishing highly accurate monitoring tools and suitable mitigation strategies.

The Ficus erecta genome aids Ceratocystis canker resistance breeding in common fig (F. carica)

Ficus erecta, a wild relative of the common fig (F. carica), is a donor of Ceratocystis canker resistance in fig breeding programmes. Interspecific hybridization followed by recurrent backcrossing is an effective method to transfer the resistance trait from wild to cultivated fig. However, this process is time consuming and labour intensive for trees, especially for gynodioecious plants such as fig. In this study, genome resources were developed for F. erecta to facilitate fig breeding programmes. The genome sequence of F. erecta was determined using single-molecule real-time sequencing technology. The resultant assembly spanned 331.6 Mb with 538 contigs and an N50 length of 1.9 Mb, from which 51 806 high-confidence genes were predicted. Pseudomolecule sequences corresponding to the chromosomes of F. erecta were established with a genetic map based on single nucleotide polymorphisms from double-digest restriction-site-associated DNA sequencing. Subsequent linkage analysis and whole-genome resequencing identified a candidate gene for the Ceratocystis canker resistance trait. Genome-wide genotyping analysis enabled the selection of female lines that possessed resistance and effective elimination of the donor genome from the progeny. The genome resources provided in this study will accelerate and enhance disease-resistance breeding programmes in fig.

Characteristics and sex dimorphism of 17β-hydroxysteroid dehydrogenase family genes in the olive flounder Paralichthys olivaceus

Sex steroid hormones play important roles in fish sex differentiation, gonadal development and secondary sexual characteristics. Olive flounder Paralichthys olivaceus is a valuable commercial marine fish species and has marked sexual dimorphism. However, the mechanisms of action of sex hormones in flounder sex are still unclear. In this study, a total of ten Hsd17b family genes, including Hsd17b3, -4, -7, -8, -9, -10, -12a, -12b, -14 and -15, were identified in the flounder, which encoded critical enzymes acting on sex steroid synthesis and metabolism. Hsd17b genes were distributed on eight chromosomes. Hsd17b12a and -12b were located on chromosomes 19 and 7, respectively. It was speculated that these two genes were just highly similar rather than different transcripts derived from the same gene. According to the results of domain and motif analyses, they all belonged to the SDR superfamily and contained conserved Hsd17b motifs TGxxxGxG, PGxxxT, NNAG and YxxxK. Analysis of amino acid sequences predicted that Hsd17b1, -4, -7, -12a and -14 were hydrophilic proteins. The stability of Hsd17b1, -3 and -12b proteins was predicted to be low. The various Hsd17b family genes differed in tissue expression pattern, and Hsd17b10, -12a and -12b were highly expressed in the flounder ovary. Moreover, throughout gonadal development, Hsd17b3 was highly expressed in the testis, and Hsd17b1, -12a and -12b were highly expressed in the ovary, suggesting that they might play an important role in testosterone synthesis in the testis or estrogen synthesis in the ovary. Activities of Hsd17b3 at stages I-V were all significantly higher in the testis than in the ovary (P < 0.05, P < 0.01). Transfection analysis in HEK293T cells showed that Hsd17b1 and -3 were located in both the cytoplasm and nucleus. Additionally, after challenging fish with tamoxifen, Hsd17b3 expression level in the testis decreased significantly (P < 0.01), and in the ovary no significant change was observed. Moreover, the expression of Hsd17b1 in the ovary was significantly upregulated after injection with flutamide (P < 0.05). These findings introduce the characteristics of the flounder Hsd17b in subfamily, which contribute to our understanding of the regulation of sex steroid hormone synthesis in fish gonadal development.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.)

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.