Autosomal gsdf acts as a male sex initiator in the fish medaka

The potential regulatory role of non-coding RNAs in mifepristone-induced masculinization in Takifugu rubripes gonads

BACKGROUND

The regulatory roles of non-coding RNAs (ncRNAs) during sex differentiation in teleosts have received widespread attention recently. Mifepristone (RU486, a progesterone antagonist), which acts as an endocrine disruptor, can affect reproduction and sex differentiation in teleosts.

RESULTS

The expression of ncRNAs in the gonads of tiger puffer (Takifugu rubripes) during RU486 (500 µg/g diet) induced masculinization process was examined by RNA-sequencing. A total of 4,381 long non-coding RNAs (lncRNAs), 309 circular RNAs (circRNAs), and 1,020 microRNAs (miRNAs) were identified. The expression of 41 differentially expressed (DE) lncRNAs and 20 DE miRNAs, which showed sexual dimorphic expression patterns in genetic female gonads in the control group (C-XX) vs. genetic male gonads in the control group (C-XY), were altered in genetic female gonads in the RU486 treated group (RU-XX). The genes targeted by DE ncRNAs were mainly enriched in sex-related pathways, such as calcium signaling, ovarian steroidogenesis, and cortisol synthesis and secretion. The results of co-expression and competing endogenous RNA (ceRNA) network analysis indicated that miRNAs (e.g., miR-205-z and fru-miR-122) and lncRNAs (including XR_003890915.1 and XR_003885862.1) may have pivotal roles, and lncRNAs (including XR_003890295.1, MSTRG.11750.1, and XR_003888827.1) may act as miRNA sponges, involved in the competition between miRNAs and sex-related genes during tiger puffer masculinization process. Dual luciferase reporter assay results identified that ovarian steroidogenesis related gene hsd17b1 is a downstream target of fru-miR-122. The expression of 4 lncRNAs, 4 circRNAs, and 6 miRNAs were validated by qPCR, indicating the accuracy and dependability of RNA-Seq.

CONCLUSIONS

This study provided the evidence that ncRNAs may participate in RU486-induced masculinization in T. rubripes, and may enhance our understanding of the regulatory network of sex differentiation in fugu.

A detailed transcriptome study uncovers the epigenetic characteristics associated with Aromatase inhibitor-induced masculinization in Takifugu rubripes larvae gonads

BACKGROUND

Takifugu rubripes is an economically valuable fish species in Asia. The implementation of all-male culture for T. rubripes is highly anticipated in aquaculture. Aromatase inhibitor (AI, letrozole) treatment was found to be an efficient method to induced masculinization in T. rubripes, as reported in our previous study. Here, to further explore the underlying regulation mechanism of AI-induced masculinization, a whole-transcriptome analysis comparing was conducted between AI-induced masculinized XX (AI-XX) gonads and control (Con) gonads in T. rubripes.

RESULTS

In Con-XX/Con-XY comparison, 1,172 differential expression (DE) mRNAs, 129 DEmiRNAs, 210 DElncRNAs, and 4 DEcircRNAs were identified. In the Con-XX/AI-XX comparison, 1,329 DEmRNAs, 174 DEmiRNAs, 6 DEcircRNAs and 280 DElncRNAs were found. Con-XX/Con-XY and Con-XX/AI-XX comparisons shared 690 DEmRNAs, 50 DEmiRNAs, 3 DEcircRNAs, and 105 DElncRNAs. The analyses of protein-protein interaction (PPI) and competitive endogenous RNA (ceRNA) network identified interactions among these shared DERNAs. Kcnh2b, trim27, cnnm2b, reln, cckbra, pkd1l2, steap4, gsg1l, hamp, and foxg1c were predicted as the top ten of hub genes. miRNAs included miRNA-27 family and miRNA-489 family showed targeting relationship with hub genes. GO and KEGG functional enrichment analysis showed that the targeted genes were mainly enriched in GO:0065008 regulation of biological quality and TGF-beta signaling pathway. qPCR validation confirmed the differential expression of selected mRNAs, and ncRNAs.

CONCLUSIONS

This research comprehensively reveals the potential regulatory effects of ncRNAs on cellular motility, fate regulation, and hormonal regulation during gonadal masculinization in T. rubripes. It may provide significant insights into the regulation mechanisms underlying sex reversal in fish.

Transient receptor potential a1b regulates primordial germ cell numbers and sex differentiation in developing zebrafish

Temperature is a leading environmental factor determining the sex ratio of some animal populations, such as fish, amphibians, and reptiles. However, the underlying mechanism by which temperature affects gender is still poorly understood. Transient receptor potential a1b (Trpa1b) belongs to the ion channel family of transient receptor potentials and exhibits dual thermosensitivity to heat and cold. In this study, we have unveiled a novel function of the trpa1b gene. Zebrafish generated through clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 with Trpa1b-null manifest a male-biased sex ratio. The quantity of primordial germ cells (PGCs) in zebrafish is closely linked to gender determination and gonadal development. Yet the role of the trpa1b gene in zebrafish reproductive development remains unexplored in the literature. Our investigation revealed a significant reduction in PGCs in Trpa1b mutant zebrafish compared to their wild-type counterparts 24-h postfertilization (hpf). Transcriptome sequencing of tissues near the reproductive crest of embryos at 1.25 days postfertilization (dpf) revealed differential changes in PGC-related marker genes and genes related to sperm cell development and differentiation. The relative expression of ddx4 and sycp3 genes was significantly downregulated, whereas amh was significantly upregulated at 20 dpf in trpa1b-/- zebrafish. The results of this study provide valuable insights and references for studying the molecular mechanism of sex determination in zebrafish. Undoubtedly, these results will further enhance our understanding of gender differentiation and gonadal development in fish and other vertebrates.

Meta-analysis of gonadal transcriptome provides novel insights into sex change mechanism across protogynous fishes

Protogyny, being capable of changing from female to male during their lifetime, is prevalent in 20 families of teleosts but is believed to have evolved within specific evolutionary lineages. Therefore, shared regulatory factors governing the sex change process are expected to be conserved across protogynous fishes. However, a comprehensive understanding of this mechanism remains elusive. To identify these factors, we conducted a meta-analysis using gonadal transcriptome data from seven species. We curated data pairs of ovarian tissue and transitional gonad, and employed ratios of expression level as a unified criterion for differential expression, enabling a meta-analysis across species. Our approach revealed that classical sex change-related genes exhibited differential expression levels between the ovary and transitional gonads, consistent with previous reports. These results validate our methodology's robustness. Additionally, we identified novel genes not previously linked to gonadal sex change in fish. Notably, changes in the expression levels of acetoacetyl-CoA synthetase and apolipoprotein Eb, which are involved in cholesterol synthesis and transport, respectively, suggest that the levels of cholesterol, a precursor of steroid hormones crucial for sex change, are decreased upon sex change onset in the gonads. This implies a potential universal influence of cholesterol dynamics on gonadal transformation in protogyny.

