Gonadal development and expression of sex-specific genes during sex differentiation in the Japanese eel

The expression profiles of cyp19a1, sf-1, esrs and gths in the brain-pituitary during gonadal sex differentiation in juvenile Japanese eels

Eels are gonochoristic species whose gonadal differentiation initiates at the yellow eel stage and is influenced by environmental factors. We revealed some sex-related genes were sex dimorphically expressed in gonads during gonadal sex differentiation of Japanese eel (Anguilla japonica); however, the expression of sex-related genes in the brain-pituitary during gonadal sex differentiation in eels is still unclear. This study aimed to investigate the sex-related gene expressions in the brain-pituitary and tried to clarify their roles in the brain and gonads during gonadal sex differentiation. Based on our previous histological study, the control eels developed as males, and estradiol-17β (E2) was used for feminization. Our results showed that during testicular differentiation, the brain cyp19a1 transcripts and aromatase proteins were increased significantly; moreover, the cyp19a1, sf-1, foxl2s, and esrs (except gperb) transcripts in the midbrain/pituitary also were increased significantly. Forebrain gnrh1 transcripts increased slightly during gonadal differentiation of both sexes, but the gnrhr1b and gnrhr2 transcripts in the midbrain/pituitary were stable during gonadal differentiation. The expression levels of gths and gh in the midbrain/pituitary were significantly increased during testicular differentiation and were much higher in males than in E2-feminized females. These results implied that endogenous estrogens might play essential roles in the brain/pituitary during testicular differentiation, sf-1, foxl2s, and esrs may have roles in cyp19a1 regulation in the midbrain/pituitary of Japanese eels. For the GnRH-GTH axis, gths, especially fshb, may be regulated by esrs and involved in regulating testicular differentiation and development in Japanese eels.

Genome-wide identification and expressional profile of the Dmrt gene family in the swimming crab (Portunus trituberculatus)

The swimming crab, Portunus trituberculatus is one of crucial aquaculture crabs with significant differences in growth and economic performance between male and female swimming crabs. Consequently, the culture of female populations presents higher economic value. The doublesex and mab-3 related transcription factor (Dmrt) gene family are known to play crucial role in gonad differentiation and development. However, there is limited information about this gene family in Portunus trituberculatus. In this study, we identified seven members of the Dmrt gene family in P. trituberculatus based on the published transcriptome and genome data and designated as Ptdmrt-1, Ptdoublesex (Ptdsx), Ptidmrt-1, Ptdmrt-11E, Ptidmrt-2, Ptdmrt-99B, and Ptdmrt-3 based on the homology analysis results, respectively. These Ptdmrt genes distributed across 6 chromosomes and were predicted to encode 283 aa, 288 aa, 529 aa, 436 aa, 523 aa, 224 aa, and 435 aa protein precursors, respectively. The expression patterns of these dmrt genes were characterized by qRT-PCR and gonad transcriptome data. The results showed that five members (Ptdmrt-99B, Ptdsx, Ptdmrt-1, Ptdmrt-3, and Ptdmrt-11E) were differentially expressed between the testis and ovary. In addition, their expression patterns from zoea 2 to juvenile 1 were characterized by published transcriptome data and the results showed that they were lowly expressed and did not exhibit notable difference except for Ptdsx during early development. Overall, majority of Ptdmrt genes were involved in gonad differentiation in the swimming crab. Current findings provide a solid foundation for further exploration of the roles of these genes in gonad development and differentiation in P. trituberculatus.

Characterization of potential spermatogonia biomarker genes in the European eel (Anguilla anguilla)

