The effect of gonadal hormones on the gene expression of brain-pituitary in protandrous black porgy, Acanthopagrus schlegelii

In black porgy (Acanthopagrus schlegelii), the brain-pituitary-testis (Gnrh-Gths-Dmrt1) axis plays a vital role in male fate determination and maintenance, and then inhibiting female development in further (puberty). However, the feedback of gonadal hormones on regulating brain signaling remains unclear. In this study, we conducted short-term sex steroid treatment and surgery of gonadectomy to evaluate the feedback regulation between the gonads and the brain. The qPCR results show that male phase had the highest gths transcripts; treatment with estradiol-17β (E2) or 17α-methyltestosterone (MT) resulted in the increased pituitary lhb transcripts. After surgery, apart from gnrh1, there is no difference in brain signaling genes between gonadectomy and sham fish. In the diencephalon/mesencephalon transcriptome, de novo assembly generated 283,528 unigenes; however, only 443 (0.16%) genes showed differentially expressed between sham and gonadectomy fish. In the present study, we found that exogenous sex steroids affect the gths transcription; this feedback control is related to the gonadal stage. Furthermore, gonadectomy may not affect gene expression of brain signaling (Gnrh-Gths axis). Our results support the communication between ovotestis and brain signaling (Gnrh-Gths-testicular Dmrt1) for the male fate.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

The Ovarian Transcriptome at the Early Stage of Testis Removal-Induced Male-To-Female Sex Change in the Protandrous Black Porgy Acanthopagrus schlegelii

Unlike gonochoristic fishes, sex is fixed after gonadal differentiation (primary sex determination), and sex can be altered in adults (secondary sex determination) of hermaphroditic fish species. The secondary sex determination of hermaphroditic fish has focused on the differences between testicular tissue and ovarian tissue during the sex change process. However, comprehensive studies analyzing ovarian tissue or testicular tissue independently have not been performed. Hermaphroditic black porgy shows a digonic gonad (ovarian tissue with testicular tissue separated by connective tissue). Protandrous black porgy has stable maleness during the first two reproductive cycles (<2 years old), and approximately 50% enter femaleness (natural sex change) during the third reproductive cycle. Precocious femaleness is rarely observed in the estradiol-17β (E2)-induced female phase (oocytes maintained at the primary oocyte stage), and a reversible female-to-male sex change is found after E2 is withdrawn in <2-year-old fish. However, precocious femaleness (oocytes entering the vitellogenic oocyte stage) is observed in testis-removed fish in <2-year-old fish. We used this characteristic to study secondary sex determination (femaleness) in ovarian tissue via transcriptomic analysis. Cell proliferation analysis showed that BrdU (5-bromo-2′-deoxyuridine)-incorporated germline cells were significantly increased in the testis-removed fish (female) compared to the control (sham) fish (male) during the nonspawning season (2 months after surgery). qPCR analysis showed that there were no differences in pituitary-releasing hormones (lhb and gtha) in pituitary and ovarian steroidogenesis-related factors (star, cyp11a1, hsd3b1, and cyp19a1a) or female-related genes (wnt4a, bmp15, gdf9, figla, and foxl2) in ovarian tissues between intact and testis-removed fish (2 months after surgery). Low expression of pituitary fshb and ovarian cyp17a1 was found after 2 months of surgery. However, we did find small numbers of genes (289 genes) showing sexual fate dimorphic expression in both groups by transcriptomic analysis (1 month after surgery). The expression profiles of these differentially expressed genes were further examined by qPCR. Our present work identified several candidate genes in ovarian tissue that may be involved in the early period of secondary sex determination (femaleness) in black porgy. The data confirmed our previous suggestion that testicular tissue plays an important role in secondary sex determination in protandrous black porgy.

Expression and Transcript Localization of star, sf-1, and dax-1 in the Early Brain of the Orange-Spotted Grouper Epinephelus coioides

We investigated the developmental expression and localization of sf-1 and dax-1 transcripts in the brain of the juvenile orange-spotted grouper in response to steroidogenic enzyme gene at various developmental ages in relation to gonadal sex differentiation. The sf-1 transcripts were significantly higher from 110-dah (day after hatching) and gradually increased up to 150-dah. The dax-1 mRNA, on the other hand, showed a decreased expression during this period, in contrast to sf-1 expression. At the same time, the early brain had increased levels of steroidogenic gene (star). sf-1 and star hybridization signals were found to be increased in the ventromedial hypothalamus at 110-dah; however, dax-1 mRNA signals decreased in the early brain toward 150-dah. Furthermore, the exogenous estradiol upregulated star and sf-1 transcripts in the early brain of the grouper. These findings suggest that sf-1 and dax-1 may have an antagonistic expression pattern in the early brain during gonadal sex differentiation. Increased expression of steroidogenic gene together with sf-1 during gonadal differentiation strongly suggests that sf-1 may play an important role in the juvenile grouper brain steroidogenesis and brain development.

