Transdifferentiation of differentiated ovary into functional testis by long-term treatment of aromatase inhibitor in Nile tilapia

Potential Implications of Acid-Sensing Ion Channels ASIC2 and ASIC4 in Gonadal Differentiation of Dicentrarchus labrax Subjected to Water Temperature Increase during Gonadal Development

In this study, the expression and implication of acid-sensing ion channels 2 and 4 (ASIC2 and ASIC4) in the gonadal sex differentiation of Dicentrarchus labrax (D. labrax), subjected to increasing water temperatures during gonadal development, were evaluated. Two groups were selected: a control group (CG), in which the average water temperature was maintained at 15 °C and increased to 20 °C in 20 days until weaning; and an experimental group (EG), in which the water temperature was retained at 15 °C for 60 days; thereafter, the temperature was increased daily by 0.5 °C until it reached 20 °C up to the weaning time. Ten fish from the CG and 13 fish from the EG were sampled randomly on the 335th day after hatching (dph). A higher percentage of gonad differentiation in ovaries rather than in testes was observed in the EG compared to the CG (p = 0.01). ASIC2 and ASIC4 were detected for the first time in D. labrax ovaries by indirect immunofluorescence. Both ASIC2 and ASIC4 were expressed in previtellogenic oocytes of ovaries and in scattered cells within some testes, and were most likely intratesticular previtellogenic oocytes in both the CG and EG groups. The CG group showed a higher expression of ASIC4 than the EG cohort (p < 0.05). The results gathered in this study revealed the capacity of water temperature to influence both gonadal differentiation and growth in this gonochoristic fish species, and suggests the possible role of ASIC2 and ASIC4 in gonad differentiation and gamete development in D. labrax.

Dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia

Sex is determined by multiple factors derived from somatic and germ cells in vertebrates. We have identified amhy, dmrt1, gsdf as male and foxl2, foxl3, cyp19a1a as female sex determination pathway genes in Nile tilapia. However, the relationship among these genes is largely unclear. Here, we found that the gonads of dmrt1;cyp19a1a double mutants developed as ovaries or underdeveloped testes with no germ cells irrespective of their genetic sex. In addition, the gonads of dmrt1;cyp19a1a;cyp19a1b triple mutants still developed as ovaries. The gonads of foxl3;cyp19a1a double mutants developed as testes, while the gonads of dmrt1;cyp19a1a;foxl3 triple mutants eventually developed as ovaries. In contrast, the gonads of amhy;cyp19a1a, gsdf;cyp19a1a, amhy;foxl2, gsdf;foxl2 double and amhy;cyp19a1a;cyp19a1b, gsdf;cyp19a1a;cyp19a1b triple mutants developed as testes with spermatogenesis via up-regulation of dmrt1 in both somatic and germ cells. The gonads of amhy;foxl3 and gsdf;foxl3 double mutants developed as ovaries but with germ cells in spermatogenesis due to up-regulation of dmrt1. Taking the respective ovary and underdeveloped testis of dmrt1;foxl3 and dmrt1;foxl2 double mutants reported previously into consideration, we demonstrated that once dmrt1 mutated, the gonad could not be rescued to functional testis by mutating any female pathway gene. The sex reversal caused by mutation of male pathway genes other than dmrt1, including its upstream amhy and downstream gsdf, could be rescued by mutating female pathway gene. Overall, our data suggested that dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia.

Disruption of Zar1 leads to arrested oogenesis by regulating polyadenylation via Cpeb1 in tilapia (Oreochromis niloticus)

Oogenesis is a complex process regulated by precise coordination of multiple factors, including maternal genes. Zygote arrest 1 (zar1) has been identified as an ovary-specific maternal gene that is vital for oocyte-to-embryo transition and oogenesis in mouse and zebrafish. However, its function in other species remains to be elucidated. In the present study, zar1 was identified with conserved C-terminal zinc finger domains in Nile tilapia. zar1 was highly expressed in the ovary and specifically expressed in phase I and II oocytes. Disruption of zar1 led to the failed transition from oogonia to phase I oocytes, with somatic cell apoptosis. Down-regulation and failed polyadenylation of figla, gdf9, bmp15 and wee2 mRNAs were observed in the ovaries of zar1-/- fish. Cpeb1, a gene essential for polyadenylation that interacts with Zar1, was down-regulated in zar1-/- fish. Moreover, decreased levels of serum estrogen and increased levels of androgen were observed in zar1-/- fish. Taken together, zar1 seems to be essential for tilapia oogenesis by regulating polyadenylation and estrogen synthesis. Our study shows that Zar1 has different molecular functions during gonadal development by the similar signaling pathway in different species.

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.

Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.

Bioactive compounds from four Indian medicinal plants have different potency to induce sex reversal in Nile tilapia: A chromatographic, molecular docking and in silico analysis

ETHNOPHARMACOLOGICAL RELEVANCE

Herbal plants such as Basella alba (Family: Basellaceae), Tribulus terrestris (Family: Zygophyllaceae), Asparagus racemosus (Family: Asparagaceae) and Mucuna pruriens (Family: Fabaceae) are mentioned in Indian traditional system of medicine Ayurveda to possess androgenic activity and increase male virility. The plants have been reported to improve testosterone level and sperm production in experimental male rodents as well.

AIM OF THE STUDY

Male Nile tilapias grow more quickly than females and hence are preferred for monosex Nile tilapia culture. Ethanol extracts of B. alba leaves (EB) and T. terrestris seeds (ET), and methanol extract of A. racemosus roots (MA) and M. pruriens seeds (MM) were found effective to induce masculinization in Nile tilapia. The present study intends to evaluate the anti-aromatase activity of EB, ET, MA and MM, to identify the androgenic bioactive compounds in the extracts, and to determine their pharmacokinetics. The study may validate the use of those plant extracts and their major bioactive phytoconstituents in the field of aquaculture and pharmaceuticals.

MATERIALS AND METHODS

The four crude plant extracts were first fractioned through column and thin layer chromatography (TLC). Three days old Nile tilapia juveniles (mean weight 0.025 ± 0.009g; mean length 12.50 ± 0.12 mm; n = 50 fish/replicate, 3 replicates/treatment) were then fed diets fortified with the obtained fractions for 30 days. After 30 days, fish were sacrificed and gonad aromatase mRNA expression, and 11-ketotestosterone (11-KT) and estradiol (E₂) levels were measured. Fractions yielding the highest male percentage for each plant were subject to gas chromatography-mass spectrometry (GC-MS) analysis. The in silico docking and SwissADME study were conducted with the components showing higher peak percentage in chromatogram.

