Characterization of gonadal transcriptomes from Nile tilapia (Oreochromis niloticus) reveals differentially expressed genes

Gonadal Development and Differentiation of Hybrid F1 Line of Ctenopharyngodon idella (♀) × Squaliobarbus curriculus (♂)

The hybrid F1 offspring of Ctenopharyngodon idella (♂) and Squaliobarbus curriculus (♀) exhibit heterosis in disease resistance and also show abnormal sex differentiation. To understand the mechanism behind gonadal differentiation in the hybrid F1, we analyzed the transcriptomes of C. idella, S. curriculus, and the hybrid F1; screened for genes related to gonad development in these samples; and measured their expression levels. Our results revealed that compared to either C. idella or S. curriculus, the gene expressions in most sub-pathways of the SNARE interactions in the vesicular transport pathway in the hypothalamus, pituitary, and gonadal tissues of their hybrid F1 offspring were significantly up-regulated. Furthermore, insufficient transcription of genes involved in oocyte meiosis may be the main reason for the insufficient reproductive ability of the hybrid F1 offspring. Through transcriptome screening, we identified key molecules involved in gonad development, including HSD3B7, HSD17B1, HSD17B3, HSD20B2, CYP17A2, CYP1B1, CYP2AA12, UGT2A1, UGT1A1, and FSHR, which showed significant differences in expression levels in the hypothalamus, pituitary, and gonads of these fish. Notably, the expression levels of UGT1A1 in the gonads of the hybrid F1 were significantly higher than those in C. idella and S. curriculus. These results provide a scientific basis for further research on the gonadal differentiation mechanism of hybrid F1 offspring.

dmrt1 Is Responsible for Androgen-Induced Masculinization in Nile Tilapia

17α-Methyltestosterone (MT) is a widely used androgen for all-male fish production in aquaculture. However, the molecular mechanism underlying MT-induced masculinization remains unclear. In this study, we aim to identify the key gene responsible for MT-induced masculinization using the Nile tilapia (Oreochromis niloticus) amhy, dmrt1, and gsdf mutants, which exhibit male-to-female sex reversal. Nile tilapia fry from these three mutant lines were treated with 50 μg/g MT from 5 to 30 days after hatching (dah). The results showed that amhy and gsdf mutants, but not dmrt1 mutants, were masculinized by the MT treatment. Gonadal transcriptome analysis revealed that genes involved in steroidogenesis and germ cell development in MT-treated dmrt1 mutants exhibited a similar expression pattern to that of the wild type (WT) XX. In addition, the dmrt1 mutants cannot be masculinized by co-treatment with MT and the aromatase inhibitor fadrozole. The MT treatment completely blocked early steroidogenic enzyme (Star2, Cyp17a2, and Cyp19a1a) expression independent of amhy, gsdf, and dmrt1. A luciferase analysis showed that MT directly suppressed basal and Sf-1-activated cyp19a1a promoter activity through ara and arb in cultured HEK293 cells. Furthermore, MT treatment inhibited germ cell proliferation in amhy and gsdf mutants but not in dmrt1 mutants. Consistently, dmrt1 expression was induced in MT-treated WT XX, -amhy, and -gsdf mutants. Taken together, these results suggest that dmrt1 is indispensable for MT-induced masculinization in Nile tilapia and that MT functions by inhibiting early steroid synthesis and activating dmrt1 to promote testis development.

Steroid hormone-deprived sex reversal in cyp11a1 mutant XX tilapia experiences an ovary-like stage at molecular level

Fish sex is largely influenced by steroid hormones, especially sex hormones. Here, we established a steroid hormone-free genetic model by mutation of cyp11a1 in Nile tilapia, which was confirmed by EIA assay. Gonadal phenotype and transcriptome analyses showed that the XX mutants displayed sex reversal from female to male but with defective spermatogenesis. Despite the sex reversal, the aromatase encoding gene cyp19a1a was continuously expressed in the gonads of the XX mutants, which might be caused by androgen deficiency. Whole-mount fluorescence in situ hybridization and transcriptome analysis showed that the gonads of the XX mutants firstly developed towards ovary but shifted to testis between 10 to 15 days after hatching. Detailed expression analysis of key sex differentiation pathway genes foxl3 and dmrt1 combined with apoptosis analysis revealed transdifferentiation of germ cells from female to male during sex reversal. Rescue experiments showed that both P5 and E2 treatment rescued the sex reversal of cyp11a1 mutant XX fish. Overall, our results revealed a transient ovary-like stage and transdifferentiation of germ cells from female to male in the early gonads of the steroid hormone-deprived cyp11a1 mutant XX fish.

CRISPR/Cas9 Technology for Enhancing Desirable Traits of Fish Species in Aquaculture

Aquaculture, the world's fastest-growing food production sector, is critical for addressing food security concerns because of its potential to deliver high-quality, nutrient-rich supplies by 2050. This review assesses the effectiveness of CRISPR/Cas9 genome editing technology in enhancing desirable traits in fish species, including growth rates, muscle quality, disease resistance, pigmentation, and more. It also focuses on the potential effectiveness of the technology in allowing precise and targeted modifications of fish DNA to improve desirable characteristics. Many studies have reported successful applications of CRISPR/Cas9, such as knocking out reproductive genes to control reproduction and sex determination, enhancing feed conversion efficiency, and reducing off-target effects. Additionally, this technology has contributed to environmental sustainability by reducing nitrogen-rich waste and improving the nutritional composition of fish. However, the acceptance of CRISPR/Cas9 modified fish by the public and consumers is hindered by concerns regarding public perception, potential ecological impacts, and regulatory frameworks. To gain public approval and consumer confidence, clear communication about the editing process, as well as data on the safety and environmental considerations of genetically modified fish, are essential. This review paper discusses these challenges, provides possible solutions, and recommends future research on the integration of CRISPR/Cas9 into sustainable aquaculture practices, focusing on the responsible management of genetically modified fish to enable the creation of growth and disease-resistant strains. In conclusion, this review highlights the transformative potential of CRISPR/Cas9 technology in improving fish traits, while also considering the challenges and ethical considerations associated with sustainable and responsible practices in aquaculture.

