Comparative Transcriptome Analysis Reveals Differentially Expressed Genes and Signaling Pathways Between Male and Female Red-Tail Catfish (Mystus wyckioides)

Gonadal Transcriptomic Analysis Reveals Novel Sex-Related Genes in Bactrocera dorsalis

Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) is one of the most devastating agricultural pests worldwide due to its high reproductive and invasive abilities. The elucidation of its gonadal developmental characteristics and the identification of sex-related genes will provide a useful genetic basis for reproductive-based pest control. Here, the gonadal transcriptome of B. dorsalis was sequenced, and novel gonad-specific expressed genes were analyzed. A total of 1338, 336, 35, and 479 differentially expressed genes (DEGs) were found in the testis (TE), ovary (OV), female accessory gland (FAG), and male accessory gland (MAG), respectively. Furthermore, 463 highly expressed gonad-specific genes were identified, with the TE having the highest number of specific highly expressed genes, at 402, followed by 51 in the OV, 9 in the MAG, and only 1 in the FAG. Strikingly, approximately half of highly expressed gonad-specific genes were uncharacterized. Then, it was found that 35, 17, 3, 2, and 1 of 202 uncharacterized highly expressed TE-specific genes encoded proteins that contained transmembrane domains, signal peptides, high-mobility group boxes, the zinc finger domain, and the BTB/POZ domain, respectively. Interestingly, approximately 40% of uncharacterized highly expressed gonad-specific genes encoding proteins were not predicted to possess functional motifs or domains. Finally, the spatiotemporal expression and sequence characterization of six novel highly expressed gonad-specific genes were analyzed. Altogether, our findings provide a valuable dataset for future functional analyses of sex-related genes and potential target sites for pest control.

Comparative transcriptome analysis of ovaries and testes reveals sex-biased genes and pathways in zebrafish

Zebrafish (Danio rerio) is a widely recognized and extensively studied model organism in scientific research. The regulatory mechanism of gonadal development and differentiation of this species has aroused considerable attention. Nonetheless, the major sex-biased genes and pathways associated with gonadal development remain elusive. Therefore, to comprehend this intricate process, gonadal transcriptome sequencing was carried out to identify differentially expressed genes (DEGs) between the testes and ovaries of adult zebrafish. The preliminary assessment yielded a total of 23,529,272 and 23,521,368 clean reads from the cDNA libraries of ovaries and testes. Afterward, a comparative analysis of the transcriptome revealed 3,604 upregulated and 11,371 downregulated DEGs in the ovaries compared to the testes. Of these genes, 428 were exclusively expressed in females, while 3,516 were exclusively expressed in males. Additionally, further assessments were conducted to explore the functions associated with these DEGs in various biological processes. The data revealed their involvement in sex-biased pathways, such as progesterone-mediated oocyte maturation, oocyte meiosis, cytokine-cytokine receptor interaction, and cardiac muscle contraction. Finally, the expression levels of 14 sex-biased DEGs (cdc20, ccnb1, ypel3, chn1, bmp15, rspo1, tnfsf10, egfra, acta2, cox8a, gsdf, dmrt1, star, and cyp17a1) associated with the enriched pathways were subjected to further validation through qRT-PCR. The data acquired from these investigations offer valuable resources to support further exploration of the mechanisms governing sexual dimorphism and gonadal development in zebrafish.

Comparative transcriptome analysis of gonads in male and female Pseudobagrus ussuriensis (Bagridae, Siluriformes)

As an important aquaculture fish in the Heilongjiang River Basin, Pseudobagrus ussuriensis has high economic value, and all-male culture is beneficial to the economic development of this fish. In this study, the transcriptomes of gonads in males and females were analyzed, and some genes related to gonad development were found. A total of 82,931 unigenes were found (average length 1504 bp, N50 1829 bp). In addition, 4689 differentially expressed genes (DEGs; including 1424 genes upregulated and 3265 genes downregulated in males) were identified. Some genes associated with testis development (such as Dmrt1 and Ropn1l) were significantly upregulated in males, while genes related to ovary development (such as Wnt2, PLC, Cyp19a, ZP3) were significantly downregulated in males, demonstrating that these genes have a crucial influence on gonad development in P. ussuriensis. Some signaling pathways related to gonad development were found, such as the Wnt pathway and oocyte meiosis. The results of RNA-seq obtained in this study provide theoretical data for elucidating the potential mechanism of gonad development of P. ussuriensis and reliable genomic data for the establishment of mono-sex breeding of P. ussuriensis.