RESULTS

After column chromatography and TLC analysis, EB, ET, MM and MA yielded 6 (EB1 - EB6), 8 (ET1- ET8), 14 (MM1-MM14) and 5 (MA1- MA5) fractions, respectively. Fish fed EB2, ET2, MA2 and MM13 fraction fortified diets showed significantly (p < 0.05) higher male percentage (92.32%-98.39%) compared to other treatment groups. EB2, ET2, MA2 and MM13 fed fish showed significantly (p < 0.05) higher 11-KT level compared to control male (+247.52 - +397.76%) and lower E₂ level compared to control female (-95.92% to -90.65%). Aromatase mRNA expression was significantly (p < 0.05) down-regulated by all these four fractions (-1.32 to -5.65 fold) with respect to control female. GC-MS analysis revealed the presence of 1-Octadecene (OD) in EB2, Phenol, 2,4-bis(1,1-dimethylethyl) (PD) in ET2 and MA2, 9,12-Octadecadienoic acid (Z,Z)- (ODDA) in MM13. In silico molecular docking indicated that PD is more effective than ODDA and OD to inhibit aromatase. In addition, PD showed better pharmacokinetics and more drug-likeness compared to OD and ODDA in SwissADME analysis.

CONCLUSION

The present results indicate that ET and MA are more potent to produce all-male tilapia by means of aromatase inhibition. PD can be an ideal compound to achieve masculinization in Nile tilapia through dietary administration, but further investigation is required.

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.

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

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

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

Involvement of IGF-1R-PI3K-AKT-mTOR pathway in increased number of GnRH3 neurons during androgen-induced sex reversal of the brain in female tilapia

The neuroplastic mechanism of sex reversal in the fish brain remains unclear due to the difficulty in identifying the key neurons involved. Mozambique tilapia show different reproductive behaviours between sexes; males build circular breeding nests while females hold and brood fertilized eggs in their mouth. In tilapia, gonadotropin-releasing hormone 3 (GnRH3) neurons, located in the terminal nerve, regulate male reproductive behaviour. Mature males have more GnRH3 neurons than mature females, and these neurons have been indicated to play a key role in the androgen-induced female-to-male sex reversal of the brain. We aimed to elucidate the signalling pathway involved in the androgen-induced increase in GnRH3 neurons in mature female tilapia. Applying inhibitors to organotypic cultures of brain slices, we showed that the insulin-like growth factor (IGF)-1 receptor (IGF-1R)/PI3K/AKT/mTOR pathway contributed to the androgen-induced increase in GnRH3 neurons. The involvement of IGF-1 and IGF-1R in 11-ketotestosterone (11-KT)-induced development of GnRH3 neurons was supported by an increase in Igf-1 mRNA shortly after 11-KT treatment, the increase of GnRH3 neurons after IGF-1 treatment and the expression of IGF-1R in GnRH3 neurons. Our findings highlight the involvement of IGF-1 and its downstream signalling pathway in the sex reversal of the tilapia brain.

Cyp17a2 is involved in testicular development and fertility in male Nile tilapia, Oreochromis niloticus

BACKGROUND

Steroid hormones play an essential role in many reproductive processes of vertebrates. Previous studies revealed that teleost-specific Cyp17a2 (cytochrome P450 family 17 subfamily a 2) might be required for the production of cortisol in the head-kidney and 17α,20β-dihydroxy-4-pregnen-3-one (DHP) in ovary during oocyte maturation. However, the role of Cyp17a2 in male reproduction remains to be largely unknown. The aim of this study was to explore the essentiality of cyp17a2 gene in male steroidogenesis, spermatogenesis, and male fertility.

METHODS

A homozygous mutation line of cyp17a2 gene was constructed in tilapia by CRISPR/Cas9 gene editing technology. The expression level of germ cell and meiosis-related genes and steroidogenic enzymes were detected by qRT-PCR, IHC, and Western blotting. EIA and LC-MS/MS assays were used to measure the steroid production levels. And sperm quality was examined by Sperm Quality Analyzer software.

RESULTS

In this study, cyp17a2 gene mutation resulted in the significant decline of serum DHP and cortisol levels. On the contrary, significant increases in intermediate products of cortisol and DHP were found in cyp17a2-/- male fish. The deficiency of cyp17a2 led to the arrest of meiotic initiation in male fish revealing as the reduction of the expression of germ cell-related genes (vasa, piwil, oct4) and meiosis-related genes (spo11 and sycp3) by 90 dah. Afterwards, spermatogenesis was gradually recovered with the development of testis in cyp17a2-/- males, but it showed a lower sperm motility and reduced fertility compared to cyp17a2+/+ XY fish. Deletion of cyp17a2 led to the abnormal upregulation of steroidogenic enzymes for cortisol production in the head-kidney. Moreover, unaltered serum androgens and estrogens, as well as unchanged related steroidogenic enzymes were found in the testis of cyp17a2-/- male fish.

CONCLUSION

This study proved that, for the fist time, Cyp17a2 is indispensable for cortisol and DHP production, and cyp17a2 deficiency associated curtailed meiotic initiation and subfertility suggesting the essentiality of DHP and cortisol in the male fertility of fish.

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.

Cyp17a1 is Required for Female Sex Determination and Male Fertility by Regulating Sex Steroid Biosynthesis in Fish

In teleost fish, sex steroids are involved in sex determination, sex differentiation, and fertility. Cyp17a1 (Cytochrome P450 family 17 subfamily A member 1) is thought to play essential roles in fish steroidogenesis. Therefore, to further understand its roles in steroidogenesis, sex determination, and fertility in fish, we constructed a cyp17a1 gene mutant in Nile tilapia (Oreochromis niloticus). In XX fish, mutation of the cyp17a1 gene led to a female-to-male sex reversal with a significant decline in 17β-estradiol (E2) and testosterone (T) production, and ectopic expression of male-biased markers (Dmrt1 and Gsdf) in gonads from the critical window of sex determination. Sex reversal was successfully rescued via T or E2 administration, and ovarian characteristics were maintained after termination of E2 supplementation in the absence of endogenous estrogen production in cyp17a1-/- XX fish. Likewise, deficiencies in T and 11-ketotestosterone (11-KT) production in both cyp17a1-/- XX sex-reversed males and cyp17a1-/- XY mutants resulted in meiotic initiation delays, vas deferens obstruction and sterility due to excessive apoptosis and abnormal mitochondrial morphology. However, 11-KT treatment successfully rescued the dysspermia to produce normal sperm in cyp17a1-/- male fish. Significant increases in gonadotropic hormone (gth) and gth receptors in cyp17a1-/- mutants may excessively upregulate steroidogenic gene expression in Leydig cells through a feedback loop. Taken together, our findings demonstrate that Cyp17a1 is indispensable for E2 production, which is fundamental for female sex determination and differentiation in XX tilapia. Additionally, Cyp17a1 is essential for T and 11-KT production, which further promotes spermatogenesis and fertility in XY males.