The cellular expression patterns of gdnfa and gdnfb in the gonads of Nile tilapia and their differential response to retinoic acid

In mammals, glial cell derived neurotrophic factor (GDNF) plays a critical role in the self-renewal and maintenance of spermatogonial stem cells (SSCs) in testis and oogenesis in ovary, whilst retinoic acid (RA), the key factor of meiosis initiation, can downregulate its expression. Unlike mammals, two Gdnf replication genes are widely present in teleost fishes, however, our understanding of them is still poor. In the present study, two paralogous gdnf from Nile tilapia (Oreochromis niloticus), namely as Ongdnfa and Ongdnfb, were characterized, and then their cellular expression profiles in testis and ovary and responsiveness to RA treatment at the tissue and cellular levels were investigated. In phylogenetic tree, the Gdnfa and Gdnfb from teleost fishes were clustered into two different subclasses, respectively, and then clustered with the homologs from cartilaginous fish and tetrapods, suggesting that OnGdnfa and OnGdnfb are orthologous to GDNF and paralogous to each other. Ongdnfa is expressed in Sertoli cells and Leydig cells in testis and oocytes in ovary. The expression pattern of Ongdnfb is similar to Ongdnfa. In the ex vivo testicular organ culture, RA down-regulated the expression of Ongdnfa, whereas up-regulated the expression of Ongdnfb (P < 0.05), suggesting that they have differential responsiveness to RA signaling. RA treatment of the cultured cells derived from adult Nile tilapia testis which have the expression of RA receptors (RAR), Ongdnfa and Ongdnfb further confirmed the above result. Collectively, our study suggests that Ongdnfa and Ongdnfb have non-germline expression patterns in testis and germline expression patterns in ovary; furthermore, they have differential responsiveness to RA signaling, implying that they might have differential biological functions. This study broadens and enriches our understanding of fish GDNF homologs and lays foundation for the study of their biological functions in the future.

Synchronously sexual maturity in hermaphrodite fish as revealed by transcriptome analysis in Plectropomus leopardus

Leopard coral grouper (Plectropomus leopardus) is a type of hermaphrodite fish, but the mechanisms of gonadal development and gametogenesis remain unclear. In the present study, we performed histological observation and transcriptomic analysis during the process of sexual differentiation in P. leopardus. According to the histological results, sexual differentiation was completed at 15 months old, developed synchronously in male and female individuals at 2 years old, and matured synchronously at 3 years old. Comparative transcriptomic analyses showed that the gonadal had differentiated by 15 months old, with enrichment of pathways associated with cell proliferation, transcriptional metabolism, and germline stem cell differentiation. Furthermore, cilium movement and fatty acid anabolism, which are associated with spermatogenesis and oocyte growth, were significantly enriched at 3 years old. In addition, key genes associated with male and female sex differentiation, such as amh, dmrt1, dmrt2a, zp4, sox3, gdf9, and gsdf, were identified by weighted gene co-expression network analysis (WGCNA). Finally, the localization and expression of the key genes amh and sox3 were observed in different cell types within the testes and ovaries, reflecting the development of the testes and ovaries, respectively. All the evidence indicates that P. leopardus is a hermaphrodite and synchronously sexually mature fish. Our study complements the gonadal development patterns of hermaphroditic fish by providing new insights into the molecular mechanisms underlying sexual differentiation and sex change in hermaphroditic groupers.

Expression characteristics of the cyp19a1b aromatase gene and its response to 17β-estradiol treatment in largemouth bass (Micropterus salmoides)

To investigate the regulatory role of the cyp19a1b aromatase gene in the sexual differentiation of largemouth bass (Micropterus salmoides, LMB), we obtained the full-length cDNA sequence of cyp19a1b using rapid amplification of cDNA ends technique. Tissue expression characteristics and feedback with 17-β-estradiol (E2) were determined using quantitative real-time PCR (qRT-PCR), while gonad development was assessed through histological section observations. The cDNA sequence of LMB cyp19a1b was found to be1950 base pairs (bp) in length, including a 5' untranslated region of 145 bp, a 3' untranslated region of 278 bp, and an open reading frame encoding a protein consisting of 1527 bp that encoded 508 amino acids. The qRT-PCR results indicated that cyp19a1b abundantly expressed in the brain, followed by the gonads, and its expression in the ovaries was significantly higher than that observed in the testes (P < 0.05). After feeding fish with E2 for 30 days, the expression of cyp19a1b in the pseudo-female gonads (XY-F) was significantly higher than that in males (XY-M) (P < 0.05), whereas expression did not differ significantly between XX-F and XY-F fish (P > 0.05). Although the expression of cyp19a1b in XY-F and XX-F fish was not significantly different after 60 days (P>0.05), both exhibited significantly higher levels than that of XY-M fish (P<0.05). Histological sections analysis showed the presence of oogonia in both XY-F and XX-F fish at 30 days, while spermatogonia were observed in XY-M fish. At 60 days, primary oocytes were abundantly observed in both XY-F and XX-F fish, while a few spermatogonia were visible in XY-M fish. At 90 days, the histological sections' results showed that a large number of oocytes were visible in XY-F and XX-F fish. Additionally, the gonads of XY-M fish contained numerous spermatocytes. These results suggest that cyp19a1b plays a pivotal role in the development of ovaries and nervous system development in LMB.

Transcriptome analysis of Schizothorax oconnori (Cypriniformes: Cyprinidae) oocytes: The role of K+ in promoting yolk globule fusion and regulating oocyte maturation

Schizothorax oconnori (S. oconnori) is an economically important fish in Tibet. Oocyte maturation is a physiological process that is of great significance to reproduction and seed production in S. oconnori, yet little is currently known regarding the molecular mechanisms of oocyte development in this species. To identify candidate genes involved in reproduction of female fish, a combination of PacBio and Illumina HiSeq technologies was employed to provide deep coverage of the oocyte transcriptome. Transcriptome analysis revealed several candidate genes that are potentially involved in the regulation of oocyte maturation in S. oconnori, including GIRK1, CHRM3, NPY2R, GABRA3, GnRH3, mGluR1α, GPER1, GDF9, HSP90, and ESR2. Genes that are significantly expressed during oocyte maturation mainly contribute to the GPCR signaling pathway and the estrogen signaling pathway. Neurotransmitter (Ach, NPY, and GABA) and peptide hormone (GnRH3) binding to G protein-coupled receptors (GPCRs) frees G-protein βγ subunits to interact with the G protein-gated inward rectifier K+ channel 1 (GIRK1). This process helps release K+ from granulosa cells to maturing oocytes, allowing yolk globule fusion. This mechanism may play an important role in oocyte maturation in S. oconnori. In conclusion, this study provides a valuable basis for deciphering the reproductive system in S. oconnori during the oocyte maturation process.

Dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia

Sex is determined by multiple factors derived from somatic and germ cells in vertebrates. We have identified amhy, dmrt1, gsdf as male and foxl2, foxl3, cyp19a1a as female sex determination pathway genes in Nile tilapia. However, the relationship among these genes is largely unclear. Here, we found that the gonads of dmrt1;cyp19a1a double mutants developed as ovaries or underdeveloped testes with no germ cells irrespective of their genetic sex. In addition, the gonads of dmrt1;cyp19a1a;cyp19a1b triple mutants still developed as ovaries. The gonads of foxl3;cyp19a1a double mutants developed as testes, while the gonads of dmrt1;cyp19a1a;foxl3 triple mutants eventually developed as ovaries. In contrast, the gonads of amhy;cyp19a1a, gsdf;cyp19a1a, amhy;foxl2, gsdf;foxl2 double and amhy;cyp19a1a;cyp19a1b, gsdf;cyp19a1a;cyp19a1b triple mutants developed as testes with spermatogenesis via up-regulation of dmrt1 in both somatic and germ cells. The gonads of amhy;foxl3 and gsdf;foxl3 double mutants developed as ovaries but with germ cells in spermatogenesis due to up-regulation of dmrt1. Taking the respective ovary and underdeveloped testis of dmrt1;foxl3 and dmrt1;foxl2 double mutants reported previously into consideration, we demonstrated that once dmrt1 mutated, the gonad could not be rescued to functional testis by mutating any female pathway gene. The sex reversal caused by mutation of male pathway genes other than dmrt1, including its upstream amhy and downstream gsdf, could be rescued by mutating female pathway gene. Overall, our data suggested that dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

The fusion gene hsf5-rnf43 in Nile tilapia: A potential regulator in the maintenance of testis function and sexual differentiation

We identified that the genes heat shock transcription factor 5 (hsf5) and ring finger protein 43 (rnf43) happened fusion in Nile tilapia (Oreochromis niloticus), called hsf5-rnf43, and provided the characteristic and functional analysis of hsf5-rnf43 gene in fish for the first time. Analysis of spatiotemporal expression showed that hsf5-rnf43 was specifically expressed in the testis and located in primary spermatocytes of adult Nile tilapia and gradually increased during testis development from 5 to 180 days after hatching. We also found DNA methylation regulated sex-biased expression of hsf5-rnf43 in the early development of Nile tilapia, and was affected by high temperature during the thermosensitive period of Nile tilapia sex differentiation. Therefore, we first reported that the fusion gene hsf5-rnf43 was sex-biased expressed in the testis regulated by DNA methylation and affected by high temperature, which may be involved in the maintenance of testis function and sex differentiation of Nile tilapia.

Molecular biological, physiological, cytological, and epigenetic mechanisms of environmental sex differentiation in teleosts: A systematic review

Human activities have been exerting widespread stress and environmental risks in aquatic ecosystems. Environmental stress, including temperature rise, acidification, hypoxia, light pollution, and crowding, had a considerable negative impact on the life histology of aquatic animals, especially on sex differentiation (SDi) and the resulting sex ratios. Understanding how the sex of fish responds to stressful environments is of great importance for understanding the origin and maintenance of sex, the dynamics of the natural population in the changing world, and the precise application of sex control in aquaculture. This review conducted an exhaustive search of the available literature on the influence of environmental stress (ES) on SDi. Evidence has shown that all types of ES can affect SDi and universally result in an increase in males or masculinization, which has been reported in 100 fish species and 121 cases. Then, this comprehensive review aimed to summarize the molecular biology, physiology, cytology, and epigenetic mechanisms through which ES contributes to male development or masculinization. The relationship between ES and fish SDi from multiple aspects was analyzed, and it was found that environmental sex differentiation (ESDi) is the result of the combined effects of genetic and epigenetic factors, self-physiological regulation, and response to environmental signals, which involves a sophisticated network of various hormones and numerous genes at multiple levels and multiple gradations in bipotential gonads. In both normal male differentiation and ES-induced masculinization, the stress pathway and epigenetic regulation play important roles; however, how they co-regulate SDi is unclear. Evidence suggests that the universal emergence or increase in males in aquatic animals is an adaptation to moderate ES. ES-induced sex reversal should be fully investigated in more fish species and extensively in the wild. The potential aquaculture applications and difficulties associated with ESDi have also been addressed. Finally, the knowledge gaps in the ESDi are presented, which will guide the priorities of future research.

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.

Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae)

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

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.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

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.

Gonadal Transcriptome Analysis and Sequence Characterization of Sex-Related Genes in Cranoglanis bouderius

In China, the Cranoglanis bouderius is classified as a national class II-protected animal. The development of C. bouderius populations has been affected by a variety of factors over the past few decades, with severe declines occurring. Considering the likelihood of continued population declines of the C. bouderius in the future, it is critical to investigate the currently unknown characteristics of gonadal differentiation and sex-related genes for C. bouderius conservation. In this study, the Illumina sequencing platform was used to sequence the gonadal transcriptome of the C. bouderius to identify the pathways and genes related to gonadal development and analyze the expression differences in the gonads. A total of 12,002 DEGs were identified, with 7220 being significantly expressed in the ovary and 4782 being significantly expressed in the testis. According to the functional enrichment results, the cell cycle, RNA transport, apoptosis, Wnt signaling pathway, p53 signaling pathway, and prolactin signaling pathway play important roles in sex development in the C. bouderius. Furthermore, the sequence characterization and evolutionary analysis revealed that AMH, DAX1, NANOS1, and AR of the C. bouderius are highly conserved. Specifically, the qRT-PCR results from various tissues showed significant differences in AMH, DAX1, NANOS1, and AR expression levels in the gonads of both sexes of C. bouderius. These analyses indicated that AMH, DAX1, NANOS1, and AR may play important roles in the differentiation and development of C. bouderius gonads. To our best knowledge, this study is the first to analyze the C. bouderius gonadal transcriptome and identify the structures of sex-related genes, laying the foundation for future research.

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

