Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

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.

Adamts9 is required for the development of primary ovarian follicles and maintenance of female sex in zebrafish†

Previous studies have suggested that adamts9 (a disintegrin and metalloprotease with thrombospondin type-1 motifs, member 9), an extracellular matrix (ECM) metalloprotease, participates in primordial germ cell (PGC) migration and is necessary for female fertility. In this study, we found that adamts9 knockout (KO) led to reduced body size, and female-to-male sex conversion in late juvenile or adult zebrafish; however, primary sex determination was not affected in early juveniles of adamts9 KO. Overfeeding and lowering the rearing density rescued growth defects in female adamts9 KO fish but did not rescue defects in ovarian development in adamts9 KO. Delayed PGC proliferation, significantly reduced number and size of Stage IB follicles (equivalent to primary follicles) in early juveniles of adamts9 KO, and arrested development at Stage IB follicles in mid- or late-juveniles of adamts9 KO are likely causes of female infertility and sex conversion. Via RNAseq, we found significant enrichment of differentially expressed genes involved in ECM organization during sexual maturation in ovaries of wildtype fish; and significant dysregulation of these genes in adamts9 KO ovaries. RNAseq analysis also showed enrichment of inflammatory transcriptomic signatures in adult ovaries of these adamts9 KO. Taken together, our results indicate that adamts9 is critical for development of primary ovarian follicles and maintenance of female sex, and loss of adamts9 leads to defects in ovarian follicle development, female infertility, and sex conversion in late juveniles and mature adults. These results show that the ECM and extracellular metalloproteases play major roles in maintaining ovarian follicle development in zebrafish.

Dissecting the dynamic cellular transcriptional atlas of adult teleost testis development throughout the annual reproductive cycle

Teleost testis development during the annual cycle involves dramatic changes in cellular compositions and molecular events. In this study, the testicular cells derived from adult black rockfish at distinct stages - regressed, regenerating and differentiating - were meticulously dissected via single-cell transcriptome sequencing. A continuous developmental trajectory of spermatogenic cells, from spermatogonia to spermatids, was delineated, elucidating the molecular events involved in spermatogenesis. Subsequently, the dynamic regulation of gene expression associated with spermatogonia proliferation and differentiation was observed across spermatogonia subgroups and developmental stages. A bioenergetic transition from glycolysis to mitochondrial respiration of spermatogonia during the annual developmental cycle was demonstrated, and a deeper level of heterogeneity and molecular characteristics was revealed by re-clustering analysis. Additionally, the developmental trajectory of Sertoli cells was delineated, alongside the divergence of Leydig cells and macrophages. Moreover, the interaction network between testicular micro-environment somatic cells and spermatogenic cells was established. Overall, our study provides detailed information on both germ and somatic cells within teleost testes during the annual reproductive cycle, which lays the foundation for spermatogenesis regulation and germplasm preservation of endangered species.

Hexafluoropropylene oxide trimer acid (HFPO-TA) disrupts sex differentiation of zebrafish (Danio rerio) via an epigenetic mechanism of DNA methylation

Hexafluoropropylene oxide trimer acid (HFPO-TA), an alternative to perfluorooctanoic acid, has been shown to have estrogenic effects. However, its potential to disrupt fish sex differentiation during gonadal development remains unknown. Therefore, this study exposed zebrafish to HFPO-TA from approximately 2 hours post fertilization (hpf) to 60 days post fertilization (dpf) to investigate its effects on sex differentiation. Results indicated that HFPO-TA disrupted steroid hormone homeostasis, delayed gonadal development in both sexes, and resulted in a female-skewed sex ratio in zebrafish. HFPO-TA exposure up-regulated gene expressions of cyp19a1a, esr1, vtg1 and foxl2, while down-regulated those of amh, sox9a and dmrt1. These suggested that HFPO-TA dysregulated the expressions of key genes related to sex differentiation of zebrafish, promoted the production and activation of estrogen, and further induced the feminization. Interestingly, we observed promoter hypomethylation of cyp19a1a and promoter hypermethylation of amh in male zebrafish, which were negatively associated with their gene expressions. These suggested that HFPO-TA dysregulated these key genes through DNA methylation in their promoters. Therefore, the HFPO-TA disrupted the sex differentiation of zebrafish through an epigenetic mechanism involving DNA methylation, ultimately skewing the sex ratio towards females. Overall, this study demonstrated adverse effects of HFPO-TA on fish sex differentiation and provided novel insights into the underlying epigenetic mechanism.

foxl2l is a germ cell-intrinsic gatekeeper of oogenesis in zebrafish

Zebrafish serve as a valuable model organism for studying germ cell biology and reproductive processes. The AB strain of zebrafish is proposed to exhibit a polygenic sex determination system, where most males initially develop juvenile ovaries before committing to male fate. In species with chromosomal sex determination, gonadal somatic cells are recognized as key determinants of germ cell fate. Notably, the loss of germ cells in zebrafish leads to masculinization, implying that germ cells harbor an intrinsic feminization signal. However, the specific signal triggering oogenesis in zebrafish remains unclear. In the present study, we identified foxl2l as an oocyte progenitor-specific gene essential for initiating oogenesis in germ cells. Results showed that foxl2l-knockout zebrafish bypassed the juvenile ovary stage and exclusively developed into fertile males. Further analysis revealed that loss of foxl2l hindered the initiation of oocyte-specific meiosis and prevented entry into oogenesis, leading to premature spermatogenesis during early gonadal development. Furthermore, while mutation of the pro-male gene dmrt1 led to fertile female differentiation, simultaneous disruption of foxl2l in dmrt1 mutants completely blocked oogenesis, with a large proportion of germ cells arrested as germline stem cells, highlighting the crucial role of foxl2l in oogenesis. Overall, this study highlights the unique function of foxl2l as a germ cell-intrinsic gatekeeper of oogenesis in zebrafish.