Genome-wide identification, evolution of histone lysine demethylases (KDM) genes and their expression during gonadal development in Nile tilapia

Histone lysine demethylases (KDM) are responsible for histone demethylation and are involved in gene expression regulation. Previous studies have shown that histone lysine demethylation plays an important role in gonadal development of vertebrates. The KDM family consists of eight subfamilies, i.e., kdm1, kdm2, kdm3, kdm4, kdm5, kdm6, kdm7 and JmjC-only subfamily. In this study, 13 to 63 KDM genes in 23 representative species were identified based on the available version of genome assembly. Phylogenetic relationships, domain architecture, and synteny of these genes were comprehensively analyzed and the results suggested KDM genes probably originated from the early diverging metazoan and significantly expanded in vertebrates with multiple whole genome duplication, especially in the third-round whole genome duplication (3R-WGD) and polyploidization of teleosts. The subfamilies of kdm2, kdm3, kdm4, kdm5, kdm6 and kdm7 were duplicated with 1R-2R events, and duplicates of kdm2a, kdm4a, kdm5b and kdm6b were resulted from 3R-WGD. Based on transcriptome data, the KDM genes were found to be dominantly expressed in the ovary and testis. More than 80% of KDM genes displayed sexual dimorphic expression, with 15 genes dominantly expressed in ovaries, and 12 genes dominantly expressed in testes. Importantly, from transcriptome data, qRT-PCR and fluorescence in situ hybridization during sex reversal, genes with higher expression in ovary than testis, such as kdm1b and two JmjC-only subfamily members hspbap1 and riox1, were downregulated, while other genes, such as kdm3c, kdm5bb, kdm6ba, kdm6bb and kdm7b, with higher expression in testis than ovary, were upregulated in ovotestis, indicating these genes play critical roles in the gonadal development and sex reversal. This study provided new insights into the evolution of the KDM genes and a fundamental clue for understanding their important roles in sex differentiation and gonadal development in teleosts.

Growth and gonadal development retardations after long-term exposure to estradiol in little yellow croaker, Larimichthys polyactis

Endocrine disrupting chemicals (EDCs) including 17β-estradiol (E2) are widely distributed in the aquatic environment and are known to negatively affect the reproductive system of many animals, including fish. EDCs leading to feminization, altered sex ratio and reduced fecundity, it is possibly posing potential risks to the ecosystems. To investigate the potentially toxic effects of E2 exposure on little yellow croaker (Larimichthys polyactis, L. poliactis) who have a unique gonadal development pattern that males undergo a hermaphroditic stage. An experiment was set up where L. poliactis were maintained in tanks and exposed to E2 concentrations of 10 μg/L or no E2 exposure (the ethanol and control groups) from 30 to 90 days post-hatching (dph). After exposure, the E2 withdrawal and continual cultured to 150 and 365 dph. The morphological and histological analyses were used to compare the changes in the fish body and gonad under E2 exposure. The results showed that E2 exposure caused three major phenotypes at 30 and 60 days after treatment (dat), including ovary, ovotestis and gonadal development retardation compared with the control groups. The average ratio of these three phenotypes is 60.6%, 11.97% and 27.43%, respectively. The body length and weight of E2 exposure groups were repressed during the E2 exposure period, while it can recover after E2 withdrawal. However, the gonadal development (Gonadosomatic Index) of E2 exposure groups testis were retarded at 60 dat and doesn't recover until 365 dph. The sex determination/differentiation-related genes erα, erβI, erβII, fshβ and cyp11b2 were significantly decreased in E2-exposure male fish. This research highlights the E2 leads to feminization, disrupts testis maturation and spermatogenesis, this effect persisted into the stage of sexual maturity. Collectively, our findings provide insights into the molecular mechanisms underlying E2 disturbance of a marine economic fish reproduction.

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

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

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

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

Dnmt3aa but Not Dnmt3ab Is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus)

Dnmt3a, a de novo methyltransferase, is essential for mammalian germ line DNA methylation. Only one Dnmt3a is identified in mammals, and homozygous mutants of Dnmt3a are lethal, while two Dnmt3a paralogs, dnmt3aa and dnmt3ab, are identified in teleosts due to the third round of genome duplication, and homozygous mutants of dnmt3aa and dnmt3ab are viable in zebrafish. The expression patterns and roles of dnmt3aa and dnmt3ab in gonadal development remain poorly understood in teleosts. In this study, we elucidated the precise expression patterns of dnmt3aa and dnmt3ab in tilapia gonads. Dnmt3aa was highly expressed in oogonia, phase I and II oocytes and granulosa cells in ovaries and spermatogonia and spermatocytes in testes, while dnmt3ab was mainly expressed in ovarian granulosa cells and testicular spermatocytes. The mutation of dnmt3aa and dnmt3ab was achieved by CRISPR/Cas9 in tilapia. Lower gonadosomatic index (GSI), increased apoptosis of oocytes and spermatocytes and significantly reduced sperm quality were observed in dnmt3aa^−/− mutants, while normal gonadal development was observed in dnmt3ab^−/− mutants. Consistently, the expression of apoptotic genes was significantly increased in dnmt3aa^−/− mutants. In addition, the 5-methylcytosine (5-mC) level in dnmt3aa^−/− gonads was decreased significantly, compared with that of dnmt3ab^−/− and wild type (WT) gonads. Taken together, our results suggest that dnmt3aa, not dnmt3ab, plays important roles in maintaining gametogenesis in teleosts.

Turning the spotlight on the C11-oxy androgens in human fetal development