Barbel regeneration and function divergence in red-tail catfish (Hemibagrus wyckioides) based on the chromosome-level genomes and comparative transcriptomes

Catfish (Siluriformes) are one of the most diverse vertebrate orders and are characterized by whisker-like barbels, which are important sensory organs in most of teleosts. However, their specific biological functions are still unclear. Red-tail catfish (Hemibagrus wyckioides) is well-known catfish species with four pairs of barbels, of which the maxillary barbels reach two-thirds of the body length. In this study, a 776.58 Mb high-quality chromosome-level genome was assembled into 29 chromosomes. Comparative genome data indicated that the barbeled regeneration gene ccl33 has expanded into 11 tandemly duplicated copies. Transcriptome data revealed the functional differentiation of different barbels and suggested that the maxillary barbel might be necessary for water temperature perception. Taste receptor genes were also characterized in teleosts with different food habits. Selection pressures were revealed to affect the sugar-based solute transport domain of the sweet taste receptor gene t1r2 in carnivorous fishes. In addition, the bitter taste receptor gene t2r200 was found to be lost from the genomes of four catfish species. Therefore, our study provides a genomic foundation for understanding the regeneration and functional differentiation of barbels in red-tail catfish and also reveals novel insights into the feeding evolution of fish species with different feeding habits.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

Advances in Reproductive Endocrinology and Neuroendocrine Research Using Catfish Models

Catfishes, belonging to the order siluriformes, represent one of the largest groups of freshwater fishes with more than 4000 species and almost 12% of teleostean population. Due to their worldwide distribution and diversity, catfishes are interesting models for ecologists and evolutionary biologists. Incidentally, catfish emerged as an excellent animal model for aquaculture research because of economic importance, availability, disease resistance, adaptability to artificial spawning, handling, culture, high fecundity, hatchability, hypoxia tolerance and their ability to acclimate to laboratory conditions. Reproductive system in catfish is orchestrated by complex network of nervous, endocrine system and environmental factors during gonadal growth as well as recrudescence. Lot of new information on the molecular mechanism of gonadal development have been obtained over several decades which are evident from significant number of scientific publications pertaining to reproductive biology and neuroendocrine research in catfish. This review aims to synthesize key findings and compile highly relevant aspects on how catfish can offer insight into fundamental mechanisms of all the areas of reproduction and its neuroendocrine regulation, from gametogenesis to spawning including seasonal reproductive cycle. In addition, the state-of-knowledge surrounding gonadal development and neuroendocrine control of gonadal sex differentiation in catfish are comprehensively summarized in comparison with other fish models.

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.

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.

Comparative Transcriptome Analysis of Gonads for the Identification of Sex-Related Genes in Giant Freshwater Prawns (Macrobrachium Rosenbergii) Using RNA Sequencing

The giant freshwater prawn (Macrobrachium rosenbergii) exhibits sex dimorphism between the male and female individuals. To date, the molecular mechanism governing gonadal development was unclear, and limited data were available on the gonad transcriptome of M. rosenbergii. Here, we conducted comprehensive gonadal transcriptomic analysis of female (ZW), super female (WW), and male (ZZ) M. rosenbergii for gene discovery. A total of 70.33 gigabases (Gb) of sequences were generated. There were 115,338 unigenes assembled with a mean size of 1196 base pair (bp) and N50 of 2195 bp. Alignment against the National Center for Biotechnology Information (NCBI) non-redundant nucleotide/protein sequence database (NR and NT), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, SwissProt database, Protein family (Pfam), Gene ontology (GO), and the eukaryotic orthologous group (KOG) database, 36,282 unigenes were annotated at least in one database. Comparative transcriptome analysis observed that 10,641, 16,903, and 3393 genes were significantly differentially expressed in ZW vs. ZZ, WW vs. ZZ, and WW vs. ZW samples, respectively. Enrichment analysis of differentially expressed genes (DEGs) resulted in 268, 153, and 42 significantly enriched GO terms, respectively, and a total of 56 significantly enriched KEGG pathways. Additionally, 23 putative sex-related genes, including Gtsf1, IR, HSP21, MRPINK, Mrr, and other potentially promising candidate genes were identified. Moreover, 56,241 simple sequence repeats (SSRs) were identified. Our findings provide a valuable archive for further functional analyses of sex-related genes and future discoveries of underlying molecular mechanisms of gonadal development and sex determination.