Loss of Function of Vasoactive-intestinal Peptide Alters Sex Ratio and Reduces Male Reproductive Fitness in Zebrafish

Vasoactive-intestinal peptide (Vip) is a pleiotropic peptide with a wide range of distribution and functions. Zebrafish possess 2 isoforms of Vip (a and b), in which Vipa is most homologous to the mammalian form. In female zebrafish, Vipa can stimulate LH secretion from the pituitary but is not essential for female reproduction, as vipa-/- females display normal reproduction. In contrast, we have found that vipa-/- males are severely subfertile and sex ratio of offspring is female-biased. By analyzing all aspects of male reproduction with wild-type (WT) males, we show that the testes of vipa-/- are underdeveloped and contain ∼70% less spermatids compared to WT counterparts. The sperm of vipa-/- males displayed reduced potency in terms of fertilization (by ∼80%) and motility span and duration (by ∼50%). In addition, vipa-/- male attraction to WT females was largely nonexistent, indicating decreased sexual motivation. We show that vipa mRNA and protein is present in Leydig cells and in developing germ cells in the testis of WT, raising the possibility that endogenous Vipa contributes to testicular function. Absence of Vipa in vipa-/- males resulted in downregulation of 3 key genes in the androgen synthesis chain in the testis, 3β-hsd, 17β-hsd1, and cyp11c1 (11β-hydrogenase), associated with a pronounced decrease in 11-ketotestosterone production and, in turn, compromised reproductive fitness. Altogether, this study establishes a crucial role for Vipa in the regulation of male reproduction in zebrafish, like in mammals, with the exception that Vipa is also expressed in zebrafish testis.

In Search of a Target Gene for a Desirable Phenotype in Aquaculture: Genome Editing of Cyprinidae and Salmonidae Species

Aquaculture supplies the world food market with a significant amount of valuable protein. Highly productive aquaculture fishes can be derived by utilizing genome-editing methods, and the main problem is to choose a target gene to obtain the desirable phenotype. This paper presents a review of the studies of genome editing for genes controlling body development, growth, pigmentation and sex determination in five key aquaculture Salmonidae and Cyprinidae species, such as rainbow trout (Onchorhynchus mykiss), Atlantic salmon (Salmo salar), common carp (Cyprinus carpio), goldfish (Carassius auratus), Gibel carp (Carassius gibelio) and the model fish zebrafish (Danio rerio). Among the genes studied, the most applicable for aquaculture are mstnba, pomc, and acvr2, the knockout of which leads to enhanced muscle growth; runx2b, mutants of which do not form bones in myoseptae; lepr, whose lack of function makes fish fast-growing; fads2, Δ6abc/5Mt, and Δ6bcMt, affecting the composition of fatty acids in fish meat; dnd mettl3, and wnt4a, mutants of which are sterile; and disease-susceptibility genes prmt7, gab3, gcJAM-A, and cxcr3.2. Schemes for obtaining common carp populations consisting of only large females are promising for use in aquaculture. The immobilized and uncolored zebrafish line is of interest for laboratory use.

Hypogonadotropic hypogonadism in male tilapia lacking a functional rln3b gene

Relaxin 3 is a neuropeptide that plays a crucial role in reproductive functions of mammals. Previous studies have confirmed that rln3a plays an important role in the male reproduction of tilapia. To further understand the significance of its paralogous gene rln3b in male fertility, we generated a homozygous mutant line of rln3b in Nile tilapia. Our findings indicated that rln3b mutation delayed spermatogenesis and led to abnormal testes structure. Knocking out rln3b gene resulted in a decrease in sperm count, sperm motility and male fish fertility. TUNEL detection revealed a small amount of apoptosis in the testes of rln3b-/- male fish at 390 days after hatching (dah). RT-qPCR analysis demonstrated that mutation of rln3b gene caused a significant downregulation of steroid synthesis-related genes such as cyp17a1, cyp11b2, germ cell marker gene, Vasa, and gonadal somatic cell marker genes of amh and amhr2. Furthermore, we found a significant down-regulation of hypothalamic-pituitary-gonadal (HPG) axis-related genes, while a significantly up-regulation of the dopamine synthetase gene in the rln3b-/- male fish. Taken together, our data strongly suggested that Rln3b played a crucial role in the fertility of XY tilapia by regulating HPG axis genes.

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.

Mutant kri1l causes abnormal retinal development via cell cycle arrest and apoptosis induction

Damage to the ribosome or an imbalance in protein biosynthesis can lead to some human diseases, such as diabetic retinopathy (DR) and other eye diseases. Here, we reported that the kri1l gene was responsible for retinal development. The kri1l gene encodes an essential component of the rRNA small subunit processome. The retinal structure was disrupted in kri1l mutants, which resulted in small eyes. The boundaries of each layer of cells in the retina were blurred, and each layer of cells was narrowed and decreased. The photoreceptor cells and Müller glia cells almost disappeared in kri1l mutants. The lack of photoreceptor cells caused a fear of light response. The development of the retina started without abnormalities, and the abnormalities began two days after fertilization. In the kri1l mutant, retinal cell differentiation was defective, resulting in the disappearance of cone cells and Müller cells. The proliferation of retinal cells was increased, while apoptosis was also enhanced in kri1l mutants. γ-H2AX upregulation indicated the accumulation of DNA damage, which resulted in cell cycle arrest and apoptosis. The kri1l mutation reduced the expression of some opsin genes and key retinal genes, which are also essential for retinal development.

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.

New insights into the all-testis differentiation in zebrafish with compromised endogenous androgen and estrogen synthesis