Research into the biosynthesis of C11-oxy C₁₉ steroids during human fetal development, specifically fetal adrenal development and during the critical period of sex differentiation, is currently lacking. Cortisol, which possesses a C11-hydroxyl moiety has, however, been firmly established in this context. Compelling questions are whether the C11-oxy C₁₉ steroids (11β-hydroxyandrostenedione, 11β-hydroxytestosterone, 11-ketoandrostenedione and 11-ketotestosterone [11KT]) and the C11-oxy C₂₁ steroids (11β-hydroxyprogesterone and 11-ketoprogesterone) are biosynthesised during gestation, and whether these hormones circulate between the placenta and the developing fetus, and between the placenta and the mother. This review will consider the role of cortisol, 11KT and 11β-hydroxysteroid dehydrogenase type 2 (11βHSD2) in determining the sex of teleost fish, while these hormones and 11βHSD2 will also be discussed with regards to murine mammals. The focus of the review will shift to highlight the potential role of C11-oxy steroids in human fetal development based on the timely expression of steroidogenic enzymes in the adrenal, testes and ovary, as well as in the placenta; summarising reported evidence of C11-oxy steroids in neonatal life.

Characterization, expression and CpG methylation analysis of Dmrt1 and its response to steroid hormone in blotched snakehead (Channa maculata)

Dmrt1 is an important transcriptional regulator that plays critical role in male gonadogenesis, testicular differentiation and development. In this study, Dmrt1 was cloned from blotched snakehead (Channa maculata), which is designated as CmDmrt1. CmDmrt1 encoded a putative protein with 293 amino acids and presented an extremely conserved DM domain. It was nearly expressed in the gonads, and the expression was more than 15 times higher in the testis than in the ovary. 1851 bp promoter sequence of CmDmrt1 was characterized and the methylation levels of the CpG sites were analyzed to detect sex-related differences. A significant negative correlation between CmDmrt1 expression and CpG methylation level of its promoter was found in the testis and ovary. During gonadal development, CmDmrt1 transcription displayed strong male-biased expression patterns, increased with the maturation of testis and reached the peak at 195 days after hatching (dah), which indicates a significant role of Dmrt1 in spermatogenesis. Steroid treatment could influence CmDmrt1 expression, and long-term 17β-estradiol (E2) treatment could induce the male-to-female secondary sex reversal (SSR), which resulted in the differentiated testis transformed to ovary or ovotestis. Meanwhile, CmDmrt1 expression was down-regulated to fairly low level in the ovary of the SSR XY fish, which was similar to that in normal XX females ovary. Our research illustrates that Dmrt1 is linked to testis differentiation and spermatogenesis in blotched snakehead, providing information for functional studies on sex differentiation and gonadal development of C. maculata, and scientific basis for the production practice of all-male snakehead breeding.

Characterization, expression, and regulatory effects of nr0b1a and nr0b1b in spotted scat (Scatophagus argus)

Nuclear receptor subfamily 0 group B member 1 (Nr0b1) belongs to the nuclear receptor (NR) superfamily. It plays critical roles in sex determination, sex differentiation, and gonadal development in mammals. In this study, the duplicated genes nr0b1a and nr0b1b were identified in spotted scat (Scatophagus argus). Phylogenetic and synteny analyses revealed that, unlike nr0b1a, nr0b1b was retained in several species of teleosts after an nr0b1 gene duplication event but was secondarily lost in other fish species, amphibians, reptiles, birds, and mammals. In a sequence analysis, only 1.5 LXXLL-related repeat motifs were identified in spotted scat Nr0b1a, Nr0b1b, and non-mammalian Nr0b1a/Nr0b1, different from the 3.5 repeat motifs in mammalian Nr0b1. By qPCR, nr0b1a and nr0b1b were highly expressed in testes from stages IV to V and in ovaries from stages II to IV, respectively. Male-to-female sex reversal was induced in XY spotted scat by the administration of exogenous E2. A qPCR analysis showed that nr0b1b mRNA expression was higher in sex-reversed XY fish than in control XY fish, with no difference in nr0b1a. A luciferase assay showed that spotted scat Nr0b1a and Nr0b1b did not individually activate cyp19a1a gene transcription. As in mammals, spotted scat Nr0b1a suppressed Nr5a1-mediated cyp19a1a expression, despite containing only 1.5 LXXLL-related repeat motifs in its N-terminal region, while Nr0b1b stimulated Nr5a1-mediated cyp19a1a transcription. These results demonstrated that nr0b1a and nr0b1b in spotted scat have distinct expression patterns and regulatory effects and further indicate that nr0b1b might be involved in ovarian development by regulating Nr5a1-mediated cyp19a1a expression.

Global expression response of genes in sex-undifferentiated Nile tilapia gonads after exposure to trace letrozole

The aromatase inhibitor letrozole can be found in rivers, effluents, and even drinking water. Studies have demonstrated that letrozole affects various metabolic pathways and may cause reproductive toxicity, especially in fish exposed during development. However, studies on the effect of a low concentration of letrozole at the whole-gonad transcriptomic level in the early stage of fish sexual development have not been investigated. The aim of our study was to explore the potential effects of a low concentration of letrozole on the gonad transcriptome of Nile tilapia at an early stage of sexual development. In this study, 9 dpf (days postfertilization) Nile tilapia were exposed to trace letrozole for 12 days. Letrozole exposure from 9 dpf to 21 dpf persistently altered phenotypic sex development and induced the male-biased sex ratio. The transcriptome results showed that 1173 differentially expressed genes (DEGs) were present in the female control vs 1.5 μg/L letrozole-treated female comparison group and that 1576 DEGs were present in the 1.5 μg/L letrozole-treated female vs male control comparison group. Differentially expressed gene enrichment analysis revealed several crucial pathways, including the drug metabolism-cytochrome P450 pathway, the ErbB-PI3K/Akt/mTOR pathway, and the calcium signalling pathway. Further analysis of these identified DEGs indicated that some key genes correlated with metabolism and epigenetic regulation were significantly affected by letrozole, such as UDP-glucuronosyltransferase (Ugt), glutathione S-transferase omega-1 (Gsto1), lysine-specific demethylase 6bb (Kdm6bb, original name is Kdm6a), jumonji and AT-rich interaction domain containing 2 (Jarid2b, original name is Jarid2), growth arrest and DNA damage inducible gamma (Gadd45g), and chromobox protein 7 (Cbx7). The qRT-PCR validation results for twelve DEGs showed that the Pearson's correlation of the log₁₀^{fold change} values between the qPCR and RNA-Seq results was 0.90, indicating the accuracy and reliability of the RNA-Seq results. Our study is the first to report the effect of letrozole on the transcriptome of gonads from fish during early-stage sexual development. These findings will be useful for understanding the toxic effects and molecular mechanisms of letrozole exposure at the early stage of gonad development on the sexual development of aquatic organisms.

