Growth differentiation factor 9 (gdf9) and bone morphogenetic protein 15 (bmp15) are potential intraovarian regulators of steroidogenesis in Japanese flounder (Paralichthys olivaceus)

Full-length transcriptome analysis of male and female gonads in Japanese Eel (Anguilla japonica)

BACKGROUND

The Japanese eel (Anguilla japonica) holds significant economic value in East Asia, but limitations in understanding its reproductive biology have hindered advancements in artificial breeding techniques. Previous research has primarily focused on conserved sex differentiation genes, offering limited insights into the broader molecular mechanisms driving gonadal development and sexual dimorphism. To address these limitations, this study aims to investigate key genes and pathways involved in gonadal development through a comprehensive transcriptomic analysis of male and female eel gonads.

RESULTS

PacBio Iso-Seq and Illumina RNA-Seq technologies were combined to conduct a full-length transcriptome analysis of male and female Japanese eel gonads at a post-differentiation, pre-maturation stage. A total of 24,661 unigenes were identified in ovaries and 15,023 in testes, along with genomic regulatory elements such as transcription factors, simple sequence repeats, and long non-coding RNAs. Additionally, 1,210 differentially expressed genes were detected. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed significant pathways involved in cell cycle regulation, metabolic processes, apoptosis, and hormone activity. Notably, several reproductive-related genes, including bambi, ccnb1, cdc20, gdf9, prlh, ccdc39, chrebp, tspo, syce3, and ngb, demonstrated significant dimorphic expression in eel gonads.

CONCLUSIONS

This study provides valuable insights into the molecular mechanisms of gonadal differentiation and sexual dimorphism in Japanese eels. The findings expand the genetic resources available for the eel breeding industry and could facilitate the development of improved artificial breeding techniques focused on reproductive development.

In Silico Analysis: Genome-Wide Identification, Characterization and Evolutionary Adaptations of Bone Morphogenetic Protein (BMP) Gene Family in Homo sapiens

We systematically analyzed BMP gene family in H. sapiens to elucidate genetic structure, phylogenetic relationships, adaptive evolution and tissue-specific expression pattern. Total of 13 BMPs genes were identified in the H. sapiens genome. Bone morphogenetic proteins (BMPs) are composed of a variable number of exons ranging from 2 to 21. They exhibit a molecular weight ranging from 31,081.81 to 82,899.61 Da. These proteins possess hydrophilic characteristics, display thermostability, and exhibit a pH range from acidic to basic. We identified four segmental and two tandem duplication events in BMP gene family of H. sapiens. All of the vertebrate species that were studied show the presence of BMPs 1, 2, 3, 4, 5, 6, 7, 8A, and 15, however only Homo sapiens demonstrated the presence of BMP9 and BMP11. The pathway and process enrichment analysis of BMPs genes showed that these were considerably enriched in positive regulation of pathway-restricted SMAD protein phosphorylation (92%) and cartilage development (77%) biological processes. These genes exhibited positive selection signals that were shown to be conserved across vertebrate lineages. The results showed that BMP2/3/5/6/8a/15 proteins underwent adaptive selection at many amino acid locations and increased positive selection was detected in TGF-β propeptide and TGF-β super family domains which were involved in dorso-ventral patterning, limb bud development. More over the expression pattern of BMP genes revealed that BMP1 and BMP5; BMP4 and BMP6 exhibited substantially identical expression patterns in all tissues while BMP10, BMP15, and BMP3 showed tissue-specific expression.

Transforming growth factor-β (TGF-β): A master signal pathway in teleost sex determination

Sex determination and differentiation in fish has always been a hot topic in genetic breeding of aquatic animals. With the advances in next-generation sequencing (NGS) in recent years, sex chromosomes and sex determining genes can be efficiently identified in teleosts. To date, master sex determination genes have been elucidated in 114 species, of which 72 species have sex determination genes belonging to TGF-β superfamily. TGF-β is the only signaling pathway that the largest proportion of components, which including ligands (amhy, gsdfy, gdf6), receptors (amhr, bmpr), and regulator (id2bby), have opportunity recognized as a sex determination gene. In this review, we focus on the recent studies about teleost sex-determination genes within TGF-β superfamily and propose several hypotheses on how these genes regulate sex determination process. Differing from other reviews, our review specifically devotes significant attention to all members of the TGF-β signal pathway, not solely the sex determination genes within the TGF-β superfamily. However, the functions of the paralogous genes of TGF superfamily are still needed ongoing research. Further studies are required to more accurately interpret the molecular mechanism of TGF-β superfamily sex determination genes.