The regulatory mechanism of gonadal sex differentiation, which is complex and regulated by multiple factors, remains poorly understood in teleosts. Recently, we have shown that compromised androgen and estrogen synthesis with increased progestin leads to all-male differentiation with proper testis development and spermatogenesis in cytochrome P450 17a1 (cyp17a1)-/- zebrafish. In the present study, the phenotypes of female-biased sex ratio were positively correlated with higher Fanconi anemia complementation group L (fancl) expression in the gonads of doublesex and mab-3 related transcription factor 1 (dmrt1)-/- and cyp17a1-/-;dmrt1-/- fish. The additional depletion of fancl in cyp17a1-/-;dmrt1-/- zebrafish reversed the gonadal sex differentiation from all-ovary to all-testis (in cyp17a1-/-;dmrt1-/-;fancl-/- fish). Luciferase assay revealed a synergistic inhibitory effect of Dmrt1 and androgen signaling on fancl transcription. Furthermore, an interaction between Fancl and the apoptotic factor Tumour protein p53 (Tp53) was found in vitro. The interaction between Fancl and Tp53 was observed via the WD repeat domain (WDR) and C-terminal domain (CTD) of Fancl and the DNA binding domain (DBD) of Tp53, leading to the K48-linked polyubiquitination degradation of Tp53 activated by the ubiquitin ligase, Fancl. Our results show that testis fate in cyp17a1-/- fish is determined by Dmrt1, which is thought to stabilize Tp53 by inhibiting fancl transcription during the critical stage of sexual fate determination in zebrafish.

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

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

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

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.

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.

Morphological and transcriptional analysis of sexual differentiation and gonadal development in a burrowing fish, the four-eyed sleeper (Bostrychus sinensis)

Four-eyed sleeper (Bostrychus sinensis) is a commercially important sea water fish, and the male individuals exhibit significant advantages in somatic growth and stress resistance, so developing sex control strategy to create all-male progeny will produce higher economic value. However, little is known about the genetic background associated with sex differentiation in this species. In this study, we investigated gonadal development and uncovered critical window stages of sexual differentiation (about 2 mph), transition from proliferation to differentiation in female germ stem cells (GSCs) (2-3 mph) and male GSCs (3-4 mph). De novo transcriptome analysis revealed candidate genes and signaling pathways associated with sexual differentiation and gonadal development in four-eyed sleeper. The results showed that sox9 and zglp1 were the earliest sex-biased transcription factors during sex differentiation. Down-regulation of chemokine, cytokines-cytokine receptors and up-regulation of cellular senescence pathway might be involved in GSC differentiation. Weighted gene correlation network analysis showed that metabolic pathway and occludin were the hub signaling and gene in ovarian development, meanwhile the MAPK signaling pathways, cellular senescence pathway and ash1l (histone H3-lysine4 N-trimethyltransferase) were the hub pathways and gene in testicular development. The present work elucidated the developmental processes of sexual differentiation and gonadal development and revealed their associated revealed genes and signaling pathways in four-eyed sleeper, providing theoretical basis for developing sex-control techniques.

Genetic evidence for differential functions of figla and nobox in zebrafish ovarian differentiation and folliculogenesis

FIGLA and NOBOX are important oocyte-specific transcription factors. Both figla-/- and nobox-/- mutants showed all-male phenotype in zebrafish due to increased dominance of the male-promoting pathway. The early diversion towards males in these mutants has precluded analysis of their roles in folliculogenesis. In this study, we attenuated the male-promoting pathway by deleting dmrt1, a key male-promoting gene, in figla-/- and nobox-/- fish, which allows a sufficient display of defects in folliculogenesis. Germ cells in figla-/-;dmrt1-/- double mutant remained in cysts without forming follicles. In contrast, follicles could form well but exhibited deficient growth in nobox-/-;dmrt1-/- double mutants. Follicles in nobox-/-;dmrt1-/- ovary could progress to previtellogenic (PV) stage but failed to enter vitellogenic growth. Such arrest at PV stage suggested a possible deficiency in estrogen signaling. This was supported by lines of evidence in nobox-/-;dmrt1-/-, including reduced expression of ovarian aromatase (cyp19a1a) and level of serum estradiol (E2), regressed genital papilla (female secondary sex characteristics), and more importantly the resumption of vitellogenic growth by E2 treatment. Expression analysis suggested Nobox might regulate cyp19a1a by controlling Gdf9 and/or Bmp15. Our discoveries indicate that Figla is essential for ovarian differentiation and follicle formation whereas Nobox is important for driving subsequent follicle development.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

Expression pattern analysis of anti-Mullerian hormone in testis development of pearlscale angelfish (Centropyge vrolikii)

In vertebrates, anti-Mullerian hormone (Amh) secreted by Sertoli cells (SC) performs a pivotal function in male sex differentiation. Compared with that of higher vertebrates, the expression pattern of Amh is more diversified in fish. In this study, the full-length complementary DNA (cDNA) of Amh in Centropyge vrolikii (Cv-Amh) was cloned and analysed, which was 2,470 bp, including a 238 bp 5'UTR, a 1,602 bp ORF and a 633 bp 3'UTR; the similarity of Amh between Cv-Amh and other fish is relatively high. The quantitative real-time PCR (qRT-PCR) results of healthy tissues and gonads at sex reversal stages in C. vrolikii showed that the expression level of Amh in the testis was significantly higher than that in other tissues (P < 0.05). Amh was weakly expressed in the vitellogenic stage ovary and perinucleolus stage ovary, but its expression significantly increased in the gonads at the hermaphroditic stage, and finally reached the highest in the pure testis after sexual reversal. The results of in situ hybridization indicated that the positive signal of Amh was strongly concentrated in SCs of testis. After Amh knockdown in the gonads, the effect on sex-related genes was tested using qRT-PCR. Among these, the expression of Dmrt1, Cyp11a, Hsd11b2, Sox8 and Sox9 significantly decreased, whereas that of Cyp19a, Sox4, Foxl2 and Sox3 increased. These results suggested that Amh could be the pivotal gene in reproductive regulation in C. vrolikii, and the data will contribute to sex-related research of C. vrolikii in the future.

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.

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.

