Characterization and expression pattern of zebrafish Anti-Müllerian hormone (Amh) relative to sox9a, sox9b, and cyp19a1a, during gonad development

Insights into the role of FoxL2 in tebuconazole-induced male- biased sex differentiation of zebrafish

Tebuconazole (TEB) is a commonly used fungicide that inhibits the aromatase Cyp19A and downregulates the transcription factor forkhead box L2 (FoxL2), leading to male-biased sex differentiation in zebrafish larvae. However, the specific mechanism by which FoxL2 functions following TEB exposure remains unclear. In this study, the phosphorylation sites and kinase-specific residues in zebrafish FoxL2 protein (zFoxL2) were predicted. Subsequently, recombinant zFoxL2 was prepared via prokaryotic expression, and a polyclonal rabbit-anti-zFoxL2 antibody was generated. Zebrafish fibroblast (ZF4) cells were exposed to 100-μM TEB alone for 8 h, after which changes in the expression of genes involved in the foxl2 regulatory pathway (akt1, pi3k, cyp19a1b, c/ebpb and sox9a) were detected. When co-exposed to 1-μM estradiol and 100-μM TEB, the expression of these key genes tended to be restored. Interestingly, TEB did not affect the expression of the foxl2 gene or protein but it significantly suppressed the phosphorylation of FoxL2 (pFoxL2) at serine 238 (decreased by 43.64 %, p = 0.009). Co-immunoprecipitation assays showed that, following exposure to 100-μM TEB, the total precipitated proteins in ZF4 cells decreased by 17.02 % (p = 0.029) and 31.39 % (p = 0.027) in the anti-zFoxL2 antibody group and anti-pFoxL2 (ser238) antibody group, respectively, indicating that TEB suppressed the capacity of the FoxL2 protein to bind to other proteins via repression of its own phosphorylation. The pull-down assay confirmed this conclusion. This study preliminarily elucidated that the foxl2 gene functions via post-translational regulation through hypophosphorylation of its encoded protein during TEB-induced male-biased sex differentiation.

Disruption of zebrafish sex differentiation by emerging contaminants hexafluoropropylene oxides at environmental concentrations via antagonizing androgen receptor pathways

As alternatives of perfluorooctanoic acid (PFOA), hexafluoropropylene oxide dimeric acid (HFPO-DA) and trimeric acid (HFPO-TA) have been detected increasingly in environmental media and even humans. They have been shown to exhibit reproductive toxicity to model species, but their effects on human remain unclear due to the knowledge gap in their mode of action. Herein, (anti-)androgenic effects of the two HFPOs and PFOA were investigated and underlying toxicological mechanism was explored by combining zebrafish test, cell assay and molecular docking simulation. Exposure of juvenile zebrafish to the chemicals during sex differentiation promoted feminization, with HFPO-TA acting at an environmental concentration of 1 μg/L. The chemicals inhibited proliferation of human prostate cells and transcriptional activity of human and zebrafish androgen receptors (AR), with HFPO-TA displaying the strongest potency. Molecular docking revealed that the chemicals bind to AR in a conformation similar to a known AR antagonist. Combined in vivo, in vitro and in silico results demonstrated that the chemicals disrupted sex differentiation likely by antagonizing AR-mediated pathways, and provided more evidence that HFPO-TA is not a safe alternative to PFOA.

Molecular and Physiological Effects of 17α-methyltestosterone on Sex Differentiation of Black Rockfish, Sebastes schlegelii

It is widely known that all-female fish production holds economic value for aquaculture. Sebastes schlegelii, a preeminent economic species, exhibits a sex dimorphism, with females surpassing males in growth. In this regard, achieving all-female black rockfish production could significantly enhance breeding profitability. In this study, we utilized the widely used male sex-regulating hormone, 17α-methyltestosterone (MT) at three different concentrations (20, 40, and 60 ppm), to produce pseudomales of S. schlegelii for subsequent all-female offspring breeding. Long-term MT administration severely inhibits the growth of S. schlegelii, while short term had no significant impact. Histological analysis confirmed sex reversal at all MT concentrations; however, both medium and higher MT concentrations impaired testis development. MT also influenced sex steroid hormone levels in pseudomales, suppressing E2 while increasing T and 11-KT levels. In addition, a transcriptome analysis revealed that MT down-regulated ovarian-related genes (cyp19a1a and foxl2) while up-regulating male-related genes (amh) in pseudomales. Furthermore, MT modulated the TGF-β signaling and steroid hormone biosynthesis pathways, indicating its crucial role in S. schlegelii sex differentiation. Therefore, the current study provides a method for achieving sexual reversal using MT in S. schlegelii and offers an initial insight into the underlying mechanism of sexual reversal in this species.

Assessing testicular morphofunctionality under Roundup WG® herbicide exposure in zebrafish

Glyphosate-based herbicides, like Roundup WG® (RWG) used for a range of crops, such as corn, soybean, coffee, sugarcane, rice, apple, and citrus, can reach aquatic ecosystems and impact non-target organisms like fish. Thus, the fish were exposed to three RWG concentrations plus one negative control, which represents the concentration allowed for inland Brazilian waters and concentrations found in surface water worldwide (0.0, 0.065, 0.65, and 6.5 mg a.i./L) for 7 and 15 days. Morphological analysis revealed significant alterations in the testicular structure, particularly in Sertoli cell extensions and cytoplasmic bridges between germ cells. Subcellular compartments also displayed alterations, including dilated mitochondria and the loss of electron density and autophagic vesicles. Gene transcript levels related to autophagy and steroidogenic regulation were upregulated in exposed fish. Germ cell quality was also affected, increasing ROS (reactive oxygen species) production and DNA fragmentation. The study highlighted the RWG reproductive toxicity, providing valuable insights into understanding the morphofunctional alterations in somatic and germ cells of Danio rerio. In conclusion, the environmental relevant concentrations used in this study were toxic to male somatic and germ cells, which raises a concern about the concentrations considered safe for human and animal use.

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.