Comparative Proteomics of the Acanthopagrus schlegelii Gonad in Different Sex Reversal

A substantial proportion of Acanthopagrus schlegelii individuals change sex from male to female during their lifetime. However, the mechanisms underlying sex change are unknown. In this research, iTRAQ analyses of proteins obtained from A.schlegelii gonads in four different stages of development were compared. In total, 4692 proteins were identified, including common sex-specific proteins, such as sperm-associated antigen 6 and cilia- and flagella-associated proteins in males, and zona pellucida sperm-binding proteins in females. Furthermore, proteins involved in the integrin signaling pathway, inflammation mediated by the chemokine and cytokine signaling pathways, pyruvate metabolism, CCKR signaling map, de novo purine biosynthesis and the ubiquitin proteasome pathway were upregulated in female gonads, whereas proteins implicated in DNA replication, the heterotrimeric G-protein signaling pathway, Gi alpha- and Gs alpha-mediated pathways, wnt signaling pathway, and hedgehog signaling pathway were upregulated in male gonads. Interestingly, cathepsins were only identified in ovaries, indicating their potential involvement in rapid ovarian development. Apoptosis-related proteins expressed in ovaries (such as MAPK and Cdc42) may protect them from cancer. This is the first report on the gonad proteome from A.schlegelii in different stages of sex reversal, and the results provide important fundamental data for studying the molecular mechanisms of sex reversal.

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.

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.

Comparison of Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

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.

Molecular and cellular regulation on sex change in hermaphroditic fish, with a special focus on protandrous black porgy, Acanthopagrus schlegelii

In teleost fish, sex can be determined by genetic factors, environmental factors, or both. Unlike in gonochoristic fish, in which sex is fixed in adults, sex can change in adults of hermaphroditic fish species. Thus, sex is generated during the initial gonadal differentiation stage (primary sex differentiation) and later during sexual fate alternation (secondary sex differentiation) in hermaphroditic fish species. Depending on the species, sex phase alternation can be induced by endogenous cues (such as individual age and body size) or by social cues (such as sex ratio or relative body size within the population). In general, the fluctuation in plasma estradiol-17β (E2) levels is correlated with the sexual fate alternation in hermaphroditic fish. Hormonal treatments can artificially induce sexual phase alternation in sequential hermaphroditic fishes, but in a transient and reversible manner. This is the case for the E2-induced female phase in protandrous black porgy and the methyltestosterone (MT)- or aromatase inhibitor (AI)-induced male phase in protogynous grouper. Recent reviews have focused on the different forms of sex change in fish who undergo sequential sex change, especially in terms of gene expression and the role of hormones. In this review, we use the protandrous black porgy, a nonsocial cue-influenced hermaphroditic species, with digonic gonads (ovarian and testis separated by a connective tissue), as a model to describe our findings and discuss the molecular and cellular regulation of sexual fate determination in hermaphroditic fish.

Activation of the brain-pituitary-gonadotropic axis in the black porgy Acanthopagrus schlegelii during gonadal differentiation and testis development and effect of estradiol treatment

Previous studies revealed an estradiol (E2)-dependent peak in brain activity, including neurosteroidogenesis and neurogenesis in the black porgy during the gonadal differentiation period. The brain-pituitary-gonadotropic axis is a key regulator of reproduction and may also be involved in gonadal differentiation, but its activity and potential role in black porgy during the gonadal differentiation period is still unknown. The present study analyzed the expression of regulatory factors involved in the gonadotropic axis at the time of gonadal differentiation (90, 120, 150 days after hatching [dah]) and subsequent testicular development (180, 210, 300 dah). In agreement with previous studies, expression of brain aromatase cyp19a1b peaked at 120 dah, and this was followed by a gradual increase during testicular development. The expression of gonadotropin subunits increased slightly but not significantly during gonadal differentiation and then increased significantly at 300 dah. In contrast, the expression of brain gnrh1 and pituitary gnrh receptor 1 (gnrhr1) exhibited a pattern with two peaks, the first at 120 dah, during the period of gonadal differentiation, and the second peak during testicular development. Gonad fshr and lhcgr increased during gonadal differentiation period with highest transcript level in prespawning season during testicular development. This suggests that the early activation of brain gnrh1, pituitary gnrhr1 and gths, and gonad gthrs might be involved in the control of gonadal differentiation. E2 treatment increased brain cyp19a1b expression at each sampling time, in agreement with previous studies in black porgy and other teleosts. E2 also significantly stimulated the expression of pituitary gonadotropin subunits at all sampling times, indicating potential E2-mediated steroid feedback. In contrast, no significant effect of E2 was observed on gnrh1. Moreover, treatment of AI or E2 had no statistically significant effect on brain gnrh1 transcription levels during gonadal differentiation. This indicated that the early peak of gnrh1 expression during the gonadal differentiation period is E2-independent and therefore not directly related to the E2-dependent peak in brain neurosteroidogenesis and neurogenesis also occurring during this period in black porgy. Both E2-independent and E2-dependent mechanisms are thus involved in the peak expression of various genes in the brain of black porgy at the time of gonadal differentiation.

