Transcriptomic Analysis Provides New Insights into Oocyte Growth and Maturation in Greater Amberjack (Seriola dumerili)

The greater amberjack (Seriola dumerili) is an emerging marine fish that is increasingly favored in aquaculture. Currently, there are few studies on the development and regulation of greater amberjack ovaries. In this study, the ovary transcriptome profiles of greater amberjack at three different stages (stage II, III, and IV) were performed, and identified the genes and pathways that may play significant roles in the processes of follicle growth and maturation. A total of 6597, and 1061 differentially expressed genes (DEGs) were detected in FII vs. FIII, FIII vs. FIV, and FII vs. FIV stages, respectively. GO and KEGG enrichment analyses revealed that these DEGS are primarily involved in steroid hormone biosynthesis (e.g., cyp11a1, cyp17a1, cyp19a1a, hsd3b1, esr1), lipid metabolism (e.g., plpp3, lpl, pld1, and fabp10a), and meiotic arrest and resumption (e.g., pgr, arb, ccnd2, adcy2, adcy9, myl9, calm1). Additionally, several signaling pathways involved in ovarian development have been identified, including the PI3K-Akt, Wnt, TGF-beta, GnRH, and immune-related signaling pathways. qPCR results of nine representative genes related to steroid hormone synthesis and cell growth verified the reliability of the generated RNA-seq data. This research contributes to our comprehension of the molecular processes underlying ovarian growth and maturation in marine fishes and provides a theoretical basis for the investigation of functional genes associated with oogenesis in greater amberjack.

Understanding Genetic Regulation of Sex Differentiation in Hermaphroditic Fish

As a fundamental taxonomic group within vertebrates, fish represent an invaluable resource for investigating the mechanisms underlying sex determination and differentiation owing to their extensive geographical distribution and rich biodiversity. Within this biological cohort, the processes of sex determination and differentiation are intricately governed by both genetic factors and the complex interplay of environmental cues. While variations in external environmental factors, particularly temperature, can exert a modulatory influence on sex differentiation in fish to a limited degree, genetic factors remain the primary determinants of sexual traits. Hermaphroditic fish display three distinct types of sexual transitions: protandry (male to female), protogyny (female-to-male), bidirectional sex change (both directions serially). These fish, characterized by their unique reproductive strategies and sexual plasticity, serve as exemplary natural models for elucidating the mechanisms of sex differentiation and sexual transitions in fish. The present review delves into the histological dynamics during gonadal development across three types of sequential hermaphroditic fish, meticulously delineating the pivotal characteristics at each stage, from the inception of primordial gonads to sexual specialization. Furthermore, it examines the regulatory genes and associated signaling pathways that orchestrate sex determination and differentiation. By systematically synthesizing these research advancements, this paper endeavors to offer a comprehensive and profound insight into the intricate mechanisms governing sex differentiation in sequential hermaphroditic fish.

Cyp19a1a Promotes Ovarian Maturation through Regulating E2 Synthesis with Estrogen Receptor 2a in Pampus argenteus (Euphrasen, 1788)

In the artificial breeding of Pampus argenteus (Euphrasen, 1788), female fish spawn before male release sperm, which indicates rapid ovarian development. In fish, aromatase is responsible for converting androgens into estrogens and estrogen plays a crucial role in ovarian development. In this study, we aimed to investigate the potential role of brain-type and ovarian-type aromatase to study the rapid ovarian development mechanism. The results showed that cyp19a1a was mainly expressed in the ovary and could be classified as the ovarian type, whereas cyp19a1b could be considered as the brain type for its expression was mainly in the brain. During ovarian development, the expression of cyp19a1a in the ovary significantly increased from stage IV to stage V and Cyp19a1a signals were present in the follicle cells, while cyp19a1b expression in the pituitary gland decreased from stage IV to stage V. To further investigate the function of Cyp19a1a, recombinant Cyp19a1a (rCyp19a1a) was produced and specific anti-Cyp19a1a antiserum was obtained. The expressions of cyp19a1a, estrogen receptors 2 alpha (esr2a), and androgen receptor alpha (arα) were significantly upregulated in the presence of rCyp19a1a. Meanwhile, cyp19a1a was expressed significantly after E2 treatment in both ovarian and testicular tissue culture. Taken together, we found two forms of aromatase in silver pomfret. The ovarian-type aromatase might play an important role in ovarian differentiation and maturation, and participate in E2 synthesis through co-regulation with esr2a. The brain-type aromatase cyp19a1b might be involved in the regulation of both brain and gonadal development.

