Identification of gonadal soma-derived factor involvement in Monopterus albus (protogynous rice field eel) sex change

Meta-analysis of gonadal transcriptome provides novel insights into sex change mechanism across protogynous fishes

Protogyny, being capable of changing from female to male during their lifetime, is prevalent in 20 families of teleosts but is believed to have evolved within specific evolutionary lineages. Therefore, shared regulatory factors governing the sex change process are expected to be conserved across protogynous fishes. However, a comprehensive understanding of this mechanism remains elusive. To identify these factors, we conducted a meta-analysis using gonadal transcriptome data from seven species. We curated data pairs of ovarian tissue and transitional gonad, and employed ratios of expression level as a unified criterion for differential expression, enabling a meta-analysis across species. Our approach revealed that classical sex change-related genes exhibited differential expression levels between the ovary and transitional gonads, consistent with previous reports. These results validate our methodology's robustness. Additionally, we identified novel genes not previously linked to gonadal sex change in fish. Notably, changes in the expression levels of acetoacetyl-CoA synthetase and apolipoprotein Eb, which are involved in cholesterol synthesis and transport, respectively, suggest that the levels of cholesterol, a precursor of steroid hormones crucial for sex change, are decreased upon sex change onset in the gonads. This implies a potential universal influence of cholesterol dynamics on gonadal transformation in protogyny.

Characterization and Functional Analysis of the 17-Beta Hydroxysteroid Dehydrogenase 2 (hsd17b2) Gene during Sex Reversal in the Ricefield Eel (Monopterus albus)

In mammals, 17-beta hydroxysteroid dehydrogenase 2 (Hsd17b2) enzyme specifically catalyzes the oxidation of the C17 hydroxyl group and efficiently regulates the activities of estrogens and androgens to prevent diseases induced by hormone disorders. However, the functions of the hsd17b2 gene involved in animal sex differentiation are still largely unclear. The ricefield eel (Monopterus albus), a protogynous hermaphroditic fish with a small genome size (2n = 24), is usually used as an ideal model to study the mechanism of sex differentiation in vertebrates. Therefore, in this study, hsd17b2 gene cDNA was cloned and its mRNA expression profiles were determined in the ricefield eel. The cloned cDNA fragment of hsd17b2 was 1230 bp, including an open reading frame of 1107 bp, encoding 368 amino acid residues with conserved catalytic subunits. Moreover, real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) analysis showed that hsd17b2 mRNA expressed strongly in the ovaries at early developmental stages, weakly in liver and intestine, and barely in testis and other tissues. In particular, hsd17b2 mRNA expression was found to peak in ovaries of young fish and ovotestis at the early stage, and eventually declined in gonads from the late ovotestis to testis. Likewise, chemical in situ hybridization results indicated that the hsd17b2 mRNA signals were primarily detected in the cytoplasm of oogonia and oocytes at stage I-II, subsequently concentrated in the granulosa cells around the oocytes at stage Ⅲ-Ⅳ, but undetectable in mature oocytes and male germ cells. Intriguingly, in ricefield eel ovaries, hsd17b2 mRNA expression could be significantly reduced by 17β-estradiol (E2) or tamoxifen (17β-estradiol inhibitor, E2I) induction at a low concentration (10 ng/mL) and increased by E2I induction at a high concentration (100 ng/mL). On the other hand, both the melatonin (MT) and flutamide (androgen inhibitor, AI) induction could significantly decrease hsd17b2 mRNA expression in the ovary of ricefield eel. This study provides a clue for demonstrating the mechanism of sexual differentiation in fish. The findings of our study imply that the hsd17b2 gene could be a key regulator in sexual differentiation and modulate sex reversal in the ricefield eel and other hermaphroditic fishes.

Smad4 and FoxH1 potentially interact to regulate cyp19a1a promoter in the ovary of ricefield eel (Monopterus albus)

BACKGROUND

Cyp19a1a is a key enzyme in the pathway that converts androgens into estrogen and is regulated by TGF-β signaling. Smad4 and FoxH1 are downstream effectors of TGF-β signaling and may play important roles in ovarian development in M. albus.

METHODS