Determination of the timing of early gonadal differentiation in silver pomfret, Pampus argenteus

Silver pomfret is a species of global significance due to its high nutritional in fisheries sector. To accurately ascertain the timing of sex differentiation mechanism and mRNA level in this species, this study examined gonad morphology and patterns of gene expression related to sex differentiation in males and females from 51 to 180 days post hatch (dph), the temperature of water was maintained at 26 ± 1 ℃. Distinct morphological differentiation of the silver pomfret ovaries, marked by the emergence of primary oocytes, became apparent from 68 dph. By 108 dph, the testes began to differentiate, as evidenced by the appearance of the efferent duct. Early oocytes exhibited a diameter ranged from 0.077 mm to 0.682 mm, with an average diameter of 0.343 ± 0.051 mm. The proportions of various types of germ cells within the testes were subjected to analysis. The localization of Vasa during the early stages of sexual differentiation was a subject to analysis as well. Vasa was predominantly localized within the cytoplasm of gonocyte, peri-nucleolus stage oocytes, primary oocytes and type A spermatogonocytes, indicating that Vasa is involved in the early gonadal differentiation of silver pomfret. The study investigated the expression patterns of dmrt1, gsdf, amh, foxl2, cyp19a1a, cyp11a, sox3 and vasa, all of which are involved in the sex differentiation of teleosts. Among these genes, amh, gsdf, sox3, foxl2, vasa were indentified as crucial contributors to the early gonadal development of silver pomfret. Significant sex-related differences were observed in the expression patterns of amh, dmrt1, gsdf, cyp11a, sox3, cyp19a1a, vasa. This study provides novel insights into the timing of physiological changes associated with the sexual differentiation of silver pomfret. Collectively, the present data indicates that the differentiation of ovaries and testes take place approximately at 68 dph in females and 108 dph in males.

Cyp19a1a Promotes Ovarian Maturation through Regulating E2 Synthesis with Estrogen Receptor 2a in Pampus argenteus (Euphrasen, 1788)

In the artificial breeding of Pampus argenteus (Euphrasen, 1788), female fish spawn before male release sperm, which indicates rapid ovarian development. In fish, aromatase is responsible for converting androgens into estrogens and estrogen plays a crucial role in ovarian development. In this study, we aimed to investigate the potential role of brain-type and ovarian-type aromatase to study the rapid ovarian development mechanism. The results showed that cyp19a1a was mainly expressed in the ovary and could be classified as the ovarian type, whereas cyp19a1b could be considered as the brain type for its expression was mainly in the brain. During ovarian development, the expression of cyp19a1a in the ovary significantly increased from stage IV to stage V and Cyp19a1a signals were present in the follicle cells, while cyp19a1b expression in the pituitary gland decreased from stage IV to stage V. To further investigate the function of Cyp19a1a, recombinant Cyp19a1a (rCyp19a1a) was produced and specific anti-Cyp19a1a antiserum was obtained. The expressions of cyp19a1a, estrogen receptors 2 alpha (esr2a), and androgen receptor alpha (arα) were significantly upregulated in the presence of rCyp19a1a. Meanwhile, cyp19a1a was expressed significantly after E2 treatment in both ovarian and testicular tissue culture. Taken together, we found two forms of aromatase in silver pomfret. The ovarian-type aromatase might play an important role in ovarian differentiation and maturation, and participate in E2 synthesis through co-regulation with esr2a. The brain-type aromatase cyp19a1b might be involved in the regulation of both brain and gonadal development.

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.

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.

The use of in silico extreme pathway (ExPa) analysis to identify conserved reproductive transcriptional-regulatory networks in humans, mice, and zebrafish

Vertebrate sex determination and differentiation are coordinated by the activations and maintenance of reproductive transcriptional-regulatory networks (TRNs). There is considerable interest in studying the conserved design principles and functions of reproductive TRNs given that their intricate regulation is susceptible to disruption by gene mutations or exposures to exogenous endocrine disrupting chemicals (or EDCs). In this manuscript, the Boolean rules describing reproductive TRNs in humans, mice, and zebrafish, were represented as a pseudo-stoichiometric matrix model. This model mathematically described the interactions of 35 transcription factors with 21 sex determination and differentiation genes across the three species. The in silico approach of Extreme Pathway (ExPa) analysis was used to predict the extent of TRN gene activations subject to the species-specific transcriptomics data, from across various developmental life-stages. A goal of this work was to identify conserved and functional reproductive TRNs across the three species. ExPa analyses predicted the sex differentiation genes, DHH, DMRT1, and AR, to be highly active in male humans, mice, and zebrafish. Whereas FOXL2 was the most active gene in female humans and mice; and CYP19A1A in female zebrafish. These results agree with the expectation that regardless of a lack of sex determination genes in zebrafish, the TRNs responsible for canalizing male vs. female sexual differentiation are conserved with mammalian taxa. ExPa analysis therefore provides a framework with which to study the TRNs that influence the development of sexual phenotypes. And the in silico predicted conservation of sex differentiation TRNs between mammals and zebrafish identifies the piscine species as an effective in vivo model to study mammalian reproductive systems under normal or perturbed pathologies.

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.

Development of epigenetic biomarkers for the identification of sex and thermal stress in fish using DNA methylation analysis and machine learning procedures