The Kiss2/GPR54 system stimulates the reproductive axis in male black porgy, Acanthopagrus schlegelii

Although it is well established that the Kiss1/GPR54 system stimulates the reproductive axis in mammals, its functional roles, especially in male reproduction of non-mammalian species, is less clear. In this study, we have isolated the full-length kiss2 and gpr54 cDNAs from black porgy (Acanthopagrus schlegelii). The Kiss2 precursor expressed from kiss2 comprises 124 amino acids and contains a highly conserved 10-amino acid sequence, Kiss2-10 (FNFNPFGLRF). GPR54 comprises 375 amino acid residues and contains distinct characteristics of G protein-coupled receptors. Real-time PCR analysis indicated that kiss2 and gpr54 were expressed highly in the brain regions. Moreover, intraperitoneal injection of porgy Kiss2-10 could stimulate genes expression of the gpr54, gnrh1, gnrh3, fshβ, lhβ, p450c17, star, and ar, and the serum testerone level in male black porgy. Our findings demonstrate that the Kisspeptin stimulates the male reproductive axis in black porgy.

A comparative transcriptomic study on developmental gonads provides novel insights into sex change in the protandrous black porgy (Acanthopagrus schlegelii)

Protandrous black porgy (Acanthopagrus schlegelii) is a popular and valuable commercial marine fish in China and East Asian countries. Controlling and managing its breeding has been an imperative step towards obtaining a sustainable supply of this fish in aquaculture production systems. Therefore, study on the molecular mechanisms of sex change in black porgy has both scientific and commercial importance. Previously, we identified some candidate genes related to sex determination and differentiation from a high-quality genome assembly of the black porgy. In the present study, transcriptome sequencing of developmental gonads (including testis, ovotestis and ovary) of black porgy was performed to further investigate the sex-change mechanisms. Our results showed that the highly expressed male-related genes (dmrt1, piwi1, piwi2, sox9, sox30 and amh) at the male phase were significantly down-regulated to a substantial degree at the intersexual stage, and the female-related genes (jnk1, vasa, wnt4, figla and foxl2) were distinctly up-regulated when the fish grows into a female adult, suggesting the potential roles of these genes in sex change of the black porgy. These data also support a previous hypothesis that the femaleness will be switched on when the testis is entering the degenerated stage through the diminished dmrt1 expression. Our transcriptome data provide a very useful genomic resource for future studies on sex change and practical aquaculture in the black porgy.

Primary males guide the femaleness through the regulation of testicular Dmrt1 and ovarian Cyp19a1a in protandrous black porgy

Controlling the development of the sexes is critically important for the broodstock management in aquaculture. Sex steroids are widely used for sex control of fish. However, hermaphroditic fish have a plastic sex, and a stable sex is difficult to maintain with sex steroids. We used the black porgy (Acanthopagrus schlegelii) as a model to understand the possible mechanism of sexual fate decision. Low exogenous estradiol (E2) induced male development. In contrast, high exogenous E2 induced the regression of the testis and the development of the ovary and resulted in an unstable expression of femaleness (passive femaleness, with ovaries containing only the primary oocytes). The removal of testicular tissue by surgery resulted in the early development of vitellogenic oocytes and active femaleness. Our data also demonstrated that the male-to-female sex change is blocked by the maintenance of male function with gonadotropin-induced dmrt1 expression in the testis. Furthermore, our data also indicated that ovarian cyp19a1a expression is regulated by the testis through epigenetic modifications. Therefore, the primary male guides the femaleness in the protandrous black porgy and the transition of sexual fate from male to female is determined by the status of the testicular tissue.