We investigated the expression pattern of the Smad4 and FoxH1 using qRT‒PCR and immunofluorescence, then tested the changes of smad4 and foxh1 by qRT‒PCR after ovary incubation with FSH in vitro, and analysed the regulation of cyp19a1a transcription by Smad4 and FoxH1 by dual-luciferase reporter assays.

RESULTS

We found that Smad4 encoded a putative protein of 449 amino acids and harbored the three conserved domains typical of this protein family. Smad4 and foxh1 exhibited similar expression patterns during ovarian development and after FSH incubation, with Pearson's coefficients of 0.873 and 0.63-0.81, respectively. Furthermore, Smad4, FoxH1 and Cyp19a1a colocalized in the granulosa cells and theca cells of ovaries during the mid-to-late vitellogenic stage. Smad4 repressed cyp19a1a activity via SBE1 (- 1372/-1364) and SBE2 (- 415/-407) in the cyp19a1a promoter, whereas mutating SBE1 or SBE2 restored cyp19a1a promoter activity. Co-overexpression of Smad4 and FoxH1 significantly reduced cyp19a1a promoter activity.

CONCLUSIONS

This study provides new insights into the potential functions of transcription factors Smad4 and FoxH1 in ovarian development and the transcriptional regulation mechanism of cyp19a1a in M. albus, which will reveal Smad4/FoxH1-mediated TGF-β signaling in reproduction and the regulation of the cyp19a1a. Aromatase, encoded by cyp19a1a, is involved in ovarian development and plays an important role in the quality of eggs, as well the sex ratio, of the teleost fish, M. albus. The research on the transcriptional regulation of cyp19a1a has contributed to the understanding of its role in ovarian development. In previous study, it was shown that FoxH1 inhibits cyp19a1a transcription. In the present study, Smad4 was confirmed as a cyp19a1a transcriptional repressor and Smad4 may also coordinate with FoxH1 to repress cyp19a1a transcription. At present, we provide a new perspective for the transcriptional regulation of cyp19a1a by transcription factors Smad4 and FoxH1 in teleost fish ovary. In the future, the regulatory networks of Smad4 and FoxH1 will be further studied and the gene editing technology will be applied to screen specific regulatory factors of cyp191a1a gene, so as to alter the female cycle and modulate the sex ratio of the eggs production.

Genome-wide identification, evolution and expression of TGF-β signaling pathway members in mandarin fish (Siniperca chuatsi)

Members of the transforming growth factor β (TGF-β) signaling pathway regulate diverse cellular biological processes in embryonic and tissue development. We took mandarin fish (Siniperca chuatsi) as the research object to identify all members of the TGF-β signaling pathway, measure their expression pattern in the key period post hatching, and further explore their possible role in the process of sex regulation. Herein, we identified eighty-three TGF-β signaling pathway members and located them on chromosomes based on the genome of mandarin fish. TGF-β signaling pathway members were highly conserved since each TGF-β subfamily clustered with orthologs from other species. Transcriptome analysis, qRT-PCR and in situ hybridization demonstrated that most mandarin fish TGF-β signaling pathway members presented stage-specific and/or sex-dimorphic expression during gonadal development, and different members of the TGF-β signaling pathway participated in different stages of gonadal development. Taken together, our results provide new insight into the role of TGF-β signaling pathway members in the sex regulation of mandarin fish.

Bisphenol A alters sexual dimorphism and gene expression in marine medaka Oryzias melastigma

Bisphenol A (BPA) is an endocrine disruptor that is present in freshwater and marine environments. However, conclusive evidence for the toxicity of chronic BPA exposure to marine fishes remains lacking. Therefore, we investigated the influence of BPA on male marine medaka (Oryzias melastigma). BPA exposure induced formation of testis-ova at 2610 µg/L, and male-type anal fins became more female type in a concentration-dependent manner. Some males with female-type anal fins had normal testes, indicating that anal fin shape is more sensitive to BPA. Gonadal soma-derived factor (gsdf) expression decreased after BPA exposure in the 746 and 2610 µg/L exposure groups, although the changes were not statistically significant. Additionally, liver vitellogenin (vtg) expression increased in a dose-dependent manner and was significantly higher in all exposure groups. vtg and gsdf are likely to be useful biomarkers for the impact of estrogenic endocrine disrupters in O. melastigma.