Transcription factors Dmrt1a, Foxl2, and Nr5a1a potentially interact to regulate cyp19a1a transcription in ovarian follicles of ricefield eel (Monopterus albus)

Aromatase (encoded by Cyp19a1) in the ovarian follicular cells catalyzes the production of estradiol from testosterone, which plays important roles in the ovarian development of vertebrates. In the present study, the interaction of Dmrt1, Foxl2, and Nr5a1a on the regulation of cyp19a1a transcription in ovarian follicles was examined in a teleost, the ricefield eel Monopterus albus. The expression of dmrt1a, foxl2, and nr5a1a was detected in ovarian follicular cells together with cyp19a1a at the mRNA and/or protein levels. Sequence analysis identified one conserved Foxo binding site in the proximal promoter region of ricefield eel cyp19a1a. Transient transfection assay showed that Foxl2 may bind to the conserved Foxo site to activate cyp19a1a transcription and act synergistically with Nr5a1a. Mutation of either the conserved Nr5a1 site or Foxo site abolished or significantly decreased the synergistic effects of Nr5a1a and Foxl2 on cyp19a1a transcription. The sequence between Region III and I-box of Nr5a1a was critical to this synergistic effect. Dmrt1a modulated the Foxl2- and Nr5a1a-induced activation of cyp19a1a transcription and their synergistic effects in a biphasic manner, with inhibitory roles observed at lower doses (10 to 50ng) but release of the inhibition or even potentiating effects observed at higher doses (100 to 200ng). Collectively, data of the present study suggest that the interaction of Dmrt1a, Foxl2, and Nr5a1a in the ovarian follicular cells may facilitate the adequate expression of cyp19a1a and the production of estradiol, and contribute to the development and maturation of ovarian follicles in ricefield eels and other vertebrates as well.

miR-27a-3p targets NR5A2 to regulate CYP19A1 expression and 17-β estradiol synthesis in ovine granulosa cells

Although 17-β estradiol (E2) synthesis is important in regulating female fertility, we know little regarding the molecular mechanism of miRNA-regulated ovine E2 synthesis. Here, our experiments with granulosa cells (GCs) from Hu sheep revealed miR-27a-3p involvement in E2 synthesis and its association with ovine litter size. First, we showed that miR-27a-3p of sheep and other mammals share a high nucleotide identity. Next, gain- and loss-of-function assays indicated that miR-27a-3p inhibits CYP19A1 expression and E2 synthesis in GCs. Moreover, we demonstrated that NR5A2 is a direct target of miR-27a-3p. Ovine miR-27a-3p suppresses E2 synthesis via the NR5A2 and CYP19A1 axes. We also identified four single nucleotide polymorphisms in the ovine miR-27a gene, and g.-13 G>A and g 0.24 T > G were significantly associated with the first and the second parity litter size, respectively (P < 0.05). In summary, our findings reveal that miR-27a-3p is a novel regulator of E2 synthesis and may predict litter size of Hu sheep, providing insight into mechanisms underlying granulosa cell function and female fertility.

Effect of liver receptor homolog-1 on cell apoptosis and steroid hormone secretion on granulosa cells from Hu sheep

The apoptosis of granulosa cells can result in follicular atresia, a key process in follicle selection and development. The related molecular mechanisms were explored to study the effects of liver receptor homolog-1 (Lrh-1) on apoptosis, sodium-titanium system, and steroid synthesis of granulosa cells from Hu sheep in vitro. After constructing, designing, and synthesizing the Lrh-1 overexpression vector, liposomes were used to transfect Hu sheep granulosa cells. Thereafter, qRT-PCR and Western blot were used to detect the mechanism of Lrh-1 on the titanium system and steroid synthesis of Hu sheep granulosa cells. The overexpression efficiency was determined by fluorescence 48 h after liposome transfection into the Hu sheep granulosa cells in vitro culture. The transfection efficiency was higher, hence providing a basis for subsequent experiments. However, the protein level and transcriptional level of Lrh-1 significantly increased after transfection with Lrh-1 overexpression vector in Hu sheep granulosa cells. We further revealed that after Lrh-1 overexpression in Hu sheep granulosa cells, the expression of the pro-apoptotic gene Bax decreased significantly, while that of the anti-apoptotic gene Bcl-2 increased significantly, as well as the sodium peptide system. Moreover, the expression levels of natriuretic peptide precursor A, B, and C (NPPA, NPPB, NPPC) all showed an upward trend, while the expression levels of steroid synthesis-related genes (P450arom and P450scc) decreased. The above results showed that Lrh-1 overexpression influenced the apoptosis, sodium peptide system, and steroid synthesis.

