A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Somatostatin signalling coordinates energy metabolism allocation to reproduction in zebrafish

BACKGROUND

Energy allocation between growth and reproduction determines puberty onset and fertility. In mammals, peripheral hormones such as leptin, insulin and ghrelin signal metabolic information to the higher centres controlling gonadotrophin-releasing hormone neurone activity. However, these observations could not be confirmed in lower vertebrates, suggesting that other factors may mediate the energetic trade-off between growth and reproduction. A bioinformatic and experimental study suggested co-regulation of the circadian clock, reproductive axis and growth-regulating genes in zebrafish. While loss-of-function of most of the identified co-regulated genes had no effect or only had mild effects on reproduction, no such information existed about the co-regulated somatostatin, well-known for its actions on growth and metabolism.

RESULTS

We show that somatostatin signalling is pivotal in regulating fecundity and metabolism. Knock-out of zebrafish somatostatin 1.1 (sst1.1) and somatostatin 1.2 (sst1.2) caused a 20-30% increase in embryonic primordial germ cells, and sst1.2-/- adults laid 40% more eggs than their wild-type siblings. The sst1.1-/- and sst1.2-/- mutants had divergent metabolic phenotypes: the former had 25% more pancreatic α-cells, were hyperglycaemic and glucose intolerant, and had increased adipocyte mass; the latter had 25% more pancreatic β-cells, improved glucose clearance and reduced adipocyte mass.

CONCLUSIONS

We conclude that somatostatin signalling regulates energy metabolism and fecundity through anti-proliferative and modulatory actions on primordial germ cells, pancreatic insulin and glucagon cells and the hypothalamus. The ancient origin of the somatostatin system suggests it could act as a switch linking metabolism and reproduction across vertebrates. The results raise the possibility of applications in human and animal fertility.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

Letrozole induced a polycystic ovary syndrome model in zebrafish by interfering with the hypothalamic-pituitary-gonadal axis

Polycystic ovary syndrome (PCOS) is a common endocrine and metabolic disease in women of childbearing age, with an incidence of 5-10%. This study compared the traits of zebrafish with three diagnostic criteria for human PCOS, and the diagnostic criteria for zebrafish PCOS were proposed: decreased fecundity, elevated testosterone (T) or 11-ketotestosterone (11-KT) levels and increased cortical-alveolar oocyte (CO) ratio, enhancing the zebrafish PCOS model's accuracy. According to the mammalian PCOS classification, the type of zebrafsh PCOS is divided into four phenotypes (A, B, C and D), but the four phenotypes of zebrafish PCOS are not fully covered in the existing studies (A and D). In this study, we successfully induced phenotype B zebrafish PCOS model using the aromatase inhibitor, letrozole (LET). That is, wild-type female zebrafish were exposed to 1000 μg/L LET for 30 days. Reproductive tests showed decreased fecundity in female zebrafish exposed to LET (Control: 132.63, 146.00, 173.00; LET: 29.20, 90.00, 82.71). Hormone analysis showed that female zebrafish exposed to LET had significantly lower 17β-estradiol/testosterone (E2/T) ratios, indicating elevated T levels. Meanwhile, levels of 11-KT in the ovaries exposed to LET were significantly up-regulated (Control: 0.0076 pg/μg; LET: 0.0138 pg/μg). Pathological sections of the ovary showed fewer CO in the LET-exposed group (Control: 16.27%; LET: 8.38%). In summary, the zebrafish PCOS model summarized and studied in this study provide a reliable and economical tool for the screening of therapeutic drugs, as well as for the etiology research and treatment strategies of PCOS.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost 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.

Gonadal transcriptome analysis of paradise fish Macropodus opercularis to reveal sex-related genes

Macropodus opercularis is an ornamental fish species endemic to China, with obvious sexual dimorphism in phenotype. To obtain the gene expression profile of the gonads of M. opercularis and explore its sex-related genes, six cDNA libraries were constructed from the sexually mature M. opercularis, and RNA-seq analysis was performed. The sequenced clean data were assembled by de novo splicing to generate 171,415 unigenes, and differentially expressed genes (DEGs) screening revealed that there were 41,638 DEGs in the gonads of M. opercularis. By comparing those DEGS in the ovary with the testis, we found 29,870 DEGs were upregulated and 11,768 DEGs were downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis showed that GO terms related to cell cycle and gamete formation were enriched, and pathway signals related to sex differences, such as FoxO signalling pathway and PI3K-Akt signalling pathway, were also detected. Reverse transcript fluorescence quantitative PCR (RT-qPCR) validation of 14 DEGs associated with sex differences showed that the RT-qPCR results were consistent with RNA-Seq analysis, and five genes, foxl2, sox3, foxo, zar1, cyp19a1, were significantly expressed in the ovaries. dmrt1, cyp11b, amh, sf1, sox9, gdf6, dmrt3, fstl1 and hsd11b2, a total of nine genes were significantly expressed in the testis. The results of this study provide a basis for the study of gonadal differentiation, developmental mechanisms and related functional genes in M. opercularis.