The ontogenic gonadal transcriptomes provide insights into sex change in the ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a freshwater protogynous hermaphrodite fish and has become an important aquaculture species in China. The sex change of ricefield eel is impeding its aquaculture practice, particularly the large-scale artificial breeding. Many studies including transcriptomes of mixed gonadal samples from different individuals have been aimed to elucidate mechanisms underlying the sex change. However, the key physiological factors involved in the initiation of sex change remain to be identified.  RESULTS: The present study performed transcriptomic analysis on gonadal samples of different sexual stages obtained through biopsy from the same fish undergoing sex change. A total of 539,764,816 high-quality reads were generated from twelve cDNA libraries of gonadal tissues at female (F), early intersexual (EI), mid-intersexual (MI), and late intersexual (LI) stages of three individual sex-changing fish. Pairwise comparisons between EI and F, MI and EI, and LI and MI identified 886, 319, and 10,767 differentially expressed genes (DEGs), respectively. Realtime quantitative PCR analysis of 12 representative DEGs showed similar expression profiles to those inferred from transcriptome data, suggesting the reliability of RNA-seq data for gene expression analysis. The expression of apoeb, csl2, and enpp2 was dramatically increased and peaked at EI while that of cyp19a1a, wnt4a, fgf16, and foxl2a significantly downregulated from F to EI and remained at very low levels during subsequent development until LI, which suggests that apoeb, csl2, enpp2, cyp19a1a, wnt4a, fgf16, and foxl2a may be closely associated with the initiation of sex change of ricefield eels.

CONCLUSIONS

Collectively, results of the present study confirmed that the down-regulation of female-related genes, such as cyp19a1a, wnt4a, fgf16, and foxl2a, is important for the sex change of ricefield eels. More importantly, some novel genes, including apoeb, csl2, and enpp2, were shown to be expressed with peak values at EI, which are potentially involved in the initiation of sex change. The present transcriptomic data may provide an important research resource for further unraveling the mechanisms underlying the sex change and testicular development in ricefield eels as well as other teleosts.

Molecular Cloning and Expression Analysis of Thyrotropin-Releasing Hormone, and Its Possible Role in Gonadal Differentiation in Rice Field eel Monopterus albus

Rice field eel (Monopterus albus), a protogynous hermaphrodite fish, is a good model for the research of sex determination and gonadal differentiation in teleosts. In this study, we cloned the full-length cDNA sequence of trh, which encoded a predicted protein with 270 amino acids. Trh mainly expressed in the brain, followed by the ovary, testis, muscle and pituitary, and had low levels in other peripheral tissues. During natural sex reversal, trh mRNA expression levels exhibited a significant increase at the late intersexual stage in the hypothalamus. In the gonad, trh mRNA expression levels showed a trend of increase followed by decrease, and only increased significantly at the middle intersexual stage. No matter static incubation or intraperitoneal (IP) injection, TRH had no significant effect on trh and thyroid-stimulating hormone βsubunit (tshβ) mRNA expression levels, and serum T3, T4 and TRH release. After static incubation of ovarian fragments by TRH, the expression of gonadal soma derived factor (gsdf) was up-regulated significantly at both the doses of 10 and 100 nM. IP injection of TRH stimulated the expression of gsdf, and inhibited the expression of ovarian aromatase gene (cyp19a1a), accompanied by the increase of serum 11-KT levels. The results indicated that TRH may play a novel role in gonadal differentiation by the regulation of gonadal differentiation-related gene expression and sex steroid hormone secretion in rice field eel.

Chromosome-Scale, Haplotype-Resolved Genome Assembly of Non-Sex-Reversal Females of Swamp Eel Using High-Fidelity Long Reads and Hi-C Data

