Disruption of dmrt1 rescues the all-male phenotype of the cyp19a1a mutant in zebrafish - a novel insight into the roles of aromatase/estrogens in gonadal differentiation and early folliculogenesis

Disruption of Thyroid Hormone Receptor Thrab Leads to Female Infertility in Zebrafish

Thyroid hormones (THs) T4 and T3 are vital for development, growth, and metabolism. Thyroid dysfunction can also cause problems in fertility, suggesting involvement of THs in reproduction. In zebrafish, there exist 2 forms of TH receptor alpha gene (thraa and thrab). Disruption of these genes by CRISPR/Cas9 showed no reproductive irregularities in the thraa mutant; however, inactivation of the thrab gene resulted in female infertility. Although young female mutants (thrabm/m) showed normal ovarian development and folliculogenesis before sexual maturation, they failed to release eggs during oviposition after sexual maturation. This spawning failure was due to oviductal blockage at the genital papilla. The obstruction of the oviduct subsequently caused an accumulation of the eggs in the ovary, resulting in severe ovarian hypertrophy, abdominal distention, and disruption of folliculogenesis. Gene expression analysis showed expression of both TH receptors and estrogen receptors in the genital papilla, suggesting a direct TH action and potential interactions between thyroid and estrogen signaling pathways in controlling genital papilla development and function. In addition to their actions in the reproductive tracts, THs may also have direct effects in the ovary, as suggested by follicle atresia and cessation of folliculogenesis in the heterozygous mutant (thrab+/m), which was normal in all aspects of female reproduction in young and sexually mature fish but exhibited premature ovarian failure in aged females. In summary, this study provides substantial evidence for roles of THs in controlling the development and functions of both reproductive tract and ovary.

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.

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

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

Tilapia, a good model for studying reproductive endocrinology

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

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.

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

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

Waterborne Tebuconazole Exposure Induces Male-Biased Sex Differentiation in Zebrafish (Danio rerio) Larvae via Aromatase Inhibition

Tebuconazole is a widely used fungicide for various crops that targets sterol 14-α-demethylase (CYP51) in fungi. However, attention has shifted to aromatase (CYP19) due to limited research indicating its reproductive impact on aquatic organisms. Herein, zebrafish were exposed to 0.5 mg/L tebuconazole at different developmental stages. The proportion of males increased significantly after long-term exposure during the sex differentiation phase (0-60, 5-60, and 19-60 days postfertilization (dpf)). Testosterone levels increased and 17β-estradiol and cyp19a1a expression levels decreased during the 5-60 dpf exposure, while the sex ratio was equally distributed on coexposure with 50 ng/L 17β-estradiol. Chemically activated luciferase gene expression bioassays determined that the male-biased sex differentiation was not caused by tebuconazole directly binding to sex hormone receptors. Protein expression and phosphorylation levels were specifically altered in the vascular endothelial growth factor signaling pathway despite excluding the possibility of tebuconazole directly interacting with kinases. Aromatase was selected for potential target analysis. Molecular docking and aromatase activity assays demonstrated the interactions between tebuconazole and aromatase, highlighting that tebuconazole poses a threat to fish populations by inducing a gender imbalance.

Rescue of bmp15 deficiency in zebrafish by mutation of inha reveals mechanisms of BMP15 regulation of folliculogenesis