The sex ratio is a key ecological demographic parameter crucial for population viability. However, the epigenetic mechanisms operating during gonadal development regulating gene expression and the sex ratio remain poorly understood. Moreover, there is interest in the development of epigenetic markers associated with a particular phenotype or as sentinels of environmental effects. Here, we profiled DNA methylation and gene expression of 10 key genes related to sex development and stress, including steroidogenic enzymes, and growth and transcription factors. We provide novel information on the sex-related differences and on the influence of elevated temperature on these genes in zebrafish, a species with mixed genetic and environmental influences on sex ratios. We identified both positive (e.g., amh, cyp11c and hsd11b2) and negative (e.g., cyp11a1 and dmrt1) correlations in unexposed males, and negative correlation (amh) in exposed females between DNA methylation and gene expression levels. Further, we combined DNA methylation analysis with machine learning procedures and found a series of informative CpGs capable not only of correctly identifying sex (based on cyp19a1a DNA methylation levels) but also of identifying whether males and females had been exposed to abnormally elevated temperature when young (based on amh and foxl2a DNA methylation levels, respectively). This was achieved in the absence of conspicuous morphological alterations of the gonads. These DNA methylation-based epigenetic biomarkers represent molecular resources that can correctly recapitulate past thermal history and pave the way for similar findings in other species to assess potential ecological effects of environmental disturbances in the context of climate change.

Amh/Amhr2 Signaling Causes Masculinization by Inhibiting Estrogen Synthesis during Gonadal Sex Differentiation in Japanese Flounder (Paralichthys olivaceus)

The anti-Müllerian hormone (Amh) is a protein belonging to the TGF-β superfamily, the function of which has been considered important for male sex differentiation in vertebrates. The Japanese flounder (Paralichthys olivaceus) is a teleost fish that has an XX/XY sex determination system and temperature-dependent sex determination. In this species, amh expression is up-regulated in genetic males and in temperature-induced masculinization during the sex differentiation period. However, to the best of our knowledge, no reports on the Amh receptor (Amhr2) in flounder have been published, and the details of Amh signaling remain unclear. In this study, we produced amhr2-deficient mutants using the CRISPR/Cas9 system and analyzed the gonadal phenotypes and sex-related genes. The results revealed that the gonads of genetically male amhr2 mutants featured typical ovaries, and the sex differentiation-related genes showed a female expression pattern. Thus, the loss of Amhr2 function causes male-to-female sex reversal in Japanese flounder. Moreover, the treatment of genetically male amhr2 mutants with an aromatase inhibitor fadrozole, which inhibits estrogen synthesis, resulted in testicular formation. These results strongly suggest that Amh/Amhr2 signaling causes masculinization by inhibiting estrogen synthesis during gonadal sex differentiation in the flounder.

cAMP signaling in ovarian physiology in teleosts: A review

Ovarian function in teleosts, like in other vertebrates, is regulated by two distinct gonadotropins, follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Gonadotropin effects are mediated by membrane-bound G protein-coupled receptors localized on the surface of follicle cells. Gonadotropin receptor activation results in increased intracellular cAMP, the most important second cellular signaling molecule. FSH stimulation induces the production of 17β-estradiol in the cells of growing follicles to promote vitellogenesis in oocytes. In contrast, in response to LH, fully grown post-vitellogenic follicles gain the ability to synthesize maturation-inducing steroids, which induce meiotic resumption and ovulation. All these events were induced downstream of cAMP. In this review, we summarize studies addressing the role of the cAMP pathway in gonadotropin-induced processes in teleost ovarian follicles. Furthermore, we discuss future problems concerning cAMP signaling in relation to teleost ovarian function and the differences and similarities in the gonadotropin-induced cAMP signaling pathways between mammals and teleosts.

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.

Becoming female: Ovarian differentiation from an evolutionary perspective

Differentiation of the bipotential gonadal primordium into ovaries and testes is a common process among vertebrate species. While vertebrate ovaries eventually share the same functions of producing oocytes and estrogens, ovarian differentiation relies on different morphogenetic, cellular, and molecular cues depending on species. The aim of this review is to highlight the conserved and divergent features of ovarian differentiation through an evolutionary perspective. From teleosts to mammals, each clade or species has a different story to tell. For this purpose, this review focuses on three specific aspects of ovarian differentiation: ovarian morphogenesis, the evolution of the role of estrogens on ovarian differentiation and the molecular pathways involved in granulosa cell determination and maintenance.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology

The vertebrate sodium-iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single-cell transcriptomics (scRNA-seq) on zebrafish sex-genotyped gonads revealed that NIS clade-expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

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.

Single-cell transcriptome reveals insights into the development and function of the zebrafish ovary

Zebrafish are an established research organism that has made many contributions to our understanding of vertebrate tissue and organ development, yet there are still significant gaps in our understanding of the genes that regulate gonad development, sex, and reproduction. Unlike the development of many organs, such as the brain and heart that form during the first few days of development, zebrafish gonads do not begin to form until the larval stage (≥5 dpf). Thus, forward genetic screens have identified very few genes required for gonad development. In addition, bulk RNA sequencing studies which identify genes expressed in the gonads do not have the resolution necessary to define minor cell populations that may play significant roles in development and function of these organs. To overcome these limitations, we have used single-cell RNA sequencing to determine the transcriptomes of cells isolated from juvenile zebrafish ovaries. This resulted in the profiles of 10,658 germ cells and 14,431 somatic cells. Our germ cell data represents all developmental stages from germline stem cells to early meiotic oocytes. Our somatic cell data represents all known somatic cell types, including follicle cells, theca cells and ovarian stromal cells. Further analysis revealed an unexpected number of cell subpopulations within these broadly defined cell types. To further define their functional significance, we determined the location of these cell subpopulations within the ovary. Finally, we used gene knockout experiments to determine the roles of foxl2l and wnt9b for oocyte development and sex determination and/or differentiation, respectively. Our results reveal novel insights into zebrafish ovarian development and function and the transcriptome profiles will provide a valuable resource for future studies.