Histological observations and transcriptome analyses reveal the dynamic changes in the gonads of the blotched snakehead (Channa maculata) during sex differentiation and gametogenesis

BACKGROUND

Blotched snakehead (Channa maculata) displays significant sexual dimorphism, with males exhibiting faster growth rates and larger body sizes compared to females. The cultivation of the all-male population of snakeheads holds substantial economic and ecological value. Nonetheless, the intricate processes governing the development of bipotential gonads into either testis or ovary in C. maculata remain inadequately elucidated. Therefore, it is necessary to determine the critical time window of sex differentiation in C. maculata, providing a theoretical basis for sex control in production practices.

METHODS

The body length and weight of male and female C. maculata were measured at different developmental stages to reveal when sexual dimorphism in growth initially appears. Histological observations and spatiotemporal comparative transcriptome analyses were performed on ovaries and testes across various developmental stages to determine the crucial time windows for sex differentiation in each sex and the sex-related genes. Additionally, qPCR and MG2C were utilized to validate and locate sex-related genes, and levels of E2 and T were quantified to understand sex steroid synthesis.

RESULTS

Sexual dimorphism in growth became evident starting from 90 dpf. Histological observations revealed that morphological sex differentiation in females and males occurred between 20 and 25 dpf or earlier and 30-35 dpf or earlier, respectively, corresponding to the appearance of the ovarian cavity or efferent duct anlage. Transcriptome analyses revealed divergent gene expression patterns in testes and ovaries after 30 dpf. The periods of 40-60 dpf and 60-90 dpf marked the initiation of molecular sex differentiation in females and males, respectively. Male-biased genes (Sox11a, Dmrt1, Amh, Amhr2, Gsdf, Ar, Cyp17a2) likely play crucial roles in male sex differentiation and spermatogenesis, while female-biased genes (Foxl2, Cyp19a1a, Bmp15, Figla, Er) could be pivotal in ovarian differentiation and development. Numerous biological pathways linked to sex differentiation and gametogenesis were also identified. Additionally, E2 and T exhibited sexual dimorphism during sex differentiation and gonadal development. Based on these results, it is hypothesized that in C. maculata, the potential male sex differentiation pathway, Sox11a-Dmrt1-Sox9b, activates downstream sex-related genes (Amh, Amhr2, Gsdf, Ar, Cyp17a2) for testicular development, while the antagonistic pathway, Foxl2/Cyp19a1a, activates downstream sex-related genes (Bmp15, Figla, Er) for ovarian development.

CONCLUSIONS

This study provides a comprehensive overview of gonadal dynamic changes during sex differentiation and gametogenesis in C. maculata, establishing a scientific foundation for sex control in this species.

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.

Theoretical Analysis and Expression Profiling of 17β-Hydroxysteroid Dehydrogenase Genes in Gonadal Development and Steroidogenesis of Leopard Coral Grouper (Plectropomus leopardus)

The differentiation and developmental trajectory of fish gonads, significantly important for fish breeding, culture, and production, has long been a focal point in the fields of fish genetics and developmental biology. However, the mechanism of gonadal differentiation in leopard coral grouper (Plectropomus leopardus) remains unclear. This study investigates the 17β-Hydroxysteroid Dehydrogenase (Hsd17b) gene family in P. leopardus, with a focus on gene characterization, expression profiling, and functional analysis. The results reveal that the P. leopardus's Hsd17b gene family comprises 11 members, all belonging to the SDR superfamily. The amino acid similarity is only 12.96%, but conserved motifs, such as TGxxxGxG and S-Y-K, are present in these genes. Hsd17b12a and Hsd17b12b are unique homologs in fish, and chromosomal localization has confirmed that they are not derived from different transcripts of the same gene, but rather are two independent genes. The Hsd17b family genes, predominantly expressed in the liver, heart, gills, kidneys, and gonads, are involved in synthesizing or metabolizing sex steroid hormones and neurotransmitters, with their expression patterns during gonadal development categorized into three distinct categories. Notably, Hsd17b4 and Hsd17b12a were highly expressed in the testis and ovary, respectively, suggesting their involvement in the development of reproductive cells in these organs. Fluorescence in situ hybridization (FISH) further indicated specific expression sites for these genes, with Hsd17b4 primarily expressed in germ stem cells and Hsd17b12a in oocytes. This comprehensive study provides foundational insights into the role of the Hsd17b gene family in gonadal development and steroidogenesis in P. leopardus, contributing to the broader understanding of fish reproductive biology and aquaculture breeding.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

Genome-wide identification, evolution and expression of TGF-β signaling pathway members in mandarin fish (Siniperca chuatsi)

Members of the transforming growth factor β (TGF-β) signaling pathway regulate diverse cellular biological processes in embryonic and tissue development. We took mandarin fish (Siniperca chuatsi) as the research object to identify all members of the TGF-β signaling pathway, measure their expression pattern in the key period post hatching, and further explore their possible role in the process of sex regulation. Herein, we identified eighty-three TGF-β signaling pathway members and located them on chromosomes based on the genome of mandarin fish. TGF-β signaling pathway members were highly conserved since each TGF-β subfamily clustered with orthologs from other species. Transcriptome analysis, qRT-PCR and in situ hybridization demonstrated that most mandarin fish TGF-β signaling pathway members presented stage-specific and/or sex-dimorphic expression during gonadal development, and different members of the TGF-β signaling pathway participated in different stages of gonadal development. Taken together, our results provide new insight into the role of TGF-β signaling pathway members in the sex regulation of mandarin 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 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.

Transcriptome characterization of gonadal sex differentiation in Pacific bluefin tuna, Thunnus orientalis (Temminck et Schlegel)

Tunas (genus Thunnus) are one of the most ecologically and commercially important fish worldwide. To establish a biological basis for reproduction in this globally essential species, we have recently studied crucial reproductive aspects of the Pacific bluefin tuna (T. orientalis; PBT), as a model of tuna species, based on our closed-cycle aquaculture technology. In this study, we clarified the global expression profile of the genes regulating gonadal sex differentiation in PBT, as this developmental process is vital to sexual reproduction. Based on the results of our comparative (RNA-sequencing) and temporal (qRT-PCR) transcriptome analyses using the updated genome dataset, we propose the molecular mechanisms of gonadal sex differentiation in PBT. In female gonads, foxl2 and cyp19a1a (coding aromatase) are expressed at the onset of sex differentiation. Active aromatase-mediated estrogen biosynthesis, which includes positive regulation of cyp19a1a expression by Foxl2, induces ovarian differentiation. By contrast, dmrt1 and gsdf are upregulated in differentiating male gonads lacking active estrogen synthesis. Dmrt1 and Gsdf would mainly promote testicular differentiation. Furthermore, androgen biosynthesis is upregulated in differentiating male gonad. Endogenous androgens may also be vital to testicular differentiation. This study provides the first comprehensive data clarifying the molecular basis for gonadal sex differentiation in tunas.