The Asian swamp eel (Monopterus albus) is an excellent model species for studying sex change and chromosome evolution. M. albus is also widely reared in East Asia and South-East Asia because of its great nutritional value. The low fecundity of this species (about 300 eggs per fish) greatly hinders fries production and breeding programs. Interestingly, about 3–5% of the eels could remain as females for 3 years and lay more than 3,000 eggs per fish, which are referred to as non-sex-reversal (NSR) females. Here, we presented a new chromosome-level genome assembly of such NSR females using Illumina, HiFi, and Hi-C sequencing technologies. The new assembly (Mal.V2_NSR) is 838.39 Mb in length, and the N50 of the contigs is 49.8 Mb. Compared with the previous assembly obtained using the continuous long-read sequencing technology (Mal.V1_CLR), we found a remarkable increase of continuity in the new assembly Mal.V2_NSR with a 20-times longer contig N50. Chromosomes 2 and 12 were assembled into a single contig, respectively. Meanwhile, two highly contiguous haplotype assemblies were also obtained, with contig N50 being 14.54 and 12.13 Mb, respectively. BUSCO and Merqury analyses indicate completeness and high accuracy of these three assemblies. A comparative genomic analysis revealed substantial structural variations (SVs) between Mal.V2_NSR and Mal.V1_CLR and two phased haplotype assemblies, as well as whole chromosome fusion events when compared with the zig-zag eel. Additionally, our newly obtained assembly provides a genomic view of sex-related genes and a complete landscape of the MHC genes. Therefore, these high-quality genome assemblies would provide great help for future breeding works of the swamp eel, and it is a valuable new reference for genetic and genomic studies of this species.

Widespread alterations to hypothalamic-pituitary-gonadal (HPG) axis signaling underlie high temperature reproductive inhibition in the eurythermal sheepshead minnow (Cyprinodon variegatus)

Fish experiencing abnormally high or prolonged elevations in temperature can exhibit impaired reproduction, even for species adapted to warm water environments. Such high temperature inhibition of reproduction has been linked to diminished gonadal steroidogenesis, but the mechanisms whereby hypothalamic-pituitary-gonadal (HPG) axis signaling is impacted by high temperature are not fully understood. Here, we characterized differences in HPG status in adult sheepshead minnow (Cyprinodon variegatus), a eurythermal salt marsh and estuarine species of eastern North America, exposed for 14 d to temperatures of 27 °C or 37 °C. Males and females at 37 °C had lower gonadosomatic index (GSI) values compared to fish at 27 °C, and females at 37 °C had fewer spawning capable eggs and lower circulating 17β-estradiol (E₂). Gene transcripts encoding gonadotropin-inhibitory hormone (gnih) and gonadotropin-releasing hormone-3 (gnrh3) were higher in relative abundance in the hypothalamus of both sexes at 37 °C. While pituitary mRNAs for the β-subunits of follicle-stimulating hormone (fshβ) and luteinizing hormone (lhβ) were lowered only in males at 37 °C, Fsh and Lh receptor mRNA levels in the gonads were at lower relative levels in both the ovary and testis of fish at 37 °C. Females at 37 °C also showed reduced ovarian mRNA levels for steroid acute regulatory protein (star), P450 side-chain cleavage enzyme (cyp11a1), 3β-hydroxysteroid dehydrogenase (3βhsd), 17β-hydroxysteroid dehydrogenase (hsd17β3), and ovarian aromatase (cyp19a1a). Females at the higher 37 °C temperature also had a lower liver expression of mRNAs encoding estrogen receptor α (esr1) and several vitellogenin and choriogenin genes, but elevated mRNA levels for hepatic sex hormone-binding globulin (shbg). Our results substantiate prior findings that exposure of fish to high temperature can inhibit gonadal steroidogenesis and oogenesis, and point to declines in reproductive performance emerging from alterations at several levels of HPG axis signaling including increased hypothalamic Gnih expression, depressed gonadal steroidogenesis, and reduced egg yolk and egg envelope protein production in the liver.

Evolution, Expression, and Function of Gonadal Somatic Cell-Derived Factor