Transcriptomic Analysis of the Cold-Pretreated Larimichthys crocea Showing Enhanced Growth Fitness in Cold Water

The large yellow croaker Larimichthys crocea is an economically important marine fish species endemic to China and East Asia. Ningde area of Fujian Province is a major L. crocea aquaculture and spawning center in China. L. crocea cultivated at the Zhoushan area appears to be popular but suffered high mortality in cold water during winter seasons. To reduce the mortality rate, we pretreated fish with cold shocks prior to shift to cold water. In this study, we show that cold-pretreated L. crocea 12 days after shift to cold water increase the viability by 5.77-fold compared to the unpretreated (live fish 75 versus 13, p value = 1.775e-06, n = 100). The highest loss of 31 out of 100 fish in the unpretreated group occurred in day 3 after temperature shift. To identify the pretreatment-induced transcriptional changes that may be attributed to cold-resistance and survival, we performed RNA-seq analysis of a total of 48 fish that were prior to and 48 h, 54 h, and 72 h after temperature shift in pretreated and unpretreated groups in sextuplicate. Transcriptomic profiling analysis indicates that pretreatment-induced transcriptional alterations of enzymes involved in FASI, β-oxidation, PUFA synthesis, oxidative phosphorylation, and molecular chaperones persisted after temperature shift, suggesting that these metabolic pathways may play a role in L. crocea cold-resistance and survival. Our study provides insights on how the pretreatment enhances the L. crocea growth fitness in cold water.

Identification and characterization of microRNAs in the gonads of Crassostrea hongkongensis using high-throughput sequencing

Crassostrea hongkongensis is one of the three most-commonly cultivated oyster species in China. Although microRNAs (miRNAs) expression in the gonads have been widely investigated, few studies of miRNAs in mollusk gonads are available, particularly in oyster. In the present study, we analyzed the miRNAs expressed in the ovaries and testes of C. hongkongensis. We obtained 14,166,409 and 15,133,900 raw reads from the ovaries and testes, respectively, yielding 13,634,997 (ovarian) and 14,494,149 (testicular) 18-35-nt sequences. We mapped these sequences to the C. hongkongensis genome reference sequence, and identified 8,771,717 (ovarian) and 9,926,014 (testicular) sequences corresponding to miRNAs in the Rfam database. After blasting the miRNA sequences against the miRBase database, we identified 50 known mature miRNAs and 53 novel miRNAs. Of these, 27 miRNAs were significantly upregulated in ovaries as compared to the testes, and 43 miRNAs were significantly upregulated in the testes as compared to the ovaries. To validate the differential expression results generated by Illumina sequencing, we used RT-real-time quantitative PCR (RT-qPCR) to characterize the expression patterns of the six most differently expressed miRNAs (lgi-miR-1990, lgi-miR-1986, lgi-miR-263b, lgi-miR-279, lgi-miR-1992, and novel_98) as well as two miRNAs associated with gonad development (lgi-miR-29 and lgi-miR-8). Most of the RT-qPCR miRNA expression patterns were similar to those recovered by high-throughput sequencing with the exceptions of novel_98 and lgi-miR-1992. Gene Ontology (GO) annotations indicated that the multi-organism cellular process GO category was enriched with the target genes of the differentially expressed miRNAs. Additionally, the target genes were enriched in several KEGG pathways, including the ECM-receptor interaction, galactose metabolism, phagosome, and notch signaling pathway. These pathways are involved in gonadal differentiation and the maintenance of gonad function. This identification and characterization of the miRNAs differentially expressed between the ovaries and testes of C. hongkongensis will increase our understanding of the role of miRNAs in gonad differentiation in the oyster.