Levonorgestrel causes feminization and dose-dependent masculinization in medaka fish (Oryzias latipes): Endocrine-disruption activity and its correlation with sex reversal

The effect of a synthetic progestin, levonorgestrel (LNG), on the sex of exposed embryos was examined in medaka fish (Oryzias latipes). The aims of this study are to clarify the dual effect of LNG on sex and the correlation with its androgenic/estrogenic potential in medaka. LNG exposure causes significant dose-dependent masculinization (0.1-100 μg/L), whereas a decrease in the masculinization ratio is observed at 100 μg/L. LNG also causes significant feminization at 1-100 μg/L, but not in a dose-dependent manner. Exposure of estrogen-responsive gene (choriogeninH-EGFP) transgenic embryos to 100 μg/L LNG produced significant fluorescent signals in hatched fry. In vitro transcriptional assays indicated that LNG at 10^-7-10^-5 M induced significant activity for estrogen receptor (ESR)2a and ESR2b, but not for ESR1. In pre-self-feeding fry at 5 days post hatching (dph), 1-100 μg/L LNG caused a significant increase in the mRNA of choriogeninH, irrespective of genetic sex. Moreover, LNG (10^-10-10^-5 M) also caused a significant increase in the transcriptional activity of androgen receptor (AR) α and ARβ in vitro, and 0.1 μg/L LNG significantly increased the mRNA levels of a testis-differentiation initiation factor, gonadal soma-derived factor (gsdf), as an androgen-upregulated and estrogen-downregulated gene, in 5 dph XX fry to levels similar to those in the control XY fry. However, 100 and 10 μg/L LNG suppressed or did not induce gsdf mRNA expression in XY and XX fry, respectively. Together, these findings show that LNG exerts estrogenic and androgenic activities in different concentration ranges, which correlate with the ratio of LNG-induced sex reversal. These results suggest for the first time, that medaka exposure to LNG can induce masculinization and feminization, based on the balance between androgenic and estrogenic activities, and the protocol applied in this study represents an alternative to the traditional animal model used to screen for endocrine-disrupting potential.

Development of Ovaries and Sex Change in Fish: Bringing Potential into Action

BACKGROUND

Encompassing about half of the 60,000 species of vertebrates, fish display the greatest diversity of sex determination mechanisms among metazoans. As such that phylum offers a unique playground to study the impressive variety of gonadal morphogenetic strategies, ranging from gonochorism, with either genetic or environmental sex determination, to unisexuality, with either simultaneous or consecutive hermaphroditism.

SUMMARY

From the two main types of gonads, the ovaries embrace the important role to produce the larger and non-motile gametes, which is the basis for the development of a future organism. The production of the egg cells is complex and involves the formation of follicular cells, which are necessary for the maturation of the oocytes and the production of feminine hormones. In this vein, our review focuses on the development of ovaries in fish with special emphasis on the germ cells, including those that transition from one sex to the other as part of their life cycle and those that are capable of transitioning to the opposite sex depending on environmental cues.

KEY MESSAGES

Clearly, establishing an individual as either a female or a male is not accomplished by the sole development of two types of gonads. In most cases, that dichotomy, be it final or transient, is accompanied by coordinated transformations across the entire organism, leading to changes in the physiological sex as a whole. These coordinated transformations require both molecular and neuroendocrine networks, but also anatomical and behavioural adjustments. Remarkably, fish managed to tame the ins and outs of sex reversal mechanisms to take the most advantages of changing sex as adaptive strategies in some situations.

Soy Isoflavones Induce Feminization of Japanese Eel (Anguilla japonica)

Under aquaculture conditions, Japanese eels (Anguilla japonica) produce a high percentage of males. However, females gain higher body weight and have better commercial value than males, and, therefore, a high female ratio is required in eel aquaculture. In this study, we examined the effects of isoflavones, genistein, and daidzein on sex differentiation and sex-specific genes of eels. To investigate the effects of these phytoestrogens on the gonadal sex, we explored the feminizing effects of soy isoflavones, genistein, and daidzein in a dose-dependent manner. The results showed that genistein induced feminization more efficiently than daidzein. To identify the molecular mechanisms of sex-specific genes, we performed a comprehensive expression analysis by quantitative real-time PCR and RNA sequencing. Phenotypic males and females were produced by feeding elvers a normal diet or an estradiol-17β- or genistein-treated diet for 45 days. The results showed that female-specific genes were up-regulated and male-specific genes were down-regulated in the gonads, suggesting that genistein induces feminization by altering the molecular pathways responsible for eel sex differentiation.

Sex-specific meiosis responses to Gsdf in medaka (Oryzias latipes)

The meiotic entry of undifferentiated germ cells is sexually specific and strictly regulated by the testicular or ovarian environment. Germline stem cells with a set of abnormal sex chromosomes and associated autosomes undergo defective meiotic processes and are eventually eliminated by yet to be defined post-transcriptional modifications. Herein, we report the role of gsdf, a member of BMP/TGFβ family uniquely found in teleost, in the regulation of meiotic entry in medaka (Oryzias latipes) via analyses of gametogenesis in gsdf-deficient XX and XY gonads in comparison with their wild-type siblings. Several differentially expressed genes, including the FKB506-binding protein 7 (fkbp7), were significantly upregulated in pubertal gsdf-deficient gonads. The increase in alternative pre-mRNA isoforms of meiotic synaptonemal complex gene sycp3 was visualized using Integrative Genomics Viewer and confirmed by real-time qPCR. Nevertheless, immunofluorescence analysis showed that Sycp3 protein products reduced significantly in gsdf-deficient XY oocytes. Transmission electron microscope observations showed that normal synchronous cysts were replaced by asynchronous cysts in gsdf-deficient testis. Breeding experiments showed that the sex ratio deviation of gsdf^-/- XY gametes in a non-Mendelian manner might be due to the non-segregation of XY chromosomes. Taken together, our results suggest that gsdf plays a role in the proper execution of cytoplasmic and nuclear events through receptor Smad phosphorylation and Sycp3 dephosphorylation to coordinate medaka gametogenesis, including sex-specific mitotic divisions and meiotic recombination.

Gonadal Soma-Derived Factor Expression is a Potential Biomarker for Predicting the Effects of Endocrine-Disrupting Chemicals on Gonadal Differentiation in Japanese Medaka (Oryzias Latipes)