Expression and regulation of 42Sp50 in spotted scat (Scatophagus argus)

42Sp50 is an isoform of the eukaryotic translation elongation factor 1 A (eEF1A) and is vital for fish ovarian development. Spotted scat (Scatophagus argus) is a popular marine cultured fish species in Southern Asia and China, and its artificial reproduction is complicated, with a relatively low success ratio in practice. In this study, the 42Sp50 gene was cloned from spotted scat. Tissue distribution analysis showed that 42Sp50 was mainly expressed in the ovary. qRT-PCR showed that 42Sp50 expression levels gradually decreased insignificantly in the ovaries from phase II to IV. Western blot analysis showed that 42Sp50 was highly expressed in the ovary, while it was almost undetectable in the testis. Immunohistochemistry analysis stained 42Sp50 mainly in the cytoplasm of the previtellogenic oocytes in ovaries of normal XX-female and sex-reversed XY-female. Aside from fish and amphibians, 42Sp50 was also identified in some reptile species using genomic database searching. Analyses of the transcriptome data from four different fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), and Hong Kong catfish (Clarias fuscus)) revealed ovaries biased expression of 42Sp50 in all, similar to spotted scat. While the neighbor genes of 42Sp50 did not show ovary biased expression in the fish species analyzed. Bisulfite Sequencing PCR (BSP) results showed that the DNA methylation level of 42Sp50 promoter was low in ovaries, testes, and muscles. The luciferase reporter assay demonstrated that Dmrt4 activated 42Sp50 expression in the presence of Sf1 or Foxh1. These results suggest that 42Sp50 may be involved in regulating the early phase oocytes development of spotted scat.

Function of Foxl2 and Dmrt1 proteins during gonadal differentiation in the olive flounder Paralichthys olivaceus

Study on fish sex differentiation is important both from academic and practical aspects. Foxl2 and Dmrt1 are important transcription factors that should be involved in fish gonadal differentiation, but there is still no direct evidence to clarify their protein functions. Olive flounder Paralichthys olivaceus, an important mariculture fish in China, Japan, and Korea, shows sex-dimorphic growth. In this study, the Foxl2 and Dmrt1 proteins were detected in granulosa cells of the ovary and Sertoli cells of the testis, respectively, showing significant sex-dimorphic expression patterns. Then, bioactive high-purity Foxl2 and Dmrt1 recombinant proteins were obtained in vitro. Furthermore, effects of the recombinant Foxl2 and Dmrt1 during gonadal differentiation period were evaluated by intraperitoneal injection in juvenile fish. Compared with the control group, the male rate in the Dmrt1 group increased from 0 % to 82 %, showing for the first time in fish that the recombinant Dmrt1 could alter the sex phenotype. In addition, transcription levels of cyp19a and its transcription factors also changed after the recombinant Foxl2 and Dmrt1 injection. These findings reveal that Foxl2 and Dmrt1 are vital regulators for fish gonadal differentiation by regulating cyp19a expression, and also provide a new approach for sex control in fish aquaculture.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

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.

Pervasive male-biased expression throughout the germline-specific regions of the sea lamprey genome supports key roles in sex differentiation and spermatogenesis

Sea lamprey undergo programmed genome rearrangement (PGR) in which ∼20% of the genome is jettisoned from somatic cells during embryogenesis. Although the role of PGR in embryonic development has been studied, the role of the germline-specific region (GSR) in gonad development is unknown. We analysed RNA-sequence data from 28 sea lamprey gonads sampled across life-history stages, generated a genome-guided de novo superTranscriptome with annotations, and identified germline-specific genes (GSGs). Overall, we identified 638 GSGs that are enriched for reproductive processes and exhibit 36x greater odds of being expressed in testes than ovaries. Next, while 55% of the GSGs have putative somatic paralogs, the somatic paralogs are not differentially expressed between sexes. Further, putative orthologs of some the male-biased GSGs have known functions in sex determination or differentiation in other vertebrates. We conclude that the GSR of sea lamprey plays an important role in testicular differentiation and potentially sex determination.

RNA-sequencing of sea lamprey gonads at different life-history stage identifies germline-specific genes which are highly expressed in males during spermatogenesis. This suggests a link between male-biased germline expression and sex differentiation in the sea lamprey.

Integrated analysis of the miRNA and mRNA expression profiles in Leiocassis longirostris at gonadal maturation

Leiocassis longirostris is a commercially important fish species that shows a sexually dimorphic growth pattern. A lack of molecular data from the gonads of this species has hindered research and selective breeding efforts. In this study, we conducted a comprehensive analysis of the expression profile of miRNA and mRNA to explore their regulatory roles in the gonadal maturation stage of L. longirostris. We identified 60 differentially expressed miRNAs and 20,752 differentially expressed genes by sequencing. A total of 90 miRNAs and 21 target genes involved in gonad development and sex determination were identified. Overall, the results of this study enhance our understanding of the molecular mechanisms underlying sex determination and differentiation and provide valuable genomic information for the selective breeding of L. longirostris.

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

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

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Identification of SF-1 and FOXL2 and Their Effect on Activating P450 Aromatase Transcription via Specific Binding to the Promoter Motifs in Sex Reversing Cheilinus undulatus