Exposure to tris(2,6-dimethylphenyl) phosphate interferes with sexual differentiation via estrogen receptors 2a and 2b in zebrafish

Tris(2,6-dimethylphenyl) phosphate (TDMPP), an emerging organophosphate flame retardant, is frequently detected in multiple environmental media. Although TDMPP has been proven as a compound with estrogenic activity, its feminizing effects on reproductive system remain unclear. This study investigated the adverse effects of TDMPP on gonadal development by exposing zebrafish for 105 days from 15 days post-fertilization. Exposure to TDMPP (0.5 and 5 μM, corresponding to about 200 and 2000 μg/L) induced ovarian formation in aromatase mutant (cyp19a1a-/-) line which normally presents all-male phenotype for deficiency of endogenous estrogen (E2), suggesting its feminizing effect on sexual differentiation. In addition, TDMPP also interfered with other aspects of reproduction by delaying puberty onset, retarding sexual maturation, impairing gametogenesis and subfertility. Molecular docking and reporter gene assay indicated that all three nuclear estrogen receptors (nERs) can be binded to and activated by TDMPP. Using a series of nERs mutant lines, we confirmed the indispensable role of esr2a and esr2b in mediating the feminizing effects of TDMPP. Further analysis revealed that the prominent effects of TDMPP on sexual differentiation correlated to upregulation of female-promoting genes and downregulation of male-promoting genes. Taken together, the present study provided unequivocal genetic evidence for estrogenic effects of TDMPP on reproductive system and its molecular mechanisms of action.

Metzincin metalloproteases in PGC migration and gonadal sex conversion

Development of a functional gonad includes migration of primordial germ cells (PGCs), differentiations of somatic and germ cells, formation of primary follicles or spermatogenic cysts with somatic gonadal cells, development and maturation of gametes, and subsequent releasing of mature germ cells. These processes require extensive cellular and tissue remodeling, as well as broad alterations of the surrounding extracellular matrix (ECM). Metalloproteases, including MMPs (matrix metalloproteases), ADAMs (a disintegrin and metalloproteinases), and ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs), are suggested to have critical roles in the remodeling of the ECM during gonad development. However, few research articles and reviews are available on the functions and mechanisms of metalloproteases in remodeling gonadal ECM, gonadal development, or gonadal differentiation. Moreover, most studies focused on the roles of transcription and growth factors in early gonad development and primary sex determination, leaving a significant knowledge gap on how differentially expressed metalloproteases exert effects on the ECM, cell migration, development, and survival of germ cells during the development and differentiation of ovaries or testes. We will review gonad development with focus on the evidence of metalloprotease involvements, and with an emphasis on zebrafish as a model for studying gonadal sex differentiation and metalloprotease functions.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Mature Gonads in Spinibarbus hollandi

Spinibarbus hollandi is an important commercial aquaculture species in southeastern China, but with long maturity period and low egg laying amount. However, there has been little study of its gonad development and reproductive regulation, which limits aquaculture production. Here, for the first time, gonadal transcriptomes of male and female S. hollandi were analyzed. A total of 167,152 unigenes were assembled, with only 48,275 annotated successfully. After comparison, a total of 21,903 differentially expressed genes were identified between male and female gonads, of which 16,395 were upregulated and 5508 were downregulated in the testis. In addition, a large number of differentially expressed genes participating in reproduction, gonad formation and differentiation, and gametogenesis were screened out and the differential expression profiles of partial genes were further validated using quantitative real-time PCR. These results will provide basic information for further research on gonad differentiation and development in S. hollandi.

Identification and expression analysis of zebrafish gnaq in the hypothalamic–Pituitary–Gonadal axis

The G proteins have emerged as essential molecular switches in a wide variety of signal transduction pathways. Gαq, encoded by G protein subunit alpha q (gnaq), is a member of the G proteins and participates in regulating important biological activities in mammals; however, its function and regulatory mechanism in teleost remain largely unclear. In the current study, we cloned the cDNA of gnaq from zebrafish (Danio rerio) and investigated the expression characteristics of Gαq/gnaq in reproductive tissues. RT-PCR and WISH analyses showed that gnaq was widely expressed in zebrafish tissues, with high expression in the brain, olfactory brain, and hypothalamus. During the embryonic development stage, the gnaq was mainly distributed in the hypothalamus after 72 h post-fertilization. In addition, immunohistochemistry analysis revealed that the Gαq protein was highly expressed in the diffuse nucleus of the inferior hypothalamic lobe (DIL), ventral zone of the periventricular hypothalamus (Hv), and caudal zone of the periventricular hypothalamus (Hc) in adult zebrafish. Furthermore, in the gonads, the Gαq protein was found in oocytes of all stages, except spermatids. Lastly, the gnaq mRNA exhibited relatively low expression in gonads on Day 4 during the reproductive cycle, while increasing drastically in the hypothalamus and pituitary afterward. Altogether, our results suggest that gnaq/Gαq might be important in fish reproduction.

Knock-down of amh transcription by antisense RNA reduces FSH and increases follicular atresia in female Oreochromis niloticus