Tefluthrin induced toxicities in zebrafish: Focusing on enantioselectivity

Tefluthrin is one of widely used chiral pyrethroid pesticides. The potential enantioselective risk posed by tefluthrin to the aquatic ecosystem is still unclear. In this study, the toxicity differences and corresponding mechanism of tefluthrin on zebrafish were investigated at the enantiomeric level. The results indicated that two tefluthrin enantiomers showed different acute toxicity, developmental toxicity and oxidative stress to zebrafish. The acute toxicity of (1R,3R)-tefluthrin was 130-176 fold as that of (1S,3S)-tefluthrin on zebrafish embryos, larvae and adults. (1R,3R)-Tefluthrin presented approximately 10, 3 and 2 times inhibition effect on the deformity rate, hatching rate and spontaneous movements on embryos as that of (1S,3S)-tefluthrin. Meanwhile, (1R,3R)-tefluthrin caused stronger oxidative stress on zebrafish embryo than (1S,3S)-tefluthrin. The molecular docking results revealed that there were stereospecific binding affinities between tefluthrin enantimers and sodium channel protein (Nav1.6), which may lead to acute toxicity differences. Transcriptome analysis showed that the two tefluthrin enantiomers markedly disturbed differential embryonic genes expression, thereby potentially causing the chronic enantioselective toxicity. The findings of the study reveal the toxicity differences and potential mechanism of tefluthrin enantiomers on zebrafish. These results also provides a foundation for a systematic evaluation of tefluthrin at enantiomer level.

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.

Cortisol safeguards oogenesis by promoting follicular cell survival

The role of glucocorticoids in oogenesis remains to be elucidated. cyp11c1 encodes the key enzyme involved in the synthesis of cortisol, the major glucocorticoid in teleosts. In our previous study, we mutated cyp11c1 in tilapia and analyzed its role in spermatogenesis. In this study, we analyzed its role in oogenesis. cyp11c1^+/- XX tilapia showed normal ovarian morphology but poor egg quality, as indicated by the mortality of embryos before 3 d post fertilization, which could be partially rescued by the supplement of exogenous cortisol to the mother fish. Transcriptome analyses revealed reduced expression of maternal genes in the eggs of the cyp11c1^+/- XX fish. The cyp11c1^-/- females showed impaired vitellogenesis and arrested oogenesis due to significantly decreased serum cortisol. Further analyses revealed decreased serum E2 level and expression of amh, an important regulator of follicular cell development, and increased follicular cell apoptosis in the ovaries of cyp11c1^-/- XX fish, which could be rescued by supplement of either exogenous cortisol or E2. Luciferase assays revealed a direct regulation of cortisol and E2 on amh transcription via GRs or ESRs. Taken together, our results demonstrate that cortisol safeguards oogenesis by promoting follicular cell survival probably via Amh signaling.

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene’s size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix–loop–helix domain that is typical of the ovary development regulator—factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

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.

Liposome-Encapsulated Rec8 and Dmrt1 Plasmids Induce Red-Spotted Grouper (Epinephelus akaara) Testis Maturation

In fish, the maturity of gonads plays an important role in the development and reproduction of the population, and it also dictates the success of captive breeding. Therefore, finding ways to promote gonadal maturation is an important goal in aquaculture. In this study, we injected recombinant dmrt1 and rec8 overexpression plasmids packaged in liposomes into the immature testis of red-spotted grouper (Epinephelus akaara) and measured the expression of Dmrt1 and Rec8 protein in vivo. Gonadosomatic index (GSI) and gonadal histology analyses showed that the testis developed from the immature to the mature state within 7 days after plasmid injection. Additionally, the spermatozoa concentration and motility in plasmid-injected fish was the same as that of naturally mature fish. These results provided evidence that delivery of dmrt1 and rec8 expression plasmids into the testis via injection induced testis maturation in vivo.

Intensive masculinization caused by chronic heat stress in juvenile Cynoglossus semilaevis: Growth performance, gonadal histology and gene responses