Comparative Transcriptome Analyses Characterize Expression Signatures Among Males, Females, Neo-Males, and Gynogenetic Females in the Yellow Drum (Nibea albiflora)

The yellow drum (Nibea albiflora) is one of the most important marine economic fish in China, and its sexually dimorphic growth makes it preferable for mono-sex culture. Although gynogenesis and neo-male induction techniques have been established, the molecular pathways and regulatory mechanisms of sex determination and maintenance in gynogenetic females and neo-males remains far from fully understood. In this study, the gene expression profiles were investigated in the gonads and brains of wild-type male, wild-type female, neo-male, and gynogenetic female yellow drum using comparative transcriptome analyses. Generally, a total of 52,999 novel transcripts were obtained in RNA-seq, of which 45,651 were isoforms of known protein-coding genes, 1,358 novel protein-coding genes, and 5,990 long non-coding RNAs. We found that the differences between wild-type males and neo-males and between wild-type females and gynogenetic females were relatively small at both the histological and transcriptomic levels, indicating that artificial gynogenesis or hormonal sex reversal may have minimal effects on normal female or male life function, respectively. In the brain, pathways such as “Oocyte meiosis”, “Cell cycle”, and “Riboflavin metabolism” were found to be significantly enriched. In the gonads, pathways such as “Prolactin signaling pathway”, “PPAR signaling pathway”, “Cholesterol metabolism”, and “Jak-STAT signaling pathway” were found to play important roles in maintaining the regular proliferation and differentiation of females and males in yellow drum. In particular, we found that zp4 might be an effective molecular marker to differentiate between gynogenetic and normal females owing to its unique expression pattern. The results of this study may help to elucidate the molecular mechanisms involved in sex maintenance in the gonads and brain and provide basic data for genetic breeding of the yellow drum.

Masculinization of Adult Gambusia holbrooki: A Case of Recapitulation of Protogyny in a Gonochorist?

Simple Summary

Although gonochoristic fish express one sex or the other in adulthood, some adults display traits similar to those of their opposite sex, suggesting plasticity similar to hermaphrodites. To fully test this potential in the eastern mosquitofish, a gonochorist, two adult stages of females were fed 17α-Methyltestosterone incorporated feed (0–200 mg/kg diet) for 50 days. The hormone (particularly at 50 mg/kg diet) stimulated the formation of complex male copulatory structures and testicular tissue as well as upregulated expression of Anti-Müllerian Hormone gene and altered the behavior from females to males. Collectively, we infer that Gambusia holbrooki retains potential for sex reversal at adulthood, similar to what occurs naturally in hermaphroditic fish and can provide an ideal system to investigate these mechanisms in a regulated fashion.

Abstract

17α-Methyltestosterone (MT) is a synthetic steroid that has been widely used to masculinize many fish species when administered early during larval development, however, reports on its efficacy on adults is limited. To this end, this study investigated the efficacy of MT in the masculinization of the eastern mosquitofish (G. holbrooki) at two adult stages (maiden and repeat gravid females). The treated females were fed control or respective MT incorporated feed (0–200 mg/kg diet) for 50 days. Effects of the hormone on secondary sexual characteristics, internal gonad morphology, expression of the Anti-Müllerian Hormone (amh) gene and sexual behavior of the treated females were investigated. The results showed that MT at the dose of 50 mg/kg feed stimulated secondary sexual character development, upregulated expression of amh, formation of testicular tissue and a shift in the behavior similar to those of normal males, prominently so in treated maiden gravid females. Post-treatment, long-term observations indicated that only two masculinized females reverted back to being females and gave birth to young. Induction of masculinizing effects in most individuals suggests that the sexual phenotype of this species appears to be highly plastic with potential to sex reverse at adulthood. This in combination with its small size and short reproductive cycle could provide an ideal system to explore the mechanisms of sequential hermaphroditism in fish and contribute to genetic control of this pest fish.

17β-Trenbolone binds to androgen receptor, decreases number of primordial germ cells, modulates expression of genes related to sexual differentiation, and affects sexual differentiation in zebrafish (Danio rerio)

Exposure to 17β-trenbolone caused a skewed sex ratio in fish. However, the molecular initiating event and key molecular event(s) remain unknown. In this study, zebrafish were exposed to 17β-trenbolone at nominal concentrations of 2 ng/L, 20 ng/L, 200 ng/L, and 2000 ng/L from fertilization to 60 days post fertilization (dpf). First, the sex ratio at 60 dpf was calculated to evaluate adverse outcomes on sexual differentiation. 17β-Trenbolone caused a skewed sex ratio toward males, with intersex individuals observed in the 20 ng/L group and all-male populations found in the 200 ng/L and 2000 ng/L groups. Then, the distribution and number of primordial germ cells, the expression of sex differentiation-related genes, and plasma vitellogenin concentrations were detected in wild-type zebrafish and the EGFP-nanos-3'UTR transgenic line using whole-mount in situ hybridization, real-time PCR, EGFP fluorescence quantification, and enzyme-linked immunosorbent assay. The results indicated that 17β-trenbolone exposure decreased the number of primordial germ cells at 1 dpf and 3 dpf, decreased expression of ovarian differentiation-related genes foxl2 and cyp19a1a at 60 dpf, increased expression of testis differentiation-related genes dmrt1, sox9a, and amh at 60 dpf, and decreased plasma vitellogenin levels at 60 dpf, revealing the key molecular events at different time points involved in affected sexual differentiation by 17β-trenbolone. Finally, molecular docking showed that 17β-trenbolone docked into ligand-binding domain of zebrafish androgen receptor with high binding energy (-3.72 kcal/mol), suggesting that binding to androgen receptor is the molecular initiating event affecting sexual differentiation by 17β-trenbolone. We found that 17β-trenbolone can bind to the zebrafish androgen receptor, decrease the number of primordial germ cells during the early embryonic stage, modulate the expression of genes related to sexual differentiation during gonadal differentiation, and eventually cause a skewed sex ratio toward males in adults.

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.