The giant wrasse Cheilinus undulatus is a protogynous socially hermaphroditic fish. However, the physiological basis of its sex reversal remains largely unknown. cyp19 is a key gender-related gene encoding P450 aromatase, which converts androgens to estrogens. cyp19 transcription regulation is currently unknown in socially sexually reversible fish. We identified NR5A1 by encoding SF-1, and FOXL2 from giant wrasse cDNA and cyp19a1a and cyp19a1b promoter regions were cloned from genomic DNA to determine the function of both genes in cyp19a1 regulation. Structural analysis showed that SF-1 contained a conserved DNA-binding domain (DBD) and a C-terminal ligand-binding domain (LBD). FOXL2 was comprised of an evolutionarily conserved Forkhead domain. In vitro transfection assays showed that SF-1 could upregulate cyp19a1 promoter activities, but FOXL2 could only enhance cyp19a1b promoter transcriptional activity in the HEK293T cell line. Furthermore, HEK293T and COS-7 cell lines showed that co-transfecting the two transcription factors significantly increased cyp19a1 promoter activity. The -120 to -112 bp (5'-CAAGGGCAC-3') and -890 to -872 bp (5'-AGAGGAGAACAAGGGGAG-3') regions of the cyp19a1a promoter were the core regulatory elements for SF-1 and FOXL2, respectively, to regulate cyp19a1b promoter transcriptional activity. Collectively, these results suggest that both FOXL2 and SF-1 are involved in giant wrasse sex reversal.

Transcriptomics integrated with metabolomics reveals the effect of Lycium barbarum polysaccharide on apoptosis in Nile tilapia (Oreochromis niloticus)

Lycium barbarum polysaccharide (LBP) is one of the main active ingredients in the fruit of L. barbarum L. It has been used as herbal medicine for thousands of years in China. In this study, Nile tilapia (Oreochromis niloticus) was taken as the research object. After feeding tilapia with 5 different doses of LBP (0 mg/kg, 500 mg/kg, 1000 mg/kg, 1500 mg/kg, 2000 mg/kg) for 55 d, it was found that LBP could promote the growth of tilapia, and this effect was the strongest at Group 1500 mg/kg. Apoptosis analysis in the liver and spleen showed that dietary supplementation with 1000 mg/kg LBP had the best protective effect on the spleen and liver in tilapia. Combined transcriptomics and metabolomics of the spleen in tilapia at Group 0 mg/kg and 1000 mg/kg showed that the differentially expressed genes (DEGs) such as NT5C2L1, pmm1, FasL and the differentially metabolites such as xanthine, dGMP, guanine and glutamate were mainly concentrated in signaling pathways such as Purine metabolism and FoxO signaling pathway. In conclusion, LBP regulates the metabolic waste levels of tilapia mainly through Purine metabolism and the FoxO signaling pathway, thereby inhibiting cell apoptosis, improving the utilization of nutrients, and promoting the growth of tilapia. This study not only provides a theoretical basis for the application of LBP in aquatic animals but also provides useful information for the healthy development of the aquaculture.

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.

Comprehensive Draft Genome Analyses of Three Rockfishes (Scorpaeniformes, Sebastiscus) via Genome Survey Sequencing

Sebastiscus species, marine rockfishes, are of essential economic value. However, the genomic data of this genus is lacking and incomplete. Here, whole genome sequencing of all species of Sebastiscus was conducted to provide fundamental genomic information. The genome sizes were estimated to be 802.49 Mb (S. albofasciatus), 786.79 Mb (S. tertius), and 776.00 Mb (S. marmoratus) by using k-mer analyses. The draft genome sequences were initially assembled, and genome-wide microsatellite motifs were identified. The heterozygosity, repeat ratios, and numbers of microsatellite motifs all suggested possibly that S. tertius is more closely related to S. albofasciatus than S. marmoratus at the genetic level. Moreover, the complete mitochondrial genome sequences were assembled from the whole genome data and the phylogenetic analyses genetically supported the validation of Sebastiscus species. This study provides an important genome resource for further studies of Sebastiscus species.

Role of sex steroids in fish sex determination and differentiation as revealed by gene editing

The involvement of sex steroids in sex determination and differentiation is relatively conserved among non-mammalian vertebrates, especially in fish. Thanks to the advances in genome sequencing and genome editing, significant progresses have been made in the understanding of steroidogenic pathway and hormonal regulation of sex determination and differentiation in fish. It seems that loss of function study of single gene challenges the traditional views that estrogen is required for ovarian differentiation and androgen is needed for testicular development, but it is not so in essence. Steroidogenic enzymes can be classified into two categories based on expression and enzyme activities in fish. One type, encoded by star2, cyp17a1 and cyp19a1a, is involved in estrogen production and exclusively expressed in the gonads. Mutation of these genes results in the up-regulation of male pathway genes and sex reversal from genetic female to male. The other type, encoded by the duplicated paralogs of the above genes, including star1, cyp11a1, cyp17a2 and cyp19a1b, as well as cyp11c1 gene, is dominantly expressed both in gonads and extra-gonadal tissues. Mutation of these genes alters the steroids (androgen, DHP and cortisol) production and spermatogenesis, fertility, secondary sexual characteristics and sexual behavior, but usually does not affect the sex differentiation. For the estrogen receptors (esr1, esr2a and esr2b), single mutation failed to, but double and triple mutation leads to sex reversal from female to male, indicating that at least Esr2a and Esr2b are required to mediate the role of estrogen in sex determination proved by gene editing experiments. Taken together, results from gene editing enrich our understanding of steroid synthesis pathways and further confirm the critical role of estrogen in female sex determination by antagonizing the male pathway in fish.

Genome-wide identification, evolution of histone lysine demethylases (KDM) genes and their expression during gonadal development in Nile tilapia

Histone lysine demethylases (KDM) are responsible for histone demethylation and are involved in gene expression regulation. Previous studies have shown that histone lysine demethylation plays an important role in gonadal development of vertebrates. The KDM family consists of eight subfamilies, i.e., kdm1, kdm2, kdm3, kdm4, kdm5, kdm6, kdm7 and JmjC-only subfamily. In this study, 13 to 63 KDM genes in 23 representative species were identified based on the available version of genome assembly. Phylogenetic relationships, domain architecture, and synteny of these genes were comprehensively analyzed and the results suggested KDM genes probably originated from the early diverging metazoan and significantly expanded in vertebrates with multiple whole genome duplication, especially in the third-round whole genome duplication (3R-WGD) and polyploidization of teleosts. The subfamilies of kdm2, kdm3, kdm4, kdm5, kdm6 and kdm7 were duplicated with 1R-2R events, and duplicates of kdm2a, kdm4a, kdm5b and kdm6b were resulted from 3R-WGD. Based on transcriptome data, the KDM genes were found to be dominantly expressed in the ovary and testis. More than 80% of KDM genes displayed sexual dimorphic expression, with 15 genes dominantly expressed in ovaries, and 12 genes dominantly expressed in testes. Importantly, from transcriptome data, qRT-PCR and fluorescence in situ hybridization during sex reversal, genes with higher expression in ovary than testis, such as kdm1b and two JmjC-only subfamily members hspbap1 and riox1, were downregulated, while other genes, such as kdm3c, kdm5bb, kdm6ba, kdm6bb and kdm7b, with higher expression in testis than ovary, were upregulated in ovotestis, indicating these genes play critical roles in the gonadal development and sex reversal. This study provided new insights into the evolution of the KDM genes and a fundamental clue for understanding their important roles in sex differentiation and gonadal development in teleosts.

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.