As an oocyte-specific growth factor, bone morphogenetic protein 15 (BMP15) plays a critical role in controlling folliculogenesis. However, the mechanism of BMP15 action remains elusive. Using zebrafish as the model, we created a bmp15 mutant using CRISPR/Cas9 and demonstrated that bmp15 deficiency caused a significant delay in follicle activation and puberty onset followed by a complete arrest of follicle development at previtellogenic (PV) stage without yolk accumulation. The mutant females eventually underwent female-to-male sex reversal to become functional males, which was accompanied by a series of changes in secondary sexual characteristics. Interestingly, the blockade of folliculogenesis and sex reversal in bmp15 mutant could be partially rescued by the loss of inhibin (inha-/-). The follicles of double mutant (bmp15-/-;inha-/-) could progress to mid-vitellogenic (MV) stage with yolk accumulation and the fish maintained their femaleness without sex reversal. Transcriptome analysis revealed up-regulation of pathways related to TGF-β signaling and endocytosis in the double mutant follicles. Interestingly, the expression of inhibin/activin βAa subunit (inhbaa) increased significantly in the double mutant ovary. Further knockout of inhbaa in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-) resulted in the loss of yolk granules again. The serum levels of estradiol (E2) and vitellogenin (Vtg) both decreased significantly in bmp15 single mutant females (bmp15-/-), returned to normal in the double mutant (bmp15-/-;inha-/-), but reduced again significantly in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-). E2 treatment could rescue the arrested follicles in bmp15-/-, and fadrozole (a nonsteroidal aromatase inhibitor) treatment blocked yolk accumulation in bmp15-/-;inha-/- fish. The loss of inhbaa also caused a reduction of Vtg receptor-like molecules (e.g., lrp1ab and lrp2a). In summary, the present study provided comprehensive genetic evidence that Bmp15 acts together with the activin-inhibin system in the follicle to control E2 production from the follicle, Vtg biosynthesis in the liver and its uptake by the developing oocytes.

Transgenerational endocrine disruptor effects of cadmium in zebrafish and contribution of standing epigenetic variation to adaptation

Evidence has emerged that environmentally-induced epigenetic changes can have long-lasting effects on gene transcription across generations. These recent findings highlight the need to investigate the transgenerational impacts of pollutants to assess their long term effects on populations. In this study, we investigated the transgenerational effect of cadmium on zebrafish across 4 generations. A first whole methylome approach carried out on fish of the first two generations led us to focus our investigations on the estradiol receptor alpha gene (esr1). We observed a sex-dependent transgenerational inheritance of Cd-induced DNA methylation changes up to the last generation. These changes were associated with single nucleotide polymorphisms (SNPs) that were themselves at the origin of the creation or deletion of methylation sites. Thus, Cd-induced genetic selection gave rise to DNA methylation changes. We also analyzed the transcription level of various sections of esr1 as well as estrogen responsive genes. While Cd triggered transgenerational disorders, Cd-induced epigenetic changes in esr1 contributed to the rapid transgenerational adaptation of fish to Cd. Our results provide insight into the processes underpinning rapid adaptation and highlight the need to maintain genetic diversity within natural populations to bolster the resilience of species faced with the global environmental changes.

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

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

Loss of growth differentiation factor 9 causes an arrest of early folliculogenesis in zebrafish-A novel insight into its action mechanism

Growth differentiation factor 9 (GDF9) was the first oocyte-specific growth factor identified; however, most information about GDF9 functions comes from studies in the mouse model. In this study, we created a mutant for Gdf9 gene (gdf9-/-) in zebrafish using TALEN approach. The loss of Gdf9 caused a complete arrest of follicle development at primary growth (PG) stage. These follicles eventually degenerated, and all mutant females gradually changed to males through sex reversal, which could be prevented by mutation of the male-promoting gene dmrt1. Interestingly, the phenotypes of gdf9-/- could be rescued by simultaneous mutation of inhibin α (inha-/-) but not estradiol treatment, suggesting a potential role for the activin-inhibin system or its signaling pathway in Gdf9 actions. In gdf9-null follicles, the expression of activin βAa (inhbaa), but not βAb (inhbab) and βB (inhbb), decreased dramatically; however, its expression rebounded in the double mutant (gdf9-/-;inha-/-). These results indicate clearly that the activation of PG follicles to enter the secondary growth (SG) requires intrinsic factors from the oocyte, such as Gdf9, which in turn works on the neighboring follicle cells to trigger follicle activation, probably involving activins. In addition, our data also support the view that estrogens are not involved in follicle activation as recently reported.

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.

First account of a transient intersex in spotted scat, Scatophagus argus: a marine gonochoristic fish