miR-133b targets tagln2 and functions in tilapia oogenesis

microRNAs (miRNAs) are important components of non-coding RNAs that participate in diverse life activities by regulating gene expression at the post transcriptional level through base complementary pairing with 3'UTRs of target mRNAs. miR-133b is a member of the miR-133 family, which play important roles in muscle differentiation and tumorigenesis. Recently, miR-133b was reported to affect estrogen synthesis by targeting foxl2 in mouse, while its role in fish reproduction remains to be elucidated. In the present study, we isolated the complete sequence of miR-133b, which was highly expressed in tilapia ovary at 30 and 90 dah (days after hatching) and subsequently decreased at 120 to 150 dah by qPCR. Interestingly, only a few oogonia were remained in the antagomir-133b treated tilapia ovary, while phase I and II oocytes were observed in the ovaries of the control group. Unexpectedly, the expression of foxl2 and cyp19a1a, as well as estradiol levels in serum were increased in the treated group. Furthermore, tagln2, an important factor for oogenesis, was predicted as the target gene of miR-133b, which was confirmed by dual luciferase reporter vector experiments. miR-133b and tagln2 were co-expressed in tilapia ovaries. Taken together, miR-133b may be involved in the early oogenesis of tilapia by regulating tagln2 expression. This study enriches the understanding of miR-133b function during oogenesis and lays a foundation for further study of the regulatory network during oogenesis.

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.

Identification, Expression and Evolution of Short-Chain Dehydrogenases/Reductases in Nile Tilapia (Oreochromis niloticus)

The short-chain dehydrogenases/reductases (SDR) superfamily is involved in multiple physiological processes. In this study, genome-wide identification and comprehensive analysis of SDR superfamily were carried out in 29 animal species based on the latest genome databases. Overall, the number of SDR genes in animals increased with whole genome duplication (WGD), suggesting the expansion of SDRs during evolution, especially in 3R-WGD and polyploidization of teleosts. Phylogenetic analysis indicated that vertebrates SDRs were clustered into five categories: classical, extended, undefined, atypical, and complex. Moreover, tandem duplication of hpgd-a, rdh8b and dhrs13 was observed in teleosts analyzed. Additionally, tandem duplications of dhrs11-a, dhrs7a, hsd11b1b, and cbr1-a were observed in all cichlids analyzed, and tandem duplication of rdh10-b was observed in tilapiines. Transcriptome analysis of adult fish revealed that 93 SDRs were expressed in more than one tissue and 5 in one tissue only. Transcriptome analysis of gonads from different developmental stages showed that expression of 17 SDRs were sexually dimorphic with 11 higher in ovary and 6 higher in testis. The sexually dimorphic expressions of these SDRs were confirmed by in situ hybridization (ISH) and qPCR, indicating their possible roles in steroidogenesis and gonadal differentiation. Taken together, the identification and the expression data obtained in this study contribute to a better understanding of SDR superfamily evolution and functions in teleosts.

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.

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

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

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

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.

Activation of stress response axis as a key process in environment-induced sex plasticity in fish

The determination of sex is an important hallmark in the life cycle of organisms, in which the fate of gonads and then the individual sex are defined. In gonochoristic teleost fish, this process is characterized by a high plasticity, considering that in spite of genotypic sex many environmental factors can cause shifts from one to another molecular pathway, resulting in organisms with mismatching genotypic and phenotypic sexes. Interestingly, in most instances, both female-to-male or male-to-female sex-reversed individuals develop functional gonads with normal gametogenesis and respective progenies with full viability. The study of these mechanisms is being spread to other non-model species or to those inhabiting more extreme environmental conditions. Although water temperature is an important mechanism involved in sex determination, there are other environmental stressors affected by the climate change which are also implicated in stress response-induced masculinization in fish. In this regard, the brain has emerged as the transducer of the environment input that can influence the gonadal fate. Furthermore, the evaluation of other environmental stressors or their synergic effect on sex determination at conditions that simulate the natural environments is growing gradually. Within such scope, the concerns related to climate change impacts rely on the fact that many of biotic and abiotic parameters reported to affect sex ratios are expected to increase concomitantly as a result of increased greenhouse gas emissions and, particularly worrying, many of them are related to male bias in the populations, such as high temperature, hypoxia, and acidity. These environmental changes can also generate epigenetic changes in sex-related genes affecting their expression, with implications on sex differentiation not only of exposed individuals but also in following generations. The co-analysis of multi-stressors with potential inter- and transgenerational effects is essential to allow researchers to perform long-term predictions on climate change impacts in wild populations and for establishing highly accurate monitoring tools and suitable mitigation strategies.

Regulation of Female Folliculogenesis by Tsp1a in Nile Tilapia (Oreochromis niloticus)

TSP1 was reported to be involved in multiple biological processes including the activation of TGF-β signaling pathways and the regulation of angiogenesis during wound repair and tumor growth, while its role in ovarian folliculogenesis remains to be elucidated. In the present study, Tsp1a was found to be expressed in the oogonia and granulosa cells of phase I to phase IV follicles in the ovaries of Nile tilapia by immunofluorescence. tsp1a homozygous mutants were generated by CRISPR/Cas9. Mutation of tsp1a resulted in increased oogonia, reduced secondary growth follicles and delayed ovary development. Expression of the cell proliferation marker PCNA was significantly up-regulated in the oogonia of the mutant ovaries. Furthermore, transcriptomic analysis revealed that expressions of DNA replication related genes were significantly up-regulated, while cAMP and MAPK signaling pathway genes which inhibit cell proliferation and promote cell differentiation were significantly down-regulated. In addition, aromatase (Cyp19a1a) expression and serum 17β-estradiol (E2) concentration were significantly decreased in the mutants. These results indicated that lacking tsp1a resulted in increased proliferation and inhibited differentiation of oogonia, which in turn, resulted in increased oogonia, reduced secondary growth follicles and decreased E2. Taken together, our results indicated that tsp1a was essential for ovarian folliculogenesis in Nile tilapia.