Dnmt3aa but Not Dnmt3ab Is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus)

Dnmt3a, a de novo methyltransferase, is essential for mammalian germ line DNA methylation. Only one Dnmt3a is identified in mammals, and homozygous mutants of Dnmt3a are lethal, while two Dnmt3a paralogs, dnmt3aa and dnmt3ab, are identified in teleosts due to the third round of genome duplication, and homozygous mutants of dnmt3aa and dnmt3ab are viable in zebrafish. The expression patterns and roles of dnmt3aa and dnmt3ab in gonadal development remain poorly understood in teleosts. In this study, we elucidated the precise expression patterns of dnmt3aa and dnmt3ab in tilapia gonads. Dnmt3aa was highly expressed in oogonia, phase I and II oocytes and granulosa cells in ovaries and spermatogonia and spermatocytes in testes, while dnmt3ab was mainly expressed in ovarian granulosa cells and testicular spermatocytes. The mutation of dnmt3aa and dnmt3ab was achieved by CRISPR/Cas9 in tilapia. Lower gonadosomatic index (GSI), increased apoptosis of oocytes and spermatocytes and significantly reduced sperm quality were observed in dnmt3aa^−/− mutants, while normal gonadal development was observed in dnmt3ab^−/− mutants. Consistently, the expression of apoptotic genes was significantly increased in dnmt3aa^−/− mutants. In addition, the 5-methylcytosine (5-mC) level in dnmt3aa^−/− gonads was decreased significantly, compared with that of dnmt3ab^−/− and wild type (WT) gonads. Taken together, our results suggest that dnmt3aa, not dnmt3ab, plays important roles in maintaining gametogenesis in teleosts.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Homozygous mutation of foxh1 arrests oogenesis causing infertility in female Nile tilapia†

Foxh1, a member of fox gene family, was first characterized as a transcriptional partner in the formation of the Smad protein complex. Recent studies have shown foxh1 is highly expressed in the cytoplasm of oocytes in both tilapia and mouse. However, its function in oogenesis remains unexplored. In the present study, foxh1-/- tilapia was created by CRISPR/Cas9. At 180 dah (days after hatching), the foxh1-/- XX fish showed oogenesis arrest and a significantly lower GSI. The transition of oocytes from phase II to phase III and follicle cells from one to two layers was blocked, resulting in infertility of the mutant. Transcriptomic analysis revealed that expression of genes involved in estrogen synthesis and oocyte growth were altered in the foxh1-/- ovaries. Loss of foxh1 resulted in significantly decreased Cyp19a1a and increased Cyp11b2 expression, consistent with significantly lower concentrations of serum estradiol-17β (E2) and higher concentrations of 11-ketotestosterone (11-KT). Moreover, administration of E2 rescued the phenotypes of foxh1-/- XX fish, as indicated by the appearance of phase III and IV oocytes and absence of Cyp11b2 expression. Taken together, these results suggest that foxh1 functions in the oocytes to regulate oogenesis by promoting cyp19a1a expression, and therefore estrogen production. Disruption of foxh1 may block the estrogen synthesis and oocyte growth, leading to the arrest of oogenesis and thus infertility in tilapia.

Transcriptome analysis of five ovarian stages reveals gonad maturation in female Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically important species of freshwater shrimp in China. Unlike other marine shrimps, the ovaries in adult female M. nipponense can mature rapidly and periodically during the reproductive period, but the resulting high stocking densities and environmental deterioration can negatively impact the harvest yield and economic benefits. To better understand ovary development in female M. nipponense, we performed systematic transcriptome sequencing of five different stages of ovarian maturation.

Results

We obtained 255,966 Gb of high quality transcriptome data from 15 samples. Of the 105,082 unigenes that were selected, 30,878 were successfully annotated. From these unigenes, we identified 17 differentially expressed genes and identified three distinct gene expression patterns related to different biological processes. We found that cathepins, legumains, and cystatin were enriched in the lysosome pathway, and they are related to vitellogenin hydrolysis. Additionally, we found that myosin heavy chain 67 participated in oocyte excretion.

Conclusions

We provide the first detailed transcriptome data relating to the ovarian maturation cycle in M. nipponense. Our results provide important reference information about the genomics, molecular biology, physiology, and population genetics of M. nipponense and other crustaceans. It is conducive to further solve the problem of M. nipponense rapid ovarian maturation from the aspects of energy supply and cell division.

Supplementary Information

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

Characterization, expression, and regulatory effects of nr0b1a and nr0b1b in spotted scat (Scatophagus argus)

Nuclear receptor subfamily 0 group B member 1 (Nr0b1) belongs to the nuclear receptor (NR) superfamily. It plays critical roles in sex determination, sex differentiation, and gonadal development in mammals. In this study, the duplicated genes nr0b1a and nr0b1b were identified in spotted scat (Scatophagus argus). Phylogenetic and synteny analyses revealed that, unlike nr0b1a, nr0b1b was retained in several species of teleosts after an nr0b1 gene duplication event but was secondarily lost in other fish species, amphibians, reptiles, birds, and mammals. In a sequence analysis, only 1.5 LXXLL-related repeat motifs were identified in spotted scat Nr0b1a, Nr0b1b, and non-mammalian Nr0b1a/Nr0b1, different from the 3.5 repeat motifs in mammalian Nr0b1. By qPCR, nr0b1a and nr0b1b were highly expressed in testes from stages IV to V and in ovaries from stages II to IV, respectively. Male-to-female sex reversal was induced in XY spotted scat by the administration of exogenous E2. A qPCR analysis showed that nr0b1b mRNA expression was higher in sex-reversed XY fish than in control XY fish, with no difference in nr0b1a. A luciferase assay showed that spotted scat Nr0b1a and Nr0b1b did not individually activate cyp19a1a gene transcription. As in mammals, spotted scat Nr0b1a suppressed Nr5a1-mediated cyp19a1a expression, despite containing only 1.5 LXXLL-related repeat motifs in its N-terminal region, while Nr0b1b stimulated Nr5a1-mediated cyp19a1a transcription. These results demonstrated that nr0b1a and nr0b1b in spotted scat have distinct expression patterns and regulatory effects and further indicate that nr0b1b might be involved in ovarian development by regulating Nr5a1-mediated cyp19a1a expression.