Loss of bmp15 function in the seasonal spawner Atlantic salmon results in ovulatory failure

Bone morphogenetic protein 15 (BMP15) is an oocyte-specific growth factor important for successful female reproduction in mammals. While mutations in BMP15/Bmp15 cause ovulatory deficiency and/or infertility in certain mammalian species, loss of bmp15 in zebrafish, a continuous spawner and the only bmp15 knockout model in fish to date, results in complete arrest of follicle development and later female-to-male sex reversal, preventing to examine effects on ovulation/fertilization. Here, we used Atlantic salmon, a seasonal spawner, and generated bmp15 mutants to investigate ovarian development and fertility. Histological and morphometric analyses revealed that in biallelic frameshift (bmp15 fs/fs) mutant ovaries, folliculogenesis started earlier, resulting in an advanced development compared to wild-type (WT) controls, accompanied by a weaker expression of the (early) oocyte-specific factor figla. This precocious ovarian development was followed in bmp15 fs/fs females by enhanced follicle atresia during vitellogenic stages. Although genes involved in steroid synthesis and signaling (star, cyp11b, cyp17a1 and esr1) were dramatically higher in late vitellogenic bmp15 fs/fs mutant ovaries, estradiol-17β plasma levels were lower than in WT counterparts, potentially reflecting compensatory changes at the level of ovarian gene expression. At spawning, bmp15 fs/fs females displayed lower gonado-somatic index values and reduced oocyte diameter, and the majority (71.4%), showed mature non-ovulating ovaries with a high degree of atresia. The remaining (28.6%) females spawned eggs but they either could not be fertilized or, upon fertilization, showed severe malformations and embryonic mortality. Our results show that Bmp15 is required for proper follicle recruitment and growth and later ovulatory success in Atlantic salmon, providing an alternative candidate target to induce sterility in farmed salmon. Moreover, since loss of bmp15 in salmon, in contrast to zebrafish, does not result in female-to-male sex change, this is the first mutant model in fish allowing further investigations on Bmp15-mediated functions in the ovulatory period.

Exploring the molecular mechanism of Chlorella vulgaris in response to androstenedione exposure based on genes continuously up-regulated in transcription analysis

Androstenedione (ADSD) is one of the widely detected androgens in diverse aquatic environments. However, there were few reports on the molecular mechanism of Chlorella vulgaris exposure to ADSD. In our previous research, we have investigated the genes associated with chlorophyll metabolism in Chlorella vulgaris response to ADSD. In this study, we focus on continuously up-regulated genes to explore the mechanism underlying Chlorella vulgaris resistance to ADSD toxicity. Chlorella vulgaris was exposed to ADSD with five concentration gradients. The continuously up-regulated genes were enriched by Series Test of Cluster (STC) analysis and verified by qRT-PCR. Microalgae Super Oxidase Dimutase (SOD) and Microalgae Malonic dialdehyde (MDA), two indicators of oxidative stress, were determined by ELISA after exposure to ADSD. The results showed that ADSD can stimulate the production of extracellular polymeric substances (EPS) and lead to enlargement in the cell body of Chlorella vulgaris. In addition, steroid biosynthesis and oxidoreductase activity processes were consistently up-regulated upon exposure to ADSD. In conclusion, our study highlighted the crucial role of phenotypic modification, hormone synthesis, and redox mechanisms in protecting Chlorella vulgaris cells from the harmful effects of ADSD contamination.

Growth factors and female reproduction in vertebrates