Genome-wide identification and transcriptome-based expression profiling of the Sox gene family in the spinyhead croaker (Collichthys lucidus)

Sox genes encode transcription factors with a high-mobility group (HMG) box, playing critical roles in the initiation and maintenance of a variety of developmental processes, such as sex determination and differentiation. In the present study, we identified 26 Sox genes in the genome of spinyhead croaker Collichthys lucidus (Richardson, 1844) with homology-based analysis of the HMG box. The transcriptome-based expression profiles revealed that the expression of the Sox gene in gonads began to differ between sexes when the body length was 2.74 ± 0.24 cm. At that time, three Sox genes (Sox11b, Sox8a and Sox19) were significantly upregulated, accompanied by the downregulation of 12 Sox genes in the ovary, and six Sox genes were temporarily significantly upregulated in the testis. Afterwards, the expression profile of Sox genes changed only with a small amplitude in both the ovary and testis. For adult tissues, huge differences were observed in the expression profiles of Sox genes between ovaries and testes, as well as small differences in somatic tissues between sexes. These results provide clues to further decipher the role of Sox genes in the processes of sex determination and differentiation in spinyhead croaker and other teleosts.

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.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

Establishment of transgenic zebrafish (Danio rerio) models expressing fluorescence proteins in the oocytes and somatic supporting cells

The interaction between gonadal somatic support cells and germ cells plays a crucial role in gonadal development. In fish, the process involves various local growth factors such as growth differentiation factor 9 (Gdf9) and gonadal soma-derived factor (Gsdf), which are both members of the transforming growth factor-β (TGF-β) superfamily. Gdf9, an oocyte-secreted factor, is a potent regulator of folliculogenesis in both mammals and fish. By contrast, Gsdf is expressed by the gonadal somatic cells (i.e., Sertoli cells in the testis and granulosa cells in the ovary) that support germ cell development. In this study, we established two transgenic zebrafish models, and demonstrated that the 2.7-kb proximal promoter region of gdf9 drove mCherry expression specifically in the oocytes, whereas the 2.1-kb proximal promoter region of gsdf drove enhanced green fluorescent protein (eGFP) expression in the Sertoli cells and granulosa cells. These proximal promoters contained sufficient information to respectively mimic the spatiotemporal expression patterns of endogenous gdf9 and gsdf in zebrafish. In the Tg(gdf9:mCherry) fish, mCherry was weakly expressed in the oocytes at primary growth stage but strongly expressed in those entering the secondary growth phase. In the Tg(gsdf:eGFP) fish, eGFP-positive Sertoli cells were distributed around spermatogenic cysts in the testis, whereas eGFP-positive granulosa cells were located at the outer side of the follicle layer in the ovary. The eGFP-positive Sertoli cells and granulosa cells seemed to have originated from the dorsal epithelium of the gonads. These Tg(gdf9:mCherry) and Tg(gsdf:eGFP) zebrafish models are suitable for studying gonadal development and function especially on the interaction between germ cells and supporting somatic cells.

Involvement of Transforming Growth Factor Beta Family Genes in Gonadal Differentiation in Japanese Eel, Anguilla japonica, According to Sex-Related Gene Expressions

The gonochoristic feature with environmental sex determination that occurs during the yellow stage in the eel provides an interesting model to investigate the mechanisms of gonadal development. We previously studied various sex-related genes during gonadal sex differentiation in Japanese eels. In the present study, the members of transforming growth factor beta (TGF-β) superfamily were investigated. Transcript levels of anti-Müllerian hormone, its receptor, gonadal soma-derived factor (amh, amhr2, and gsdf, respectively) measured by real-time polymerase chain reaction (qPCR) showed a strong sexual dimorphism. Transcripts were dominantly expressed in the testis, and their levels significantly increased with testicular differentiation. In contrast, the expressions of amh, amhr2, and gsdf transcripts were low in the ovary of E2-feminized female eels. In situ hybridization detected gsdf (but not amh) transcript signals in undifferentiated gonads. amh and gsdf signals were localized to Sertoli cells and had increased significantly with testicular differentiation. Weak gsdf and no amh signals were detected in early ovaries of E2-feminized female eels. Transcript levels of amh and gsdf (not amhr2) decreased during human chorionic gonadotropin (HCG)-induced spermatogenesis in males. This study suggests that amh, amhr2, and especially gsdf might be involved in the gene pathway regulating testicular differentiation of Japanese eels.

Meiotic Chromosome Dynamics in Zebrafish

Recent studies in zebrafish have revealed key features of meiotic chromosome dynamics, including clustering of telomeres in the bouquet configuration, biogenesis of chromosome axis structures, and the assembly and disassembly of the synaptonemal complex that aligns homologs end-to-end. The telomere bouquet stage is especially pronounced in zebrafish meiosis and sub-telomeric regions play key roles in mediating pairing and homologous recombination. In this review, we discuss the temporal progression of these events in meiosis prophase I and highlight the roles of proteins associated with meiotic chromosome architecture in homologous recombination. Finally, we discuss the interplay between meiotic mutants and gonadal sex differentiation and future research directions to study meiosis in living cells, including cell culture.

Heterozygous loss-of-function variants significantly expand the phenotypes associated with loss of GDF11

PURPOSE