The sea temperature has been observed to chronically increase during the past decades, leaving unpredictable influences to the marine biological resources. Thus, it is of vital significance to study the biological responses of ocean inhabited organisms with the artificially stimulated heat stress environment. Cynoglossus semilaevis provides us with an ideal model to study the influence of chronic heat stress on the sexual differentiation in marine teleosts for its genetic sex determination (GSD) + environmental effected (EE) sex determination system. In this study, the comparative experiment was conducted employing heated seawater (HT group) and ambient seawater (CT group) to cultivate juvenile C. semilaevis respectively. Significant differences were exhibited in growth performance and a delayed germ cell development effect was found in pseudomales formed under chronic heat stress. Using transcriptome analysis, the transcription profile of 55 days post fertilization (dpf) and 100 dpf juveniles' gonads were studied. A total of 47 libraries were constructed with an average mapping rate of 94.63% after assembling. GO and KEGG enrichment were proceeded using DEGs screened out between (1) pseudomale gonads at 55 dpf and 100 dpf in HT and CT group (2) pseudomale and female gonads at 55 dpf and 100 dpf in HT and CT group. Terms and pathways involved in steroid stimulation, reproduction ability, germ cell proliferation et al. were shed light on. The expression pattern of 29 DEGs including amh, hsp90b1, pgr et al. were also provided to supplement the results of functional enrichment. Weighted gene co-expression networks analysis (WGCNA) was constructed and hspb8-like, histone H2A.V were exhibited to play vital roles in the heat-induced masculinization. Our findings facilitate the understanding for transcriptional variations in intensive masculinization cause by chronic heat stress of C. semilaevis and provide referable study of the influences on the teleosts in elevated sea temperature.

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.

Evidence-based hormonal, mutational, and endocrine-disrupting chemical-induced zebrafish as an alternative model to study PCOS condition similar to mammalian PCOS model

Polycystic ovarian syndrome (PCOS) causes swollen ovaries in women at reproductive age due to hormonal disorder with small cysts on the outer edges. The cause of the disorder is still yet to be found. Multiple factors have increased PCOS prevalence, hyperandrogenism, oxidative stress, inflammation, and insulin resistance. Various animal PCOS models have been developed to imitate the pathophysiology of PCOS in humans. Zebrafish is one of the most versatile animal experimental models because of the transparency of the embryos, small size, and rapid growth. The zebrafish similarity to higher vertebrates made it a useful non-mammalian model for PCOS drug testing and screening. This review provides an insight into the usage of zebrafish, a non-mammalian model for PCOS, as an opportunity for evaluating future initiatives in such a research domain.

Tissue-Specific Expression Pattern in Ancherythroculter nigrocauda, a Sexually Size Dimorphic Fish

Certain members of the Actinopterygii class are known to exhibit sexual dimorphism (SD) that results in major phenotypic differences between male and female fishes of a species. One of the most common differences between the two sexes is in body weight, a factor with a high economic value in aquaculture. In this study, we used RNA sequencing (RNA-seq) to study the liver and brain transcriptomes of Ancherythroculter nigrocauda, a fish exhibiting SD. Females attain about fourfold body weight of males at sexual maturity. Sample clustering showed that both sexes were grouped well with their sex phenotypes. In addition, 2,395 and 457 differentially expressed genes (DEGs) were identified in the liver and brain tissues, respectively. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses predicted the association of PPAR signaling, cytochrome P450, and steroid hormone biosynthesis to the differences in sexual size. In addition, weighted gene co-expression network analyses (WGCNA) were conducted, and the green module was identified to be significantly correlated with sexual size dimorphism (SSD). Altogether, these results improve our understanding of the molecular mechanism underlying SSD in A. nigrocauda.

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.

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.

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.

A sex-linked SNP mutation in amhr2 is responsible for male differentiation in obscure puffer (Takifugu obscurus)

BACKGROUND

Anti-Mullerian hormone receptor type II (Amhr2) is a key receptor of Amh signaling in regulating gonad development. The amhr2 gene has been identified in numerous species, including a few teleost fishes. However, the roles of Amhr2 in Amh signaling in fish are poorly studied.

METHODS AND RESULTS

In this study, an amhr2 homolog from obscure puffer (Takifugu obscurus) was identified, and its molecular characteristics were systematically analyzed. Expression analysis revealed that amhr2 was highly expressed in the gonads of adult pufferfish and significantly upregulated during sex differentiation. Significantly, a sex-linked SNP site was verified in obscure puffer amhr2. Females exhibited a homozygous genotype (C/C), while males possessed a heterozygous genotype (C/G), resulting in an amino acid variation (His/Asp384) in the kinase domain of Amhr2. Then, the functions of the different Amhr2 genotypes were further investigated. The male genotype protein (Amhr2^D384) showed an enhanced ability to interact with the type I receptor (Bmpr1a) compared to the female genotype (Amhr2^H384). The phosphorylation levels of Smads and activity of the target gene (id3) induced by the male genotype were also much higher than those induced by the female genotype. These results confirmed that the male genotype had an enhanced effect on the Amh signaling pathway compared with the female genotype.