This study presents the first incidence of intersex associated with testis-ova in spotted scat (Scatophagus argus) reared in a controlled environment. The testis-ova is associated with the abnormal occurrence of primary oocytes (POs) in some male testis and is referred to as ectopic primary oocytes (Ecto-PO), whiles individuals with Ecto-PO are called "Ecto-PO gonad/individuals." We investigated gonads of 129 male spotted scat aged 4-12 and 18 months after hatch (mah) by histological studies for the presence of female sexual characteristics. A total of 20 out of 88 gonads representing 22.7% of male fish aged 6-12, or 15.5% of all male fish sampled, were found to have visible Ecto-PO. At least, the Ecto-PO had an average of 7 oocytes per gonadal section, indicating high severity. The Ecto-PO appears after sex differentiation and degenerates during sexual maturation. The Ecto-PO did not significantly influence the expression pattern of male and female sex-related genes performed using qPCR. Immunofluorescence of 42sp50 specifically stained the Ecto-PO without influence from the surrounding testicular tissues. In addition, temperature did not correlate with the proliferation of the Ecto-PO, but rather gonad developmental strategy. The results show that the naturally occurring Ecto-PO might not be detrimental to testis development and could be considered a frequent-high-level incidence of natural aberration. This study highlights the intricacy of fish sex differentiation and provides a new research chapter to ascertain the mystery behind the occurrence of Ecto-PO.

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

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

Sex-specific differences in zebrafish brains

In this systematic review, we highlight the differences between the male and female zebrafish brains to understand their differentiation and their use in studying sex-specific neurological diseases. Male and female brains display subtle differences at the cellular level which may be important in driving sex-specific signaling. Sex differences in the brain have been observed in humans as well as in non-human species. However, the molecular mechanisms of brain sex differentiation remain unclear. The classical model of brain sex differentiation suggests that the steroid hormones derived from the gonads are the primary determinants in establishing male and female neural networks. Recent studies indicate that the developing brain shows sex-specific differences in gene expression prior to gonadal hormone action. Hence, genetic differences may also be responsible for differentiating the brain into male and female types. Understanding the signaling mechanisms involved in brain sex differentiation could help further elucidate the sex-specific incidences of certain neurological diseases. The zebrafish model could be appropriate for enhancing our understanding of brain sex differentiation and the signaling involved in neurological diseases. Zebrafish brains show sex-specific differences at the hormonal level, and recent advances in RNA sequencing have highlighted critical sex-specific differences at the transcript level. The differences are also evident at the cellular and metabolite levels, which could be important in organizing sex-specific neuronal signaling. Furthermore, in addition to having one ortholog for 70% of the human gene, zebrafish also shares brain structural similarities with other higher eukaryotes, including mammals. Hence, deciphering brain sex differentiation in zebrafish will help further enhance the diagnostic and pharmacological intervention of neurological diseases.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Identification and functional analysis of foxl2 and nodal in sea cucumber, Apostichopus japonicus

Sea cucumber (Apostichopus japonicus) is an important mariculture species in China. To date, the mechanisms of sex determination and differentiation in sea cucumber remain unclear. Identifying sex-specific molecular markers is an effective method for revealing the genetic basis of sex determination and sex differentiation. In this study, foxl2 and nodal homologous genes were identified in A. japonicus. Foxl2 exhibited dynamic and sexually dimorphic expression patterns in the gonads, with prominent expression in the ovaries and minimal expression in the testis according to real-time quantitative PCR (RT-qPCR) study. As nodal was specifically expressed in the ovary, it could serve as an ovary-specific marker in sea cucumber. Additionally, knockdown of foxl2 or nodal using RNA interference (RNAi) led to the down-regulation of piwi, germ cell-less, and dmrt1, suggesting that foxl2 and nodal may play important roles in gonad maintenance of sea cucumber. Overall, this study adds to our understanding of the roles of foxl2 and nodal in the gonadal development of A. japonicus, which provides further insight into the mechanisms of sea cucumber sex determination and differentiation.

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.

Functional analysis of co-expression networks of zebrafish ace2 reveals enrichment of pathways associated with development and disease

Human Angiotensin I Converting Enzyme 2 (ACE2) plays an essential role in blood pressure regulation and SARS-CoV-2 entry. ACE2 has a highly conserved, one-to-one ortholog (ace2) in zebrafish, which is an important model for human diseases. However, the zebrafish ace2 expression profile has not yet been studied during early development, between genders, across different genotypes, or in disease. Moreover, a network-based meta-analysis for the extraction of functionally enriched pathways associated with differential ace2 expression is lacking in the literature. Herein, we first identified significant development-, tissue-, genotype-, and gender-specific modulations in ace2 expression via meta-analysis of zebrafish Affymetrix transcriptomics datasets (n_datasets = 107); and the correlation analysis of ace2 meta-differential expression profile revealed distinct positively and negatively correlated local functionally enriched gene networks. Moreover, we demonstrated that ace2 expression was significantly modulated under different physiological and pathological conditions related to development, tissue, gender, diet, infection, and inflammation using additional RNA-seq datasets. Our findings implicate a novel translational role for zebrafish ace2 in organ differentiation and pathologies observed in the intestines and liver.