Rapid Induction of Female-to-Male Sex Change in Adult Zebrafish by Injection of an Aromatase Inhibitor

Previously, we examined whether aromatase inhibitor (AI) treatment induces a sex change in adult female zebrafish. A 5-month AI treatment regime resulted in the retraction of the ovaries and testis formation. Eight weeks after changing the diet to AI-free food, a large number of normal sperm were obtained. Artificial fertilization using sperm from the sex-changed females was successful. These results demonstrated that sex plasticity remains in the mature ovaries of zebrafish. However, >7 months of treatment was necessary; thus, pairing was unsuccessful. In this study, we tried to induce sex change through the injection of an AI to shorten the time course of sex change. When the AI solution was directly injected into the abdomen of zebrafish, retraction of the ovary was induced within 2 months. The natural mating of sex-changed females with normal females was successful at 3 months. Although the fertilization rate was low, juveniles resulting from these matings developed normally. We succeeded in establishing a method for inducing sex changes in adult zebrafish within 3 months. The procedure will support the study of how sexual plasticity persists in adult zebrafish following sex differentiation and the identification of undifferentiated stem cells.

Genome-Wide Identification, Evolution and Expression of the Complete Set of Cytoplasmic Ribosomal Protein Genes in Nile Tilapia

Ribosomal proteins (RPs) are indispensable in ribosome biogenesis and protein synthesis, and play a crucial role in diverse developmental processes. In the present study, we carried out a comprehensive analysis of RPs in chordates and examined the expression profiles of the complete set of 92 cytoplasmic RP genes in Nile tilapia. The RP genes were randomly distributed throughout the tilapia genome. Phylogenetic and syntenic analyses revealed the existence of duplicated RP genes from 2R (RPL3, RPL7, RPL22 and RPS27) and 3R (RPL5, RPL19, RPL22, RPL41, RPLP2, RPS17, RPS19 and RPS27) in tilapia and even more from 4R in common carp and Atlantic salmon. The RP genes were found to be expressed in all tissues examined, but their expression levels differed among different tissues. Gonadal transcriptome analysis revealed that almost all RP genes were highly expressed, and their expression levels were highly variable between ovaries and testes at different developmental stages in tilapia. No sex- and stage-specific RP genes were found. Eleven RP genes displayed sexually dimorphic expression with nine higher in XY gonad and two higher in XX gonad at all stages examined, which were proved to be phenotypic sex dependent. Quantitative real-time PCR and immunohistochemistry ofRPL5b and RPL24 were performed to validate the transcriptome data. The genomic resources and expression data obtained in this study will contribute to a better understanding of RPs evolution and functions in chordates.

Chronic simultaneous exposure of common carp (Cyprinus carpio) from embryonic to juvenile stage to drospirenone and gestodene at low ng/L level caused intersex

Synthetic progestins are emerging contaminants of the aquatic environment with endocrine disrupting potential. The main aim of the present study was to investigate the effects of the synthetic progestins gestodene, and drospirenone on sex differentiation in common carp (Cyprinus carpio) by histological analysis. To gain insights into the mechanisms behind the observations from the in vivo experiment on sex differentiation, we analyzed expression of genes involved in hypothalamus-pituitary-gonad (HPG) and hypothalamus-pituitary-thyroid (HPT) axes, histology of hepatopancreas, and in vitro bioassays. Carp were continuously exposed to concentrations of 2 ng/L of single progestins (gestodene or drospirenone) or to their mixture at concentration 2 ng/L of each. The exposure started 24 h after fertilization of eggs and concluded 160 days post-hatching. Our results showed that exposure of common carp to a binary mixture of drospirenone and gestodene caused increased incidence of intersex (32%) when compared to clean water and solvent control groups (both 3%). Intersex most probably was induced by a combination of multiple modes of action of the studied substances, namely anti-gonadotropic activity, interference with androgen receptor, and potentially also with HPT axis or estrogen receptor.

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

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

A conceptual framework for understanding sexual differentiation of the teleost brain

Vertebrate brains are sexually differentiated, giving rise to differences in various physiological and behavioral phenotypes between the sexes. In developing mammals and birds, the neural substrate underlying sex-dependent physiology and behavior undergoes an irreversible process of sexual differentiation due to the effects of perinatal gonadal steroids and sex chromosome complement. The differentiated neural substrate is then activated in the adult by the sex-specific steroid milieu to facilitate the expression of sex-typical phenotypes. However, this well-established concept does not hold for teleost fish, whose sexual phenotypes (behavioral or otherwise) are highly labile throughout life and can be reversed even in adulthood. Indeed, the available evidence suggests that, in teleosts, neither gonadal steroids early in development nor the sex chromosome complement contribute much to brain sexual differentiation; instead, steroids in adulthood serve to both differentiate the neural substrate and activate it to elicit sex-typical phenotypes in a transient and reversible manner. Evidence further suggests that marked sexual dimorphisms and adult steroid-dependent lability in the neural expression of sex steroid receptors constitute the primary molecular basis for sexual differentiation and lability of the teleost brain. The consequent sexually dimorphic but reversible steroid sensitivity in response to the adult steroid milieu may enable the teleost brain to maintain lifelong sexual lability and to undergo phenotypic sex reversal.

Transcriptional effects in the estuarine guppy Poecilia vivipara exposed to sanitary sewage in laboratory and in situ

The urban growth has increased sanitary sewage discharges in coastal ecosystems, negatively affecting the aquatic biota. Mangroves, one of the most human-affected coastal biomes, are areas for reproduction and nursing of several species. In order to evaluate the effects of sanitary sewage effluents in mangrove species, this study assessed the hepatic transcriptional responses of guppy fish Poecilia vivipara exposed to sanitary sewage 33% (v:v), using suppressive subtraction hybridization (SSH), high throughput sequencing of RNA (Ion-proton) and quantification of transcript levels by qPCR of some identified genes in fish kept in a sewage-contaminated environment. Genes identified are related predominantly to xenobiotic biotransformation, immune system and sexual differentiation. The qPCR results confirmed the induction of cytochrome P450 1A (CYP1A), glutathione S transferase A-like (GST A-like) methyltransferase (MET) and UDP glycosyltransferase 1A (UDPGT1A), and repression of complement component C3 (C3), doublesex and mab-3 related transcription factor 1 (DMRT1), and transferrin (TF) in the laboratory experiment. In the field exposure, the transcript levels of CYP1A, DMRT1, MET, GST A-like and UDPGT1A were higher in fishes exposed at the contaminated sites compared to the reference site. Chemical analysis in fish from the laboratory and in situ experiments, and surface sediment from the sewage-contaminated sites revealed relevant levels of polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyl (PCBs) and linear alkylbenzenes (LABs). These data reinforce the use of P. vivipara as a sentinel for monitoring environmental contamination in coastal regions.

