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

The potential regulatory role of the non-coding RNAs in regulating the exogenous estrogen-induced feminization in Takifugu rubripes gonad

Estrogen plays a pivotal role in the early stage of sex differentiation in teleost. However, the underlying mechanisms of estrogen-induced feminization process are still needed for further clarification. Here, the comparative analysis of whole-transcriptome RNA sequencing was conducted between 17beta-Estradiol induced feminized XY (E-XY) gonads and control gonads (C) in Takifugu rubripes. A total of 57 miRNAs, 65 lncRNAs, and 4 circRNAs were found to be expressed at lower levels in control-XY (C-XY) than that in control-XX (C-XX), and were up-regulated in XY during E2-induced feminization process. The expression levels of 24 miRNAs, and 55 lncRNAs were higher in C-XY than that in C-XX, and were down-regulated in E2-treated XY. Furthermore, a correlation analysis was performed between miRNA-seq and mRNA-seq data. In C-XX/C-XY, 114 differential expression (DE) miRNAs were predicted to target to 904 differential expression genes (DEGs), while in C-XY/E-XY, 226 DEmiRNAs were predicted to target to 2,048 DEGs. In C-XX/C-XY, and C-XY/E-XY, KEGG pathway enrichment analysis showed that those targeted genes were mainly enriched in MAPK signaling, calcium signaling, steroid hormone biosynthesis and ovarian steroidogenesis pathway. Additionally, the competitive endogenous RNA (ceRNA) regulatory network was constructed by 24 miRNAs, 21 lncRNAs, 4 circRNAs and 5 key sex-related genes. These findings suggested that the expression of critical genes in sex differentiation were altered in E2-treated XY T. rubripes may via the lncRNA-miRNA-mRNA regulation network to facilitate the differentiation and maintenance of ovaries. Our results provide a new insight into the comprehensive understanding of the effects of estrogen signaling pathways on sex differentiation in teleost gonads.

Comparable to 17α- methyl testosterone, dietary supplements of Tribulus terrestris and Mucuna pruriens promote the development of mono-sex, all-male tilapia fry, growth, survival rate and sex-related genes (Amh, Sox9, Foxl2, Dmrt1)

To evaluate Tribulus terrestris and Mucuna pruriens for inducing all-male tilapia, mixed-sex Nile tilapia, Oreochromis niloticus, (mean weight 0.025 ± 0.009 g; mean length 1.25 ± 0.012 cm), were given a meal supplemented with either T. terrestris powder (commercial fish feed, 40% crude protein) (TT group), M. pruriens seed extract (MP group), MP + TT (mixed group), 17α-methyl testosterone (MT, control positive), or without supplements (control negative). The MP extracts significantly increased (P < 0.05) the final weight, weight gain, weight gain rate, and specific growth rate while feed conversion ratio was significantly decreased (P < 0.05). Plant extracts markedly improved (P < 0.05) the survival rate, proportion of males, and total testosterone compared to control and MT. Estrogen levels were lower in groups with plant extract than other groups. Fifteen days post-feeding, the Amh gene was expressed in the brain of O. niloticus fries with higher levels in MP, TT, and MT groups. Additionally, the expression of the Sox9 and Dmrt1 genes as a male related genes in fish fry gonads revealed significantly (P < 0.05) higher levels in groups fed on MP, TT, and MT compared to control after 30-day post-feeding, whereas; Foxl2 gene expression as a female related gene was significantly (P < 0.05) lower in fish fed on MP, TT, and MT compared to other groups after 30 days post feeding. Histologically, MT, MP, TT, and the mixture all exhibited solely male reproductive traits without noticeable abnormalities. This study concluded that each of the TT or MP extracts can induce sex reversal in tilapia while having no negative health impact compared to MT as the growth and survival rate in the treated groups with TT and MP were higher than control and group treated with MT.

Sex-biased gene expression precedes sexual dimorphism in the agonadal annelid Platynereis dumerilii