Comparison of differential expression genes in ovaries and testes of Pearlscale angelfish Centropyge vrolikii based on RNA-Seq analysis

Pearlscale angelfish Centropyge vrolikii is a kind of protogynous hermaphrodite fish with a natural sexual reversion. Under appropriate social conditions, a female fish can transform into a male fish spontaneously. It is an important prerequisite for artificial breeding to understand the process of its gonadal development and sexual reversion. Gonadal development is regulated by many sex-related genes. In this study, we used unreferenced RNA-Seq technology to sequence the ovary at the perinucleolus stage (O_II), ovary at the yolk vesicle stage (O_IV),IV and testis (T), respectively; screened the gonadal differential expression genes (DEGs); and analyzed the expression of these genes in different developmental stages of ovary and different sex gonads. The results showed that a total of 142,589 all-unigene samples were assembled, and gene annotation was performed by COG, GO, KEGG, KOG, Pfam, Swissprot, eggNOG, and NR functional database. Comparative analysis revealed that there were 1919 genes that were up-regulated and 1289 genes were down-regulated in comparison to O_IV vs O_II, while there were 3653 genes that were up-regulated and 2874 genes were down-regulated in comparison of O_IV vs T, there were 3345 genes that were up-regulated and 2995 genes were down-regulated in comparison of the O_II vs the T. At the same time, the results verified by RT-qPCR were consistent with the variation trend of transcriptome data. Among the results, amh, sox9b, dmrt1, dmrt2, cyp11a, cyp17a, and cyp19a were significantly expressed in the testes, while sox3, sox4, sox11, sox17, and hsd3b7 were significantly expressed in the ovaries. And, the expression of the amh, sox9b, dmrt2, and dmrt1 were low in the O_II and O_IV, while significantly increased during the ovotestis in the hermaphroditic period (O_T), and finally reached the highest level in pure testis after sex reversal. The expression of sox3, sox4, hsd3b7, sox11, and sox17 was significantly reduced during the hermaphroditic period (O_T). These results suggested that these genes may play an important role in the process of sex reversal. This study is helpful to further understand the molecular regulation mechanism of gonadal development and sexual reversion in Pearlscale angelfish and also provide important clues for future studies.

Sex Determination and Differentiation in Teleost: Roles of Genetics, Environment, and Brain

The fish reproductive system is a complex biological system. Nonetheless, reproductive organ development is conserved, which starts with sex determination and then sex differentiation. The sex of a teleost is determined and differentiated from bipotential primordium by genetics, environmental factors, or both. These two processes are species-specific. There are several prominent genes and environmental factors involved during sex determination and differentiation. At the cellular level, most of the sex-determining genes suppress the female pathway. For environmental factors, there are temperature, density, hypoxia, pH, and social interaction. Once the sexual fate is determined, sex differentiation takes over the gonadal developmental process. Environmental factors involve activation and suppression of various male and female pathways depending on the sexual fate. Alongside these factors, the role of the brain during sex determination and differentiation remains elusive. Nonetheless, GnRH III knockout has promoted a male sex-biased population, which shows brain involvement during sex determination. During sex differentiation, LH and FSH might not affect the gonadal differentiation, but are required for regulating sex differentiation. This review discusses the role of prominent genes, environmental factors, and the brain in sex determination and differentiation across a few teleost species.

The application of genome editing technology in fish

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

A transgenic system for targeted ablation of reproductive and maternal-effect genes