Robust gdf9 and bmp15 expression in the oocytes of ovotestes through the Figla-independent pathway in the hermaphroditic black porgy, Acanthopagrus schlegelii

More than 1,500 fish species are hermaphroditic, but no hermaphroditic lineage appears to be evolutionarily ancient in fishes. Thus, whether more than one sex at a time was present during the evolutionary shift from gonochorism to hermaphroditism in fishes is an intriguing question. Ectopic oocytes were created in the ovotestes of protandrous black porgy via the withdrawal of estradiol (E2) administration. These ectopic oocytes reprogrammed the surrounding cells, which changed from Sertoli cells to follicle-like cells. We observed that gdf9 and bmp15 expression was localized in the primary oocytes and gradually decreased after oocytes entered a secondary oocyte stage. Robust expression of gdf9 and bmp15 in ectopic oocytes was associated with the surrounding Sertoli cells. However, blocking Cyp19a1a activity and increasing androgen levels did not stimulate the expression of gdf9 and bmp15. Thus, the robust gdf9 and bmp15 expression was not related to the inappropriate male microenvironment. Furthermore, in vitro data demonstrated that gdf9 and bmp15 were not downstream genes of Figla signaling. Therefore, our results suggest that there are two independent mechanisms, a Figla-dependent pathway and a Figla-independent pathway, by which oocyte-surrounding cells are altered from a male somatic fate to a female somatic fate. This functional switch might clarify how oocytes created an appropriate microenvironment during the transition from the ancient gonochorism to the present hermaphroditism.

Expression profile of amh/Amh during bi-directional sex change in the protogynous orange-spotted grouper Epinephelus coioides

Gonadal differentiation is tightly regulated by the initial sex determining gene and the downstream sex-related genes in vertebrates. However, sex change in fish can alter the sexual fate from one sex to the other. Chemical-induced maleness in the protogynous orange-spotted grouper is transient, and a reversible sex change occurs after the chemical treatment is withdrawn. We used these characteristics to study Amh signaling during bi-directional sex change in the grouper. We successfully induced the female-to-male sex change by chemical (aromatase inhibitor, AI, or methyltestosterone, MT) treatment. A dormant gonad (a low proliferation rate of early germ cells and no characteristics of both sexes) was found during the transient phase of reversible male-to-female sex change after the withdrawal of chemical administration. Our results showed that amh (anti-mullerian hormone) and its receptor amhr2 (anti-mullerian hormone receptor type 2) were significantly increased in the gonads during the process of female-to-male sex change. Amh is expressed in the Sertoli cells surrounding the type A spermatogonia in the female-to-male grouper. Male-related gene (dmrt1 and sox9) expression was immediately decreased in MT-terminated males during the reversible male-to-female sex change. However, Amh expression was found in the surrounding cells of type A spermatogonia-like cells during the transient phase of reversible male-to-female sex change. This phenomenon is correlated with the dormancy of type A spermatogonia-like cells. Thus, Amh signaling is suggested to play roles in regulating male differentiation during the female-to-male sex change and in inhibiting type-A spermatogonia-like cell proliferation/differentiation during the reversible male-to-female sex change. We suggest that Amh signaling might play dual roles during bi-directional sex change in grouper.

Transcriptome analysis of three critical periods of ovarian development in Yellow River carp (Cyprinus carpio)

Ovary development is a complex process involving numerous genes; the molecular mechanism underlying the ovary development of carp is still unknown. Here we used Illumina HiSeq™ 2500 to explore the transcriptome of undifferentiated gland (PG), juvenile ovary (OJ) and adult ovary (OA) of Yellow River carp (Cyprinus carpio). A total of 58,749 unigenes were obtained, comprising 45,707 known genes and 13,042 new genes. We identified differentially-expressed genes (DEGs) during development and characterized the functional properties of DEGs by comparison with the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes databases. qRT-PCR was used to analyze the expression of 22 DEGs and the results corresponded with those of RNA-Seq. Among DEGs between PG and OJ, some upstream regulators of gonad development were upregulated in PG, such as cyp19a and sox9, while some oocyte-specific genes were upregulated in OJ, such as nobox, bmp15 and zp2. Among DEGs between OJ and OA, many oocyte physiological function-related genes were upregulated in OA, such as fem-1 and foxl2. GO analysis showed a higher number of DEGs from PG-OJ analysis were assigned to reproduction terms. Furthermore, our investigation has also revealed DEGs identified from PG-OJ analysis were enriched in several important functional pathways, such as Fanconi anemia and the notch signal pathway. These data suggested a dynamic shift in gene expression during ovary development, and DEGs between PG and OJ provided crucial candidate gene data for the study of ovarian differentiation. Additionally, a total of 1,776,769 single nucleotide polymorphisms and 157,279 INDEs were revealed from transcriptome data. This result will contribute to knowledge of ovary differentiation of Yellow River carp.