Fish gonads develop in very diverse ways different from mammalian gonads. This diversity is contributed by species-specific factors. Gonadal somatic cell-derived factor (Gsdf) is one such factor. The gsdf gene exists mostly in teleosts and is absent in many tetrapods, probably as a result of two gene losses during evolution. The gsdf transcript is expressed mainly in gonadal somatic cells, including Sertoli cell in testis and granulosa cells in ovary; however, these gonadal somatic cells can surround many types of germ cells at different developmental stages depending on the fish species. The function of gsdf is also variable. It is involved in germ cell proliferation, testicular formation, ovarian development and even male sex determination. Here, we summarize the common and diverse expression, regulation and functions of gsdf among different fish species with aspect of evolution.

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.

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.

Metabonomic Insights into the Sperm Activation Mechanisms in Ricefield Eel (Monopterus albus)

In fish, sperm motility activation is one of the most essential procedures for fertilization. Previous studies have mainly focused on the external environmental effects and intracellular signals in sperm activation; however, little is known about the metabolic process of sperm motility activation in fish. In the present study, using ricefield eel (Monopterus albus) sperm as a model, metabonomics was used to analyze the metabolic mechanism of the sperm motility activation in fish. Firstly, 529 metabolites were identified in the sperm of ricefield eel, which were clustered into the organic acids, amino acids, nucleotides, benzene, and carbohydrates, respectively. Among them, the most abundant metabolites in sperm were L-phenylalanine, DL-leucine, L-leucine, lysolecithin choline 18:0, L-tryptophan, adenine, hypoxanthine, 7-Methylguanine, shikimic acid, and L-tyrosine. Secondly, compared to pre-activated sperm, the level of S-sulfo-L-cysteine and L-asparagine were both increased in the post-activated sperm. Ninety-two metabolites were decreased in the post-activated sperm, including quinic acid, acetylsalicylic acid, 7,8-dihydro L-biopterin, citric acid, glycylphenylalanine, and dihydrotachysterol (DHT). Finally, basing on the pathway analysis, we found that the changed metabolites in sperm motility activation were mainly clustered into energy metabolism and anti-oxidative stress. Fish sperm motility activation would be accompanied by the release of a large amount of energy, which might damage the genetic material of sperm. Thus, the anti-oxidative stress function is a critical process to maintain the normal physiological function of sperm.

Gonadal Transcriptome Analysis of Sex-Related Genes in the Protandrous Yellowfin Seabream (Acanthopagrus latus)

Yellowfin seabream (Acanthopagrus latus), a protandrous hermaphroditic fish, is a good model for studying the mechanism of sex reversal. However, limited knowledge is known about the genetic information related to reproduction and sex differentiation in this species. Here, we performed de novo transcriptome sequencing analysis of the testis, ovotestis, and ovary to identify sex-related genes in yellowfin seabream. The results assembled 71,765 unigenes in which 16,126 and 17,560 unigenes were differentially expressed in the ovotestis and ovary compared to the testis, respectively. The most differentially expressed gene (DEG)-enriched Kyoto Encyclopedia of Genes and Genomes and GO pathways were closely associated with the synthesis of sex steroid hormones. Functional analyses identified 55 important sex-related DEGs, including 32 testis-biased DEGs (dmrt1, amh, and sox9, etc.), 20 ovary-biased DEGs (cyp19a, foxl2, and wnt4, etc.), and 3 ovotestis-biased DEGs (lhb, dmrt2, and foxh1). Furthermore, the testis-specific expression of dmrt1 and the brain-pituitary-ovary axis expression of foxl2 were characterized, suggesting that they might play important roles in sex differentiation in yellowfin seabream. Our present work provided an important molecular basis for elucidating the mechanisms underlying sexual transition and reproductional regulation in yellowfin seabream.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Transcriptomic analysis of dead end knockout testis reveals germ cell and gonadal somatic factors in Atlantic salmon

Background

Sustainability challenges are currently hampering an increase in salmon production. Using sterile salmon can solve problems with precocious puberty and genetic introgression from farmed escapees to wild populations. Recently sterile salmon was produced by knocking out the germ cell-specific dead end (dnd). Several approaches may be applied to inhibit Dnd function, including gene knockout, knockdown or immunization. Since it is challenging to develop a successful treatment against a gene product already existing in the body, alternative targets are being explored. Germ cells are surrounded by, and dependent on, gonadal somatic cells. Targeting genes essential for the survival of gonadal somatic cells may be good alternative targets for sterility treatments. Our aim was to identify and characterize novel germ cell and gonadal somatic factors in Atlantic salmon.