Identification of specific molecular markers for spermatogonial stem cells in teleost is crucial for enhancing the efficacy of reproductive biotechnologies in aquaculture, such as transplantation and surrogate production in fishes. Since it is not yet possible to distinguish spermatogonial stem cells of European eel (Anguilla anguilla) using specific molecular markers, we isolated spermatogonial cells from immature European eels to find these potential markers. We attempted this by studying three candidate genes: vasa, nanos2, and dnd1. Two vasa (vasa1 and vasa2) genes, nanos2, and dnd1 were identified, characterized, and studied in the muscle, testis, and isolated spermatogonia. Our results showed that vasa1 and vasa2 had the highest levels of expression when measured by qPCR. In situ hybridization and immunochemistry assays showed that the four genes were localized explicitly in type A spermatogonia. However, vasa1 and vasa2 exhibited stronger signals in the immature testicular tissue than the other two potential markers. According to this, vasa1 and vasa2 were found to be the most effective markers for spermatogonial cells in the European eel.

Molecular Characterization, Expression Pattern, DNA Methylation and Gene Disruption of Figla in Blotched Snakehead (Channa maculata)

Figla is one of the earliest expressed genes in the oocyte during ovarian development. In this study, Figla was characterized in C. maculata, one of the main aquaculture species in China, and designated as CmFigla. The length of CmFigla cDNA was 1303 bp, encoding 197 amino acids that contained a conserved bHLH domain. CmFigla revealed a female-biased expression patterns in the gonads of adult fish, and CmFigla expression was far higher in ovaries than that in testes at all gonadal development stages, especially at 60~180 days post-fertilization (dpf). Furthermore, a noteworthy inverse relationship was observed between CmFigla expression and the methylation of its promoter in the adult gonads. Gonads at 90 dpf were used for in situ hybridization (ISH), and CmFigla transcripts were mainly concentrated in oogonia and the primary oocytes in ovaries, but undetectable in the testes. These results indicated that Figla would play vital roles in the ovarian development in C. maculata. Additionally, the frame-shift mutations of CmFigla were successfully constructed through the CRISPR/Cas9 system, which established a positive foundation for further investigation on the role of Figla in the ovarian development of C. maculata. Our study provides valuable clues for exploring the regulatory mechanism of Figla in the fish ovarian development and maintenance, which would be useful for the sex control and reproduction of fish in aquaculture.

Characterization of the transcription factor Sox3 regulating the gonadal development of pearlscale angelfish (Centropyge vrolikii)

As a member of the Sox gene family, Sox3 plays a vital role in gonadal development and gametogenesis. Nevertheless, the exact expression pattern of this gene in fish is still unknown. Here, we identified the Sox3 gene of Centropyge vrolikii, namely, Cv-Sox3. The Cv-Sox3 mRNA expression in the ovary and testis was detected by reverse transcription-polymerase chain reaction (RT-PCR) analysis, and the mRNA expression level of Cv-Sox3 in the ovary in the resting stage was significantly higher than that in other tissues. The phylogenetic tree and alignment of multiple sequences were constructed to analyze the evolutionary relationships of Cv-Sox3. Cv-Sox3 was relatively conserved in the evolution of teleost fish, indicating the importance and similarity of its function. The in situ hybridization results demonstrate that Cv-Sox3 was present in the follicle cells and cytoplasm of oocytes in the ovary of different stages, and the positive signals occurred in germ cells of the testis. After interfering with Cv-Sox3, the growth rate of ovarian cells in culture became slow, and the expression of ovary-bias-related genes Cyp19a and Foxl2 significantly increased. Meanwhile, the expression of testis-bias-related genes Dmrt1, Sox9, Cyp11a, Amh, and Sox8 significantly decreased. These results suggest that Cv-Sox3 gene might be expressed in the germ cells of male and female gonads during gonadal development. This study provides a precise expression pattern of Cv-Sox3 and demonstrates that Cv-Sox3 might play a significant role in the reproductive regulation of C. vrolikii. In this study, Sox3 of C. vrolikii (Cv-Sox3) was cloned to understand the expression pattern in the gonadal development, which is expressed in germ cells, involved in the process of gonadal development. The results demonstrated that Cv-Sox3 may play a significant role in the reproductive regulation of C. vrolikii.