Co-expression analysis and identification of fecundity-related long non-coding RNAs in sheep ovaries

Small Tail Han sheep, including the FecB^BFecB^B (Han BB) and FecB⁺ FecB⁺ (Han++) genotypes, and Dorset sheep exhibit different fecundities. To identify novel long non-coding RNAs (lncRNAs) associated with sheep fecundity to better understand their molecular mechanisms, a genome-wide analysis of mRNAs and lncRNAs from Han BB, Han++ and Dorset sheep was performed. After the identification of differentially expressed mRNAs and lncRNAs, 16 significant modules were explored by using weighted gene coexpression network analysis (WGCNA) followed by functional enrichment analysis of the genes and lncRNAs in significant modules. Among these selected modules, the yellow and brown modules were significantly related to sheep fecundity. lncRNAs (e.g., NR0B1, XLOC_041882, and MYH15) in the yellow module were mainly involved in the TGF-β signalling pathway, and NYAP1 and BCORL1 were significantly associated with the oxytocin signalling pathway, which regulates several genes in the coexpression network of the brown module. Overall, we identified several gene modules associated with sheep fecundity, as well as networks consisting of hub genes and lncRNAs that may contribute to sheep prolificacy by regulating the target mRNAs related to the TGF-β and oxytocin signalling pathways. This study provides an alternative strategy for the identification of potential candidate regulatory lncRNAs.

Sexual Fate Reprogramming in the Steroid-Induced Bi-Directional Sex Change in the Protogynous Orange-Spotted Grouper, Epinephelus coioides

Androgen administration has been widely used for masculinization in fish. The mechanism of the sex change in sexual fate regulation is not clear. Oral administration or pellet implantation was applied. We orally applied an aromatase inhibitor (AI, to decrease estrogen levels) and 17α-methyltestosterone (MT, to increase androgen levels) to induce masculinization to clarify the mechanism of the sex change in the protogynous orange-spotted grouper. After 3 mo of AI/MT administration, male characteristics were observed in the female-to-male sex change fish. These male characteristics included increased plasma 11-ketotestosterone (11-KT), decreased estradiol (E2) levels, increased male-related gene (dmrt1, sox9, and cyp11b2) expression, and decreased female-related gene (figla, foxl2, and cyp19a1a) expression. However, the reduced male characteristics and male-to-female sex change occurred after AI/MT-termination in the AI- and MT-induced maleness. Furthermore, the MT-induced oocyte-depleted follicle cells (from MT-implantation) had increased proliferating activity, and the sexual fate in a portion of female gonadal soma cells was altered to male function during the female-to-male sex change. In contrast, the gonadal soma cells were not proliferative during the early process of the male-to-female sex change. Additionally, the male gonadal soma cells did not alter to female function during the male-to-female sex change in the AI/MT-terminated fish. After MT termination in the male-to-female sex-changed fish, the differentiated male germ cells showed increased proliferating activities together with dormancy and did not show characteristics of both sexes in the early germ cells. In conclusion, these findings indicate for the first time in a single species that the mechanism involved in the replacement of soma cells is different between the female-to-male and male-to-female sex change processes in grouper. These results also demonstrate that sexual fate determination (secondary sex determination) is regulated by endogenous sex steroid levels.

Molecular cloning and characterization of amh and dax1 genes and their expression during sex inversion in rice-field eel Monopterus albus

The full-length cDNAs of amh and dax1 in the hermaphrodite, rice-field eel (Monopterus albus), were cloned and characterized in this study. Multiple sequence alignment revealed Dax1 was well conserved among vertebrates, whereas Amh had a low degree of similarity between different vertebrates. Their expression profiles in gonads during the course of sex inversion and tissues were investigated. The tissue distribution indicated amh was expressed mostly in gonads and was scarcely detectable in other tissues, whereas the expression of dax1 was widespread among the different tissues, especially liver and gonads. amh was scarcely detectable in ovaries whereas it was abundantly expressed in both ovotestis and testis. By contrast, dax1 was highly expressed in ovaries, especially in ♀IV (ovaries in IV stage), but it was decreased significantly in ♀/♂I (ovotestis in I stage). Its expression was increased again in ♀/♂III (ovotestis in III stage), and then decreased to a low level in testis. These significant different expression patterns of amh and dax1 suggest the increase of amh expression and the decline of dax1 expression are important for the activation of testis development, and the high level of amh and a low level of dax1 expression are necessary for maintenance of testis function.