Gametogenesis is the process by which germ cells differentiate into mature sperm and oocytes, cells essential for sexual reproduction. The sex-specific molecular programs that drive spermatogenesis and oogenesis can also serve as sex identification markers. Platynereis dumerilii is a research organism that has been studied in many areas of developmental biology. However investigations often disregard sex, as P. dumerilii juveniles lack sexual dimorphism. The molecular mechanisms of gametogenesis in the segmented worm P. dumerilii are also largely unknown. In this study, we used RNA sequencing to investigate the transcriptomic profiles of gametogenesis in P. dumerilii juveniles. Our analysis revealed that sex-biased gene expression becomes increasingly pronounced during the advanced developmental stages, particularly during the meiotic phases of gametogenesis. We identified conserved genes associated with spermatogenesis, such as dmrt1, and a novel gene psmt, that is associated with oogenesis. Additionally, putative long non-coding RNAs were upregulated in both male and female gametogenic programs. This study provides a foundational resource for germ cell research in P. dumerilii, markers for sex identification, and offers comparative data to enhance our understanding of the evolution of gametogenesis mechanisms across species.

Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.

Influence of Bisphenol A and 17β-Trenbolone Exposure in Oryzias Congeners

Japanese medaka is specified as a model fish in the test guidelines of the Organisation for Economic Co-operation and Development. Recently, populations of Japanese medaka in Japan were divided into two species, the northern Oryzias sakaizumii and the southern O. latipes. Previously, we reported that induction concentrations for sex reversal by exposure to 17α-methyltestosterone differed significantly between these two species, indicating that they respond differently to endocrine-disrupting chemica. In the present study, we examined the effects of exposure to two more endocrine-disrupting chemicals (bisphenol A and 17β-trenbolone) in O. sakaizumii, and compared the results with those previously reported for O. latipes. Exposure to both bisphenol A and 17β-trenbolone induced testis-ova formation or sex reversal in O. sakaizumii. Exposure to 17β-trenbolone also increased expression of gonadal soma-derived factor (gsdf). Least-observed-effect concentrations for gonadal sex differentiation and gsdf expression were lower for O. latipes than for O. sakaizumii after exposure to bisphenol A, and were lower for O. sakaizumii than for O. latipes after exposure to 17β-trenbolone. These results demonstrate that O. sakaizumii and O. latipes respond differently to androgenic and estrogenic endocrine-disrupting chemicals. Environ Toxicol Chem 2023;42:673-678. © 2022 SETAC.

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

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

Comparative gonad transcriptome analysis in cobia (Rachycentron canadum)

Background: Cobia (Rachycentron canadum) is a species of fish with high commercial potential particularly due to fast growth rates. The evidence of sexual size dimorphism favoring females indicate potential benefits in having a monosex culture. However, the involvement of genetic factors responsible for sexual development and gonadal maintenance that produces phenotypic sex in cobia is largely unknown. Methods: In the present study, we performed transcriptome sequencing of cobia to identify sex-biased significantly differentially expressed genes (DEGs) in testes and ovaries. The reliability of the gonad transcriptome data was validated by qPCR analysis of eight selected significantly differential expressed sex-related candidate genes. Results: This comparative gonad transcriptomic analysis revealed that 7,120 and 4,628 DEGs are up-regulated in testes or ovaries, respectively. Further functional annotation analyses identified 76 important candidate genes involved in sex determination cascades or sex differentiation, including 42 known testis-biased DEGs (dmrt1, amh and sox9 etc.), and 34 known ovary-biased DEGs (foxl2, sox3 and cyp19a etc.). Moreover, eleven significantly enriched pathways functionally related to sex determination and sex differentiation were identified, including Wnt signaling pathway, oocyte meiosis, the TGF-beta signaling pathway and MAPK signaling pathway. Conclusion: This work represents the first comparative gonad transcriptome study in cobia. The putative sex-associated DEGs and pathways provide an important molecular basis for further investigation of cobia's sex determination, gonadal development as well as potential control breeding of monosex female populations for a possible aquaculture setting.

The evolution of the testis transcriptome in pregnant male pipefishes and seahorses