CONCLUSIONS

This study provides direct experimental evidence for the roles of different Amhr2 genotypes in pufferfish and suggests that amhr2 is responsible for male sex differentiation in obscure puffer.

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.

Molecular characterization and expression analysis of anti-Müllerian hormone in common carp (Cyprinus carpio) adult testes

Anti-Müllerian hormone (Amh) is a member of the transforming growth factor-β (Tgf-β) superfamily required in the regression of Müllerian ducts during gonadal sex differentiation of higher vertebrates. Teleost fish lack Müllerian ducts, but identified Amh orthologs have been shown to exert crucial functions during sex determination and differentiation of several species of teleosts. However, the function of Amh during gametogenesis in adult fish remains poorly investigated. Therefore, to expand present knowledge on the role of Amh in teleosts, the present study aimed to isolate and clone full-length amh cDNA in the common carp, Cyprinus carpio, and examine its expression levels throughout the male reproductive cycle and in response to different hormone treatments of testicular explants. Molecular cloning and characterization showed that the common carp Amh precursor amino acid sequence shared common features to other fish Amh precursors, including a conserved C-terminus (Tgf-β domain) and a double proteolytic cleavage site (R-X-X-R-X-X-R) upstream to the Tgf-β domain. Expression analysis showed amh dimorphic expression in the adult gonads with higher expression in the testes than ovaries. In testes, amh mRNA was detected in Sertoli cells contacting different types of germ cells, although the expression was greatest in Sertoli cells associated with type A undifferentiated spermatogonia. Expression analysis during the reproductive cycle showed that amh transcripts were down-regulated during the developing phase, which is characterized by an increased proliferation of type A undifferentiated spermatogonia and Sertoli cells and appearance of spermatocytes (meiosis) in the testes. Furthermore, ex vivo experiments showed that a 7 day exposure to Fsh or estrogens was required to decrease amh mRNA levels in common carp testicular explants. In summary, this study provided information on the molecular characterization and transcript abundance of amh in common carp adult testes. Altogether, these data will be useful for further investigations on sex determination and differentiation in this species, and also to improved strategies for improved carp aquaculture, such as inhibiting precocious maturation of males.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

Glycan expression as a tool for a deeper understanding of a reproductive gland in a skate of economic importance

Atlantoraja platana is an oviparous skate endemic to the south-west Atlantic Ocean, and is one of the skate species most exploited by local industrial bottom trawl fisheries. Oviparous elasmobranchs encapsulate their eggs in complex egg cases produced by the oviductal gland (OG). This organ is exclusively present in these fishes and comprises four distinct zones: club, baffle, papillary and terminal. The relative size and structural complexity of these zones correlate with mode of reproduction. Glycans are known to play major roles in reproduction so their distribution in each zone of the OG could explain the functional multiplicity of the gland in skates, but this topic has not been previously investigated. In this study, morphological, histochemical and lectin-histochemical analysis revealed various novel aspects of A. platana's OG. The club, papillary and terminal zones positively stained for periodic acid Schiff's reagent (PAS) and Alcian Blue (AB), indicating the presence of neutral and acid mucopolysaccharides. However, the buffle zone was negative for PAS and AB stains, but was positive for all the lectins used. Each zone of the OG had a characteristic pattern of glycan expression. Finally, we confirmed the presence of sperm but not sperm storage. This is the first lectin-histochemical study of the OG in chondrichtyan fish and it has proven to be an important tool to understand some of the mechanisms of fertility and reproductive success in economic important species such as A. platana.

Knockout of myoc Provides Evidence for the Role of Myocilin in Zebrafish Sex Determination Associated with Wnt Signalling Downregulation