Anti-Müllerian hormone (Amh) plays an important role in regulating gonad development in teleosts. However, little is known about the effects of Amh on follicle development. In this study, we transfected the vector containing antisense RNA fragments of the amh gene to produce Nile tilapia, Oreochromis niloticus, with knocked-down Amh function in vivo. The results confirmed that the antisense RNA effectively inhibited amh transcription and Amh protein expression in female tilapia ovarian tissue. At 180 days of age, compared with control fish, female tilapia with knocked-down Amh function showed significantly increased growth and significantly decreased ovary weight and gonadosomatic index (P < 0.05). Female fish in the control group had ruddy-colored external genitalia, eggs extruded from the abdomen when gently squeezed, and most oocytes were developmental stage V. In contrast, the external genitalia of female fish with knocked-down Amh function did not have the ruddy color, no eggs extruded from the abdomen when squeezed, most oocytes were at developmental stages II and III, and considerable follicular atresia was apparent. At 180 days of age, the transcript levels of amhrII, cyp19a1a, foxl2 and sox9b in ovarian tissue, and the titers of luteinizing hormone, follicle stimulating hormone, and estradiol in the serum, were significantly lower in fish with knocked-down Amh function than in control fish (P < 0.05). We concluded that decreased serum hormone levels and an abnormal AMH signal delayed development and caused follicular degeneration in Nile tilapia with knocked-down Amh function. These findings show that antisense RNA is a feasible approach for gene silencing in fish, and represents an accurate and effective strategy to study gene function.

ceRNA crosstalk mediated by ncRNAs is a novel regulatory mechanism in fish sex determination and differentiation

Competing endogenous RNAs (ceRNAs) are vital regulators of gene networks in mammals. The involvement of noncoding RNAs (ncRNAs) as ceRNA in genotypic sex determination (GSD) and environmental sex determination (ESD) in fish is unknown. The Chinese tongue sole, which has both GSD and ESD mechanisms, was used to map the dynamic expression pattern of ncRNAs and mRNA in gonads during sex determination and differentiation. Transcript expression patterns shift during the sex differentiation phase, and ceRNA modulation occurs through crosstalk of differentially expressed long ncRNAs (lncRNAs), circular RNAs (circRNAs), microRNAs (miRNAs), and sex-related genes in fish. Of note was the significant up-regulation of a circRNA from the sex-determining gene dmrt1 (circular RNA dmrt1) and a lncRNA, called AMSDT (which stands for associated with male sex differentiation of tongue sole) in Chinese tongue sole testis. These two ncRNAs both share the same miRNA response elements with gsdf, which has an up-regulated expression when they bind to miRNA cse-miR-196 and concurrent down-regulated female sex-related genes to facilitate testis differentiation. This is the first demonstration in fish that ceRNA crosstalk mediated by ncRNAs modulates sexual development and unveils a novel regulatory mechanism for sex determination and differentiation.

Ovarian transcriptome analysis of diploid and triploid rainbow trout revealed new pathways related to gonadal development and fertility

Triploidisation represents several advantages (e.g. sterility) and therefore is routinely applied in aquaculture of several commercially important fish species, including rainbow trout. The comparative transcriptomic analysis of ovaries of triploid (3N) and diploid (2N) female rainbow trout revealed a total of 9 075 differentially expressed genes (DEGs; 4 105 genes upregulated in 2N and 4 970 genes upregulated in 3N ovaries, respectively). Identified clusters for DEGs upregulated in 3N and 2N ovaries were different, including carbohydrate and lipid metabolic process and transport, protein modification, signalling (related to folliculogenesis) and response to stimulus for DEGs upregulated in 2N, and developmental process, signalling (related to apoptosis, cellular senescence and adherence junctions) and regulation of RNA metabolic process for DEGs upregulated in 3N. The enrichment of processes involved in carbohydrate and lipid metabolism in 2N ovaries indicated high metabolism of ovarian tissue and the energy reservoir generation indispensable during the earliest stages of development. Our results highlight the importance of oocyte hydration along with oestrogen, insulin, leptin, fibroblast growth factor, and Notch signalling and pathways related to the regulation of cyclic adenosine monophosphate (cAMP) levels in proper oocyte meiotic maturation prior to ovulation in 2N ovaries. Conversely, triploidisation may lead to an increase in ovarian cellular senescence and apoptosis, which in turn can result in abnormal gonadal morphology and fibrosis. The downregulation of genes responsible for the precise regulation of meiosis and proper chromosome segregation during meiosis probably affects meiotic maturation via irregular meiotic division of chromosomes. The induction of triploidy of the rainbow trout genome resulted in enhanced expression of male-specific genes, genes responsible for re-establishing the transcriptional balance after genome reorganisation and genes involved in regulatory mechanisms, including gene silencing and DNA methylation. To the best of our knowledge, this is the first genome-wide investigation providing in-depth comprehensive and comparative gene expression patterns in the ovary from 2N and 3N rainbow trout females helping in elucidating the molecular mechanisms leading to impaired gonadal development and sterility of female triploids.

Function of Foxl2 and Dmrt1 proteins during gonadal differentiation in the olive flounder Paralichthys olivaceus

Study on fish sex differentiation is important both from academic and practical aspects. Foxl2 and Dmrt1 are important transcription factors that should be involved in fish gonadal differentiation, but there is still no direct evidence to clarify their protein functions. Olive flounder Paralichthys olivaceus, an important mariculture fish in China, Japan, and Korea, shows sex-dimorphic growth. In this study, the Foxl2 and Dmrt1 proteins were detected in granulosa cells of the ovary and Sertoli cells of the testis, respectively, showing significant sex-dimorphic expression patterns. Then, bioactive high-purity Foxl2 and Dmrt1 recombinant proteins were obtained in vitro. Furthermore, effects of the recombinant Foxl2 and Dmrt1 during gonadal differentiation period were evaluated by intraperitoneal injection in juvenile fish. Compared with the control group, the male rate in the Dmrt1 group increased from 0 % to 82 %, showing for the first time in fish that the recombinant Dmrt1 could alter the sex phenotype. In addition, transcription levels of cyp19a and its transcription factors also changed after the recombinant Foxl2 and Dmrt1 injection. These findings reveal that Foxl2 and Dmrt1 are vital regulators for fish gonadal differentiation by regulating cyp19a expression, and also provide a new approach for sex control in fish aquaculture.

Application of CRISPR/Cas Technology in Spermatogenesis Research and Male Infertility Treatment