Chemicals with androgenic or estrogenic activity induce the sex reversal and/or intersex condition in various teleost fish species. Previously, we reported that exposure to 17α-methyltestosterone, bisphenol A, or 4-nonylphenol induces changes in expression of the gonadal soma-derived factor (gsdf) gene accompanied by disruption of gonadal differentiation in Japanese medaka (Oryzias latipes). These findings suggest that gsdf expression might be a useful biomarker for predicting the potential effect of chemicals on gonadal differentiation. We examined the gsdf expression in Japanese medaka exposed to chemicals with estrogenic or androgenic activity. Exposure to the androgenic steroid 17β-trenbolone at 0.5-22.1 μg/L induced the development of ovotestis (presence of ovarian tissue with testicular tissue) and female-to-male sex reversal in XX embryos, and exposure at 6.32 and 22.1 μg/L significantly increased gsdf expression in XX embryos compared with controls at developmental stage 38 (1 day before hatching). In the present study, no statistically significant difference in gsdf mRNA expression was observed after exposure to 17β-estradiol, 17α-ethinylestradiol, and 4-t-octylphenol, which have estrogenic activity. In addition, antiandrogenic chemicals or chemicals without endocrine-disrupting activity did not induce changes in gsdf expression in XX or XY embryos. Thus, an increase in gsdf expression after androgen exposure was observed in XX embryos. Together, these findings indicate that gsdf expression might be useful for predicting the adverse effect of chemicals on gonadal differentiation. Environ Toxicol Chem 2022;41:1875-1884. © 2022 SETAC.

A Potential Role of the Translation Elongation Factor eef1a1 in Gonadal High-Temperature Perception in Chinese Tongue Sole (Cynoglossus semilaevis)

Simple Summary

The eukaryotic elongation factor 1 alpha (eef1a) gene is vital for protein translation by delivering aminoacylated tRNAs to the A/P site of the ribosome via the GTP-dependent reaction. Here, the Chinese tongue sole (Cynoglossus semilaevis) eef1a1 gene was identified, and its potential role in gonadal high-temperature perception was assessed. The full-length sequence of eef1a1 cDNA was 1809 base pair (bp) encoding a putative protein of 461 amino acids. The expression levels of eef1a1 in the ovary were significantly higher than that in the testis from 6 mpf to 3 ypf. Under high-temperature induction during sex differentiation, eef1a1 was significantly down-regulated in males, while the difference was not detected in females. Furthermore, the rapid response of eef1a1 to environmental high temperature was assessed in vitro. Our findings suggest that C. semilaevis eef1a1 might be essential for the molecular response regulatory network of external temperature affecting internal sex differentiation.

Abstract

The eukaryotic translation elongation factor 1 alpha (eef1a) gene has a well-defined role in protein synthesis. However, its role in external temperature perception and internal sex differentiation and development is still unclear. In this study, eef1a1 was identified and functionally analyzed in Chinese tongue sole (Cynoglossus semilaevis). The eef1a1 cDNA, 1809 bp in length, had a 1386 bp open reading frame (ORF) that encoded a 461 amino acid polypeptide containing one EF-1_alpha domain. eef1a1 expression levels were investigated across different tissues and during gonadal development. In the gonad, eef1a1 showed a sexually dimorphic expression pattern with a statistically higher expression level in the ovary than in the testis from 6 months postfertilization to 3 years postfertilization. Under high temperature (28 °C) treatment during C. semilaevis sex differentiation (from 30 days postfertilization to 3 months postfertilization), eef1a1 was statistically down-regulated in males, while the difference was not detected in females. In addition, the dual-luciferase assay exhibited that eef1a1 can respond to high temperature rapidly. Based on these results, C. semilaevis eef1a1 might have a dual role in the perception of external temperature changes and sex differentiation regulation.

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.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

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 miR-200 Family Targeting amh Affects the Gonadal Development of Japanese Flounder

Four members of the miR-200 family in Japanese flounder (Paralichthys olivaceus) have sex-biased expression patterns, but their target genes and how they work in the development of the gonads are rarely known. Anti-Müllerian hormone (AMH) can inhibit the development of Muller’s duct in female mammals and regulate the formation of gametes after sexual maturity. There is no Muller’s duct in teleosts, but the amh gene still exists. Knockout of amh results in sex reversal from male to female. Therefore, it is essential to explore the relationship between the miR-200 family and amh to clarify what role miR-200 plays in the development of the gonads. In Japanese flounder, the two binding sites for the miR-200 family in the 3′UTR of amh were found through bioinformatic prediction. Double luciferase and green fluorescent protein reporter experiments demonstrated amh to be directly targeted by miR-200a and miR-200b. Moreover, miR-200a and miR-200b reduced the expression of amh through site 1 rather than site 2. To explore the regulatory role of miR-200a in gonadal development, we further overexpressed miR-200a in the primary Sertoli cells of the testis. With the overexpression of miR-200a, the expression of amh decreased, while the expression of the other two male sex-related genes, dmrt1 (doublesex and mab-3 related transcription factor 1) and gsdf (diagonal soma driven factor), increased significantly. This result indicates that the miR-200 family regulates the gonadal differentiation and development by targeting amh in Japanese flounder.

Developmental changes in gene expression and gonad morphology during sex differentiation in Atlantic salmon (Salmo salar)

The Atlantic salmon (Salmo salar) is a globally important species for its value in fisheries and aquaculture, and as a research model. In order to characterise aspects of sex differentiation at the morphological and mRNA level in this species, the present study examined developmental changes in gonad morphology and gene expression in males and females between 0 and 79 days post hatch (dph). Morphological differentiation of the ovary (indicated by the formation of germ cell cysts) became apparent from 52 dph. By 79 dph, ovarian phenotype was evident in 100% of genotypic females. Testes remained in an undifferentiated-like state throughout the experiment, containing germ cells dispersed singularly within the gonadal region distal to the mesentery. There were no significant sex-related differences in gonad cross-section size, germ cell number or germ cell diameter during the experiment. The expression of genes involved in teleost sex differentiation (anti-müllerian hormone (amh), cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a), forkhead box L2a (foxl2a), gonadal soma-derived factor (gsdf), r-spondin 1 (rspo1), sexually dimorphic on the Y chromosome (sdY)), retinoic acid-signalling (aldehyde dehydrogenase 1a2 (aldh1a2), cytochrome P450 family 26 a1 (cyp26a1), cytochrome P450 family 26 b1 (cyp26b1), t-box transcription factor 1 (tbx1a)) and neuroestrogen production (cytochrome P450, family 19, subfamily A, polypeptide 1b (cyp19a1b)) was investigated. Significant sex-related differences were observed only for the expression of amh, cyp19a1a, gsdf and sdY. In males, amh, gsdf and sdY were upregulated from 34, 59 and 44 dph respectively. In females, cyp19a1a was upregulated from 66 dph. Independent of sex, foxl2a expression was highest at 0 dph and had reduced ∼ 47-fold by the time of morphological sex differentiation at 52 dph. This study provides new insights into the timing and sequence of some physiological changes associated with sex differentiation in Atlantic salmon. These findings also reveal that some aspects of the mRNA sex differentiation pathways in Atlantic salmon are unique compared to other teleost fishes, including other salmonids.