Female reproductive efficiency is influenced by the outcomes of various processes, including folliculogenesis, apoptosis, response to gonadotropin signaling, oocyte maturation, and ovulation. The role of hormones in regulating these processes and other reproductive activities has been well established. It is becoming increasingly evident that in addition to well-characterized hormones, growth factors play vital roles in regulating some of these reproductive activities. Growth factors and their receptors are widely distributed in vertebrate ovaries at different stages of ovarian development, indicating their involvement in intraovarian reproductive functions. In the ovary, cell surface receptors allow growth factors to regulate intraovarian reproductive activities. Understanding these actions in the reproductive axis would provide a tool to target growth factors and/or their receptors to yield desirable reproductive outcomes. These include enrichment of in vitro maturation and fertilization culture media, and management of infertility. This review discusses some widely characterized growth factors belonging to the TGF, EGF, IGF, FGF, and BDNF family of peptides and their role in female reproduction in vertebrates, with a focus on mammals.

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

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

Rescue of bmp15 deficiency in zebrafish by mutation of inha reveals mechanisms of BMP15 regulation of folliculogenesis

As an oocyte-specific growth factor, bone morphogenetic protein 15 (BMP15) plays a critical role in controlling folliculogenesis. However, the mechanism of BMP15 action remains elusive. Using zebrafish as the model, we created a bmp15 mutant using CRISPR/Cas9 and demonstrated that bmp15 deficiency caused a significant delay in follicle activation and puberty onset followed by a complete arrest of follicle development at previtellogenic (PV) stage without yolk accumulation. The mutant females eventually underwent female-to-male sex reversal to become functional males, which was accompanied by a series of changes in secondary sexual characteristics. Interestingly, the blockade of folliculogenesis and sex reversal in bmp15 mutant could be partially rescued by the loss of inhibin (inha-/-). The follicles of double mutant (bmp15-/-;inha-/-) could progress to mid-vitellogenic (MV) stage with yolk accumulation and the fish maintained their femaleness without sex reversal. Transcriptome analysis revealed up-regulation of pathways related to TGF-β signaling and endocytosis in the double mutant follicles. Interestingly, the expression of inhibin/activin βAa subunit (inhbaa) increased significantly in the double mutant ovary. Further knockout of inhbaa in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-) resulted in the loss of yolk granules again. The serum levels of estradiol (E2) and vitellogenin (Vtg) both decreased significantly in bmp15 single mutant females (bmp15-/-), returned to normal in the double mutant (bmp15-/-;inha-/-), but reduced again significantly in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-). E2 treatment could rescue the arrested follicles in bmp15-/-, and fadrozole (a nonsteroidal aromatase inhibitor) treatment blocked yolk accumulation in bmp15-/-;inha-/- fish. The loss of inhbaa also caused a reduction of Vtg receptor-like molecules (e.g., lrp1ab and lrp2a). In summary, the present study provided comprehensive genetic evidence that Bmp15 acts together with the activin-inhibin system in the follicle to control E2 production from the follicle, Vtg biosynthesis in the liver and its uptake by the developing oocytes.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

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.

Genome-wide identification, evolution and expression analysis of bone morphogenetic protein (BMP) gene family in chinese soft-shell turtle (Pelodiscus sinensis)

Introduction: Bone morphogenetic proteins (BMPs) play a crucial role in bone formation and differentiation. Recent RNA-Seq results suggest that BMPs may be involved in the sex differentiation of P. sinensis, yet more relevant studies about BMPs in P. sinensis are lacking. Methods: Herein, we identified BMP gene family members, analyzed the phylogeny, collinear relationship, scaffold localization, gene structures, protein structures, transcription factors and dimorphic expression by using bioinformatic methods based on genomic and transcriptomic data of P. sinensis. Meanwhile, qRT-PCR was used to verify the RNA-Seq results and initially explore the function of the BMPs in the sex differentiation of P. sinensis. Results: A total of 11 BMP genes were identified, 10 of which were localized to their respective genomic scaffolds. Phylogenetic analysis revealed that BMP genes were divided into eight subfamilies and shared similar motifs ("WII", "FPL", "TNHA", "CCVP", and "CGC") and domain (TGF-β superfamily). The results of the sexually dimorphic expression profile and qRT-PCR showed that Bmp2, Bmp3, Bmp15l, Bmp5, Bmp6 and Bmp8a were significantly upregulated in ovaries, while Bmp2lb, Bmp7, Bmp2bl and Bmp10 were remarkable upregulated in testes, suggesting that these genes may play a role in sex differentiation of P. sinensis. Discussion: Collectively, our comprehensive results enrich the basic date for studying the evolution and functions of BMP genes in P. sinensis.