Activation of brain steroidogenesis and neurogenesis during the gonadal differentiation in protandrous black porgy, Acanthopagrus schlegelii

The early brain development, at the time of gonadal differentiation was investigated using a protandrous teleost, black porgy. This natural model of monosex juvenile fish avoids the potential complexity of sexual dimorphism. Brain neurogenesis was evaluated by histological analyses of the diencephalon, at the time of testicular differentiation (in fish between 90 and 150 days after hatching). Increases in the number of both Nissl-stained total brain cells, and Pcna-immunostained proliferative brain cells were observed in specific area of the diencephalon, such as ventromedialis thalami and posterior preoptic area, revealing brain cell proliferation. qPCR analyses showed significantly higher expression of the radial glial cell marker blbp and neuron marker bdnf. Strong immunohistochemical staining of Blbp and extended cellular projections were observed. A peak expression of aromatase (cyp19a1b), as well as an increase in estradiol (E2 ) content were also detected in the early brain. These data demonstrate that during gonadal differentiation, the early brain exhibits increased E2 synthesis, cell proliferation, and neurogenesis. To investigate the role of E2 in early brain, undifferentiated fish were treated with E2 or aromatase inhibitor (AI). E2 treatment upregulated brain cyp19a1b and blbp expression, and enhanced brain cell proliferation. Conversely, AI reduced brain cell proliferation. Castration experiment did not influence the brain gene expression patterns and the brain cell number. Our data clearly support E2 biosynthesis in the early brain, and that brain E2 induces neurogenesis. These peak activity patterns in the early brain occur at the time of gonad differentiation but are independent of the gonads.

Genes encoding aromatases in teleosts: evolution and expression regulation

Cytochrome P450 aromatases, encoded by cyp19a1 genes, catalyzes the conversion of androgens to estrogens and plays important roles in the reproduction of vertebrates. Vertebrate cyp19a1 genes showed high synteny in chromosomal locations and conservation in sequences during evolution. However, amphioxus cyp19a1 does not show synteny to vertebrate cyp19a1. Teleost fish possess two copies of the cyp19a1 gene, which were postulated to result from a fish-specific genome duplication. The duplicated copies of fish cyp19a1 genes evolved into the brain and ovarian forms of cytochrome P450 aromatase genes, cyp19a1a and cyp19a1b, respectively, with different regulatory mechanisms of expression, through subfunctionalization under long-term selective pressure. In addition to the estradiol (E2) auto-regulatory loop, there may be other mechanisms responsible for the high expression of aromatase in the teleost brain. The study of the two cyp19a1 copies in teleost fish will shed light on the general evolution, function, and regulation of vertebrate cyp19a1.

The sex-specific transcriptome of the hermaphrodite sparid sharpsnout seabream (Diplodus puntazzo)

Background

Teleosts are characterized by a remarkable breadth of sexual mechanisms including various forms of hermaphroditism. Sparidae is a fish family exhibiting gonochorism or hermaphroditism even in closely related species. The sparid Diplodus puntazzo (sharpsnout seabream), exhibits rudimentary hermaphroditism characterized by intersexual immature gonads but single-sex mature ones. Apart from the intriguing reproductive biology, it is economically important with a continuously growing aquaculture in the Mediterranean Sea, but limited available genetic resources. Our aim was to characterize the expressed transcriptome of gonads and brains through RNA-Sequencing and explore the properties of genes that exhibit sex-biased expression profiles.

Results

Through RNA-Sequencing we obtained an assembled transcriptome of 82,331 loci. The expression analysis uncovered remarkable differences between male and female gonads, while male and female brains were almost identical. Focused search for known targets of sex determination and differentiation in vertebrates built the sex-specific expression profile of sharpsnout seabream. Finally, a thorough genetic marker discovery pipeline led to the retrieval of 85,189 SNPs and 29,076 microsatellites enriching the available genetic markers for this species.

Conclusions

We obtained a nearly complete source of transcriptomic sequence as well as marker information for sharpsnout seabream, laying the ground for understanding the complex process of sex differentiation of this economically valuable species. The genes involved include known candidates from other vertebrate species, suggesting a conservation of the toolkit between gonochorists and hermaphrodites.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-655) contains supplementary material, which is available to authorized users.

The switch of secondary sex determination in protandrous black porgy, Acanthopagrus schlegeli