Transcriptomes of testis and pituitary from male Nile tilapia (O. niloticus L.) in the context of social status

African cichlids are well established models for studying social hierarchies in teleosts and elucidating the effects social dominance has on gene expression. Ascension in the social hierarchy has been found to increase plasma levels of steroid hormones, follicle stimulating hormone (Fsh) and luteinizing hormone (Lh) as well as gonadosomatic index (GSI). Furthermore, the expression of genes related to gonadotropins and steroidogenesis and signaling along the brain-pituitary-gonad axis (BPG-axis) is affected by changes of an animal’s social status. In this study, we use RNA-sequencing to obtain an in-depth look at the transcriptomes of testes and pituitaries from dominant and subordinate male Nile tilapia living in long-term stable social hierarchies. This allows us to draw conclusions about factors along the brain-pituitary-gonad axis that are involved in maintaining dominance over weeks or even months. We identify a number of genes that are differentially regulated between dominant and subordinate males and show that in high-ranking fish this subset of genes is generally upregulated. Genes differentially expressed between the two social groups comprise growth factors, related binding proteins and receptors, components of Wnt-, Tgfβ- and retinoic acid-signaling pathway, gonadotropin signaling and steroidogenesis pathways. The latter is backed up by elevated levels of 11-ketotestosterone, testosterone and estradiol in dominant males. Luteinizing hormone (Lh) is found in higher concentration in the plasma of long-term dominant males than in subordinate animals. Our results both strengthen the existing models and propose new candidates for functional studies to expand our understanding of social phenomena 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.

Full-Length Transcriptome Sequences Provide Insight Into Hermaphroditism of Freshwater Pearl Mussel Hyriopsis schlegelii

The freshwater mussel Hyriopsis schlegelii is a cultured bivalve in China, and the quality of the pearls produced is affected by the type of gonads. However, because of the lack of a published genome and the complexity of sex determination, research on sex reversal and development of this species is limited. In this study, Illumina RNA-seq and PacBio Isoform Sequencing (Iso-Seq) were combined to analyze the gonads of H. schlegelii. A total of 201,481 high-quality transcripts were generated. The study identified 7,922 differentially expressed genes in three comparison group (females versus males, hermaphrodites versus females, and hermaphrodites versus males). Twenty-four genes were identified as potential sex-related genes, including sox9 and wnt4 involved in sex determination, and vtg, cyp17a1 and 17β-hsd2 involved in gonadal development. We also speculated a possible pathways for the formation of hermaphroditism in H. schlegelii. The data provide a clear view of the transcriptome for H. schlegelii gonads and will be valuable in elucidating the mechanisms of gonad development.

Homozygous Mutation of gsdf Causes Infertility in Female Nile Tilapia (Oreochromis niloticus)

Gonadal somatic cell-derived factor (Gsdf) is a member of the TGF-β superfamily, which exists mainly in fishes. Homozygous gsdf mutations in Japanese medaka and zebrafish resulted in infertile females, and the reasons for their infertility remain unknown. This study presents functional studies of Gsdf in ovary development using CRISPR/Cas9 in Nile tilapia (Oreochromis niloticus). The XX wild type (WT) female fish regularly reproduced from 12 months after hatching (mah), while the XX gsdf-/- female fish never reproduced and were infertile. Histological observation showed that at 24 mah, number of phase IV oocyte in the XX gsdf-/- female fish was significantly lower than that of the WT fish, although their gonadosomatic index (GSI) was similar. However, the GSI of the XX gsdf-/- female at 6 mah was higher than that of the WT. The mutated ovaries were hyperplastic with more phase I oocytes. Transcriptome analysis identified 344 and 51 up- and down-regulated genes in mutants compared with the WT ovaries at 6 mah. Some TGF-β signaling genes that are critical for ovary development in fish were differentially expressed. Genes such as amh and amhr2 were up-regulated, while inhbb and acvr2a were down-regulated in mutant ovaries. The cyp19a1a, the key gene for estrogen synthesis, was not differentially expressed. Moreover, the serum 17β-estradiol (E2) concentrations between XX gsdf-/- and WT were similar at 6 and 24 mah. Results from real-time PCR and immunofluorescence experiments were similar and validated the transcriptome data. Furthermore, Yeast-two-hybrid assays showed that Gsdf interacts with TGF-β type II receptors (Amhr2 and Bmpr2a). Altogether, these results suggest that Gsdf functions together with TGF-β signaling pathway to control ovary development and fertility. This study contributes to knowledge on the function of Gsdf in fish oogenesis.

Oxidative Stress Causes Masculinization of Genetically Female Medaka Without Elevating Cortisol

Medaka (Oryzias latipes) is a teleost fish with an XX/XY sex determination system. Sex reversal from female-to-male (masculinization of XX fish) can be induced through cortisol elevation from exposure to environmental stress such as high temperature during sexual differentiation. However, the effects of oxidative stress, generated via metabolic reactions and biological defense mechanisms, on the sexual differentiation of medaka are unclear. Here, we investigated the effect of oxidative stress on medaka sexual differentiation using hydrogen peroxide (H₂O₂), which induces oxidative stress in vertebrates. H₂O₂ treatment from 0 to 5 days post-hatching induced masculinization of wild-type XX medaka, but not of gonadal soma-derived growth factor (gsdf) or peroxisome proliferator-activated receptor alpha-a (pparaa) knockout XX fish. Co-treatment with an oxidative stress inhibitor caused masculinization recovery but co-treatment with a cortisol synthesis inhibitor did not. H₂O₂ treatment significantly upregulated gsdf and pparaa expression in XX medaka. However, H₂O₂ did not elevate cortisol levels in medaka larvae during sexual differentiation. These results strongly indicate that oxidative stress induces masculinization of XX medaka without causing elevation of cortisol.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

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.

The application of genome editing technology in fish

The advent and development of genome editing technology has opened up the possibility of directly targeting and modifying genomic sequences in the field of life sciences with rapid developments occurring in the last decade. As a powerful tool to decipher genome data at the molecular biology level, genome editing technology has made important contributions to elucidating many biological problems. Currently, the three most widely used genome editing technologies include: zinc finger nucleases (ZFN), transcription activator like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats (CRISPR). Researchers are still striving to create simpler, more efficient, and accurate techniques, such as engineered base editors and new CRISPR/Cas systems, to improve editing efficiency and reduce off-target rate, as well as a near-PAMless SpCas9 variants to expand the scope of genome editing. As one of the important animal protein sources, fish has significant economic value in aquaculture. In addition, fish is indispensable for research as it serves as the evolutionary link between invertebrates and higher vertebrates. Consequently, genome editing technologies were applied extensively in various fish species for basic functional studies as well as applied research in aquaculture. In this review, we focus on the application of genome editing technologies in fish species detailing growth, gender, and pigmentation traits. In addition, we have focused on the construction of a zebrafish (Danio rerio) disease model and high-throughput screening of functional genes. Finally, we provide some of the future perspectives of this technology.