Comparative transcriptome provides insights into differentially expressed genes between testis and ovary of Onychostoma macrolepis in reproduction period

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered fishery species inhabiting the river of Qinling Mountains and some flowing freshwaters in China. The declining population of O. macrolepis caused by asynchrony of male and female development prompted us to focus on genetic regulation of its reproduction. In this study, high-throughput RNA-sequencing technology was applied to assemble and annotate the transcriptome of O. macrolepis testis and ovary. The results showed that a number of 338089335 (ovary:163216500, testis:174872835) raw sequences were obtained. After non-redundant analysis, a number of 207826065 (ovary:102334008, testis:105492057) high quality reads were obtained and predicted as unigenes, in which 201,038,682 unigenes were annotated with multiple databases. Taking the ovarian transcriptome as a control, comparative transcriptome analysis showed that 9918 differentially expressed genes (DEGs) up-regulated in the testis and 13,095 DEGs down-regulated. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, gonadal development, steroid biosynthesis pathways, MAPK signaling pathway and Wnt signaling pathway. Finally, the expression patterns of 19 unigenes were validated by using quantitative real-time polymerase chain reaction (qRT-PCR). This study illustrates a potential molecular mechanism on the unsynchronized male and female development of the O. macrolepis during the reproduction period in June and provides a theoretical basis for future artificial reproduction.

Tandem Mass Tag-Based Quantitative Proteomics Analysis of Gonads Reveals New Insight into Sexual Reversal Mechanism in Chinese Soft-Shelled Turtles

Chinese soft-shelled turtles display obvious sex dimorphism. The exogenous application of hormones (estradiol and methyltestosterone) can change the direction of gonadal differentiation of P. sinensis to produce sex reversed individuals. However, the molecular mechanism remains unclear. In this study, TMT-based quantitative proteomics analysis of four types of P. sinensis (female, male, pseudo-female, and pseudo-male) gonads were compared. Quantitative analysis of 6107 labeled proteins in the four types of P. sinensis gonads was performed. We identified 440 downregulated and 423 upregulated proteins between pseudo-females and males, as well as 394 downregulated and 959 upregulated proteins between pseudo-males and females. In the two comparisons, the differentially expressed proteins, including K7FKG1, K7GIQ2, COL4A6, K7F2U2, and K7FF80, were enriched in some important pathways, such as focal adhesion, endocytosis, apoptosis, extracellular matrix-receptor interaction, and the regulation of actin cytoskeleton, which were upregulated in pseudo-female vs. male and downregulated in pseudo-male vs. female. In pathways such as ribosome and spliceosome, the levels of RPL28, SRSF3, SNRNP40, and HNRNPK were increased from male to pseudo-female, while they decreased from female to pseudo-male. All differentially expressed proteins after sexual reversal were divided into six clusters, according to their altered levels in the four types of P. sinensis, and associated with cellular processes, such as embryonic development and catabolic process, that were closely related to sexual reversal. These data will provide clues for the sexual reversal mechanism in P. sinensis.

Morphological and Molecular Gonadal Sex Differentiation in the Wild Japanese eel Anguilla japonica

Most cultured Japanese eels (Anguilla japonica) show male sex differentiation; however, natural gonadal sex differentiation has not been evaluated. In this study, this process was characterized in wild eels. Differentiated ovaries and testes were observed after the eels grew to 320 and 300 mm in total length, respectively. The youngest ovary and testis appeared at 3 and 4 years old, respectively; however, undifferentiated gonads were found up to 7 years, suggesting that sex differentiation was triggered by growth rather than aging. gsdf, amh, foxl2b and foxl3b were highly expressed in the testes, whereas figla, sox3, foxn5, zar1, and zp3 were highly expressed in the ovaries. The expression of cyp19a1a and foxl2a did not differ significantly between the testis and ovary. In the ovaries, the cyp19a1a and foxl2a levels were highest in the early stages, suggesting that their function is limited to early ovarian differentiation. The foxn5, zar1 and zp3 levels tended to increase in the later stages, suggesting that they function after the initiation of ovarian differentiation. In undifferentiated gonads, dimorphic gene expression was not observed, suggesting that the molecular sex differentiation phase is short and difficult to detect. These findings provide the first demonstration of the whole course of natural gonadal sex differentiation in eels at molecular and morphological levels.

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

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