Hermaphrodites have both sexes during their life, including an initial primary sex determination and in later stage maintenance one of the sexual fates (secondary sex determination). Sex change (secondary sex determination) occurs in animals, but it is lost in amphibians through, mammals in vertebrates. Teleosts have various strategies and mechanisms of sex determination including genetic and environmental cues. However, the mechanisms by which the cues guide sex change are complicated in fish. This manuscript reviews our understanding of these processes in protandrous black porgy at the gonadal and neuroendocrine levels. Our studies addressed the process of sex change through brain-pituitary-gonad axis, and then secondary sex determination was switched by the fate of testis.

Sex differentiation and sex change in the protandrous black porgy, Acanthopagrus schlegeli

Protandrous black porgy fish, Acanthopagrus schlegeli, have a striking life cycle with a male sex differentiation at the juvenile stage and male-to-female sex change at 3 years of age. We had characterized the sex differentiation and sex change in this species by the integrative approaches of histology, endocrine and molecular genetics. The fish differentiated in gonad at the age around 4-months and the gonad further developed with a bisexual gonad for almost for 3 years and sex change at 3 year of age. An antagonistic relationship in the testicular and ovarian tissues was found during the development of the gonadal tissue. Male- (such as sf-1, dmrt1, dax-1 and amh) and female- (such as wnt4, foxl2 and cyp19a1a) promoting genes were associated with testicular and ovarian development, respectively. During gonadal sex differentiation, steroidogenic pathway and estrogen signaling were also highly expressed in the brain. The increased expression of sf-1 and wnt4, cyp19a1a in ovarian tissue and decreased expression of dax-1 in the ovarian tissue may play important roles in sex change from a male-to-female. Endocrine factors such as estradiol and luteinizing hormone may also involve in the natural sex change. Estradiol induced the expression of female-promoting genes and resulted in the precocious sex change in black porgy. Our series of studies shed light on the sex differentiation and sex change in protandrous black porgy and other animals.

Ovarian aromatase and estrogens: a pivotal role for gonadal sex differentiation and sex change in fish

The present review focuses on the roles of estrogens and aromatase (Cyp19a1a), the enzyme needed for their synthesis, in fish gonadal sex differentiation. Based on the recent literature, we extend the already well accepted hypothesis of an implication of estrogens and Cyp19a1a in ovarian differentiation to a broader hypothesis that would place estrogens and Cyp19a1a in a pivotal position to control not only ovarian, but also testicular differentiation, in both gonochoristic and hermaphrodite fish species. This working hypothesis states that cyp19a1a up-regulation is needed not only for triggering but also for maintaining ovarian differentiation and that cyp19a1a down-regulation is the only necessary step for inducing a testicular differentiation pathway. When considering arguments for and against, most of the information available for fish supports this hypothesis since either suppression of cyp19a1a gene expression, inhibition of Cyp19a1a enzymatic activity, or blockage of estrogen receptivity are invariably associated with masculinization. This is also consistent with reports on normal gonadal differentiation, and steroid-modulated masculinization with either androgens, aromatase inhibitors or estrogen receptor antagonists, temperature-induced masculinization and protogynous sex change in hermaphrodite species. Concerning the regulation of fish cyp19a1a during gonadal differentiation, the transcription factor foxl2 has been characterized as an ovarian specific upstream regulator of a cyp19a1a promoter that would co-activate cyp19a1a expression, along with some additional partners such as nr5a1 (sf1) or cAMP. In contrast, upstream factors potentially down-regulating cyp19a1a during testicular differentiation are still hypothetical, such as the dmrt1 gene, but their definitive characterization as testicular repressors of cyp19a1a would strongly strengthen the hypothesis that early testicular differentiation would need active repression of cyp19a1a expression.

Molecular characterization of cyp11a1 and cyp11b1 and their gene expression profile in pejerrey (Odontesthes bonariensis) during early gonadal development

Sex steroids are known to be involved in gonadal differentiation in fish, but whether androgens are early mediators of testis differentiation remains unclear. We studied the sex-related developmental variations in the gene expression of two key enzymes involved in steroids and androgen synthesis (cyp11a1 and cyp11b1) in trunks and isolated gonads of pejerrey (Odontesthes bonariensis) larvae during and after the sex determination period. Also, and in order to have a better characterization of this process we studied the expression of Sertoli (dmrt1, amh, sox9) and Leydig (nr5a1 or sf-1) cell markers as well as a gene with higher expression in females (cyp19a1a). No clear differences were observed in the expression of cyp11a1 and cyp11b1 during the temperature-sensitive window in the trunk of pejerrey larvae. Nevertheless, a clear increase of cyp11b1 was observed in isolated gonads taken from fish reared at the male producing temperature. In these gonads we also confirmed the trends of genes with higher expression in males (dmrt1, amh) and females (cyp19a1a) as previously described in larval trunks of pejerrey. Our results showed that the expression of cyp11b1 was positively associated with the morphological differentiation of the testis. Nevertheless the involvement of 11-oxygenated androgens during the temperature-sensitive window could not be clearly established.