Global expression response of genes in sex-undifferentiated Nile tilapia gonads after exposure to trace letrozole

The aromatase inhibitor letrozole can be found in rivers, effluents, and even drinking water. Studies have demonstrated that letrozole affects various metabolic pathways and may cause reproductive toxicity, especially in fish exposed during development. However, studies on the effect of a low concentration of letrozole at the whole-gonad transcriptomic level in the early stage of fish sexual development have not been investigated. The aim of our study was to explore the potential effects of a low concentration of letrozole on the gonad transcriptome of Nile tilapia at an early stage of sexual development. In this study, 9 dpf (days postfertilization) Nile tilapia were exposed to trace letrozole for 12 days. Letrozole exposure from 9 dpf to 21 dpf persistently altered phenotypic sex development and induced the male-biased sex ratio. The transcriptome results showed that 1173 differentially expressed genes (DEGs) were present in the female control vs 1.5 μg/L letrozole-treated female comparison group and that 1576 DEGs were present in the 1.5 μg/L letrozole-treated female vs male control comparison group. Differentially expressed gene enrichment analysis revealed several crucial pathways, including the drug metabolism-cytochrome P450 pathway, the ErbB-PI3K/Akt/mTOR pathway, and the calcium signalling pathway. Further analysis of these identified DEGs indicated that some key genes correlated with metabolism and epigenetic regulation were significantly affected by letrozole, such as UDP-glucuronosyltransferase (Ugt), glutathione S-transferase omega-1 (Gsto1), lysine-specific demethylase 6bb (Kdm6bb, original name is Kdm6a), jumonji and AT-rich interaction domain containing 2 (Jarid2b, original name is Jarid2), growth arrest and DNA damage inducible gamma (Gadd45g), and chromobox protein 7 (Cbx7). The qRT-PCR validation results for twelve DEGs showed that the Pearson's correlation of the log₁₀^{fold change} values between the qPCR and RNA-Seq results was 0.90, indicating the accuracy and reliability of the RNA-Seq results. Our study is the first to report the effect of letrozole on the transcriptome of gonads from fish during early-stage sexual development. These findings will be useful for understanding the toxic effects and molecular mechanisms of letrozole exposure at the early stage of gonad development on the sexual development of aquatic organisms.

Germline sexual fate is determined by the antagonistic action of dmrt1 and foxl3/foxl2 in tilapia

Germline sexual fate has long been believed to be determined by the somatic environment, but this idea is challenged by recent studies of foxl3 mutants in medaka. Here, we demonstrate that the sexual fate of tilapia germline is determined by the antagonistic interaction of dmrt1 and foxl3, which are transcriptionally repressed in male and female germ cells, respectively. Loss of dmrt1 rescued the germ cell sex reversal in foxl3Δ7/Δ7 XX fish, and loss of foxl3 partially rescued germ cell sex reversal but not somatic cell fate in dmrt1Δ5/Δ5 XY fish. Interestingly, germ cells lost sexual plasticity in dmrt1Δ5/Δ5 XY and foxl3Δ7/Δ7 XX single mutants, as aromatase inhibitor (AI) and estrogen treatment failed to rescue the respective phenotypes. However, recovery of germ cell sexual plasticity was observed in dmrt1/foxl3 double mutants. Importantly, mutation of somatic cell-specific foxl2 resulted in testicular development in foxl3Δ7/Δ7 or dmrt1Δ5/Δ5 mutants. Our findings demonstrate that sexual plasticity of germ cells relies on the presence of both dmrt1 and foxl3. The existence of dmrt1 and foxl3 allows environmental factors to influence the sex fate decision in vertebrates.

Identification, Expression and Evolution of Short-Chain Dehydrogenases/Reductases in Nile Tilapia (Oreochromis niloticus)

The short-chain dehydrogenases/reductases (SDR) superfamily is involved in multiple physiological processes. In this study, genome-wide identification and comprehensive analysis of SDR superfamily were carried out in 29 animal species based on the latest genome databases. Overall, the number of SDR genes in animals increased with whole genome duplication (WGD), suggesting the expansion of SDRs during evolution, especially in 3R-WGD and polyploidization of teleosts. Phylogenetic analysis indicated that vertebrates SDRs were clustered into five categories: classical, extended, undefined, atypical, and complex. Moreover, tandem duplication of hpgd-a, rdh8b and dhrs13 was observed in teleosts analyzed. Additionally, tandem duplications of dhrs11-a, dhrs7a, hsd11b1b, and cbr1-a were observed in all cichlids analyzed, and tandem duplication of rdh10-b was observed in tilapiines. Transcriptome analysis of adult fish revealed that 93 SDRs were expressed in more than one tissue and 5 in one tissue only. Transcriptome analysis of gonads from different developmental stages showed that expression of 17 SDRs were sexually dimorphic with 11 higher in ovary and 6 higher in testis. The sexually dimorphic expressions of these SDRs were confirmed by in situ hybridization (ISH) and qPCR, indicating their possible roles in steroidogenesis and gonadal differentiation. Taken together, the identification and the expression data obtained in this study contribute to a better understanding of SDR superfamily evolution and functions in teleosts.

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.

Pomc Plays an Important Role in Sexual Size Dimorphism in Tilapia