As the basis of animal reproductive activity, normal spermatogenesis directly determines the efficiency of livestock production. An in-depth understanding of spermatogenesis will greatly facilitate animal breeding efforts and male infertility treatment. With the continuous development and application of gene editing technologies, they have become valuable tools to study the mechanism of spermatogenesis. Gene editing technologies have provided us with a better understanding of the functions and potential mechanisms of action of factors that regulate spermatogenesis. This review summarizes the applications of gene editing technologies, especially CRISPR/Cas9, in deepening our understanding of the function of spermatogenesis-related genes and disease treatment. The problems of gene editing technologies in the field of spermatogenesis research are also discussed.

Comprehensive Transcriptome Analysis of Gonadal and Somatic Tissues for Identification of Sex-Related Genes in the Largemouth Bass Micropterus salmoides

Largemouth bass (Micropterus salmoides) is an economically important fish. It can spawn many times during a breeding season, and there are no obvious morphological characteristics to distinguish male and female juvenile fish. So far, little is known about the genes regulating their sexual development in this species. Here, we performed RNA sequencing (RNA-Seq) analysis of the testis, ovary, and somatic tissue to identify sex-related genes in the largemouth bass. A total of 51,672 unigenes were obtained via the transcriptome analysis, and 5900 differential expression genes (DEGs), including 3028 up-regulated and 2872 down-regulated DEGs, were obtained in the somatic tissue, testis, and ovary. DEGs were retrieved by making comparisons: somatic tissue vs testis (1733-up and 1382-down), testis vs ovary (841-up and 807-down), and ovary vs somatic tissue (454-up and 683-down). Finally, functional annotation identified 22 key sex-related DEGs, including 13 testis-biased DEGs (dmrt1, cyp11b1, sox9, spata4, spata22, spata17, fshr, fem-1a, wt1, daz1, amh, vasa, and piwi1) and 9 ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, bmp15, fem-1b, fig. la, and piwi2). This result was further confirmed by the tissue expression detection via RT-PCR and RT-qPCR. Protein-protein interacting (PPI) network analysis revealed that the testis-specific dmrt1 interacts directly with the testis-biased DEGs (cyp11b1 and spata4) and the ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, and bmp15), suggesting that the dmrt1 as a sex-determining gene can play a dual role through inducing the testis-biased DEGs and inhibiting the ovary-biased DEGs during the testicular development. Our present results provide useful molecular data for a better understanding of sexual development in the largemouth bass.

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.

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.

Gdnf Acts as a Germ Cell-Derived Growth Factor and Regulates the Zebrafish Germ Stem Cell Niche in Autocrine- and Paracrine-Dependent Manners

Glial cell line-derived neurotrophic factor (GDNF) and its receptor (GDNF Family Receptor α1-GFRα1) are well known to mediate spermatogonial stem cell (SSC) proliferation and survival in mammalian testes. In nonmammalian species, Gdnf and Gfrα1 orthologs have been found but their functions remain poorly investigated in the testes. Considering this background, this study aimed to understand the roles of the Gdnf-Gfrα1 signaling pathway in zebrafish testes by combining in vivo, in silico and ex vivo approaches. Our analysis showed that zebrafish exhibit two paralogs for Gndf (gdnfa and gdnfb) and its receptor, Gfrα1 (gfrα1a and gfrα1b), in accordance with a teleost-specific third round of whole genome duplication. Expression analysis further revealed that both ligands and receptors were expressed in zebrafish adult testes. Subsequently, we demonstrated that gdnfa is expressed in the germ cells, while Gfrα1a/Gfrα1b was detected in early spermatogonia (mainly in types A_und and A_diff) and Sertoli cells. Functional ex vivo analysis showed that Gdnf promoted the creation of new available niches by stimulating the proliferation of both type A_und spermatogonia and their surrounding Sertoli cells but without changing pou5f3 mRNA levels. Strikingly, Gdnf also inhibited late spermatogonial differentiation, as shown by the decrease in type B spermatogonia and down-regulation of dazl in a co-treatment with Fsh. Altogether, our data revealed that a germ cell-derived factor is involved in maintaining germ cell stemness through the creation of new available niches, supporting the development of spermatogonial cysts and inhibiting late spermatogonial differentiation in autocrine- and paracrine-dependent manners.

Single Cell Transcriptome Sequencing of Zebrafish Testis Revealed Novel Spermatogenesis Marker Genes and Stronger Leydig-Germ Cell Paracrine Interactions

Spermatogenesis in testis is an important process for sexual reproduction, and worldwide about 10-15 percent of couples suffer from infertility. It is of importance to study spermatogenesis at single cell level in both of human and model organisms. Currently, single-cell RNA sequencing technologies (scRNA-seq) had been extensively applied to the study of cellular components and its gene regulations in the testes of different species, including human, monkey, mouse, and fly, but not in zebrafish. Zebrafish was a widely used model organism in biology and had been extensively used for the study of spermatogenesis in the previous studies. Therefore, it is also important to profile the transcriptome of zebrafish testis at single cell level. In this study, the transcriptomes of 14, 315 single cells from adult male zebrafish testes were profiled by scRNA-seq, and 10 cell populations were revealed, including Leydig cell, Sertoli cell, spermatogonia cell (SPG), spermatocyte, and spermatids. Notably, thousands of cell-type specific novel marker genes were identified, including sumo3b for SPG, krt18a.1 for Sertoli cells, larp1b and edrf1 for spermatids, which were also validated by RNA in situ hybridization experiments. Interestingly, through Ligand-Receptor (LR) analyses, zebrafish Leydig cells demonstrated stronger paracrine influence on germ cells than Sertoli cells. Overall, this study could be an important resource for the study of spermatogenesis in zebrafish and might also facilitate the study of the genes associated with human infertility through using zebrafish as a model organism.