In many animals, sperm competition and sexual conflict are thought to drive the rapid evolution of male-specific genes, especially those expressed in the testes. A potential exception occurs in the male pregnant pipefishes, where females transfer eggs to the males, eliminating testes from participating in these processes. Here, we show that testis-related genes differ dramatically in their rates of molecular evolution and expression patterns in pipefishes and seahorses (Syngnathidae) compared to other fish. Genes involved in testis or sperm function within syngnathids experience weaker selection in comparison to their orthologs in spawning and livebearing fishes. An assessment of gene turnover and expression in the testis transcriptome suggests that syngnathids have lost (or significantly reduced expression of) important classes of genes from their testis transcriptomes compared to other fish. Our results indicate that more than 50 million years of male pregnancy have removed syngnathid testes from the molecular arms race that drives the rapid evolution of male reproductive genes in other taxa.

Gonadal Soma-Derived Factor Expression is a Potential Biomarker for Predicting the Effects of Endocrine-Disrupting Chemicals on Gonadal Differentiation in Japanese Medaka (Oryzias Latipes)

Chemicals with androgenic or estrogenic activity induce the sex reversal and/or intersex condition in various teleost fish species. Previously, we reported that exposure to 17α-methyltestosterone, bisphenol A, or 4-nonylphenol induces changes in expression of the gonadal soma-derived factor (gsdf) gene accompanied by disruption of gonadal differentiation in Japanese medaka (Oryzias latipes). These findings suggest that gsdf expression might be a useful biomarker for predicting the potential effect of chemicals on gonadal differentiation. We examined the gsdf expression in Japanese medaka exposed to chemicals with estrogenic or androgenic activity. Exposure to the androgenic steroid 17β-trenbolone at 0.5-22.1 μg/L induced the development of ovotestis (presence of ovarian tissue with testicular tissue) and female-to-male sex reversal in XX embryos, and exposure at 6.32 and 22.1 μg/L significantly increased gsdf expression in XX embryos compared with controls at developmental stage 38 (1 day before hatching). In the present study, no statistically significant difference in gsdf mRNA expression was observed after exposure to 17β-estradiol, 17α-ethinylestradiol, and 4-t-octylphenol, which have estrogenic activity. In addition, antiandrogenic chemicals or chemicals without endocrine-disrupting activity did not induce changes in gsdf expression in XX or XY embryos. Thus, an increase in gsdf expression after androgen exposure was observed in XX embryos. Together, these findings indicate that gsdf expression might be useful for predicting the adverse effect of chemicals on gonadal differentiation. Environ Toxicol Chem 2022;41:1875-1884. © 2022 SETAC.

Developmental changes in gene expression and gonad morphology during sex differentiation in Atlantic salmon (Salmo salar)

The Atlantic salmon (Salmo salar) is a globally important species for its value in fisheries and aquaculture, and as a research model. In order to characterise aspects of sex differentiation at the morphological and mRNA level in this species, the present study examined developmental changes in gonad morphology and gene expression in males and females between 0 and 79 days post hatch (dph). Morphological differentiation of the ovary (indicated by the formation of germ cell cysts) became apparent from 52 dph. By 79 dph, ovarian phenotype was evident in 100% of genotypic females. Testes remained in an undifferentiated-like state throughout the experiment, containing germ cells dispersed singularly within the gonadal region distal to the mesentery. There were no significant sex-related differences in gonad cross-section size, germ cell number or germ cell diameter during the experiment. The expression of genes involved in teleost sex differentiation (anti-müllerian hormone (amh), cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a), forkhead box L2a (foxl2a), gonadal soma-derived factor (gsdf), r-spondin 1 (rspo1), sexually dimorphic on the Y chromosome (sdY)), retinoic acid-signalling (aldehyde dehydrogenase 1a2 (aldh1a2), cytochrome P450 family 26 a1 (cyp26a1), cytochrome P450 family 26 b1 (cyp26b1), t-box transcription factor 1 (tbx1a)) and neuroestrogen production (cytochrome P450, family 19, subfamily A, polypeptide 1b (cyp19a1b)) was investigated. Significant sex-related differences were observed only for the expression of amh, cyp19a1a, gsdf and sdY. In males, amh, gsdf and sdY were upregulated from 34, 59 and 44 dph respectively. In females, cyp19a1a was upregulated from 66 dph. Independent of sex, foxl2a expression was highest at 0 dph and had reduced ∼ 47-fold by the time of morphological sex differentiation at 52 dph. This study provides new insights into the timing and sequence of some physiological changes associated with sex differentiation in Atlantic salmon. These findings also reveal that some aspects of the mRNA sex differentiation pathways in Atlantic salmon are unique compared to other teleost fishes, including other salmonids.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