Results

We have for the first time analysed RNA-sequencing data from germ cell-free (GCF)/dnd knockout and wild type (WT) salmon testis and searched for genes preferentially expressed in either germ cells or gonadal somatic cells. To exclude genes with extra-gonadal expression, our dataset was merged with available multi-tissue transcriptome data. We identified 389 gonad specific genes, of which 194 were preferentially expressed within germ cells, and 11 were confined to gonadal somatic cells. Interestingly, 5 of the 11 gonadal somatic transcripts represented genes encoding secreted TGF-β factors; gsdf, inha, nodal and two bmp6-like genes, all representative vaccine targets. Of these, gsdf and inha had the highest transcript levels. Expression of gsdf and inha was further confirmed to be gonad specific, and their spatial expression was restricted to granulosa and Sertoli cells of the ovary and testis, respectively. Finally, we show that inha expression increases with puberty in both ovary and testis tissue, while gsdf expression does not change or decreases during puberty in ovary and testis tissue, respectively.

Conclusions

This study contributes with transcriptome data on salmon testis tissue with and without germ cells. We provide a list of novel and known germ cell- and gonad somatic specific transcripts, and show that the expression of two highly active gonadal somatic secreted TGF-β factors, gsdf and inha, are located within granulosa and Sertoli cells.

Loss of Gsdf leads to a dysregulation of Igf2bp3-mediated oocyte development in medaka

Gonadal soma-derived factor (Gsdf) is a unique TGF-β factor essential for both ovarian and testicular development in Hd-rR medaka (Oryzias latipes). However, the downstream genes regulated by Gsdf signaling remain unknown. Using a high-throughput proteomic approach, we identified a significant increase in the expression of the RNA-binding protein Igf2bp3 in gsdf-deficient ovaries. We verified this difference in transcription and protein expression against normal gonads using real-time PCR quantification and Western blotting. The genomic structure of igf2bp3 and the syntenic flanking segments are highly conserved across fish and mammals. igf2bp3 expression was correlated with oocyte development, which is consistent with the expression of the igf2bp3 ortholog Vg1-RBP/Vera in Xenopus. In contrast to the normal ovary, cysts of H3K27me3- and Igf2bp3-positive germ cells were dramatically increased in the one-month-old gsdf-deficient ovary, indicating that the gsdf depletion led to a dysregulation of Igf2bp3-mediated oocyte development. Our results provide novel insights into the Gsdf-Igf2bp3 signaling mechanisms that underlie the fundamental process of gametogenesis; these mechanisms may be well conserved across phyla.

Expression and transcriptional regulation of gsdf in spotted scat (Scatophagus argus)

Gonadal soma-derived factor (Gsdf) is critical for testicular differentiation and early germ cell development in teleosts. The spotted scat (Scatophagus argus), with a stable XX-XY sex-determination system and the candidate sex determination gene dmrt1, provides a good model for understanding the mechanism of sex determination and differentiation in teleosts. In this study, we analyzed spotted scat gsdf tissue distribution and gene expression patterns in gonads, as well as further analysis of transcriptional regulation. Tissue distribution analysis showed that gsdf was only expressed in testis and ovary. Real-time PCR showed that both gsdf and dmrt1 were expressed significantly higher in testes at different phases (phase III, IV and V) compared to ovaries at phase II, III and IV, while gsdf was expressed significantly higher in phase II ovaries than those of phase III and IV. Western blot analysis also showed that Gsdf was more highly expressed in the testis than ovary. Immunohistochemistry analysis showed that Gsdf was expressed in Sertoli cells surrounding spermatogonia in the testis, while it was expressed in the somatic cells surrounding the oogonia of the ovary. Approximately 2.7 kb of the 5' upstream region of gsdf was cloned from the spotted scat genomic DNA and in silico promoter analysis revealed the putative transcription factor binding sites of Dmrt1 and Sf1. The luciferase reporter assay, using the human embryonic kidney cells, demonstrated that Dmrt1 activated gsdf expression in a dose-dependent manner in the presence of Sf1 in spotted scat. These results suggest that Gsdf could play a role in regulating the development of spermatogonia and oogonia, and also participate in male sex differentiation by acting as a downstream gene of Dmrt1 in spotted scat.

