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

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.

Whole-genome assembly and annotation for the little yellow croaker (Larimichthys polyactis) provide insights into the evolution of hermaphroditism and gonochorism

The evolutionary direction of gonochorism and hermaphroditism is an intriguing mystery to be solved. The special transient hermaphroditic stage makes the little yellow croaker (Larimichthys polyactis) an appealing model for studying hermaphrodite formation. However, the origin and evolutionary relationship between of L. polyactis and Larimichthys crocea, the most famous commercial fish species in East Asia, remain unclear. Here, we report the sequence of the L. polyactis genome, which we found is ~706 Mb long (contig N50 = 1.21 Mb and scaffold N50 = 4.52 Mb) and contains 25,233 protein-coding genes. Phylogenomic analysis suggested that L. polyactis diverged from the common ancestor, L. crocea, approximately 25.4 million years ago. Our high-quality genome assembly enabled comparative genomic analysis, which revealed several within-chromosome rearrangements and translocations, without major chromosome fission or fusion events between the two species. The dmrt1 gene was identified as the male-specific gene in L. polyactis. Transcriptome analysis showed that the expression of dmrt1 and its upstream regulatory gene (rnf183) were both sexually dimorphic. Rnf183, unlike its two paralogues rnf223 and rnf225, is only present in Larimichthys and Lates but not in other teleost species, suggesting that it originated from lineage-specific duplication or was lost in other teleosts. Phylogenetic analysis shows that the hermaphrodite stage in male L. polyactis may be explained by the sequence evolution of dmrt1. Decoding the L. polyactis genome not only provides insight into the genetic underpinnings of hermaphrodite evolution, but also provides valuable information for enhancing fish aquaculture.

Development of epigenetic biomarkers for the identification of sex and thermal stress in fish using DNA methylation analysis and machine learning procedures

The sex ratio is a key ecological demographic parameter crucial for population viability. However, the epigenetic mechanisms operating during gonadal development regulating gene expression and the sex ratio remain poorly understood. Moreover, there is interest in the development of epigenetic markers associated with a particular phenotype or as sentinels of environmental effects. Here, we profiled DNA methylation and gene expression of 10 key genes related to sex development and stress, including steroidogenic enzymes, and growth and transcription factors. We provide novel information on the sex-related differences and on the influence of elevated temperature on these genes in zebrafish, a species with mixed genetic and environmental influences on sex ratios. We identified both positive (e.g., amh, cyp11c and hsd11b2) and negative (e.g., cyp11a1 and dmrt1) correlations in unexposed males, and negative correlation (amh) in exposed females between DNA methylation and gene expression levels. Further, we combined DNA methylation analysis with machine learning procedures and found a series of informative CpGs capable not only of correctly identifying sex (based on cyp19a1a DNA methylation levels) but also of identifying whether males and females had been exposed to abnormally elevated temperature when young (based on amh and foxl2a DNA methylation levels, respectively). This was achieved in the absence of conspicuous morphological alterations of the gonads. These DNA methylation-based epigenetic biomarkers represent molecular resources that can correctly recapitulate past thermal history and pave the way for similar findings in other species to assess potential ecological effects of environmental disturbances in the context of climate change.

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.

Neurohypophysial Hormones Associated with Osmotic Challenges in the Brain and Pituitary of the Euryhaline Black Porgy, Acanthopagrus schlegelii

Our study showed differential expression of the arginine vasotocin (avt)/isotocin (it) in the brain and pituitary gland of the euryhaline black porgy (Acanthopagrus schlegelii) during osmotic stress. A decrease in serum osmolality and increased cortisol levels were observed after acute transfer from seawater (SW) to freshwater (FW). The increased expressions of avt, avt receptor (avtr: v1a), and isotocin receptor (itr: itr1) transcripts on day 1 and it and itr transcripts on days 7 and 30 were found in the brains and pituitary glands of FW fish. Increased levels of avt mRNA in the diencephalon and avtr mRNA in the pituitary together with serum cortisol on day 1 of FW exposure indicated activation of the hypothalamic-pituitary-interrenal (HPI) axis. The expression levels of avtr and itr after FW transfer were increased in the pituitary on days 7 and 30. Furthermore, in situ hybridization demonstrated spatially differential expression of avt and itr transcripts in nucleus preopticus parvocellularis of pars gigantocellularis (PMgc), magnocellularis (PMmc), and parvocellularis (PMpc) of the preoptic area (POA). Positive signals for avt and it were highly abundant in PMpc after FW exposure. The data suggest involvement of neurohypophysial hormones in the brain (telencephalon and diencephalon) and pituitary for osmotic stress.

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

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

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.

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.

Comparative Genomics Studies on the dmrt Gene Family in Fish