Profile of gene expression changes during estrodiol-17β-induced feminization in the Takifugu rubripes brain

Background

As the critical tissue of the central nervous system, the brain has been found to be involved in gonad development. Previous studies have suggested that gonadal fate may be affected by the brain. Identifying brain-specific molecular changes that occur during estrodiol-17β (E₂) -induced feminization is crucial to our understanding of the molecular control of sex differentiation by the brains of fish.

Results

In this study, the differential transcriptomic responses of the Takifugu rubripes larvae brain were compared after E₂ treatment for 55 days. Our results showed that 514 genes were differentially expressed between E₂-treated-XX (E-XX) and Control-XX (C-XX) T. rubripes, while 362 genes were differentially expressed between E₂-treated-XY (E-XY) and Control-XY (C-XY). For example, the expression of cyp19a1b, gnrh1 and pgr was significantly up-regulated, while st, sl, tshβ, prl and pit-1, which belong to the growth hormone/prolactin family, were significantly down-regulated after E₂ treatment, in both sexes. The arntl1, bhlbe, nr1d2, per1b, per3, cry1, cipc and ciart genes, which are involved in the circadian rhythm, were also found to be altered. Differentially expressed genes (DEGs), which were identified between E-XX and C-XX, were significantly enriched in neuroactive ligand-receptor interaction, arachidonic acid metabolism, cytokine-cytokine receptor interaction and the calcium signaling pathway. The DEGs that were identified between E-XY and C-XY were significantly enriched in tyrosine metabolism, phenylalanine metabolism, arachidonic acid metabolism and linoleic acid metabolism.

Conclusion

A number of genes and pathways were identified in the brain of E₂-treated T. rubripes larvae by RNA-seq. It provided the opportunity for further study on the possible involvement of networks in the brain-pituitary-gonadal axis in sex differentiation in T. rubripes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08158-0.

Comparison of differential expression genes in ovaries and testes of Pearlscale angelfish Centropyge vrolikii based on RNA-Seq analysis

Pearlscale angelfish Centropyge vrolikii is a kind of protogynous hermaphrodite fish with a natural sexual reversion. Under appropriate social conditions, a female fish can transform into a male fish spontaneously. It is an important prerequisite for artificial breeding to understand the process of its gonadal development and sexual reversion. Gonadal development is regulated by many sex-related genes. In this study, we used unreferenced RNA-Seq technology to sequence the ovary at the perinucleolus stage (O_II), ovary at the yolk vesicle stage (O_IV),IV and testis (T), respectively; screened the gonadal differential expression genes (DEGs); and analyzed the expression of these genes in different developmental stages of ovary and different sex gonads. The results showed that a total of 142,589 all-unigene samples were assembled, and gene annotation was performed by COG, GO, KEGG, KOG, Pfam, Swissprot, eggNOG, and NR functional database. Comparative analysis revealed that there were 1919 genes that were up-regulated and 1289 genes were down-regulated in comparison to O_IV vs O_II, while there were 3653 genes that were up-regulated and 2874 genes were down-regulated in comparison of O_IV vs T, there were 3345 genes that were up-regulated and 2995 genes were down-regulated in comparison of the O_II vs the T. At the same time, the results verified by RT-qPCR were consistent with the variation trend of transcriptome data. Among the results, amh, sox9b, dmrt1, dmrt2, cyp11a, cyp17a, and cyp19a were significantly expressed in the testes, while sox3, sox4, sox11, sox17, and hsd3b7 were significantly expressed in the ovaries. And, the expression of the amh, sox9b, dmrt2, and dmrt1 were low in the O_II and O_IV, while significantly increased during the ovotestis in the hermaphroditic period (O_T), and finally reached the highest level in pure testis after sex reversal. The expression of sox3, sox4, hsd3b7, sox11, and sox17 was significantly reduced during the hermaphroditic period (O_T). These results suggested that these genes may play an important role in the process of sex reversal. This study is helpful to further understand the molecular regulation mechanism of gonadal development and sexual reversion in Pearlscale angelfish and also provide important clues for future studies.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

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

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