Decoding genome recombination and sex reversal

Over the past 440 years since the discovery of the medicinal value of swamp eels, much progress has been made in the study of their biology. The fish is emerging as an important model animal in sexual development, in addition to economic and pharmaceutical implications. Tracing genomic history that shapes speciation of the fish has led to discovery of the whole genome-wide chromosome fission/fusion events. Natural intersex differentiation is a compelling feature for sexual development research. Notably, identification of progenitors of germline stem cells that have bipotential to differentiate into either male or female germline stem cells provides new insight into sex reversal. Here, we review these advances that have propelled the field forwards and present unsolved issues that will guide future investigations to finally elucidate vertebrate sexual development using the new model.

Augmentation of progestin signaling rescues testis organization and spermatogenesis in zebrafish with the depletion of androgen signaling

Disruption of androgen signaling is known to cause testicular malformation and defective spermatogenesis in zebrafish. However, knockout of cyp17a1, a key enzyme responsible for the androgen synthesis, in ar-/- male zebrafish paradoxically causes testicular hypertrophy and enhanced spermatogenesis. Because Cyp17a1 plays key roles in hydroxylation of pregnenolone and progesterone (P4), and converts 17α-hydroxypregnenolone to dehydroepiandrosterone and 17α-hydroxyprogesterone to androstenedione, we hypothesize that the unexpected phenotype in cyp17a1-/-;androgen receptor (ar)-/- zebrafish may be mediated through an augmentation of progestin/nuclear progestin receptor (nPgr) signaling. In support of this hypothesis, we show that knockout of cyp17a1 leads to accumulation of 17α,20β-dihydroxy-4-pregnen-3-one (DHP) and P4. Further, administration of progestin, a synthetic DHP mimetic, is sufficient to rescue testicular development and spermatogenesis in ar-/- zebrafish, whereas knockout of npgr abolishes the rescue effect of cyp17a1-/- in the cyp17a1-/-;ar-/- double mutant. Analyses of the transcriptomes among the mutants with defective testicular organization and spermatogenesis (ar-/-, ar-/-;npgr-/- and cyp17a-/-;ar-/-;npgr-/-), those with normal phenotype (control and cyp17a1-/-), and rescued phenotype (cyp17a1-/-;ar-/-) reveal a common link between a downregulated expression of insl3 and its related downstream genes in cyp17a-/-;ar-/-;npgr-/- zebrafish. Taken together, our data suggest that genetic or pharmacological augmentation of the progestin/nPgr pathway is sufficient to restore testis organization and spermatogenesis in zebrafish with the depletion of androgen signaling.

Loss of Nobox prevents ovarian differentiation from juvenile ovaries in zebrafish

As a species without master sex-determining genes, zebrafish displays high plasticity in sex differentiation, making it an excellent model for studying the regulatory mechanisms underlying gonadal differentiation and gametogenesis. Despite being a gonochorist, zebrafish is a juvenile hermaphrodite that undergoes a special phase of juvenile ovary before further differentiation into functional testis and ovary. The mechanisms underlying juvenile ovary formation and subsequent gonadal differentiation remain largely unknown. In a recent study, we demonstrated an important role for Figla (factor in the germline alpha) in zebrafish oogenesis. In this study, we explored the role of Nobox/nobox (new born ovary homeobox protein), another oocyte-specific transcription factor in females, in early zebrafish gonadogenesis using CRISPR/Cas9 technology. As in mammals, nobox is specifically expressed in zebrafish gonads with a dimorphic pattern at juvenile stage. In contrast to the mutant of figla (another oocyte-specific transcription factor), the nobox mutants showed formation of typical perinucleolar (PN) follicles at primary growth (PG) stage in juvenile gonads, suggesting occurrence of follicle assembly from cystic oocytes (chromatin nucleolar stage, CN). These follicles, however, failed to develop further to form functional ovaries, resulting in all-male phenotype. Despite its expression in adult testis, the loss of nobox did not seem to affect testis development, spermatogenesis and male spawning. In summary, our results indicate an important role for Nobox in zebrafish ovarian differentiation and early folliculogenesis.

Cyp11a2 Is Essential for Oocyte Development and Spermatogonial Stem Cell Differentiation in Zebrafish

Cytochrome P45011A1, encoded by Cyp11a1, converts cholesterol to pregnenolone (P5), the first and rate-limiting step in steroidogenesis. In zebrafish, cyp11a1 is maternally expressed and cyp11a2 is considered the ortholog of Cyp11a1 in mammals. A recent study has shown that depletion of cyp11a2 resulted in steroidogenic deficiencies and the mutants developed into males with feminized secondary sexual characteristics. Here, we independently generated cyp11a2 mutants in zebrafish and showed that the mutants can develop into males and females in the juvenile stage, but finally into infertile males with defective mating behavior in the adult stage. In the developing ovaries, the cyp11a2 mutation led to stage I oocyte apoptosis and final sex reversal, which could be partially rescued by treatment with P5 but not estradiol. In the developing testes, depletion of cyp11a2 resulted in dysfunction of Sertoli cells and lack of functional Leydig cells. Spermatogonial stem cells (SSCs) in the mutant testes underwent active self-renewal but no differentiation, resulting in a high abundance of SSCs in the testis, as revealed by immunofluorescence staining with Nanos2 antibody. The high abundance and differentiation competence of SSCs in the mutant testes were verified by a novel testicular cell transplantation method developed in this study, by transplanting mutant testicular cells into germline-depleted wild-type (WT) fish. The transplanted mutant SSCs efficiently differentiated into functional spermatids in WT hosts. Overall, our study demonstrates the functional importance of cyp11a2 in early oogenesis and differentiation of SSCs.

Roles of anti-Müllerian hormone and its duplicates in sex determination and germ cell proliferation of Nile tilapia