Rxrs and their partner receptor genes inducing masculinization plausibly mediated by endocrine disruption in Paralichthys olivaceus

Retinoid X receptors (RXRs) can form homo- or heterodimers with orphan receptors involved in multiple intertwined signaling pathways. However, there is limited study on the formation of sex phenotypes and the regulation of steroidogenesis by RXRs in fish. Here, in Paralichthys olivaceus, we first indicated that PPARγ::RXRα was predictably a transcription factor for steroidogenesis genes, and Foxl2 and Dmrt1 were also transcription factors for rxrs and their partner receptor genes. When the flounder fry were exposed to LG100268 (LG, RXRs agonist, 50 mg/kg diet), the percentage of males increased from 50% to 71.4%. Before histological differentiation of the flounder ovary (3 cm TL) and testis (6 cm TL), the transcripts of rar β and rar γ (P < 0.05) were activated, and the steroidogenesis gene Hsd3b1 was down-regulated (P < 0.05). The ratios of testosterone (T)/17β-estradiol (E2) were all greatly increased (P < 0.05), and the ratio of 11-ketotestosterone (11-KT)/E2 was elevated at 3 cm TL. Moreover, LG was used to treat the cultured gonads in vitro (10 μM) and the fish with intraperitoneal injection in vivo (12 mg/kg body weight), respectively. LG was able to up-regulate rxr γ, rar γ, and ppar δ, and Hsd3b1 was significantly up-regulated (P < 0.05). The ratios of 11-KT/E2 in the culture medium and in the ovaries of the fish were decreased. Furthermore, the recombinant flounder Foxl2 protein was able to significantly down-regulate ppar γ (P < 0.05) and tr β (P < 0.01) in the ovaries in vitro, and the result of the Dmrt1 in the testes was opposite to that of the Foxl2, probably indicating a feedback loop between RXRs' partner receptors and Foxl2/Dmrt1. These findings introduce for the first time the mode of action of RXRs on the flounder steroidogenesis and provide important data to learn the potential function of RXRs in fish sex differentiation and the potential role of RXRs in aquatic animals in the presence of water pollutants.

Dynamic changes of urotensin II and its receptor during ovarian development of olive flounder Paralichthys olivaceus

Urotensin II (UII) is a kind of fish somatostatins cyclic peptide, which was originally extracted from the caudal neurosecretory system (CNSS). The system of UII and UII receptor (UIIR) has been reported to have multiple physiological regulatory functions, such as cardiovascular control, osmoregulation, and lipid metabolism. However, the effect of UII and UIIR on the ovarian development has not been covered. This study investigated the expression pattern of UII and UIIR in the ovarian follicles and explored their impact on ovarian development in olive flounder Paralichthys olivaceus. The results showed that the highest UII and UIIR mRNA levels were observed at stage II and stage III follicles during ovarian development, respectively. In situ hybridization revealed that a strong signal of UII was expressed in the oocyte nuclei of stage II follicles, however, UIIR was found in the follicle cells and oocyte cytoplasm of stage II and stage III follicles. Similarly, immunohistochemistry found positive signal of UII was detected in the oocyte nuclei of stage II follicles. The results from in vitro culture of olive flounder follicles suggested the expression of UII and UIIR mRNA levels significantly increased by 10 IU/ml human chorionic gonadotropin (hCG) for 9 h. Furthermore, the transcriptional expression of UII and UIIR was not statistically significantly changed by 17α, 20β-dihydroxy-4-pregnen-3-one (DHP). These results firstly suggested that UII and UII receptor may play vital roles in regulating ovarian growth in olive flounder.

Loss of growth differentiation factor 9 causes an arrest of early folliculogenesis in zebrafish-A novel insight into its action mechanism