Systematic identification and expression analysis of the Sox gene family in spotted sea bass (Lateolabrax maculatus)

The Sox gene family encodes a set of transcription factors characterized by a conserved Sry-related high mobility group (HMG)-box domain, which performs a series of essential biological functions in diverse tissues and developmental processes. In this study, the Sox gene family was systematically characterized in spotted sea bass (Lateolabrax maculatus). A total of 26 Sox genes were identified and classified into eight subfamilies, namely, SoxB1, SoxB2, SoxC, SoxD, SoxE, SoxF, SoxH and SoxK. The phylogenetic relationship, exon-intron and domain structure analyses supported their annotation and classification. Comparison of gene copy numbers and chromosome locations among different species indicated that except tandem duplicated paralogs of Sox17/Sox32, duplicated Sox genes in spotted sea bass were generated from teleost-specific whole genome duplication during evolution. In addition, qRT-PCR was performed to detect the expression profiles of Sox genes during development and adulthood. The results showed that the expression of 16 out of 26 Sox genes was induced dramatically at different starting points after the multicellular stage, which is consistent with embryogenesis. At the early stage of sex differentiation, 9 Sox genes exhibited sexually dimorphic expression patterns, among which Sox3, Sox19 and Sox6b showed the most significant ovary-biased expression. Moreover, the distinct expression pattern of Sox genes was observed in different adult tissues. Our results provide a fundamental resource for further investigating the functions of Sox genes in embryonic processes, sex determination and differentiation as well as controlling the homeostasis of adult tissues in spotted sea bass.

Transcriptional changes of mouse ovary during follicle initial or cyclic recruitment mediated by extra hormone treatment

AIMS

Folliculogenesis contains gonadotropin-independent and -dependent stage. Disruption in any of this process would induce failure in retrieving capable oocytes during clinical treatment. However, there is still limited understanding of the molecular components specifically regulating this process.

MATERIAL AND METHODS

Ovaries of P3, P20 and exogenous gonadotropin-treated P22 mice were sampled and underwent RNA-seq to investigate the transcriptome variance during mouse folliculogenesis.

KEY FINDINGS

In our dataset, 1883 and 626 DEGs were captured for each stage respectively, which were further clustered into eight expression patterns. Pathway enrichment analysis identified distinct biological processes enriched in two stages, with the most prominent being the pathways related to metabolism, gene expression, cell cycle, immune system and DNA methylation. Transcriptional regulator inference yielded eight master transcription factors (i.e. Runx1, Stat3, Sox3, Pou5f1, Gata4, Foxl2, Cebpb, and Esr1) driving folliculogenesis.

SIGNIFICANCE

Our study revealed the temporal transcriptional reprogramming and gene expression dynamics during folliculogenesis mediated by extra hormone treatment, which could provide novel insights to controlled ovarian stimulation in future infertility treatment.

Profiles of 5α-Reduced Androgens in Humans and Eels: 5α-Dihydrotestosterone and 11-Ketodihydrotestosterone Are Active Androgens Produced in Eel Gonads