Transcriptome Analysis of Ovary Development in Nile Tilapia Under Different Photoperiod Regimes

This study investigated the molecular mechanisms involved in ovarian transcriptomic responses in Nile tilapia under different photoperiod regimes. Histological analysis indicated that ovarian development was significantly affected by photoperiod. The photoperiods tested were as follows: LD (12 h light:12 h dark), LL (24 h light:0 h dark), and DD (0 h light:24 h dark). The longer photoperiod (LL) was shown to induce ovary development earlier than LD and DD. Next, ovary transcriptome levels were sequenced and analyzed. These data indicated that 988, 992, and 1,036 differentially expressed genes (DEGs) were detected by comparing LD–LL, LD–DD, and LL–DD. A number of genes that may be involved in photoperiod-specific regulation of ovarian development were observed. These findings may be useful for investigating the molecular mechanisms underlying light-induced ovarian development in Nile tilapia.

Effects of long term antiprogestine mifepristone (RU486) exposure on sexually dimorphic lncRNA expression and gonadal masculinization in Nile tilapia (Oreochromis niloticus)

Mifepristone (RU486), a clinical abortion agent and potential endocrine disruptor, binds to progestin and glucocorticoid receptors and has multiple functional importance in reproductive physiology. A long-term exposure of RU486 resulted in masculinization of female fish, however, the epigenetic landscape remains elusive. Recent studies demonstrated that long non-coding RNAs (lncRNAs) might play potential roles in epigenetic modulation of sex differentiation, ovarian cancer and germline stem cell survival. To further understand the influence of RU486 exposure on epigenetic regulation, we performed a comparative investigation on sex-biased gonadal lncRNAs profiles using control XX/XY and RU486-induced sex reversed XX Nile tilapia (Oreochromis niloticus) by RNA-seq. In total, 962 sexually differentially expressed lncRNAs and their target genes were screened from the gonads of control and sex reversed fish. In comparison with the control XX group, sex reversal induced by RU486 treatment led to significant up-regulation of 757 lncRNAs and down-regulation of 221 lncRNAs. Hierarchical clustering analysis revealed that global lncRNA expression profiles in RU486-treated XX group clustered into the same branch with the control XY, whereas XX control group formed a separate branch. The KEGG pathway enrichment analysis showed that the cis-target genes between RU486-XX and control-XX were concentrated in NOD - like receptor signaling pathway, Cell adhesion molecules (CAMs) and Biosynthesis of amino acids. Real-time PCR and in situ hybridization experiments demonstrate that lncRNAs showing intense fluctuation during RU486 treatment are also sexually dimorphic during early sex differentiation, which further proves the intimate relationship between lncRNAs and sex differentiation and sexual transdifferentiation. Taken together, our data strongly indicates that a long-term exposure of RU486 resulted in sex reversal of XX female fish and the altered expression of sexually dimorphic lncRNAs might partially account for the sex reversal via epigenetic modification.

Establishment of three estrogen receptors (esr1, esr2a, esr2b) knockout lines for functional study in Nile tilapia

Estrogens play fundamental roles in regulating reproductive activities and they act through estrogen receptors (ESRs) in all vertebrates. To date, distinct roles of estrogen receptors have been characterized only in human and model organisms, including mouse, rat, zebrafish and medaka. Physiological role of estrogen/receptor signaling in reproduction remains poorly defined in non-model organisms. In the present study, we successfully generated esr1, esr2a and esr2b mutant lines in tilapia by CRISPR/Cas9 and examined their phenotypes. Surprisingly, the esr1 mutants showed no phenotypes of reproductive development and function in both females and males. The esr2a mutant females showed significantly delayed ovarian development and follicle growth at 90 and 180 dah, and the development caught up later at 360 dah. The esr2a mutant males showed no phenotypes at 90 dah, and displayed smaller gonads and efferent ducts, less spermatogonia and more abnormal sperms at 180 dah. In contrast, the esr2b mutants displayed abnormal development of ovarian ducts and efferent ducts which failed to connect to the genital orifice, and which in turn, resulted in infertility in female and male, respectively, although they produced gametes in their gonads. Taken together, our study provides evidence for differential functions of esr1, esr2a and esr2b in fish reproduction.

The germline-specific expression of Foxl3a and its paralogous Foxl3b are associated with male gonadal differentiation in the Japanese eel, Anguilla japonica

Unlike its paralog Foxl2, which is well known for its role in ovarian development in vertebrates, the function of Foxl3 is still unclear. Foxl3 is an ancient duplicated copy of Foxl2. It is present as a single copy in ray-finned fish. But, due to repeated losses, it is absent in most tetrapods. Our transcriptomic data, however, show that two Foxl3s (Foxl3a and its paralog Foxl3b) are present in Japanese eel. Foxl3a is predominantly expressed in the pituitary, and Foxl3b is predominantly expressed in the gills. Both Foxl3s show a sex-dimorphic expression, being higher expression in testes than in ovaries. Moreover, Foxl3a and Foxl3b were exclusively expressed during gonadal differentiation in control eels (100% male). Conversely, Foxl3a and Foxl3b significantly decreased after gonadal differentiation in E2-treated eels (100% female). Furthermore, in accordance the difference in adhesive ability between somatic cells and germline cells in testes, Foxl3s showed a high expression in suspension cells (putative germline cells) and low expression in adhesive cells (putative somatic cells). In situ hybridization further showed that Foxl3a and Foxl3b were expressed in the testicular germline cells. In addition, Foxl3s expression was not changed by sex steroids in in vitro testes culture. Taken together, our results suggest that the teleost-specific Foxl3 paralog was repeatedly lost in most fish after the third round of whole genome duplication. The two germline-expressed Foxl3s had higher expression levels in males than in females during gonadal differentiation in Japanese eel. These results demonstrated that Foxl3s might play an important role in germline sexual fate determination from ancient fish to modern fish.

Roles of estrogens in fish sexual plasticity and sex differentiation