Formation of autotriploid Carassius auratus and its fertility-related genes analysis

Background

Formation of triploid organism is useful in genetics and breeding. In this study, autotriploid Carassius auratus (3nRR, 3n = 150) was generated from Carassius auratus red var. (RCC, 2n = 100) (♀) and autotetraploid Carassius auratus (4nRR, 4n = 200) (♂). The female 3nRR produced haploid, diploid and triploid eggs, whereas the male 3nRR was infertile. The aim of the present study was to explore fertility of potential candidate genes of 3nRR.

Results

Gonadal transcriptome profiling of four groups (3 females RCC (FRCC), 3 males 4nRR (M4nRR), 3 males 3nRR (M3nRR) and 3 females 3nRR (F3nRR)) was performed using RNA-SEq. A total of 78.90 Gb of clean short reads and 24,262 differentially expressed transcripts (DETs), including 20,155 in F3nRR vs. FRCC and 4,107 in M3nRR vs. M4nRR were identified. A total of 106 enriched pathways were identified through KEGG enrichment analysis. Out of the enriched pathways, 44 and 62 signalling pathways were identified in F3nRR vs. FRCC and M3nRR vs. M4nRR, respectively. A total of 80 and 25 potential candidate genes for fertility-related in F3nRR and M3nRR were identified, respectively, through GO, KEGG analyses and the published literature. Moreover, protein-protein interaction (PPI) network construction of these fertility-associated genes were performed. Analysis of the PPI networks showed that 6 hub genes (MYC, SOX2, BMP4, GATA4, PTEN and BMP2) were involved in female fertility of F3nRR, and 2 hub genes (TP53 and FGF2) were involved in male sterility of M3nRR.

Conclusions

Establishment of autotriploid fish offers an ideal model to study reproductive traits of triploid fish. RNA-Seq data revealed 6 genes, namely, MYC, SOX2, BMP4, GATA4, PTEN and BMP2, involved in the female fertility of the F3nRR. Moreover, 2 genes, namely, TP53 and FGF2, were related to the male sterility of the M3nRR. These findings provide information on reproduction and breeding in triploid fish.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07753-5.

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.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Effects of the Thyroid Endocrine System on Gonadal Sex Ratios and Sex-Related Gene Expression in the Pufferfish Takifugu rubripes

To examine the effect and mechanism of thyroid hormone on gonadal sex differentiation, Takifugu rubripes larvae were treated with goitrogen (methimazole, MET, 1000 g/g), and thyroxine (T4, 2nM) from 25 to 80 days after hatching (dah). Gonadal histology and sex ratios of fish were then determined at 80 dah. MET treatment induced masculinization, but T4 treatment did not induce feminization in T. rubripes larvae. Transcriptomic analysis of gonads at 80 dah was then conducted. Among the large number of differentially expressed genes between the groups, the expression of foxl2, cyp19a1a, and dmrt1 was altered. The expression of foxl2, cyp19a1a, dmrt1 and gsdf at 25, 40, 55 days after treatment (dat) was further analyzed by qPCR. MET treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 in genetic females (p < 0.05). Additionally, T4 treatment induced an increase in the expression of cyp19a1a in genetic XY gonads only at 25 dat. However, the increase in cyp19a1a expression did not continue to 40 and 55 dat. This may explain why feminization of larvae was not found in the T4-treated group. Thus, the present study provides the first evidence that MET treatment causes masculinization in teleost fish. The effects of MET-induced masculinization in T. rubripes may act primarily via suppression of the expression of foxl2 and cyp19a1a, and stimulation of the expression of dmrt1. Moreover, the effects of higher concentrations of T4 or different concentrations of T3, on sex differentiation require further testing.