Doublesex and mab-3-related transcription factor (dmrt) genes are widely distributed across various biological groups and play critical roles in sex determination and neural development. Here, we applied bioinformatics methods to exam cross-species changes in the dmrt family members and evolutionary relationships of the dmrt genes based on genomes of 17 fish species. All the examined fish species have dmrt1–5 while only five species contained dmrt6. Most fish harbored two dmrt2 paralogs (dmrt2a and dmrt2b), with dmrt2b being unique to fish. In the phylogenetic tree, 147 DMRT are categorized into eight groups (DMRT1–DMRT8) and then clustered in three main groups. Selective evolutionary pressure analysis indicated purifying selections on dmrt1–3 genes and the dmrt1–3–2(2a) gene cluster. Similar genomic conservation patterns of the dmrt1–dmrt3–dmrt2(2a) gene cluster with 20-kb upstream/downstream regions in fish with various sex-determination systems were observed except for three regions with remarkable diversity. Synteny analysis revealed that dmrt1, dmrt2a, dmrt2b, and dmrt3–5 were relatively conserved in fish during the evolutionary process. While dmrt6 was lost in most species during evolution. The high conservation of the dmrt1–dmrt3–dmrt2(2a) gene cluster in various fish genomes suggests their crucial biological functions while various dmrt family members and sequences across fish species suggest different biological roles during evolution. This study provides a molecular basis for fish dmrt functional analysis and may serve as a reference for in-depth phylogenomics.

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.

An estradiol-17β/miRNA-26a/cyp19a1a regulatory feedback loop in the protogynous hermaphroditic fish, Epinephelus coioides

Cyp19a1a is a key gene responsible for the production of estradiol-17β (E₂), the main functional estrogen and a major downstream regulator of reproduction in teleost fish. It is widely known that CYP19 gene expression, aromatase activity, and E₂ production can influence gonadal differentiation and sex reversal in teleost fish, but the feedback mechanisms whereby E₂ regulates cyp19a1a remain poorly understood, especially regarding the potential roles of endogenous small RNA molecules (miRNAs). Here, we identified miR-26a-5p as a regulatory factor of its predicted target gene (cyp19a1a). In vitro and in vivo studies showed that miR-26a-5p can decrease cyp19a1a expression. Furthermore, high doses of E₂ act as a repressor of miR-26a-5p. This study proposes a regulatory feedback loop whereby E₂ regulates cyp19a1a through miR-26a-5p, and suggests that this positive feedback is an important aspect of the control of E₂ production.

Initiation of sex change and gonadal gene expression in black sea bass (Centropristis striata) exposed to exemestane, an aromatase inhibitor

Many teleost fishes exhibit sequential hermaphroditism, where male or female gonads develop first and later undergo sex change. Model sex change species are characterized by social hierarchies and coloration changes, which enable experimental manipulations to better understand these processes. However, other species such as the protogynous black sea bass (Centropristis striata) do not exhibit these characteristics and instead receive research attention due to their importance in fisheries or aquaculture. Black sea bass social structure is unknown, which makes sex change sampling difficult, and few molecular resources are available. The purpose of the present study was to induce sex change using exemestane, an aromatase inhibitor, and assess gonadal gene expression using sex markers (amh, zpc2) and genes involved in steroidogenesis (cyp19a1a, cyp11b), estrogen signaling (esr1, esr2b), and apoptosis or atresia (aen, casp9, fabp11, parg, pdcd4, rif1). Overall, dietary exemestane treatment was effective, and most exposed females exhibited early histological signs of sex change and significantly higher rates of ovarian atresia relative to control females. Genes associated with atresia did not reflect this, however, as expression patterns in sex changing gonads were overall similar to those of ovaries, likely due to a whole ovary dilution effect of the RNA. Still, small but insignificant expression decreases during early sex change were detected for ovary-related genes (aen, casp9, fabp11, zpc2) and anti-apoptotic factors (parg, rif1). Exemestane treatment did not impact spermatogenesis or testicular gene expression, but testes were generally characterized by elevated steroidogenic enzyme and estrogen receptor mRNAs. Further research will be needed to understand these processes in black sea bass, using isolated ovarian follicles and multiple stages of sex change.

Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

DNA methyltransferases (dnmts) are responsible for DNA methylation and play important roles in organism development. In this study, seven dnmts genes (dnmt1, dnmt2, dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1, dnmt3bb.2) were identified in Nile tilapia. Comprehensive analyses of dnmts were performed using available genome databases from representative animal species. Phylogenetic analysis revealed that the dnmts family were highly conserved in teleosts. Based on transcriptome data from eight adult tilapia tissues, the dnmts were found to be dominantly expressed in the head kidney, testis and ovary. Analyses of the gonadal transcriptome data in different developmental stages revealed that all dnmts were expressed in both ovary and testis, and four de novo dnmts (dnmt3aa, dnmt3ab, dnmt3bb.1, dnmt3bb.2) showed higher expression in the testis than in the ovary. Furthermore, during sex reversal induced by Fadrozole, the expression of these four de novo dnmts increased significantly in treated group compared to female control group. By in situ hybridization, the seven dnmts were found to be expressed mainly in phase I and II oocytes of the ovary and spermatocytes of the testis. When gonads were incubated with a methyltransferase inhibitor (5-AzaCdR) in vitro, the expression of dnmts genes were down-regulated significantly, while the expression of cyp19a1a (a key gene in female pathway) and dmrt1 (a key gene in male pathway) increased significantly. Our results revealed the conservation of dnmts during evolution and indicated a potential role of dnmts in epigenetic regulation of gonadal development.

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.