Sexual dimorphism is common across the animal kingdom. Knowledge of the mechanisms of sexual size dimorphism is limited although it is important in biology and aquaculture. Tilapia is the common name for ~ 100 species of cichlid fish. Some are important aquaculture species and males outgrow females. To gain novel insights into the mechanisms underlying sexual size dimorphism, we analyzed the differences of brain transcriptomes between males and females in Mozambique tilapia and studied the function of the pro-opiomelanocortin (Pomc) gene in tilapia and zebrafish. The transcriptome analysis identified 123, 55, and 2706 sex-biased genes at 5, 30, and 90 dph (days post-hatch), respectively, indicating sexual dimorphism of gene expressions in the brain. The expression of Pomc in the tilapia brain was a female-biased at 30, 90, and 120 dph. An analysis of the DNA sequence located upstream of the tilapia Pomc transcriptional start site identified two estrogenic response elements. In vitro luciferase assay of the two elements revealed that β-estradiol significantly enhanced the expression of luciferase activity, suggesting that the expression of Pomc is mediated by estrogen. We knocked out Pomc in zebrafish using Crispr/Cas-9. The Pomc-knockout zebrafish showed faster growth and higher sensitivity to feeding as compared to the wild-type fish. Taken together, our results indicate that Pomc contributes to sexual size dimorphism and suggest that the high estrogen level in females promotes the expression of Pomc and suppresses feeding in female tilapias, which leads to the slower growth of female tilapias.

Flow cytometry in the analysis of hematological parameters of tilapias: applications in environmental aquatic toxicology

Blood tissue has been used to assess animal health and the environment in which they live. This tissue is easily acquired and has the ability to respond to various adverse conditions. Several techniques have been employed in the detection of xenobiotic-induced cell damage in blood cells. In general, traditionally used technologies, such as cellular analysis in blood smears, are time-consuming and require great analytical capacity. The present study proposes flow cytometry as a method to detect changes in blood cell populations. Tilapia (Oreochromis niloticus) was selected as a model for plotting the profile of fish blood cell populations after exposure to xenobiotics without euthanizing animals or using cell markers. Populations of erythrocytes and lymphocytes were detected only by combining the techniques of FACSAria cell sorting and light microscopy. Systemic deleterious effects were found through blood analysis, such as an increased lymphocyte-rich population at 48 h of exposure followed by a subsequent decrease. Moreover, the time-dependent expression of Nrf2 suggests its participation in increased membrane disruption, indicating it has a central role in erythrocyte lifespan. The present results shed light on the viability of using flow cytometry for blood analysis of living fish.

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.

Comparative transcriptome analysis of three gonadal development stages reveals potential genes involved in gametogenesis of the fluted giant clam (Tridacna squamosa)

BACKGROUND

Gonad development and differentiation is an essential function for all sexually reproducing species, and many aspects of these developmental processes are highly conserved among the metazoa. However, the mechanisms underlying gonad development and gametogenesis remain unclear in Tridacna squamosa, a large-size bivalve of great ecological value. They are protandrous simultaneous hermaphrodites, with the male gonad maturing first, eventually followed by the female gonads. In this study, nine gonad libraries representing resting, male and hermaphrodite stages in T. squamosa were performed to identify the molecular mechanisms.

RESULTS

Sixteen thousand four hundred ninety-one unigenes were annotated in the NCBI non-redundant protein database. Among the annotated unigenes, 5091 and 7328 unigenes were assigned to Gene Ontology categories and the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway database, respectively. A total of 4763 differentially expressed genes (DEGs) were identified by comparing male to resting gonads, consisting of 3499 which were comparatively upregulated in males and 1264 which were downregulated in males. Six hundred-ninteen DEGs between male and hermaphroditic gonads were identified, with 518 DEGs more strongly expressed in hermaphrodites and 101 more strongly expressed in males. GO (Gene Ontology) and KEGG pathway analyses revealed that various biological functions and processes, including functions related to the endocrine system, oocyte meiosis, carbon metabolism, and the cell cycle, were involved in regulating gonadal development and gametogenesis in T. squamosa. Testis-specific serine/threonine kinases 1 (TSSK1), TSSK4, TSSK5, Doublesex- and mab-3-related transcription factor 1 (DMRT1), SOX, Sperm surface protein 17 (SP17) and other genes were involved in male gonadal development in Tridacna squamosal. Both spermatogenesis- (TSSK4, spermatogenesis-associated protein 17, spermatogenesis-associated protein 8, sperm motility kinase X, SP17) and oogenesis-related genes (zona pellucida protein, Forkhead Box L2, Vitellogenin, Vitellogenin receptor, 5-hydroxytryptamine, 5-hydroxytryptamine receptor) were simultaneously highly expressed in the hermaphroditic gonad to maintain the hermaphroditism of T. squamosa.

CONCLUSION

All these results from our study will facilitate better understanding of the molecular mechanisms underlying giant clam gonad development and gametogenesis, which can provided a base on obtaining excellent gametes during the seed production process for giant clams.

Effects of gonadotropin-releasing hormone analog (GnRHa) immunization on the gonadal transcriptome and proteome of tilapia (Oreochromis niloticus)

Gonadotropin releasing hormone (GnRH) plays an important role in the regulation of vertebrate reproduction. Studies have shown that immunization against GnRHa can induce sexually sterile tilapia. To explore the mechanism behind this, in this study, RNA-seq and data-independent acquisition (DIA) techniques were used to study the transcriptome and proteome of the gonad of tilapia immunized with GnRHa. 644 differentially expressed genes (80 upregulated and 564 downregulated) and 1150 differentially expressed proteins (351 upregulated and 799 downregulated) were identified. There were 209 genes with consistent differential expression patterns in the transcriptomic and proteomic analyses, of which 9 were upregulated and 200 downregulated, indicating that the gonad gene expression was inhibited by GnRHa immunization. The downregulated genes were particularly involved in the functions of single-organism process, binding, cellular process, metabolic process and catalytic activity, and associated with the pathways including ECM-receptor interaction, focal adhesion, cardiac muscle contraction and oxidative phosphorylation. The expression of six differentially expressed genes involved in the GnRH signaling pathway was all downregulated. In addition, several important functional genes related to gonadal development after GnRHa immunization were screened. This study confirmed the expression of corresponding genes was affected by GnRHa on the gonad development in tilapia at the molecular level, and laid a foundation for elucidating the mechanism of GnRHa immunization.