Growth differentiation factor 11 (GDF11) is a key signaling protein required for proper development of many organ systems. Only one prior study has associated an inherited GDF11 variant with a dominant human disease in a family with variable craniofacial and vertebral abnormalities. Here, we expand the phenotypic spectrum associated with GDF11 variants and document the nature of the variants.

METHODS

We present a cohort of six probands with de novo and inherited nonsense/frameshift (4/6 patients) and missense (2/6) variants in GDF11. We generated gdf11 mutant zebrafish to model loss of gdf11 phenotypes and used an overexpression screen in Drosophila to test variant functionality.

RESULTS

Patients with variants in GDF11 presented with craniofacial (5/6), vertebral (5/6), neurological (6/6), visual (4/6), cardiac (3/6), auditory (3/6), and connective tissue abnormalities (3/6). gdf11 mutant zebrafish show craniofacial abnormalities and body segmentation defects that match some patient phenotypes. Expression of the patients' variants in the fly showed that one nonsense variant in GDF11 is a severe loss-of-function (LOF) allele whereas the missense variants in our cohort are partial LOF variants.

CONCLUSION

GDF11 is needed for human development, particularly neuronal development, and LOF GDF11 alleles can affect the development of numerous organs and tissues.

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.

Surrogate production of genome-edited sperm from a different subfamily by spermatogonial stem cell transplantation

The surrogate reproduction technique, such as inter-specific spermatogonial stem cells (SSCs) transplantation (SSCT), provides a powerful tool for production of gametes derived from endangered species or those with desirable traits. However, generation of genome-edited gametes from a different species or production of gametes from a phylogenetically distant species such as from a different subfamily, by SSCT, has not succeeded. Here, using two small cyprinid fishes from different subfamilies, Chinese rare minnow (gobiocypris rarus, for brief: Gr) and zebrafish (danio rerio), we successfully obtained Gr-derived genome-edited sperm in zebrafish by an optimized SSCT procedure. The transplanted Gr SSCs supported the host gonadal development and underwent normal spermatogenesis, resulting in a reconstructed fertile testis containing Gr spermatids and zebrafish testicular somatic cells. Interestingly, the surrogate spermatozoa resembled those of host zebrafish but not donor Gr in morphology and swimming behavior. When pou5f3 and chd knockout Gr SSCs were transplanted, Gr-derived genome-edited sperm was successfully produced in zebrafish. This is the first report demonstrating surrogate production of gametes from a different subfamily by SSCT, and surrogate production of genome-edited gametes from another species as well. This method is feasible to be applied to future breeding of commercial fish and livestock.

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.

The emerging role of microRNAs in fish ovary: A mini review

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily at the post-transcriptional levels. It is now well established that miRNAs are crucial regulators of many developmental and physiological processes, including reproduction. In teleosts, expression profiling studies have shown that miRNAs are expressed in the fish ovary and their levels are regulated during follicle development and by hormones. Using CRISPR/Cas9 mediated gene knockout strategies, several recent studies have provided strong evidence that miR-202 and miR-200s on chromosome 23 play critical roles in regulating ovarian development, oogenesis, and ovulation. In this mini review, we provide a brief overview of canonical miRNA biogenesis and functions; summarize miRNAs that are expressed in fish ovary; and discuss the emerging role of miRNAs in regulating fish ovarian functions.

High temperature induced masculinization of zebrafish by down-regulation of sox9b and esr1 via DNA methylation

Elevated temperature could influence the sex differentiation by altering the expression of sex-related genes in fish. However, the underlying mechanisms by which the gene expression is altered remain poorly understood. Here, we aimed to explore the role of DNA methylation in sex differentiation of zebrafish (Danio rerio) in response to elevated temperature. The results showed that high temperature (33°C) exposure of fish from 20 to 30 days post fertilization (dpf), compared to normal temperature (28°C), resulted in male-biased sex ratio and decreased expression of female-related genes including cyp19a1a, sox9b and esr1. Meanwhile, the expressions of DNA methyltransferases dnmt3a1 and dnmt3a2, and the DNA methylation levels in sox9b and esr1 promoter were significantly increased by high temperature, strongly implying that DNA methylation is involved in high temperature-induced masculinization of zebrafish. Co-treatment with 5-aza-2'-deoxycytidine (a DNA methylation inhibitor) attenuated the high temperature-induced masculinizing effect, recovered the expression of esr1 and sox9b, suppressed the transcription of dnmt3a1 and dnmt3a2, and decreased the methylation of esr1 and sox9b promoter, further confirming that DNA methylation plays an important role in high temperature-induced masculinization of zebrafish. Furthermore, the methylation of sox9b promoter decreased the enrichment of transcription factor CREB (cAMP-responsive element binding proteins). Overall, these findings suggest that high temperature induce masculinization of zebrafish by down-regulation of female-related genes via DNA methylation, providing a new insight in understanding the epigenetic mechanism of thermal-mediated sex differentiation in fish.

Immunohistochemical expression of aromatase cyp19a1a and cyp19a1b in the ovary and brain of zebrafish (Danio rerio) exposed to different concentrations of bisphenol A

Bisphenol A (BPA) is used to produce plastic and plastic derived products in multitude of daily utensils, being one of the industrial compounds most widely used. This endocrine disrupting chemical (EDCs) is a well-known environmental pollutant released into the aquatic environment from industrial wastewater, sewage sludge or landfill leachate. Aromatases are considered potential targets of EDCs with characteristics that make them suitable biomarkers of exposure to their effects. The main objective of our study was to evaluate the expression of cyp19a aromatase as a toxicological endpoint after BPA exposure through the identification and assessment of alterations of the main cells responsible for cyp19a1a and cyp19a1b expression in the zebrafish ovary and brain using different concentrations of BPA in water. Immunohistochemistry was used to analyze the expression of these enzymes in female zebrafish exposed and not exposed to different concentrations of BPA (1, 10, 100 and 1000 μg / L) in water (n = 6/group) for 14 days. The results obtained in this study showed that the cyp19a aromatase system, involved in the synthesis of steroid compounds, is specially located in distinct oocyte stages in the ovary (cyp19a1a) and in radial glial cells of the brain (cyp19a1b). An overexpression of these aromatases was observed after BPA exposure in zebrafish, peaking from a concentration of 10 µg/L and showing to be good biomarkers of exposure to identify the early effects of low BPA concentrations. To our knowledge, this study is the first to localize and quantify the expression of cyp19a1a and cyp19a1b in the cells of brain and ovary after fish exposure to different BPA concentrations in water.

