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.