Growth differentiation factor 9 (GDF9) was the first oocyte-specific growth factor identified; however, most information about GDF9 functions comes from studies in the mouse model. In this study, we created a mutant for Gdf9 gene (gdf9-/-) in zebrafish using TALEN approach. The loss of Gdf9 caused a complete arrest of follicle development at primary growth (PG) stage. These follicles eventually degenerated, and all mutant females gradually changed to males through sex reversal, which could be prevented by mutation of the male-promoting gene dmrt1. Interestingly, the phenotypes of gdf9-/- could be rescued by simultaneous mutation of inhibin α (inha-/-) but not estradiol treatment, suggesting a potential role for the activin-inhibin system or its signaling pathway in Gdf9 actions. In gdf9-null follicles, the expression of activin βAa (inhbaa), but not βAb (inhbab) and βB (inhbb), decreased dramatically; however, its expression rebounded in the double mutant (gdf9-/-;inha-/-). These results indicate clearly that the activation of PG follicles to enter the secondary growth (SG) requires intrinsic factors from the oocyte, such as Gdf9, which in turn works on the neighboring follicle cells to trigger follicle activation, probably involving activins. In addition, our data also support the view that estrogens are not involved in follicle activation as recently reported.

Expression Patterns and Gonadotropin Regulation of the TGF-β II Receptor (Bmpr2) during Ovarian Development in the Ricefield Eel Monopterus albus

Bmpr2 plays a central role in the regulation of reproductive development in mammals, but its role during ovarian development in fish is still unclear. To ascertain the function of bmpr2 in ovarian development in the ricefield eel, we isolated and characterized the bmpr2 cDNA sequence; the localization of Bmpr2 protein was determined by immunohistochemical staining; and the expression patterns of bmpr2 in ovarian tissue incubated with FSH and hCG in vitro were analyzed. The full-length bmpr2 cDNA was 3311 bp, with 1061 amino acids encoded. Compared to other tissues, bmpr2 was abundantly expressed in the ovary and highly expressed in the early yolk accumulation (EV) stages of the ovary. In addition, a positive signal for Bmpr2 was detected in the cytoplasm of oocytes in primary growth (PG) and EV stages. In vitro, the expression level of gdf9, the ligand of bmpr2, in the 10 ng/mL FSH treatment group was significantly higher after incubation for 4 h than after incubation for different durations. However, bmpr2 expression in the 10 ng/mL FSH treatment group at 2 h, 4 h and 10 h was significantly lower. Importantly, the expression level of bmpr2 and gdf9 in the 100 IU/mL hCG group had similar changes that were significantly decreased at 4 h and 10 h. In summary, Bmpr2 might play a pivotal role in ovarian growth in the ricefield eel, and these results provide a better understanding of the function of bmpr2 in ovarian development and the basic data for further exploration of the regulatory mechanism of gdf9 in oocyte development.

Proteomic analysis implicates that postovulatory aging leads to aberrant gene expression, biosynthesis, RNA metabolism and cell cycle in mouse oocytes

Background

In mammals, oocytes display compromised quality after experiencing a process of postovulatory aging. However, the mechanisms underlying are not yet fully understood. Here, we portrayed a protein expression profile of fresh and aging metaphase II (MII) mouse oocytes by means of four-dimensional label-free quantification mass spectrometry (4D-LFQ).

Results

The analysis of 4D-LFQ data illustrated that there were seventy-six differentially expressed proteins (DEPs) between two groups of MII stage oocytes. Fifty-three DEPs were up-regulated while twenty-three DEPs were down-regulated in the MII oocytes of the aging group, and Gene Ontology (GO) analysis revealed that these DEPs were mainly enriched in regulation of gene expression, biosynthesis, RNA metabolism and cell cycle. Our detailed analysis revealed that the expression of proteins that related to gene expression processes such as transcription, translation, post-translational modifications and epigenome was changed; the relative protein expression of RNA metabolic processes, such as RNA alternative splicing, RNA export from nucleus and negative regulation of transcription from RNA polymerase II promoter was also altered.

Conclusion

In conclusion, we identified considerable DEPs and discussed how they agreed with previous researches illustrating altered protein expression associated with the quality of oocytes. Our research provided a new perspective on the mechanisms of postovulatory aging and established a theoretical support for practical methods to control and reverse postovulatory aging.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-022-01045-6.