A duplicated amh is the master sex-determining gene for Sebastes rockfish in the Northwest Pacific

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

In Silico Prediction of Transcription Factor Collaborations Underlying Phenotypic Sexual Dimorphism in Zebrafish (Danio rerio)

The transcriptional regulation of gene expression in higher organisms is essential for different cellular and biological processes. These processes are controlled by transcription factors and their combinatorial interplay, which are crucial for complex genetic programs and transcriptional machinery. The regulation of sex-biased gene expression plays a major role in phenotypic sexual dimorphism in many species, causing dimorphic gene expression patterns between two different sexes. The role of transcription factor (TF) in gene regulatory mechanisms so far has not been studied for sex determination and sex-associated colour patterning in zebrafish with respect to phenotypic sexual dimorphism. To address this open biological issue, we applied bioinformatics approaches for identifying the predicted TF pairs based on their binding sites for sex and colour genes in zebrafish. In this study, we identified 25 (e.g., STAT6-GATA4; JUN-GATA4; SOX9-JUN) and 14 (e.g., IRF-STAT6; SOX9-JUN; STAT6-GATA4) potentially cooperating TFs based on their binding patterns in promoter regions for sex determination and colour pattern genes in zebrafish, respectively. The comparison between identified TFs for sex and colour genes revealed several predicted TF pairs (e.g., STAT6-GATA4; JUN-SOX9) are common for both phenotypes, which may play a pivotal role in phenotypic sexual dimorphism in zebrafish.

Sycp1 Is Not Required for Subtelomeric DNA Double-Strand Breaks but Is Required for Homologous Alignment in Zebrafish Spermatocytes

In meiotic prophase I, homologous chromosomes are bound together by the synaptonemal complex, in which two axial elements are connected by transverse filaments and central element proteins. In human and zebrafish spermatocytes, homologous recombination and assembly of the synaptonemal complex initiate predominantly near telomeres. In mice, synapsis is not required for meiotic double-strand breaks (DSBs) and homolog alignment but is required for DSB repair; however, the interplay of these meiotic events in the context of peritelomeric bias remains unclear. In this study, we identified a premature stop mutation in the zebrafish gene encoding the transverse filament protein Sycp1. In sycp1 mutant zebrafish spermatocytes, axial elements were formed and paired at chromosome ends between homologs during early to mid-zygonema. However, they did not synapse, and their associations were mostly lost in late zygotene- or pachytene-like stages. In sycp1 mutant spermatocytes, γH2AX signals were observed, and Dmc1/Rad51 and RPA signals appeared predominantly near telomeres, resembling wild-type phenotypes. We observed persistent localization of Hormad1 along the axis in sycp1 mutant spermatocytes, while the majority of Iho1 signals appeared and disappeared with kinetics similar to those in wild-type spermatocytes. Notably, persistent Iho1 foci were observed in spo11 mutant spermatocytes, suggesting that Iho1 dissociation from axes occurs in a DSB-dependent manner. Our results demonstrated that Sycp1 is not required for peritelomeric DSB formation but is necessary for complete pairing of homologs in zebrafish meiosis.

Diverse Regulation but Conserved Function: SOX9 in Vertebrate Sex Determination

Sex determination occurs early during embryogenesis among vertebrates. It involves the differentiation of the bipotential gonad to ovaries or testes by a fascinating diversity of molecular switches. In most mammals, the switch is SRY (sex determining region Y); in other vertebrates it could be one of a variety of genes including Dmrt1 or dmy. Downstream of the switch gene, SOX9 upregulation is a central event in testes development, controlled by gonad-specific enhancers across the 2 Mb SOX9 locus. SOX9 is a ‘hub’ gene of gonadal development, regulated positively in males and negatively in females. Despite this diversity, SOX9 protein sequence and function among vertebrates remains highly conserved. This article explores the cellular, morphological, and genetic mechanisms initiated by SOX9 for male gonad differentiation.

The SARS-CoV-2 receptor and other key components of the Renin-Angiotensin-Aldosterone System related to COVID-19 are expressed in enterocytes in larval zebrafish

People with underlying conditions, including hypertension, obesity, and diabetes, are especially susceptible to negative outcomes after infection with coronavirus SARS-CoV-2, which causes COVID-19. Hypertension and respiratory inflammation are exacerbated by the Renin-Angiotensin-Aldosterone System (RAAS), which normally protects from rapidly dropping blood pressure via Angiotensin II (Ang II) produced by the enzyme Ace. The Ace paralog Ace2 degrades Ang II, counteracting its chronic effects, and serves as the SARS-CoV-2 receptor. Ace, the coronavirus, and COVID-19 comorbidities all regulate Ace2, but we do not yet understand how. To exploit zebrafish (Danio rerio) to help understand the relationship of the RAAS to COVID-19, we must identify zebrafish orthologs and co-orthologs of human RAAS genes and understand their expression patterns. To achieve these goals, we conducted genomic and phylogenetic analyses and investigated single cell transcriptomes. Results showed that most human RAAS genes have one or more zebrafish orthologs or co-orthologs. Results identified a specific type of enterocyte as the specific site of expression of zebrafish orthologs of key RAAS components, including Ace, Ace2, Slc6a19 (SARS-CoV-2 co-receptor), and the Angiotensin-related peptide cleaving enzymes Anpep (receptor for the common cold coronavirus HCoV-229E), and Dpp4 (receptor for the Middle East Respiratory Syndrome virus, MERS-CoV). Results identified specific vascular cell subtypes expressing Ang II receptors, apelin, and apelin receptor genes. These results identify genes and cell types to exploit zebrafish as a disease model for understanding mechanisms of COVID-19.

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.