Salinity effects on the expression of osmoregulatory genes in the euryhaline black porgy Acanthopagrus schlegeli

Black porgy is a marine euryhaline species with a capacity to cope with demands in a wide range of salinities and thus is a perfect model-fish to study osmoregulatory responses to salinity-acclimated processes and their hormonal control. The present study was performed to understand the regulatory changes in hormone, hormone receptors and important osmoregulatory genes in pituitary, gill, intestine and kidney in response to acute salinity stress. Transcript levels were analyzed by quantitative real-time PCR following acute salinity challenge by direct transfer of seawater (SW) acclimatized fish to fresh water (FWBP) and vice versa (SWBP). SW acclimation significantly increased plasma osmolality and intestine Na+/K+-ATPase (NKA) activity while FW acclimation increased plasma cortisol and branchial NKA activity. Plasma osmolality and chloride concentration decreased in FWBP whereas GH levels remained unchanged in both FWBP and SWBP. Comparative analysis of gene profiles between FWBP and SWBP showed that pituitary prolactin transcript increased significantly in FWBP. Prolactin receptor (PRLR) transcripts increased in gill of FWBP while it decreased in gill and kidney of SWBP. NKA transcripts increased in gill of both FWBP and SWBP, while it decreased in intestine of FWBP and increased in intestine and kidney of SWBP. Glucocorticoid receptor (GR) transcripts decreased in intestine and kidney of FWBP while it increased in gill and intestine of SWBP. No significant changes were observed in growth hormone receptor (GHR) transcripts of both FWBP and SWBP in pituitary, gill, intestine and kidney. Our current data demonstrated the correlation between PRLR gene expression in relation to FW adaptation, and GR gene expression in relation to SW adaptation in euryhaline black porgy. The results indicate that black porgy has an excellent osmoregulatory capacity and is capable of withstanding large variations in salinity.

Dual roles of cyp19a1a in gonadal sex differentiation and development in the protandrous black porgy, Acanthopagrus schlegeli

Protandrous black porgy fish, Acanthopagrus schlegeli, have a striking life cycle, with male sex differentiation at the juvenile stage, a bisexual gonad during first 2 yr of life, and a male-to-female sex change (with vitellogenic oocytes) at 3 yr of age. The present study investigated the role of aromatase (cyp19a1a/Cyp19a1a) in gonadal development in this species, especially in relation to sexual differentiation and sex change. Fish of various ages were treated with estradiol (E2) or aromatase inhibitor (AI) to determine whether manipulation of the hormonal environment has an impact on these processes. We report an integrative immunohistochemical, cellular, and molecular data set describing these interesting phenomena. During male sex differentiation, high levels of cyp19a1a/Cyp19a1a expression were observed in the undifferentiated gonad (4 mo of age), in marked contrast to the low cyp19a1a/Cyp19a1a levels detected in the differentiated testis at the age of 5-6 mo. A low dose of E2 (0.25 mg/kg feed) stimulated testicular growth and function in sexually differentiated fish, whereas a high dose of E2 (6 mg/kg feed) induced female development. Furthermore, administration of AI suppressed male development and promoted female sexual differentiation. An increased number of figla transcripts (an oocyte-specific gene) were observed prior to cyp19a1a expression, concomitant with the development of oogonia and early primary oocytes in the ovaries of both E2- and AI-treated groups. Immunohistochemical Pcna staining showed that the regression of testicular tissue occurred prior to the development of ovarian tissue in both E2- and AI-induced females. The importance of cyp19a1a in female development was further demonstrated by the increase in cyp19a1a transcripts during the naturally occurring sex change. Transcripts of foxl2 increased in the gonads of 2- to 3-yr-old black porgy during the early stages of the natural sex change, followed by a gradual elevation of cyp19a1a levels. The levels of both genes peaked in the resulting ovarian tissue. Thus, cyp19a1a/Cyp19a1a plays dual roles in the gonadal development, namely, in testicular development during the initial period of sexual differentiation and later in ovarian development during the natural sex change.