Although 11-ketotestosterone (11KT) and testosterone (T) are major androgens in both teleosts and humans, their 5α-reduced derivatives produced by steroid 5α-reductase (SRD5A/srd5a), i.e., 11-ketodihydrotestosterone (11KDHT) and 5α-dihydrotestosterone (DHT), remains poorly characterized, especially in teleosts. In this study, we compared the presence and production of DHT and 11KDHT in Japanese eels and humans. Plasma 11KT concentrations were similar in both male and female eels, whereas T levels were much higher in females. In accordance with the levels of their precursors, 11KDHT levels did not show sexual dimorphism, whereas DHT levels were much higher in females. It is noteworthy that plasma DHT levels in female eels were higher than those in men. In addition, plasma 11KDHT was undetectable in both sexes in humans, despite the presence of 11KT. Three srd5a genes (srd5a1, srd5a2a and srd5a2b) were cloned from eel gonads. All three srd5a genes were expressed in the ovary, whereas only both srd5a2 genes were expressed in the testis. Human SRD5A1 was expressed in testis, ovary and adrenal, whereas SRD5A2 was expressed only in testis. Human SRD5A1, SRD5A2 and both eel srd5a2 isoforms catalyzed the conversion of T and 11KT into DHT and 11KDHT, respectively, whereas only eel srd5a1 converted T into DHT. DHT and 11KDHT activated eel androgen receptor (ar)α-mediated transactivation as similar fashion to T and 11KT. In contrast, human AR and eel arβ were activated by DHT and11KDHT more strongly than T and 11KT. These results indicate that in teleosts, DHT and 11KDHT may be important 5α-reduced androgens produced in the gonads. In contrast, DHT is the only major 5α-reduced androgens in healthy humans.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

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.

Dmrt1 (doublesex and mab-3-related transcription factor 1) expression during gonadal development and spermatogenesis in the Japanese eel

Dmrt1, doublesex- and mab-3-related transcription factor-1, has been suggested to play critical roles in male gonadogenesis, testicular differentiation and development, including spermatogenesis, among different vertebrates. Vasa is a putative molecular marker of germ cells in vertebrates. In this study, we cloned the full-length dmrt1 cDNA from Japanese eel, and the protein comprised 290 amino acids and presented an extremely conserved Doublesex and Mab-3 (DM) domain. Vasa proteins were expressed in gonadal germ cells in a stage-specific manner, and were expressed at high levels in PGC and spermatogonia, low levels in spermatocytes, and were absent in spermatids and spermatozoa of Japanese eels. Dmrt1 proteins were abundantly expressed in spermatogonia B cells, spermatocytes, spermatids, but not in spermatozoa, spermatogonia A and Sertoli cells. To our knowledge, this study is the first to show a restricted expression pattern for the Dmrt1 protein in spermatogonia B cells, but not spermatogonia A cells, of teleosts. Therefore, Dmrt1 might play vital roles at the specific stages during spermatogenesis from spermatogonia B cells to spermatids in the Japanese eel. Moreover, the Dmrt1 protein exhibited a restricted localization in differentiating oogonia in the early differentiating gonad (ovary-like structure) of male Japanese eels and in E2-induced feminized Japanese eels. We proposed that dmrt1 may be not only required for spermatogenesis but might also play a role in oogenesis in the Japanese eel.

Medroxyprogesterone acetate affects sex differentiation and spermatogenesis in zebrafish

Medroxyprogesterone acetate (MPA) is a widely used synthetic progestin and it has been frequently detected in aquatic environments. However, its effects on aquatic organisms remain largely unknown. Here we investigated the chronic effects of MPA on sex differentiation and gonad development in zebrafish. Zebrafish larvae at 20 days post fertilization (dpf) were exposed to 4.32, 42.0, and 424 ng L^-1 of MPA until they reached 140 dpf. The results showed that chronic exposure to 42.0 ng L^-1 of MPA caused 60% proportion of males as well as significant up-regulation of dmrt1 (˜1.79 fold) and hsd17b3 (˜1.92 fold). Histological analysis showed MPA significantly increased the frequency of immature spermatocytes accompanied with the increased transcription of dmrt1 (˜2.06 fold) and ar (˜1.73 fold) in the testes. Meanwhile, MPA exposure significantly increased the transcription of lhb at all exposure concentrations in the males. In contrast, it significantly suppressed the transcription of lhb (˜-8.06-fold) and fshb (˜-6.35-fold) at 42.0 ng L^-1 in the females. Collectively our results demonstrated that MPA had androgenic activity, and could affect sex differentiation and spermatogenesis in zebrafish at environmentally relevant concentrations. The findings from this study suggest that MPA in the aquatic environment may pose potential androgenic risks to fish populations.

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.