Nr5a homologues in the ricefield eel Monopterus albus: Alternative splicing, tissue-specific expression, and differential roles on the activation of cyp19a1a promoter in vitro

Nr5a (Fushi tarazu factor 1, Ftz-F1) homologues belong to the nuclear receptor superfamily, and are involved in the regulation of reproduction in vertebrates. Four genes encoding Nr5a homologues were present in the genome of ricefield eel, which are designated as nr5a1a, nr5a1b, nr5a2, and nr5a5 in the present study. Alternatively spliced transcripts were identified for nr5a1a and nr5a1b genes. Sequence analysis indicated that nr5a5 is possibly a paralog of nr5a2, and nr5a1b is lost during evolution in some teleosts including tilapia and medaka. Ricefield eel nr5a genes exhibit tissue-specific expression patterns, with nr5a1a and nr5a1b resembling that of the SF-1/Ad4BP (NR5A1) subfamily, and nr5a2 and nr5a5 resembling that of the NR5A2/LRH/FTF subfamily. Transcriptomic analysis revealed parallel expression profiles of nr5a1a, foxl2, and cyp19a1a in ovarian follicles during vitellogenesis, with peak values at the late vitellogenic stage. Real-time PCR indicated that the expression levels of nr5a1a and foxl2 in gonads were decreased significantly during the sexual transition from female to the late intersexual stage. In vitro transient transfection assay showed that Nr5a1a up-regulated ricefield eel cyp19a1a promoter activities synergistically with Foxl2. However, Nr5a1b, Nr5a2, and Nr5a5 could neither activate ricefield eel cyp19a1a promoter alone nor enhance the stimulatory effects of Foxl2 on cyp19a1a promoter activities. Collectively, the above data suggest that Nr5a homologues may have diverse and differential roles in the tissues of ricefield eels. The up-regulation of gonadal nr5a1a and foxl2 during vitellogenesis may be important for the ovarian development whereas their down-regulation during the sexual transition period may be important for the sex change process of ricefield eels, possibly through the regulation of cyp19a1a gene expression.

Potential role of DNA methylation of cyp19a1a promoter during sex change in protogynous orange-spotted grouper, Epinephelus coioides