Characterization of two myostatin genes in pufferfish Takifugu bimaculatus: sequence, genomic structure, and expression

Myostatin (MSTN) is a negative regulator of muscle growth, which restrains the proliferation and differentiation of myoblasts. To understand the role of two mstn genes of Takifugu bimaculatus, the full-length cDNAs of 1131 bp Tbmstn1 and 1,080 bp Tbmstn2 were obtained from the T. bimaculatus' genomic database, which encodes 376 and 359 amino acids, respectively. The results of qRT-PCR showed that Tbmstn1 was expressed in the eye, kidney, spleen, skeletal muscle, gill, and brain, and the expression level in the skeletal muscle was extremely significantly higher than in other examined tissues. Tbmstn2 was expressed in the skin, skeletal muscle, gill, and brain, and had the highest expression in the skeletal muscle, followed by expression in the brain. Meanwhile, in different stages of embryonic development, the expression of Tbmstn1 started from the gastrula stage. Its expression in the eye-pigment formation stage and hatching stage was significantly higher than that in other stages. The Tbmstn2 was expressed in all examined embryonic stages with different levels, and the highest expression was detected in the eye-pigment formation stage. These results suggested that Tbmstn1 and Tbmstn2 may involve in the development of skeletal muscle, and Tbmstn2 may be related to the formation of nervous system.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Transcriptome analysis of genes related to gonad differentiation and development in Muscovy ducks

BACKGROUND

Sex-related genes play a crucial role in gonadal differentiation into testes or ovaries. However, the genetic control of gonadal differentiation in Muscovy ducks remains unknown. Therefore, the objective of our study was to screen new candidate genes associated with ovarian and testicular development.

RESULTS

In this study, 24 males before gonadal differentiation (MB), 24 females before gonadal differentiation (FB), 24 males after gonadal differentiation (MA) and 24 females after gonadal differentiation (FA) were selected from Putian Muscovy ducks, forming 4 groups. RNA-Seq revealed 101.76 Gb of clean reads and 2800 differentially expressed genes (DEGs), including 46 in MB vs FB, 609 in MA vs FA, 1027 in FA vs FB, and 1118 in MA vs MB. A total of 146 signalling pathways were enriched by KEGG analysis, among which 20, 108, 108 and 116 signalling pathways were obtained in MB vs FB, MA vs MB, MA vs FA and FA vs FB, respectively. In further GO and KEGG analyses, a total of 21 candidate genes related to gonad differentiation and development in Muscovy ducks were screened. Among these, 9 genes were involved in the differentiation and development of the testes, and 12 genes were involved in the differentiation and development of the ovaries. In addition, RNA-Seq data revealed 2744 novel genes.

CONCLUSIONS

RNA-Seq data revealed 21 genes related to gonadal differentiation and development in Muscovy ducks. We further identified 12 genes, namely, WNT5B, HTRA3, RSPO3, BMP3, HNRNPK, NIPBL, CREB3L4, DKK3, UBE2R2, UBPL3KCMF1, ANXA2, and OSR1, involved in the differentiation and development of ovaries. Moreover, 9 genes, namely, TTN, ATP5A1, DMRT1, DMRT3, AMH, MAP3K1, PIK3R1, AGT and ADAMTSL1, were related to the differentiation and development of testes. Moreover, after gonadal differentiation, DMRT3, AMH, PIK3R1, ADAMTSL1, AGT and TTN were specifically highly expressed in males. WNT5B, ANXA2 and OSR1 were specifically highly expressed in females. These results provide valuable information for studies on the sex control of Muscovy ducks and reveal novel candidate genes for the differentiation and development of testes and ovaries.

Comparative transcriptome analysis of male and female gonads reveals sex-biased genes in spotted scat (Scatophagus argus)