Characterization of the European Sea Bass (Dicentrarchus labrax) Gonadal Transcriptome During Sexual Development

The European sea bass is one of the most important cultured fish in Europe and has a marked sexual growth dimorphism in favor of females. It is a gonochoristic species with polygenic sex determination, where a combination between still undifferentiated genetic factors and environmental temperature determines sex ratios. The molecular mechanisms responsible for gonadal sex differentiation are still unknown. Here, we sampled fish during the gonadal developmental period (110 to 350 days post fertilization, dpf), and performed a comprehensive transcriptomic study by using a species-specific microarray. This analysis uncovered sex-specific gonadal transcriptomic profiles at each stage of development, identifying larger number of differentially expressed genes in ovaries when compared to testis. The expression patterns of 54 reproduction-related genes were analyzed. We found that hsd17β10 is a reliable marker of early ovarian differentiation. Further, three genes, pdgfb, snx1, and nfy, not previously related to fish sex differentiation, were tightly associated with testis development in the sea bass. Regarding signaling pathways, lysine degradation, bladder cancer, and NOD-like receptor signaling were enriched for ovarian development while eight pathways including basal transcription factors and steroid biosynthesis were enriched for testis development. Analysis of the transcription factor abundance showed an earlier increase in females than in males. Our results show that, although many players in the sex differentiation pathways are conserved among species, there are peculiarities in gene expression worth exploring. The genes identified in this study illustrate the diversity of players involved in fish sex differentiation and can become potential biomarkers for the management of sex ratios in the European sea bass and perhaps other cultured species.

Evolutionarily conserved vasa identifies embryonic and gonadal germ cells in spinyhead croaker Collichthys lucidus

In this study, a 2198 bp full-length cDNA of spinyhead croaker Collichthys lucidus vasa gene encoding 616 amino-acid residues was obtained. Multiple alignment revealed that C. lucidus vasa has eight conserved characteristic motifs of the DEAD box protein family and has the highest identity to large yellow croaker Larimichthys croceas. Reverse-transcription (RT)-PCR and Western blot analyses indicated that the vasa messenger (m)RNA and Vasa protein are specifically expressed in the gonads in both sexes. In situ hybridisation (ISH) demonstrated that vasa RNA is exclusively detected in the germ cells in C. lucidus gonads and its temporospatial expression reveals a dynamic pattern during oogenesis. Surprisingly, C. lucidus vasa 3'UTR can direct stable and specific GFP expression in the primordial germ cells (PGC) of medaka Oryzias latipes embryos. Taken together, these results suggest that because C. lucidus vasa expression delineates critical stages of oogenesis, it may be a useful molecular marker for the identification of gonadal germ cells, facilitating the isolation and utilization of germ cells in future study.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.

Loss-of-function of sox3 causes follicle development retardation and reduces fecundity in zebrafish

Folliculogenesis is essential for production of female gametes in vertebrates. However, the molecular mechanisms underlying follicle development, particularly apoptosis regulation in ovary, remain elusive. Here, we generated sox3 knockout zebrafish lines using CRISPR/Cas9. sox3 knockout led to follicle development retardation and a reduced fecundity in females. Comparative analysis of transcriptome between sox3^−/− and wild-type ovaries revealed that Sox3 was involved in pathways of ovarian steroidogenesis and apoptosis. Knockout of sox3 promoted follicle apoptosis and obvious apoptosis signals were detected in somatic cells of stages III and IV follicles of sox3^−/− ovaries. Moreover, Sox3 can bind to and activate the promoter of cyp19a1a. Up-regulation of Cyp19a1a expression promoted 17β-estradiol synthesis, which inhibited apoptosis in follicle development. Thus, Sox3 functions as a regulator of Cyp19a1a expression, via 17β-E2 linking apoptosis suppression, which is implicated in improving female fecundity.