High temperature increases the gsdf expression in masculinization of genetically female Japanese flounder (Paralichthys olivaceus)

In teleosts, sex is plastic and is influenced by environmental factors. Elevated temperatures have masculinizing effects on the phenotypic sex of certain sensitive species. In this study, we reared genetic XX Japanese flounder at a high temperature (27.5 ± 0.5 °C) and obtained a population of sex-reversal XX males (male ratio, 95.24%). We comparatively analyzed the dynamic characteristics of germ cells and gsdf (gonadal soma-derived factor) expression during sexual differentiation for the experimental (27.5 ± 0.5 °C) and control (18 °C ± 0.5 °C) groups. The results revealed that the germ cell proliferation inhibited and gsdf expression up-regulated in the experimental group, and the gsdf mRNA and proteins expressed in somatic cells that had direct contact with germline stem cells (with Nanos 2 protein expression) including spermatogonia and oogonia by ISH (in situ hybridization) and IHC (immunohistochemistry). In addition, we also overexpressed the gsdf in XX flounders, and the germ cell number of XX flounders bearing gsdf gene significantly decreased and sometimes disappeared completely, which was consistent with the results from high-temperature induction. Therefore, based on all the results, we speculated that the high expression of gsdf might inhibit germ cell proliferation during sex differentiation, and eventually cause sex reversal in the high-temperature induced masculinization of XX Japanese flounder.

The autosomal Gsdf gene plays a role in male gonad development in Chinese tongue sole (Cynoglossus semilaevis)

Gsdf is a key gene for testicular differentiation in teleost. However, little is known about the function of Gsdf in Chinese tongue sole (Cynoglossus semilaevis). In this study, we obtained the full-length Gsdf gene (CS-Gsdf), and functional characterization revealed its potential participation during germ cell differentiation in testes. CS-Gsdf transcription was predominantly detected in gonads, while the levels in testes were significantly higher than those in ovaries. During the different developmental stages in male gonads, the mRNA level was significantly upregulated at 86 dph, and a peak appeared at 120 dph; then, the level decreased at 1 and 2 yph. In situ hybridization revealed that CS-Gsdf mRNA was mainly localized in the Sertoli cells, spermatogonia, and spermatids in mature testes. After CS-Gsdf knockdown in the male testes cell line by RNA interference, a series of sex-related genes was influenced, including several sex differentiation genes, CS-Wnt4a, CS-Cyp19a1a and CS-Star. Based on these data, we speculated that CS-Gsdf may play a positive role in germ differentiation and proliferation via influencing genes related to sex differentiation.

A Comparative View on Sex Differentiation and Gametogenesis Genes in Lungfish and Coelacanths

Gonadal sex differentiation and reproduction are the keys to the perpetuation of favorable gene combinations and positively selected traits. In vertebrates, several gonad development features that differentiate tetrapods and fishes are likely to be, at least in part, related to the water-to-land transition. The collection of information from basal sarcopterygians, coelacanths, and lungfishes, is crucial to improve our understanding of the molecular evolution of pathways involved in reproductive functions, since these organisms are generally regarded as “living fossils” and as the direct ancestors of tetrapods. Here, we report for the first time the characterization of >50 genes related to sex differentiation and gametogenesis in Latimeria menadoensis and Protopterus annectens. Although the expression profiles of most genes is consistent with the intermediate position of basal sarcopterygians between actinopterygian fish and tetrapods, their phylogenetic placement and presence/absence patterns often reveal a closer affinity to the tetrapod orthologs. On the other hand, particular genes, for example, the male gonad factor gsdf (Gonadal Soma-Derived Factor), provide examples of ancestral traits shared with actinopterygians, which disappeared in the tetrapod lineage.