Application of quantitative transcriptomics in evaluating the ex vivo effects of per- and polyfluoroalkyl substances on Atlantic cod (Gadus morhua) ovarian physiology

Because of their global consumption and persistence, per- and polyfluoroalkyl substances (PFASs), are ubiquitously distributed in the environment, as well as in wildlife and humans. In the present study, we have employed an ex vivo organ culture technique, based on the floating agarose method, of Atlantic cod ovarian tissue to investigate the effects of three different concentrations of PFOS, PFOA (1, 5 and 25 μM) and PFNA (0.5, 5 and 50 μM), used singly and in also in combination (1×, 20× and 100×). In the 1× exposure mixture, concentrations were decided based on their proportional levels (in molar equivalents) relative to PFOS, which is the most abundant PFAS in cod liver from a 2013 screening project. To investigate the detailed underlying mechanisms and biological processes, transcriptome sequencing was performed on exposed ovarian tissue. The number of differentially expressed genes (DEGs) having at least 0.75 log₂-fold change was elevated in high, compared to low and medium concentration exposures. The highest PFNA, PFOA and PFOS concentrations, and the highest (100×) mixture exposure, showed 40, 68, 1295, and 802 DEGs, respectively. The latter two exposure groups shared a maximum of 438 DEGs. In addition, they both shared the majority of functionally enriched pathways belonging to biological processes such as cellular signaling, cell adhesion, lipid metabolism, immunological responses, cancer, reproduction and metabolism. Shortlisted DEGs that were specifically annotated to reproduction associated gene ontology (GO) terms were observed only in the highest PFOS and mixture exposure groups. These transcripts contributed to ovarian key events such as steroidogenesis (star, cyp19a1a), oocyte growth (amh), maturation (igfbp5b, tgfβ2, tgfβ3), and ovulation (pgr, mmp2). Contrary to other PFAS congeners, the highest PFOS concentration showed almost similar transcript expression patterns compared to the highest mixture exposure group. This indicates that PFOS is the active component of the mixture that significantly altered the normal functioning of female gonads, and possibly leading to serious reproductive consequences in teleosts.

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.

Rbm46, a novel germ cell-specific factor, modulates meiotic progression and spermatogenesis

It has been suggested that many novel RNA-binding proteins (RBPs) are required for gametogenesis, but the necessity of few of these proteins has been functionally verified. Here, we identified one RBP, Rbm46, and investigated its expression pattern and role in zebrafish reproduction. We found that rbm46 is maternally provided and specifically expressed in the germ cells of gonadal tissues using in situ hybridization, reverse transcription-PCR, and quantitative real-time polymerase chain reaction (qRT-PCR). Two independent rbm46 mutant zebrafish lines were generated via the transcription activator-like effector nuclease technique. Specific disruption of rbm46 resulted in masculinization and infertility in the mutants. Although the spermatogonia appeared grossly normal in the mutants, spermatogenesis was impaired, and meiosis events were not observed. The introduction of a tp53M214K mutation could not rescue the female-to-male sex-reversal phenotype, indicating that rbm46 acts independently of the p53-dependent apoptotic pathway. RNA sequencing and qRT-PCR subsequently indicated that Rbm46 might be involved in the posttranscriptional regulation of functional genes essential for germ cell development, such as nanos3, dazl, and sycp3, during gametogenesis. Together, our results reveal for the first time the crucial role of rbm46 in regulating germ cell development in vivo through promotion of germ cell progression through meiosis prophase I.

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.

Igf3 is essential for ovary differentiation in zebrafish†

Zebrafish gonadal sexual differentiation is an important but poorly understood subject. Previously, we have identified a novel insulin-like growth factor (Igf) named insulin-like growth factor 3 (Igf3) in teleosts. The importance of Igf3 in oocyte maturation and ovulation has been recently demonstrated by us in zebrafish. In this study, we have further found the essential role of Igf3 in gonadal sexual differentiation of zebrafish. A differential expression pattern of igf3 between ovary and testis during sex differentiation (higher level in ovary than in testis) was found in zebrafish. An igf3 knockout zebrafish line was established using TALENs-mediated gene knockout technique. Intriguingly, all igf3 homozygous mutants were males due to the female-to-male sex reversal occurred during sex differentiation. Further analysis showed that Igf3 did not seem to affect the formation of so-called juvenile ovary and oocyte-like germ cells. Oocyte development was arrested at primary growth stage, and the ovary was gradually sex-reversed to testis before 60 day post fertilization (dpf). Such sex reversal was likely due to decreased germ cell proliferation by suppressing PI3K/Akt pathway in early ovaries of igf3 mutants. Estrogen is considered as a master regulator in fish sex differentiation. Here, we found that igf3 expression could be upregulated by estrogen in early stages of ovarian follicles as evidenced in in vitro treatment assays and cyp19a1a mutant zebrafish, and E2 failed to rescue the defects of igf3 mutants in ovarian development, suggesting that Igf3 may serve as a downstream factor of estrogen signaling in sex differentiation. Taken together, we demonstrated that Igf3 is essential for ovary differentiation in zebrafish.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

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

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