Transcriptome analysis revealed gene expression feminization of testis after exogenous tetrodotoxin administration in pufferfish Takifugu flavidus

Tetrodotoxin (TTX) is a deadly neurotoxin and usually accumulates in large amounts in the ovaries but is non-toxic or low toxic in the testis of pufferfish. The molecular mechanism underlying sexual dimorphism accumulation of TTX in ovary and testis, and the relationship between TTX accumulation with sex related genes expression remain largely unknown. The present study investigated the effects of exogenous TTX treatment on Takifugu flavidus. The results demonstrated that exogenous TTX administration significantly incresed level of TTX concentration in kidney, cholecyst, skin, liver, heart, muscle, ovary and testis of the treatment group (TG) than that of the control group (CG). Transcriptome sequencing and analysis were performed to study differential expression profiles of mRNA and piRNA after TTX administration of the ovary and testis. The results showed that compared with female control group (FCG) and male control group (MCG), TTX administration resulted in 80 and 23 piRNAs, 126 and 223 genes up and down regulated expression in female TTX-treated group (FTG), meanwhile, 286 and 223 piRNAs, 2 and 443 genes up and down regulated expression in male TTX-treated group (MTG). The female dominant genes cyp19a1, gdf9 and foxl2 were found to be up-regulated in MTG. The cyp19a1, whose corresponding target piRNA uniq_554482 was identified as down-regulated in the MTG, indicating the gene expression feminization in testis after exogenous TTX administration. The KEGG enrichment analysis revealed that differentially expressed genes (DEGs) and piRNAs (DEpiRNAs) in MTG vs MCG group were more enriched in metabolism pathways, indicating that the testis produced more metabolic pathways in response to exogenous TTX, which might be a reason for the sexual dimorphism of TTX distribution in gonads. In addition, TdT-mediated dUTP-biotin nick end labeling staining showed that significant apoptosis was detected in the MTG testis, and the role of the cell apoptotic pathways was further confirmed. Overall, our research revealed that the response of the ovary and testis to TTX administration was largely different, the ovary is more tolerant whereas the testis is more sensitive to TTX. These data will deepen our understanding on the accumulation of TTX sexual dimorphism in Takifugu.

Weighted gene Co-expression network analysis (WGCNA) reveals a set of hub genes related to chlorophyll metabolism process in chlorella (Chlorella vulgaris) response androstenedione

Androstenedione (ADSD) was the main androgen detected in wastewaters. Chlorella was the most widely used plant in biological wastewater treatment process. In order to understand the toxicological response of chlorella to ADSD contamination, we used the weighted gene co-expression network analysis (WGCNA) method to systematically analyze the gene regulatory networks of chlorella after ADSD treatments. Total of 25 modules was identified from gene co-expression networks, and the turquoise module were selected for GO and KEGG enrichment analysis. Results showed that most hub genes were associated with chloroplast organizations or photosystems processes. Among them, the expressions profiles of hcar, nol, pao and sgr genes were highly correlated to the content fluctuations of chlorophylls after different ADSD treatments. All these results demonstrated that chlorophylls play a key role in preventing cell damage of chlorella caused by ADSD contamination. Besides, we proposed a possible chlorophyll metabolism pathway in chlorella response to ADSD contamination.

Functional differentiation of two lhx8 paralogs and possible regulatory role of lhx8a in Japanese flounder (Paralichthys olivaceus)

Lhx8, belonging to the LIM-Homebox family, is involved in the tooth, nervous system, and primordial follicles development in mammals. However, little is known about the regulatory roles of lhx8 in teleosts. In this study, two lhx8 duplicates were identified in Paralichthys olivaceus, termed Polhx8a and Polhx8b, respectively. Bioinformatic analysis showed that Polhx8a was more likely to be a teleost-specific paralog. According to expression analysis, Polhx8a transcripts were almost exclusively concentrated in the oocytes, while Polhx8b was weakly expressed in the spleen, gill, and some facial organs, indicating sub-functionalization of this gene pair during evolution. Furthermore, Polhx8a mRNA level elevated from perinucleolar oocyte (PNO) stage to vitellogenic oocyte (VO) stage transition and changed after exogenous hormone stimulation, proving that Polhx8a was involved in the oocyte development and could be regulated by sex hormones. Yeast two-hybrid, bimolecular fluorescence complementation (BiFC) and co-immunoprecipitation (co-IP) experiments captured the positive protein interactions between PoLhx8a and the other two oocyte-specific transcription factors: PoFigla and PoNobox. After knocking down lhx8a in embryos or adult ovaries in vivo, the expression of oocyte-associated genes was significantly down-regulated (P < 0.05). Our findings suggest the evolution and functional differentiation of lhx8 genes, and shed light on the potential role of lhx8a in protein interactions and gene regulation in teleosts.