Myocilin is a secreted glycoprotein with a poorly understood biological function and it is mainly known as the first glaucoma gene. To explore the normal role of this protein in vivo we developed a myoc knockout (KO) zebrafish line using CRISPR/Cas9 genome editing. This line carries a homozygous variant (c.236_239delinsAAAGGGGAAGGGGA) that is predicted to result in a loss-of-function of the protein because of a premature termination codon p.(V75EfsX60) that resulted in a significant reduction of myoc mRNA levels. Immunohistochemistry showed the presence of myocilin in wild-type embryonic (96 h post-fertilization) anterior segment eye structures and caudal muscles. The protein was also detected in different adult ocular and non-ocular tissues. No gross macroscopic or microscopic alterations were identified in the KO zebrafish, but, remarkably, we observed absence of females among the adult KO animals and apoptosis in the immature juvenile gonad (28 dpf) of these animals, which is characteristic of male development. Transcriptomic analysis showed that adult KO males overexpressed key genes involved in male sex determination and presented differentially expressed Wnt signalling genes. These results show that myocilin is required for ovary differentiation in zebrafish and provides in vivo support for the role of myocilin as a Wnt signalling pathway modulator. In summary, this myoc KO zebrafish line can be useful to investigate the elusive function of this protein, and it provides evidence for the unexpected function of myocilin as a key factor in zebrafish sex determination.

Genetic evidence for Amh modulation of gonadotropin actions to control gonadal homeostasis and gametogenesis in zebrafish and its noncanonical signaling through Bmpr2a receptor

Anti-Müllerian hormone (Amh) plays an important role in gonadal function. Amh deficiency causes severe gonadal dysgenesis and dysfunction in zebrafish, with gonadal hypertrophy in both sexes. However, its mechanism of action remains unknown. Intriguingly, the Amh cognate type II receptor (Amhr2) is missing in the zebrafish genome, in sharp contrast to other species. Using a series of zebrafish mutants (amh, fshb, fshr and lhcgr), we provided unequivocal evidence for actions of Amh, via modulation of gonadotropin signaling, on both germ cell proliferation and differentiation. The gonadal hypertrophy in amh mutants was abolished in the absence of Fshr in females or Fshr/Lhcgr in males. Furthermore, we demonstrated that knockout of bmpr2a, but not bmpr2b, phenocopied all phenotypes of the amh mutant in both sexes, including gonadal hypertrophy, hyperproliferation of germ cells, retarded gametogenesis and reduced fshb expression. In summary, the present study provided comprehensive genetic evidence for an intimate interaction of gonadotropin and Amh pathways in gonadal homeostasis and gametogenesis and for Bmpr2a as the possible missing link for Amh signaling in zebrafish.

Anti-Müllerian hormone (Amh/amh) plays dual roles in maintaining gonadal homeostasis and gametogenesis in zebrafish

Anti-Müllerian hormone (AMH/Amh) plays a role in gonadal differentiation and function across vertebrates. In zebrafish we demonstrated that Amh deficiency caused severe gonadal dysgenesis and dysfunction. The mutant gonads showed extreme hypertrophy with accumulation of early germ cells in both sexes, namely spermatogonia in the testis and primary growth oocytes in the ovary. In amh mutant females, the folliculogenesis was normal in young fish but receded progressively in adults, which was accompanied by progressive decrease in follicle-stimulating hormone (fshb) expression. Interestingly the expression of fshb increased in the pituitary of juvenile amh mutant males but decreased in adults. The upregulation of fshb in mutant male juveniles was likely one of the mechanisms for triggering gonadal hypergrowth, whereas the downregulation of fshb in adults might involve a negative feedback by gonadal inhibin. Further analysis using mutants of fshb and growth differentiation factor 9 (gdf9) provided evidence for a role of FSH in triggering ovarian hypertrophy in young female amh mutant as well. In summary, the present study provided comprehensive genetic evidence for dual roles of Amh in controlling zebrafish gonadal homeostasis and gametogenesis in both sexes. Amh suppresses proliferation or accumulation of early germ cells (spermatogonia in testis and primary growth oocytes in ovary) while promoting their exit to advanced stages, and its action may involve both endocrine and paracrine pathways.

Zebrafish as a model for studying ovarian development: Recent advances from targeted gene knockout studies

Ovarian development is a complex process controlled by precise coordination of multiple factors. The targeted gene knockout technique is a powerful tool to study the functions of these factors. The successful application of this technique in mice in the past three decades has significantly enhanced our understanding on the molecular mechanism of ovarian development. Recently, with the advent of genome editing techniques, targeted gene knockout research can be carried out in many species. Zebrafish has emerged as an excellent model system to study the control of ovarian development. Dozens of genes related to ovarian development have been knocked out in zebrafish in recent years. Much new information and perspectives on the molecular mechanism of ovarian development have been obtained from these mutant zebrafish. Some findings have challenged conventional views. Several genes have been identified for the first time in vertebrates to control ovarian development. Focusing on ovarian development, the purpose of this review is to briefly summarize recent findings using these gene knockout zebrafish models, and compare these findings with mammalian models. These established mutants and rapid development of gene knockout techniques have prompted zebrafish as an ideal animal model for studying ovarian development.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.