Estrogen has a pivotal role in early female differentiation and further ovarian development. Aromatase (Cyp19a) is responsible for the conversion of androgens to estrogens in vertebrates. In teleosts, cyp19a1a and it paralog cyp19a1b are mainly expressed in the ovary and hypothalamus, respectively. Decreased plasma estrogen levels and lower cyp19a1a expression are associated with the initiation of female-to-male sex change in protogynous grouper. However, an 17α-methyltestosterone (MT)-induced the sex change from a female to a precocious male is a transient phase, and a reversible sex change (induced male-to-female) occurs after chemical withdrawal. Thus, we used this characteristic to study the epigenetic modification of cyp19a1a promoter in orange-spotted grouper. CpG-rich region with a CpG island is located on the putative regulatory region of distal cyp19a1a promoter. Our results showed that cyp19a1a promoter exhibited tissue-specific methylation status. Low methylation levels of distal cyp19a1a promoter and hypomethylated (0-40%) clones of cyp19a1a promoter region were widely observed in the ovary but not shown in testis and other tissues. In femaleness, higher numbers of hypomethylated clones of cyp19a1a promoter region were observed in the vitellogenic oocyte stage compared to the primary oocyte stage. Furthermore, decreased numbers of hypomethylated clones of cyp19a1a promoter region were associated with the maleness during the female-to-male sex change. DNA methylation inhibitor (5-aza-2'-deoxycytidine) delayed the spermatogenesis process (according to germ cell stage and numbers: by decrease of sperm and increase of spermatocytes) but did not influence the reversed sex change in MT-induced bi-directional sex change. These results suggest that epigenetic modification of cyp19a1a promoter may play an important role during the sex change in orange-spotted grouper.

Nr5a1b promotes and Nr5a2 inhibits transcription of lhb in the orange-spotted grouper, Epinephelus coioides†

Nr5a1 (Sf-1) up-regulates lhb expression across vertebrates; however, its regulatory roles on fshb remain to be defined. Moreover, the involvement of Nr5a2 in the regulation of gonadotropin expression is not clear either. In the present study, the involvement of Nr5a1b (a homologue of Nr5a1) and Nr5a2 in the regulation of lhb and fshb expression in the orange-spotted grouper was examined. Dual fluorescent immunohistochemistry using homologous antisera showed that in the pituitary of orange-spotted groupers, Lh cells contain both immunoreactive Nr5a1b and Nr5a2 signals, whereas Fsh cells contain neither of them. In LβT2 cells, Nr5a1b up-regulated basal activities of lhb and fshb promoters possibly via Nr5a sites, and synergistically (on lhb promoter) or additively (on fshb promoter) with forskolin. Surprisingly, Nr5a2 inhibited basal activities of lhb promoter possibly via Nr5a sites and attenuated the stimulatory effects of both forskolin and Nr5a1b. In contrast, Nr5a2 had no effects on fshb promoter. Chromatin immunoprecipitation analysis showed that both Nr5a1b and Nr5a2 bound to lhb promoter, but not fshb promoter in the pituitary of the orange-spotted grouper. The abundance of Nr5a1b bound to lhb promoter was significantly higher at the vitellogenic stage than the pre-vitellogenic stage, whereas that of Nr5a2 exhibited an opposite trend. Taken together, data of the present study demonstrated antagonistic effects of Nr5a1b and Nr5a2 on lhb transcription in the orange-spotted grouper and revealed novel regulatory mechanisms of differential expression of lhb and fshb genes through Nr5a homologues in vertebrates.

miR-1275 controls granulosa cell apoptosis and estradiol synthesis by impairing LRH-1/CYP19A1 axis

miR-1275 is one of the microRNAs (miRNAs) that are differentially expressed during follicular atresia in pig ovaries, as identified by a miRNA microarray assay in our previous study [1]. However, its functions in follicular atresia remain unknown. In this study, we showed that miR-1275 promotes early apoptosis of porcine granulosa cells (pGCs) and the initiation of follicular atresia, and inhibits E2 release and expression of CYP19A1, the key gene in E2 production. Bioinformatics and luciferase reporter assays revealed that liver receptor homolog (LRH)-1, not CYP19A1, is a direct functional target of miR-1275. In vitro overexpression and knockdown experiments showed that LRH-1 significantly repressed apoptosis and induced E2 secretion and CYP19A1 expression in pGCs. LRH-1, whose expression was regulated by miR-1275, prevented apoptosis in pGCs. Furthermore, luciferase and chromatin immunoprecipitation assays demonstrated that LRH-1 protein bound to the CYP19A1 promoter and increased its activity. Our findings suggest that miR-1275 attenuates LRH-1 expression by directly binding to its 3'UTR. This prevents the interaction of LRH-1 protein with the CYP19A1 promoter, represses E2 synthesis, promotes pGC apoptosis, and initiates follicular atresia in porcine ovaries.

A homologue of Nr5a1 activates cyp19a1a transcription additively with Nr5a2 in ovarian follicular cells of the orange-spotted grouper

Transcription factors of nuclear receptor 5A (Nr5a) subfamily play pivotal roles in regulation of steroidogenic enzymes in vertebrates including teleosts. In the orange-spotted grouper, the expression of Nr5a1a was only detectable in the ovary, spleen, and head kidney in the female. The immunoreactive Nr5a1a was present in ovarian follicular and germ cells. In the ovarian follicular cells surrounding vitellogenic oocytes, Nr5a1a was detected both in the nucleus and cytoplasm, and co-localized with Cyp19a1a and Nr5a2. In the ovarian follicular cells surrounding fully grown oocytes, Nr5a1a was localized almost exclusively to the cytoplasm together with Nr5a2. Nr5a1a could up-regulate cyp19a1a promoter activities through Nr5a sites, and further increase the responses elicited by Nr5a2 at sub-maximal doses. Chromatin immunoprecipitation analysis showed that Nr5a1a bound to cyp19a1a promoter in the vitellogenic but not fully grown ovary. Taken together, Nr5a1a up-regulates cyp19a1a additively with Nr5a2 during vitellogenesis, and its cytoplasmic sequestration may also contribute to the down-regulation of cyp19a1a in the fully grown ovary.

Transcriptomic analysis of the differentiating ovary of the protogynous ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a protogynous hermaphroditic Synbranchiform species that changes sex naturally from female to male, which offers an interesting model for studying gonadal (particularly ovarian) differentiation in vertebrates. In the present study, transcriptome sequencing of the gonad of ricefield eel larvae was performed to explore the molecular mechanisms underlying the ovarian differentiation and development.

RESULTS

A total of 301,267,988 clean reads were generated from cDNA libraries of gonadal tissues of ricefield eel larvae at 6, 9, 12, and 20 days post hatching (dph), which contained undifferentiated gonads, differentiating ovaries, ovaries with oogonia, and ovaries with meiotic oocytes, respectively. De-novo assembly of all the clean reads generated a total of 265,896 unigenes with a mean size of 720 bp and a N50 of 1107 bp. RT-qPCR analysis of the developmental expression of 13 gonadal development-related functional genes indicated that RNA-seq data are reliable. Transcriptome data suggest that high expression of female development-related genes and low expression of male development-related genes in the early gonads of ricefield eel larvae participate in the cascade of sex differentiation leading to the final female phenotype. The contrasting expression patterns of genes involved in retinoid acid (RA) synthesis and degradation might result in peak production of RA at 12 dph in the gonad of ricefield eel larvae, and induce molecular events responsible for the initiation of meiosis before the meiotic signs could be observed at 20 dph. In addition, only stra6 but not stra8 could be identified in gonadal transcriptome data of ricefield eel larvae, and the expression pattern of stra6 paralleled those of genes involved in RA synthesis, suggesting that stra6 may be a downstream target of RA and play a role in RA metabolism and/or meiotic initiation in the gonad of ricefield eel larvae.

CONCLUSIONS

The present study depicted the first large-scale RNA sequencing of the gonad of ricefield eel larvae, and identified many important functional genes, GO terms and KEGG pathways involved in gonadal development and germ cell meiosis. Results of the present study will facilitate future study on the ovarian differentiation of ricefield eels and other teleosts as well.

Foxo3b but not Foxo3a activates cyp19a1a in Epinephelus coioides

FOXO3 has been shown to be a critical transcription factor for folliculogenesis in mammals, while the information on its roles in reproduction of nonmammalian vertebrates remains scarce. In this study, two foxo3 homologs, namely foxo3a and foxo3b, were identified in a teleost, the orange-spotted grouper Epinephelus coioides. foxo3a was mainly expressed in the central nervous system, ovary, and gut whereas foxo3b was expressed ubiquitously in tissues examined. In contrast to the dominant expression of mammalian FOXO3 in germ cells but barely detectable in ovarian follicular cells, immunoreactive Foxo3a and Foxo3b were identified both in the ovarian germ cells and follicular cells. The immunointensities of both Foxo3a and Foxo3b in ovarian follicular cells during vitellogenesis were significantly increased stage-dependently, and co-localized with Cyp19a1a. In the nucleus of ovarian follicular cells, both Foxo3a and Foxo3b immunostaining could be detected at the vitellogenic stages. Transient transfection and EMSA showed that Foxo3a and Foxo3b upregulated cyp19a1a promoter activities in vitro through a conserved Foxo-binding site, with the latter being a more potent activator. However, ChIP analysis showed that only Foxo3b binds to cyp19a1a proximal promoter region containing the conserved Foxo-binding site in the vitellogenic ovary. Taken together, these results suggested that Foxo3a and Foxo3b are involved in the ovarian development possibly through regulating the ovarian germ cells as well as follicular cells, and Foxo3b but not Foxo3a may activate cyp19a1a in the ovarian follicular cells during vitellogenesis in the orange-spotted grouper.

Changes in ovarian gene expression profiles and plasma hormone levels in maturing European eel (Anguilla anguilla); Biomarkers for broodstock selection

Complete sexual maturation of European eels (Anguilla anguilla) in captivity can only be achieved via injections with gonadotropins. For female eels this procedure takes 4-6months and the response ranges from "unresponsive" to final maturation and ovulation. Reproductive success could be significantly increased via early selection of responders based on predictive markers and minimally invasive sampling methods. To get a better understanding of the genetic background of ovarian maturation of the European eel we performed a pilot deep-sequencing transcriptome analysis of ovarian tissue derived from a yellow eel, a prepubertal silver eel and a post-spawning matured eel. Two key players in steroidogenesis were strongly correlated with advanced sexual maturation, namely P450c17 and liver receptor homolog-1, suggesting that blood plasma steroids might qualify as minimally invasive markers for early detection of responders. Since the predictive value of plasma sex steroid levels for final maturation of the European eel had not yet been carefully examined, we performed an extensive artificial maturation trial. Farmed silver eels were treated with pituitary extracts and sampled at multiple time intervals. Expression of steroidogenesis-related genes in ovarian tissue of responding and non-responding eels after four weekly injections with pituitary extract was compared using a custom-built microarray and RNAseq. Increased expression of 17β-hsd1 was strongly linked to sexual maturation. Blood plasma levels of sex steroids were measured using ELISAs. We show that a 2.5-fold increase in blood-plasma estradiol level after 4 weekly pituitary extract injections is a strong predictor of final sexual maturation of female European eel.