Duplicates of amh are crucial for fish sex determination and differentiation. In Nile tilapia, unlike in other teleosts, amh is located on X chromosome. The Y chromosome amh (amhΔ-y) is mutated with 5 bp insertion and 233 bp deletion in the coding sequence, and tandem duplicate of amh on Y chromosome (amhy) has been identified as the sex determiner. However, the expression of amh, amhΔ-y, and amhy, their roles in germ cell proliferation and the molecular mechanism of how amhy determines sex is still unclear. In this study, expression and functions of each duplicate were analyzed. Sex reversal occurred only when amhy was mutated as revealed by single, double, and triple mutation of the 3 duplicates in XY fish. Homozygous mutation of amhy in YY fish also resulted in sex reversal. Earlier and higher expression of amhy/Amhy was observed in XY gonads compared with amh/Amh during sex determination. Amhy could inhibit the transcription of cyp19a1a through Amhr2/Smads signaling. Loss of cyp19a1a rescued the sex reversal phenotype in XY fish with amhy mutation. Interestingly, mutation of both amh and amhy in XY fish or homozygous mutation of amhy in YY fish resulted in infertile females with significantly increased germ cell proliferation. Taken together, these results indicated that up-regulation of amhy during the critical period of sex determination makes it the sex-determining gene, and it functions through repressing cyp19a1a expression via Amhr2/Smads signaling pathway. Amh retained its function in controlling germ cell proliferation as reported in other teleosts, while amhΔ-y was nonfunctionalized.

Sexual determination in zebrafish

Zebrafish have emerged as a major model organism to study vertebrate reproduction due to their high fecundity and external development of eggs and embryos. The mechanisms through which zebrafish determine their sex have come under extensive investigation, as they lack a definite sex-determining chromosome and appear to have a highly complex method of sex determination. Single-gene mutagenesis has been employed to isolate the function of genes that determine zebrafish sex and regulate sex-specific differentiation, and to explore the interactions of genes that promote female or male sexual fate. In this review, we focus on recent advances in understanding of the mechanisms, including genetic and environmental factors, governing zebrafish sex development with comparisons to gene functions in other species to highlight conserved and potentially species-specific mechanisms for specifying and maintaining sexual fate.

Reproductive health and endocrine disruption in smallmouth bass (Micropterus dolomieu) from the Lake Erie drainage, Pennsylvania, USA

Smallmouth bass Micropterus dolomieu were sampled from three sites within the Lake Erie drainage (Elk Creek, Twentymile Creek, and Misery Bay, an embayment in Presque Isle Bay). Plasma, tissues for histopathological analyses, and liver and testes preserved in RNALater® were sampled from 30 smallmouth bass (of both sexes) at each site. Liver and testes samples were analyzed for transcript abundance with Nanostring nCounter® technology. Evidence of estrogenic endocrine disruption was assessed by the presence and severity of intersex (testicular oocytes; TO) and concentrations of plasma vitellogenin in male fish. Abundance of 17 liver transcripts associated with reproductive function, endocrine activity, and contaminant detoxification pathways and 40 testes transcripts associated with male and female reproductive function, germ cell development, and steroid biosynthesis were also measured. Males with a high rate of TO (87-100%) and plasma vitellogenin were noted at all sites; however, TO severity was greatest at the site with the highest agricultural land cover. Numerous transcripts were differentially regulated among the sites and patterns of transcript abundance were used to better understand potential risk factors for estrogenic endocrine disruption. The results of this study suggest endocrine disruption is prevalent in this region and further research would benefit to identify the types of contaminants that may be associated with the observed biological effects.

Functional Activity of Recombinant Forms of Amh and Synergistic Action with Fsh in European Sea Bass Ovary

Although anti-Müllerian hormone (AMH) has classically been correlated with the regression of Müllerian ducts in male mammals, involvement of this growth factor in other reproductive processes only recently come to light. Teleost is the only gnathostomes that lack Müllerian ducts despite having amh orthologous genes. In adult teleost gonads, Amh exerts a role in the early stages of germ cell development in both males and females. Mechanisms involving the interaction of Amh with gonadotropin- and growth factor-induced functions have been proposed, but our overall knowledge regarding Amh function in fish gonads remains modest. In this study, we report on Amh actions in the European sea bass ovary. Amh and type 2 Amh receptor (Amhr2) are present in granulosa and theca cells of both early and late-vitellogenic follicles and cannot be detected in previtellogenic ovaries. Using the Pichia pastoris system a recombinant sea bass Amh has been produced that is endogenously processed to generate a 12–15 kDa bioactive mature protein. Contrary to previous evidence in lower vertebrates, in explants of previtellogenic sea bass ovaries, mature Amh has a synergistic effect on steroidogenesis induced by the follicle-stimulating hormone (Fsh), increasing E2 and cyp19a1a levels.

Sex chromosome evolution among amniotes: is the origin of sex chromosomes non-random?

Sex chromosomes are a great example of a convergent evolution at the genomic level, having evolved dozens of times just within amniotes. An intriguing question is whether this repeated evolution was random, or whether some ancestral syntenic blocks have significantly higher chance to be co-opted for the role of sex chromosomes owing to their gene content related to gonad development. Here, we summarize current knowledge on the evolutionary history of sex determination and sex chromosomes in amniotes and evaluate the hypothesis of non-random emergence of sex chromosomes. The current data on the origin of sex chromosomes in amniotes suggest that their evolution is indeed non-random. However, this non-random pattern is not very strong, and many syntenic blocks representing putatively independently evolved sex chromosomes are unique. Still, repeatedly co-opted chromosomes are an excellent model system, as independent co-option of the same genomic region for the role of sex chromosome offers a great opportunity for testing evolutionary scenarios on the sex chromosome evolution under the explicit control for the genomic background and gene identity. Future studies should use these systems more to explore the convergent/divergent evolution of sex chromosomes. 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 II)'.