Scatophagus argus is a new emerging aquaculture fish in East and Southeast Asia. To date, research on reproductive development and regulation in S. argus is lacking. Additionally, genetic and genomic information about reproduction, such as gonadal transcriptome data, is also lacking. Herein, we report the first gonadal transcriptomes of S. argus and identify genes potentially involved in reproduction and gonadal development. A total of 136,561 unigenes were obtained by sequencing of testes (n = 3) and ovaries (n = 3) at stage III. Genes upregulated in males and females known to be involved in gonadal development and gametogenesis were identified, including male-biased dmrt1, amh, gsdf, wt1a, sox9b, and nanos2, and female-biased foxl2, gdf9, bmp15, sox3, zar1, and figla. Serum estradiol-17β and 11-ketotestosterone levels were biased in female and male fish, respectively. Sexual dimorphism of serum steroid hormone levels were interpreted after expression analysis of 20 steroidogenesis-related genes, including cyp19a1a and cyp11b2. This gonadal transcript dataset will help investigate functional genes related to reproduction in S. argus.

A T > G Mutation in the NR5A2 Gene Is Associated With Litter Size in Hu Sheep Through Upregulation of Promoter Activity by Transcription Factor MTF-1

Nuclear receptor subfamily 5 group A member 2 (NR5A2), also referred to as LRH-1 or FTF, is an orphan nuclear hormone receptor that is involved in regulating embryonic development, ovarian granulosa cell differentiation, gonadal sex differentiation, and steroidogenesis in mammals. However, little is known about how NR5A2 regulates reproduction in sheep. In this study, we amplified the promoter sequence of NR5A2 and determined that its core promoter region ranged from -721 nt to -281 nt. A T > G polymorphism at -700 nt was detected in the core promoter region. Association analysis found that the litter sizes of Hu ewes at their second and average parities with genotype GG (2.20 ± 0.20 and 1.97 ± 0.06, respectively) were significantly higher than those of ewes with genotype TG (1.68 ± 0.10 and 1.74 ± 0.05, respectively) (p < 0.05) and TT (1.67 ± 0.10 and 1.62 ± 0.06, respectively) (p < 0.05). The litter size of Hu ewes at their third parity with genotype GG (2.10 ± 0.10) was significantly higher than that of ewes with genotype TT (1.56 ± 0.12) (p < 0.05). A luciferase assay showed that the -700G allele increased the luciferase activity relative to the -700T allele. Furthermore, the -700T > G polymorphism created a novel binding site for metal-regulatory transcription factor 1 (MTF-1). A competitive electrophoretic mobility shift assay confirmed that MTF-1 specifically bound with the G-type promoter of NR5A2. An overexpression experiment demonstrated that MTF-1 was involved in the alteration of NR5A2 transcription activity and further increased NR5A2 gene mRNA expression. Our findings revealed that the -700T > G polymorphism promoted NR5A2 expression due to the positive effects on NR5A2 gene transcription activity by MTF-1 and thereby increased fecundity in Hu sheep.

Phylogenetic analysis and embryonic expression of panarthropod Dmrt genes

Background

One set of the developmentally important Doublesex and Male-abnormal-3 Related Transcription factors (Dmrt) is subject of intense research, because of their role in sex-determination and sexual differentiation. This likely non-monophyletic group of Dmrt genes is represented by the Drosophila melanogaster gene Doublesex (Dsx), the Caenorhabditis elegans Male-abnormal-3 (Mab-3) gene, and vertebrate Dmrt1 genes. However, other members of the Dmrt family are much less well studied, and in arthropods, including the model organism Drosophila melanogaster, data on these genes are virtually absent with respect to their embryonic expression and function.

Results

Here we investigate the complete set of Dmrt genes in members of all main groups of Arthropoda and a member of Onychophora, extending our data to Panarthropoda as a whole. We confirm the presence of at least four families of Dmrt genes (including Dsx-like genes) in Panarthropoda and study their expression profiles during embryogenesis. Our work shows that the expression patterns of Dmrt11E, Dmrt93B, and Dmrt99B orthologs are highly conserved among panarthropods. Embryonic expression of Dsx-like genes, however, is more derived, likely as a result of neo-functionalization after duplication.

Conclusions

Our data suggest deep homology of most of the panarthropod Dmrt genes with respect to their function that likely dates back to their last common ancestor. The function of Dsx and Dsx-like genes which are critical for sexual differentiation in animals, however, appears to be much less conserved.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0322-0) contains supplementary material, which is available to authorized users.