Maternally provided gene products regulate the earliest events of embryonic life, including formation of the oocyte that will develop into an egg, and eventually into an embryo. Forward genetic screens have provided invaluable insights into the molecular regulation of embryonic development, including the essential contributions of some genes whose products must be provided to the transcriptionally silent early embryo for normal embryogenesis, called maternal-effect genes. However, other maternal-effect genes are not accessible due to their essential zygotic functions during embryonic development. Identifying these regulators is essential to fill the large gaps in our understanding of the mechanisms and molecular pathways contributing to fertility and to maternally regulated developmental processes. To identify these maternal factors, it is necessary to bypass the earlier requirement for these genes so that their potential later functions can be investigated. Here, we report reverse genetic systems to identify genes with essential roles in zebrafish reproductive and maternal-effect processes. As proof of principle and to assess the efficiency and robustness of mutagenesis, we used these transgenic systems to disrupt two genes with known maternal-effect functions: kif5ba and bucky ball.

RNase1 can modulate gut microbiota and metabolome after Aeromonas hydrophila infection in blunt snout bream

Pancreatic ribonuclease (RNase1) of Megalobrama amblycephala exhibits both antimicrobial and digestive activity. The gut microbiome improve the digestion and metabolic capacity and enhance the functioning of the immune system of the host against pathogenic bacteria. In this study, we aimed to assess the protective effect of RNase1 on Aeromonas hydrophila-induced inflammation and intestinal microbial metabolism. Megalobrama amblycephala were randomly divided into three groups: control (injected PBS), infection (A. hydrophila-injected), and treatment group (RNase1 pretreatment 24 h before the A. hydrophila injection). The morphological symptoms were significantly alleviated by RNase1. RNase1 reshaped the perturbed gut microbiota by upregulating Proteobacteria and Vibrio richness and downregulating Firmicutes, Chlamydiae, Bacillus, and Gemmobacter richness. The lysophosphatidylcholine, (±) 17 HETE, D- (+) -cellobiose, and PC (20:5) in the treatment group were restored by RNase 1 protein treatment to the level of the control group. In the treatment group, phospholipid metabolism, fatty acid metabolism, glucose metabolism and lipid metabolism were different from the control and infection groups. The proinflammatory factors concentration in intestinal samples significantly increased after A. hydrophila infection. Our results revealed that RNase1 plays an important role in resistance to pathogen invasion, reducing inflammation, and improving intestinal function, thus inhibiting the occurrence of disease.

Functional Divergence of Multiple Duplicated Foxl2 Homeologs and Alleles in a Recurrent Polyploid Fish

Evolutionary fates of duplicated genes have been widely investigated in many polyploid plants and animals, but research is scarce in recurrent polyploids. In this study, we focused on foxl2, a central player in ovary, and elaborated the functional divergence in gibel carp (Carassius gibelio), a recurrent auto-allo-hexaploid fish. First, we identified three divergent foxl2 homeologs (Cgfoxl2a-B, Cgfoxl2b-A, and Cgfoxl2b-B), each of them possessing three highly conserved alleles and revealed their biased retention/loss. Then, their abundant sexual dimorphism and biased expression were uncovered in hypothalamic-pituitary-gonadal axis. Significantly, granulosa cells and three subpopulations of thecal cells were distinguished by cellular localization of CgFoxl2a and CgFoxl2b, and the functional roles and the involved process were traced in folliculogenesis. Finally, we successfully edited multiple foxl2 homeologs and/or alleles by using CRISPR/Cas9. Cgfoxl2a-B deficiency led to ovary development arrest or complete sex reversal, whereas complete disruption of Cgfoxl2b-A and Cgfoxl2b-B resulted in the depletion of germ cells. Taken together, the detailed cellular localization and functional differences indicate that Cgfoxl2a and Cgfoxl2b have subfunctionalized and cooperated to regulate folliculogenesis and gonad differentiation, and Cgfoxl2b has evolved a new function in oogenesis. Therefore, the current study provides a typical case of homeolog/allele diversification, retention/loss, biased expression, and sub-/neofunctionalization in the evolution of duplicated genes driven by polyploidy and subsequent diploidization from the recurrent polyploid fish.

Germline sexual fate is determined by the antagonistic action of dmrt1 and foxl3/foxl2 in tilapia