Sex-Related Embryotoxicity of Pulp Mill Effluent Extracts in Medaka (Oryzias latipes) Female Leucophore-free FLFII Strain

The aim of the present study was to evaluate the effects of exposure to Chilean pulp mill effluent extracts on developing postfertilized medaka embryos before and after sex definition relative to sex steroids (testosterone and 17beta-estradiol) and a wood phytoestrogen (beta-sitosterol). Our study included 2 waterborne semichronic exposure experiments, using a 24-h post fertilization (hpf) unknown-sex FLFII (female leucophore free) group and a second 72-hpf FLFII phenotypic sex-identified group (male autofluorescence leucophore) strain of medaka embryos. Chronic exposure of both FLFII strain embryo groups showed similar delay in time to hatch and decreased hatchability. Teratogenic responses such as vertebral malformation (fusion, incomplete formation, and lack of vertebral formation process) and pericardial edema were observed in both experiments, with a high percentage related to FLFII fluorescent leucophore-identified males. In addition, high mortality associated with severe malformations was observed in male and female embryos exposed to testosterone. Our research has demonstrated that exposure to Chilean mill effluent extracts caused severe male medaka embryotoxicity (in postfertilized embryos) before and after sex definition and, irrespective of the experimental group and effluent treatment, suggests partial removal following secondary treatment. Furthermore, differences in the severity and type of teratogenic effects with previous experiments (d-rR medaka strain), are associated with the unique phenotypes of this medaka mutant strain. Environ Toxicol Chem 2021;40:2297-2305. © 2021 SETAC.

The replaceable master of sex determination: bottom-up hypothesis revisited

Different group of vertebrates and invertebrates demonstrate an amazing diversity of gene regulations not only at the top but also at the bottom of the sex determination genetic network. As early as 1995, based on emerging findings in Drosophila melanogaster and Caenorhabditis elegans, Wilkins suggested that the evolution of the sex determination pathway evolved from the bottom to the top of the hierarchy. Based on our current knowledge, this review revisits the 'bottom-up' hypothesis and applies its logic to vertebrates. The basic operation of the determination network is through the dynamics of the opposing male and female pathways together with a persistent need to maintain the sexual identity of the cells of the gonad up to the reproductive stage in adults. The sex-determining trigger circumstantially acts from outside the genetic network, but the regulatory network is not built around it as a main node, thus maintaining the genetic structure of the network. New sex-promoting genes arise either through allelic diversification or gene duplication and act specially at the sex-determination period, without integration into the complete network. Due to this peripheral position the new regulator is not an indispensable component of the sex-determining network and can be easily replaced. 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)'.

Evolution, Expression, and Function of Gonadal Somatic Cell-Derived Factor

Fish gonads develop in very diverse ways different from mammalian gonads. This diversity is contributed by species-specific factors. Gonadal somatic cell-derived factor (Gsdf) is one such factor. The gsdf gene exists mostly in teleosts and is absent in many tetrapods, probably as a result of two gene losses during evolution. The gsdf transcript is expressed mainly in gonadal somatic cells, including Sertoli cell in testis and granulosa cells in ovary; however, these gonadal somatic cells can surround many types of germ cells at different developmental stages depending on the fish species. The function of gsdf is also variable. It is involved in germ cell proliferation, testicular formation, ovarian development and even male sex determination. Here, we summarize the common and diverse expression, regulation and functions of gsdf among different fish species with aspect of evolution.

Effects of synthetic sex steroid hormone exposures on gonadal sex differentiation and dynamics of a male-related gene, Gonadal soma-derived factor (Gsdf) and an estrogen up-regulated gene, Choriogenine-H (ChgH) gene expression in the euryhaline Javafish medaka, Oryzias javanicus, based on genetic sexes

To clarify the basal aspects of sexual development in Javafish medaka, Oryzias javanicus (ZZ/ZW), a model marine species for ecotoxicity testing, we examined the details of gonadal sex differentiation and exogenous sex hormone-dependent sex reversals using genetic sex-linked DNA markers. Sex differences in germ cell numbers were observed at 5 days post hatching (dph), in which there was a significant increase in the germ cells of ZW. In ZW, diplotene oocytes and the ovarian cavity appeared at approximately 10, and 30 dph, respectively. In ZZ, spermatogonial proliferation was observed at approximately 20 dph. A ZZ-dominant expression of Gonadal soma-derived factor (Gsdf) mRNA was detected before hatching. The exposure of embryos to 17α-ethinylestradiol (EE2; 0.1, 1, 10 ng/mL) did not cause sex reversals in most cases. However, EE2 exposures led to significant Choriogenin-H (ChgH) mRNA expression, an estrogen up-regulated gene, in all fry; these exposures did not suppress Gsdf expression in ZZ fry. The exposure of embryos to 17α-methyltestosterone (MT; 0.1, 1, 10 ng/mL) caused sex reversals but only at low frequencies in ZW and ZZ fish. Although the 10 ng/mL MT exposure was accompanied by induction of significant Gsdf expression in ZW fry, induction of ChgH expression was also observed in several fry. Together, the present study indicates for the first time that male-dominant sexual dimorphic expression of Gsdf precedes the first morphological sex difference, i.e., the sex difference in germ cell number, and results strongly suggest that exogenous sex hormone-dependent sex reversal is not induced easily in O. javanicus.

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.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy

Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus

The genetic regulatory network governing vertebrate gonadal differentiation appears less conserved than previously thought. Here, we investigated the gonadal development of Astyanax mexicanus Pachón cavefish by looking at primordial germ cells (PGCs) migration and proliferation, gonad histology, and gene expression patterns. We showed that PGCs are first detected at the 80% epiboly stage and then reach the gonadal primordium at 1 day post-fertilization (dpf). However, in contrast to the generally described absence of PGCs proliferation during their migration phase, PGCs number in cavefish doubles between early neurula and 8-9 somites stages. Combining both gonadal histology and vasa (germ cell marker) expression patterns, we observed that ovarian and testicular differentiation occurs around 65 dpf in females and 90 dpf in males, respectively, with an important inter-individual variability. The expression patterns of dmrt1, gsdf, and amh revealed a conserved predominant male expression during cavefish gonadal development, but none of the ovarian differentiation genes, i. e., foxl2a, cyp19a1a, and wnt4b displayed an early sexually dimorphic expression, and surprisingly all these genes exhibited predominant expression in adult testes. Altogether, our results lay the foundation for further research on sex determination and differentiation in A. mexicanus and contribute to the emerging picture that the vertebrate sex differentiation downstream regulatory network is less conserved than previously thought, at least in teleost fishes.