Fish sex could be reversed at the undifferentiated stage of gonad by administration of exogenous estrogen (E2) or blockade of endogenous estrogen synthesis with aromatase inhibitors, which is designated as primary sex reversal (PSR). Recent studies have well demonstrated that gonochoristic fish maintain their sexual plasticity after sex determination/differentiation. The differentiated ovary could be transdifferentiated into functional testis, and vice versa, the differentiated testis could be transdifferentiated into ovary. By analyzing these two secondary sex reversal (SSR) models, it was found that induction of male-to-female sex reversal initiates from dorsal (near the blood vessel) to the ventral, while induction of female-to-male sex reversal initiates from the ventral to dorsal. Down regulation of endogenous estrogen is the prerequisite for the ovarian transdifferentiation. However, exogenous estrogen alone is not sufficient for inducing differentiated testis to ovary. Administration of E2 and simultaneous blockage of androgen synthesis could induce testicular transdifferentiation. Therefore, endogenous estrogen is critical for the ovarian differentiation/maintenance and androgen is critical for testicular maintenance. Recently, genetic studies with genome editing technologies also showed that disruption of Cyp19a1a induced testicular development, indicating that cyp19a1a is the key gene essential for estrogen synthesis and ovary differentiation/maintenance. Knockout of male pathway genes or overexpression of female pathway genes could up-regulate cyp19a1a expression and increase estrogen level so as to promote ovary. Conversely, knockout of female pathway genes or overexpression of male pathway genes could down-regulate cyp19a1a expression and decrease estrogen level so as to promote testis (transgenic or knockout sex reversal, TSR). Epigenetic regulation of cyp19a1a play a critical role in natural sex reversal (NSR), but its relation with PSR, SSR and TSR needs further detailed investigations. In all, these studies further highlighted the important roles of endogenous estrogens in fish sex differentiation/maintenance.

Toxic effects of bisphenol A on goldfish gonad development and the possible pathway of BPA disturbance in female and male fish reproduction

Bisphenol A (BPA) is an abundant endocrine-disrupting compound that is found in the aquatic environment and has adverse effects on fish reproduction; however, the exact pathway of these impacts is unclear. In this study, the different effects of BPA on ovarian and testis development in goldfish (Carassius auratus) and the different mechanisms underlying these effects were investigated. The gonadosomatic index (GSI) and gonadal histology demonstrated that BPA diminished ovarian maturation in goldfish, which recovered after BPA treatment withdrawal. In males, BPA disrupted testis maturation, but this disruption could not be recovered after BPA treatment withdrawal. The hypothalamic-pituitary-gonad (HPG) axis-related genes sgnrh, fshβ and lhβ were significantly decreased in BPA-treated female fish, while no changes in sex steroid hormone levels and no TUNEL and PCNA staining were found in the ovary, suggesting that BPA may reduce ovarian maturation through the HPG axis. In male fish, TUNEL staining was found in 1 μg L^-1 BPA-exposed germ cells and 50 and 500 μg L^-1 BPA-exposed Leydig cells. Decreases in 11-KT levels were also found in 50 and 500 μg L^-1 BPA-exposed fish, but BPA did not affect genes associated with the HPG axes. This result shows that BPA disrupts testis maturation through apoptosis of germ cells and Leydig cells, thus inducing decreases in 11-KT levels that disrupt spermatogenesis. Collectively, our findings provide insights into the molecular and cellular mechanisms underlying BPA disturbance of goldfish reproduction.

Sexually Dimorphic Neuropeptide B Neurons in Medaka Exhibit Activated Cellular Phenotypes Dependent on Estrogen

Brain and behavior of teleosts are highly sexually plastic throughout life, yet the underlying neural mechanisms are largely unknown. On examining brain morphology in the teleost medaka (Oryzias latipes), we identified distinctively large neurons in the magnocellular preoptic nucleus that occurred much more abundantly in females than in males. Examination of sex-reversed medaka showed that the sexually dimorphic abundance of these neurons is dependent on gonadal phenotype, but independent of sex chromosome complement. Most of these neurons in females, but none in males, produced neuropeptide B (Npb), whose expression is known to be estrogen-dependent and associated with female sexual receptivity. In phenotypic analysis, the female-specific Npb neurons had a large euchromatic nucleus with an abundant cytoplasm containing plentiful rough endoplasmic reticulum, exhibited increased overall transcriptional activity, and typically displayed a spontaneous regular firing pattern. These phenotypes, which are probably indicative of cellular activation, were attenuated by ovariectomy and restored by estrogen replacement. Furthermore, the population of Npb-expressing neurons emerged in adult males treated with estrogen, not through frequently occurring neurogenesis in the adult teleost brain, but through the activation of preexisting, quiescent male counterpart neurons. Collectively, our results demonstrate that the morphological, transcriptional, and electrophysiological phenotypes of sexually dimorphic preoptic Npb neurons are highly dependent on estrogen and can be switched between female and male patterns. These properties of the preoptic Npb neurons presumably underpin the neural mechanism for sexual differentiation and plasticity of brain and behavior in teleosts.

New Insights Into the Role of Follicle-Stimulating Hormone in Sex Differentiation of the Protogynous Orange-Spotted Grouper, Epinephelus coioides

Follicle-stimulating hormone (FSH) signaling is considered to be essential for early gametogenesis in teleosts, but its functional roles during sex differentiation are largely unknown. In this study, we investigated the effects of long-term and short-term FSH injection on sex differentiation in the protogynous orange-spotted grouper (Epinephelus coioides). Long-term FSH treatment initially promoted the formation of ovaries but subsequently induced a male fate. The expression of female pathway genes was initially increased but then decreased, whereas the expression of male pathway genes was up-regulated only during long-term FSH treatment. The genes related to the synthesis of sex steroid hormones, as well as serum 11-ketotestosterone and estradiol, were also up-regulated during long-term FSH treatment. Short-term FSH treatment activated genes in the female pathway (especially cyp19a1a) at low doses but caused inhibition at high doses. Genes in the male pathway were up-regulated by high concentrations of FSH over the short term. Finally, we found that low, but not high, concentrations of FSH treatment activated cyp19a1a promoter activities in human embryonic kidney (HEK) 293 cells. Overall, our data suggested that FSH may induce ovarian differentiation or a change to a male sex fate in the protogynous orange-spotted grouper, and that these processes occurred in an FSH concentration-dependent manner.