Germline sexual fate has long been believed to be determined by the somatic environment, but this idea is challenged by recent studies of foxl3 mutants in medaka. Here, we demonstrate that the sexual fate of tilapia germline is determined by the antagonistic interaction of dmrt1 and foxl3, which are transcriptionally repressed in male and female germ cells, respectively. Loss of dmrt1 rescued the germ cell sex reversal in foxl3Δ7/Δ7 XX fish, and loss of foxl3 partially rescued germ cell sex reversal but not somatic cell fate in dmrt1Δ5/Δ5 XY fish. Interestingly, germ cells lost sexual plasticity in dmrt1Δ5/Δ5 XY and foxl3Δ7/Δ7 XX single mutants, as aromatase inhibitor (AI) and estrogen treatment failed to rescue the respective phenotypes. However, recovery of germ cell sexual plasticity was observed in dmrt1/foxl3 double mutants. Importantly, mutation of somatic cell-specific foxl2 resulted in testicular development in foxl3Δ7/Δ7 or dmrt1Δ5/Δ5 mutants. Our findings demonstrate that sexual plasticity of germ cells relies on the presence of both dmrt1 and foxl3. The existence of dmrt1 and foxl3 allows environmental factors to influence the sex fate decision in vertebrates.

A Transient Hermaphroditic Stage in Early Male Gonadal Development in Little Yellow Croaker, Larimichthys polyactis

Animal taxa show remarkable variability in sexual reproduction, where separate sexes, or gonochorism, is thought to have evolved from hermaphroditism for most cases. Hermaphroditism accounts for 5% in animals, and sequential hermaphroditism has been found in teleost. In this study, we characterized a novel form of the transient hermaphroditic stage in little yellow croaker (Larimichthys polyactis) during early gonadal development. The ovary and testis were indistinguishable from 7 to 40 days post-hatching (dph). Morphological and histological examinations revealed an intersex stage of male gonads between 43 and 80 dph, which consist of germ cells, somatic cells, efferent duct, and early primary oocytes (EPOs). These EPOs in testis degenerate completely by 90 dph through apoptosis yet can be rescued by exogenous 17-β-estradiol. Male germ cells enter the mitotic flourishing stage before meiosis is initiated at 180 dph, and they undergo normal spermatogenesis to produce functional sperms. This transient hermaphroditic stage is male-specific, and the ovary development appears to be normal in females. This developmental pattern is not found in the sister species Larimichthys crocea or any other closely related species. Further examinations of serum hormone levels indicate that the absence of 11-ketotestosterone and elevated levels of 17-β-estradiol delineate the male intersex gonad stage, providing mechanistic insights on this unique phenomenon. Our research is the first report on male-specific transient hermaphroditism and will advance the current understanding of fish reproductive biology. This unique gonadal development pattern can serve as a useful model for studying the evolutionary relationship between hermaphroditism and gonochorism, as well as teleost sex determination and differentiation strategies.

Transcriptomic and Epigenetic Preservation of Genetic Sex Identity in Estrogen-feminized Male Chicken Embryonic Gonads

Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within 1 year. To identify the feminization-resistant "memory" of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (eg, DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex.

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.

Loss of dmrt1 restores zebrafish female fates in the absence of cyp19a1a but not rbpms2a/b

Sex determination and differentiation is a complex process regulated by multiple factors, including factors from the germline or surrounding somatic tissue. In zebrafish, sex-determination involves establishment of a bipotential ovary that undergoes sex-specific differentiation and maintenance to form the functional adult gonad. However, the relationships among these factors are not fully understood. Here, we identify potential Rbpms2 targets and apply genetic epistasis experiments to decipher the genetic hierarchy of regulators of sex-specific differentiation. We provide genetic evidence that the crucial female factor rbpms2 is epistatic to the male factor dmrt1 in terms of adult sex. Moreover, the role of Rbpms2 in promoting female fates extends beyond repression of Dmrt1, as Rbpms2 is essential for female differentiation even in the absence of Dmrt1. In contrast, female fates can be restored in mutants lacking both cyp19a1a and dmrt1, and prolonged in bmp15 mutants in the absence of dmrt1. Taken together, this work indicates that cyp19a1a-mediated suppression of dmrt1 establishes a bipotential ovary and initiates female fate acquisition. Then, after female fate specification, Cyp19a1a regulates subsequent oocyte maturation and sustains female fates independently of Dmrt1 repression.

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.