Identification of genes involved in gonadal sex differentiation and the dimorphic expression pattern in Takifugu rubripes gonad at the early stage of sex differentiation

Quantifying the expression of mRNAs in the gonads at the critical stage of molecular sex differentiation stage might help to clarify the regulatory network during early sex differentiation and provide new information on the role of sex-related genes in gonadal function. In this study, transcriptomic analysis of sex-related genes expression profiles in fugu gonads at 60 and 90 days after hatching (dah) was conducted firstly, and a total of 112,504,991 clean reads, encompassing 28.35 Gb of sequences were retrieved. Twenty-three thousand eight hundred ten genes were found to be expressed in juvenile fugu gonads, and we mainly focused on the differentially expressed genes that have the potential to be involved in the gonadal sex differentiation. For 60-dah juveniles, we identified 1014 genes that were upregulated in the ovary and 1570 that were upregulated in the testis. For 90-dah juveniles, we identified 1287 genes that were upregulated in the ovary and 1500 that were upregulated in the testis. The dimorphic expression patterns of 15 genes in gonads at 30 and 40 dah were further investigate using qPCR. Cyp11b and star were expressed at higher levels in XY than in XX, while cyp11a1 and cyp19a1a were expressed at higher levels in XX than in XY at 30 dah. At 40 dah, the levels of gsdf, dmrt1, dmrt3, cyp11c1, star, and hsd3b expression were higher in XY, while the levels of foxl2, cyp19a1a, wnt9b, and foxD4 expression were higher in XX. Sox9, cyp11a1, cyp17a1, cyp17a2, and nr5a2 were expressed at similar levels in XX and XY at 40 dah. This is the first report of gonadal transcriptome of fugu at early sex differentiation stage, and our results provide an archive for further study on molecular mechanism underlying sex differentiation in this species.

Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish

BACKGROUND

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

RESULTS

Zebrafish expressed gsdf in somatic cells of bipotential gonads and expression continued in ovarian granulosa cells and testicular Sertoli cells. Homozygous gsdf knockout mutants delayed leaving the bipotential gonad state, but then became a male or a female. Mutant females ovulated a few oocytes, then became sterile, accumulating immature follicles. Female mutants stored excess lipid and down-regulated aromatase, gata4, insulin receptor, estrogen receptor, and genes for lipid metabolism, vitellogenin, and steroid biosynthesis. Mutant females contained less estrogen and more androgen than wild-types. Mutant males were fertile. Genomic analysis suggests that Gsdf, Bmp15, and Gdf9, originated as paralogs in vertebrate genome duplication events.

CONCLUSIONS

In zebrafish, gsdf regulates ovarian follicle maturation and expression of genes for steroid biosynthesis, obesity, diabetes, and female fertility, leading to ovarian and extra-ovarian phenotypes that mimic human polycystic ovarian syndrome (PCOS), suggesting a role for a related TGFB signaling molecule in the etiology of PCOS. Developmental Dynamics 246:925-945, 2017. © 2017 Wiley Periodicals, Inc.

Developmental tracing of oocyte development in gonadal soma-derived factor deficiency medaka (Oryzias latipes) using a transgenic approach

Gonadal soma-derived factor (gsdf) is reported to be a male initiator in medaka based on loss- and gain- of function via targeted disruption, or transgenic over-expression. However, little is known about how gsdf promotes undifferentiated gonad entry into male pathways or prevents entry into the female pathway. We utilized a visible folliculogenesis system with a reporter cassette of dual-color fluorescence expression to identify difference between oocyte development from wildtype and gsdf deficiency medaka. A red fluorescent protein (RFP) is driven by a major component of the synaptonemal complex (SYCP3) promoter which enables RFP expression solely in oocytes after the onset of meiosis, and a histone 2b-EGFP fused protein (H2BEGFP) under the control of an elongation factor (EF1α) promoter, wildly used as a mitotic reporter of cell cycle. This mitosis-meiosis visible switch revealed that early meiotic oocytes present in gsdf deficiency were more than those in wildtype ovaries, corresponding to the decrease of inhibin expression detected by real-time qPCR analysis, suggesting gsdf is tightly involved in the process of medaka oocyte development at early stage.