A Potential Role for the Gsdf-eEF1α Complex in Inhibiting Germ Cell Proliferation: A Protein-Interaction Analysis in Medaka (Oryzias latipes) From a Proteomics Perspective

Gonadal soma-derived factor (gsdf) has been demonstrated to be essential for testicular differentiation in medaka (Oryzias latipes). To understand the protein dynamics of Gsdf in spermatogenesis regulation, we used a His-tag "pull-down" assay coupled with shotgun LC-MS/MS to identify a group of potential interacting partners for Gsdf, which included cytoplasmic dynein light chain 2, eukaryotic polypeptide elongation factor 1 alpha (eEF1α), and actin filaments in the mature medaka testis. As for the interaction with transforming growth factor β-dynein being critical for spermatogonial division in Drosophila melanogaster, the physical interactions of Gsdf-dynein and Gsdf-eEF1α were identified through a yeast 2-hybrid screening of an adult testis cDNA library using Gsdf as bait, which were verified by a paired yeast 2-hybrid assay. Coimmunoprecipitation of Gsdf and eEF1α was defined in adult testes as supporting the requirement of a Gsdf and eEF1α interaction in testis development. Proteomics analysis (data are available via ProteomeXchange with identifier PXD022153) and ultrastructural observations showed that Gsdf deficiency activated eEF1α-mediated protein synthesis and ribosomal biogenesis, which in turn led to the differentiation of undifferentiated germ cells. Thus, our results provide a framework and new insight into the coordination of a Gsdf (transforming growth factor β) and eEF1α complex in the basic processes of germ cell proliferation, transcriptional and translational control of sexual RNA, which may be fundamentally conserved across the phyla during sexual differentiation.

Effects of the Thyroid Endocrine System on Gonadal Sex Ratios and Sex-Related Gene Expression in the Pufferfish Takifugu rubripes

To examine the effect and mechanism of thyroid hormone on gonadal sex differentiation, Takifugu rubripes larvae were treated with goitrogen (methimazole, MET, 1000 g/g), and thyroxine (T4, 2nM) from 25 to 80 days after hatching (dah). Gonadal histology and sex ratios of fish were then determined at 80 dah. MET treatment induced masculinization, but T4 treatment did not induce feminization in T. rubripes larvae. Transcriptomic analysis of gonads at 80 dah was then conducted. Among the large number of differentially expressed genes between the groups, the expression of foxl2, cyp19a1a, and dmrt1 was altered. The expression of foxl2, cyp19a1a, dmrt1 and gsdf at 25, 40, 55 days after treatment (dat) was further analyzed by qPCR. MET treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 in genetic females (p < 0.05). Additionally, T4 treatment induced an increase in the expression of cyp19a1a in genetic XY gonads only at 25 dat. However, the increase in cyp19a1a expression did not continue to 40 and 55 dat. This may explain why feminization of larvae was not found in the T4-treated group. Thus, the present study provides the first evidence that MET treatment causes masculinization in teleost fish. The effects of MET-induced masculinization in T. rubripes may act primarily via suppression of the expression of foxl2 and cyp19a1a, and stimulation of the expression of dmrt1. Moreover, the effects of higher concentrations of T4 or different concentrations of T3, on sex differentiation require further testing.

High temperature stress response is not sexually dimorphic at the whole-body level and is dependent on androgens to induce sex reversal

The understanding of the molecular and endocrine mechanisms behind environmentally-induced sex reversal in fish is of great importance in the context of predicting the potential effects of climate change, especially increasing temperature. Here, we demonstrate the global effects of high temperature on genome-wide transcription in medaka (Oryzias latipes) during early development. Interestingly, data analysis did not show sexual dimorphic changes, demonstrating that thermal stress is not dependent on genotypic sex. Additionally, our results revealed significant changes in several pathways under high temperature, such as stress response from brain, steroid biosynthesis, epigenetic mechanisms, and thyroid hormone biosynthesis, among others. These microarray data raised the question of what the exact molecular and hormonal mechanisms of action are for female-to-male sex reversal under high temperatures in fish. Complementary gene expression analysis revealed that androgen-related genes increase in females (XX) experiencing high water temperature. To test the involvement of androgens in thermal-induced sex reversal, an androgen antagonist was used to treat XX medaka under a high-temperature setup. Data clearly demonstrated failure of female-to-male sex reversal when androgen action is inhibited, corroborating the importance of androgens in environmentally-induced sex reversal.

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.

Molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with XX/XY sex determination

Phenotypic sex of an organism is determined by molecular changes in the gonads, so-called molecular sex differentiation, which should precede the rise of cellular or anatomical sex-distinguishing features. This study characterized molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with established XX/XY genotypic sex determination. Next generation sequencing was conducted on sablefish ovarian and testicular mRNAs to obtain sequences for transcripts associated with vertebrate sex determination and differentiation and early reproductive development. Gene-specific PCRs were developed to determine the distribution and ontogenetic gonadal expression of transcription, growth, steroidogenic and germline factors, as well as gonadotropin and steroid receptors. Molecular changes associated with sex differentiation were first apparent in both XY- and XX-genotype sablefish at ~ 60 mm in body length and prior to histological signs of sex differentiation. The earliest and most robust markers of testicular differentiation were gsdf, amh, dmrt1, cyp11b, star, sox9a, and fshr. Markedly elevated mRNA levels of several steroidogenesis-related genes and ar2 in differentiating testes suggested that androgens play a role in sablefish testicular differentiation. The earliest markers of ovarian differentiation were cyp19a1a, lhcgr, foxl2, nr0b1, and igf3. Other transcripts such as figla, zp3, and pou5f3 were expressed predominantly in XX-genotype fish and significantly increased with the first appearance and subsequent development of primary oocytes. This study provides valuable insight to the developmental sequence of events associated with gonadal sex differentiation in marine teleosts with XX/XY sex determination. It also implicates particular genes in processes of male and female development and establishes robust molecular markers for phenotypic sex in sablefish, useful for ongoing work related to sex control and reproductive sterilization.

Transcriptome Profiling and Analysis of Genes Associated with High Temperature-Induced Masculinization in Sex-Undifferentiated Nile Tilapia Gonad