The emerging role of microRNAs in fish ovary: A mini review

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily at the post-transcriptional levels. It is now well established that miRNAs are crucial regulators of many developmental and physiological processes, including reproduction. In teleosts, expression profiling studies have shown that miRNAs are expressed in the fish ovary and their levels are regulated during follicle development and by hormones. Using CRISPR/Cas9 mediated gene knockout strategies, several recent studies have provided strong evidence that miR-202 and miR-200s on chromosome 23 play critical roles in regulating ovarian development, oogenesis, and ovulation. In this mini review, we provide a brief overview of canonical miRNA biogenesis and functions; summarize miRNAs that are expressed in fish ovary; and discuss the emerging role of miRNAs in regulating fish ovarian functions.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy

Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

Exosomes as a Potential Tool for Supporting Canine Oocyte Development

The canine oviduct is a unique reproductive organ where the ovulated immature oocytes complete their maturation, while the other mammals ovulate matured gametes. Due to their peculiar reproductive characteristics, the in vitro maturation of dog oocytes is still not wellestablished compared with other mammals. Investigations of the microenvironment conditions in the oviductal canal are required to establish a reliable in vitro maturation system in the dog. Previous studies have suggested that the oviduct and its derivatives play a key role in improving fertilization as well as embryo development. In particular, the biological function of oviduct-derived exosomes on sperm and early embryo development has been investigated in porcine, bovine, and murine species. However, the information about their functions on canine cumulus-oocyte complexes is still elusive. Recent canine reproductive studies demonstrated how oviduct-derived extracellular vesicles such as microvesicles and exosomes interact with oocyte-cumulus complexes and how they can play roles in regulating canine cumulus/oocyte communications. In this review, we summarize the physiological characteristics of canine oviduct-derived exosomes and their potential effects on cumulus cells development as well as oocyte in vitro maturation via molecular signaling pathways.

Spatiotemporal expression of activin receptor-like kinase-5 and bone morphogenetic protein receptor type II in the ovary of shortfinned eel, Anguilla australis

In the eel ovary, the expression of growth differentiation factor-9 (Gdf9) appears to be largely confined to the germ cell in early stages of oogenesis. However, both the target tissue and the function of Gdf9 in fish remain unknown. This study aimed to describe the abundance and localization of activin receptor-like kinase-5 (Alk5) and bone morphogenetic protein receptor type II (Bmpr2), which together mediate the Gdf9 signal, in the ovary of a basal teleost, the shortfinned eel, Anguilla australis, during early folliculogenesis. The cDNA encoding eel alk5 and bmpr2 genes were cloned, characterized and the transcript abundances of these receptors quantified by quantitative real-time PCR. Ovarian transcript abundance for both receptors, along with that of gdf9 and of its paralogue bmp15, increased from the previtellogenic to early vitellogenic stage. Localization of receptor mRNAs by in situ hybridization revealed that these receptors are located in the somatic cells surrounding the oocyte. Furthermore, tissue distribution analysis showed that the expression of alk5 and bmpr2 were highest in ovary and thyroid, respectively. Unexpectedly, however, bmpr2 mRNA levels were lower in the ovary than in any of the other 17 tissues examined, and indeed, lower than ovarian gdf9 transcript abundance. These findings, together with the ovarian expression pattern of Gdf9, suggest that Gdf9, and conceivably, Bmp15, from the oocyte can signal through receptors that are located on the somatic cells surrounding the oocyte; this, in turn, facilitates elucidation of the function of these growth factors during oogenesis in teleost fish.