Artificially high temperatures during critical thermosensitive periods (TSPs) can induce the sex reversal of Nile tilapia (Oreochromis niloticus) females into pseudomales; Nile tilapia is a GSD + TE (genotypic plus temperature effects) fish species. Previous studies have shown that water temperature affects the expression levels of many genes in the gonad or brain in various teleost species. However, few studies on the effect of temperature at the whole-gonad transcriptomic level in the early stage of sex differentiation have been reported in fish species exhibiting GSD + TE. In this study, RNA-Seq was performed to characterize the transcriptomic profile and identify genes exhibiting temperature- and sex-biased expressions in the Nile tilapia gonad at 21 dpf. A total of 42 genes were found to be associated with both high-temperature treatment and sex development, as the expression levels of these genes differed in both FC (female control) vs MC (male control) and FC vs FT (high temperature-treated females in the TSP). Among these genes, the transcriptional alterations of many male sex determination and differentiation genes, such as Dmrt1, Gsdf, and the DNA damage-inducible protein GADD45 alpha, suggested that the male pathway is initiated after high-temperature treatment and that its initiation may play a role in high temperature-induced masculinization in Nile tilapia. The qRT-PCR validation results for thirteen differentially expressed genes showed that the Pearson's correlation of the log₁₀ fold change values between the qPCR and RNA-Seq results was 0.70 (p < 0.001), indicating the accuracy and reliability of the RNA-Seq results. Our study provides insights into how high-temperature treatment induces the sex reversal of Nile tilapia females.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Integrative analysis reveals pathways associated with sex reversal in Cynoglossus semilaevis

Sex reversal is a complex biological phenomenon exhibited by Cynoglossus semilaevis. Some genetic females may irreversibly convert to pseudomales, thus increasing aquaculture costs because males grow much more slowly than females. In this study, an integrative analysis of transcriptome and proteome was performed to compare differences in gene and protein expression in females and pseudomales after gonad differentiation in C. semilaevis. Based on RNA-Seq results, 1893 genes showed differences in expression at the transcript level between females and pseudomales. Of these differentially expressed genes (DEGs), zona pellucida sperm-binding protein 4-like (LOC103393374 , ZP4), zona pellucida sperm-binding protein 4-like (LOC103396071, ZP4) and forkhead box L2 (foxl2) were highly expressed in females and doublesex and mab-3 related transcription factor 1(dmrt1) and doublesex and mab-3 related transcription factor 3 (dmrt3) were highly expressed in pseudomales. GO enrichment analysis results indicate that wnt signaling pathways and oocyte maturation are two terms enriched in female. At the protein level, Tandem Mass Tags analysis revealed that 324 proteins differed in their relative abundance between pseudomales and females. KEGG analysis found that pseudo-highly expressed proteins were enriched in the ubiquitin mediated proteolysis pathway. For integrative analysis, the Spearman correlation coefficient between the transcriptome and proteome was 0.59. Among 52 related genes, 46 DEGs (88%) were well matched in their levels of change in protein abundance. These findings reveal major active pathways in female and pseudomale gonads after sex reversal and provide new insights into molecular mechanisms associated with sex reversal regulatory network.

Transcriptomic analysis of dead end knockout testis reveals germ cell and gonadal somatic factors in Atlantic salmon

Background

Sustainability challenges are currently hampering an increase in salmon production. Using sterile salmon can solve problems with precocious puberty and genetic introgression from farmed escapees to wild populations. Recently sterile salmon was produced by knocking out the germ cell-specific dead end (dnd). Several approaches may be applied to inhibit Dnd function, including gene knockout, knockdown or immunization. Since it is challenging to develop a successful treatment against a gene product already existing in the body, alternative targets are being explored. Germ cells are surrounded by, and dependent on, gonadal somatic cells. Targeting genes essential for the survival of gonadal somatic cells may be good alternative targets for sterility treatments. Our aim was to identify and characterize novel germ cell and gonadal somatic factors in Atlantic salmon.

Results

We have for the first time analysed RNA-sequencing data from germ cell-free (GCF)/dnd knockout and wild type (WT) salmon testis and searched for genes preferentially expressed in either germ cells or gonadal somatic cells. To exclude genes with extra-gonadal expression, our dataset was merged with available multi-tissue transcriptome data. We identified 389 gonad specific genes, of which 194 were preferentially expressed within germ cells, and 11 were confined to gonadal somatic cells. Interestingly, 5 of the 11 gonadal somatic transcripts represented genes encoding secreted TGF-β factors; gsdf, inha, nodal and two bmp6-like genes, all representative vaccine targets. Of these, gsdf and inha had the highest transcript levels. Expression of gsdf and inha was further confirmed to be gonad specific, and their spatial expression was restricted to granulosa and Sertoli cells of the ovary and testis, respectively. Finally, we show that inha expression increases with puberty in both ovary and testis tissue, while gsdf expression does not change or decreases during puberty in ovary and testis tissue, respectively.

Conclusions

This study contributes with transcriptome data on salmon testis tissue with and without germ cells. We provide a list of novel and known germ cell- and gonad somatic specific transcripts, and show that the expression of two highly active gonadal somatic secreted TGF-β factors, gsdf and inha, are located within granulosa and Sertoli cells.

Genotoxic biomarkers and histological changes in marine medaka (Oryzias melastigma) exposed to 17α-ethynylestradiol and 17β-trenbolone

Endocrine-disrupting pollutants in marine environments have aroused great concern for their adverse effects on the reproduction of marine organisms. This study aimed to seek promising biomarkers for estrogenic/androgenic chemicals. First, two possible male-specific genes, SRY-box containing gene 9a2 (sox9a2) and gonadal soma-derived factor (gsdf), were cloned from marine medaka (Oryzias melastigma). Then the responses of sox9a2, gsdf, choriogenin (chgH and chgL), vitellogenin (vtg1 and vtg2), and cytochrome P450 aromatase (cyp19a and cyp19b) were investigated after exposure to 17α-ethynylestradiol (EE₂) and 17β-trenbolone (TB) at 2, 10, and 50 ng/L. The results showed that gsdf was specifically expressed in the testes and easily induced in the ovaries after TB exposure, indicating that gsdf was a potential biomarker of environmental androgens. ChgL was a useful biomarker of weak estrogen pollution for its high sensitivity to low levels of EE₂. In addition, both EE₂ and TB exposure damaged gonadal structures and inhibited gonadal development.

Comparative transcriptome analysis of male and female gonads reveals sex-biased genes in spotted scat (Scatophagus argus)

Scatophagus argus is a new emerging aquaculture fish in East and Southeast Asia. To date, research on reproductive development and regulation in S. argus is lacking. Additionally, genetic and genomic information about reproduction, such as gonadal transcriptome data, is also lacking. Herein, we report the first gonadal transcriptomes of S. argus and identify genes potentially involved in reproduction and gonadal development. A total of 136,561 unigenes were obtained by sequencing of testes (n = 3) and ovaries (n = 3) at stage III. Genes upregulated in males and females known to be involved in gonadal development and gametogenesis were identified, including male-biased dmrt1, amh, gsdf, wt1a, sox9b, and nanos2, and female-biased foxl2, gdf9, bmp15, sox3, zar1, and figla. Serum estradiol-17β and 11-ketotestosterone levels were biased in female and male fish, respectively. Sexual dimorphism of serum steroid hormone levels were interpreted after expression analysis of 20 steroidogenesis-related genes, including cyp19a1a and cyp11b2. This gonadal transcript dataset will help investigate functional genes related to reproduction in S. argus.