Whole Transcriptome Analysis of Male and Female Northern Pike (Esox lucius)

The northern pike (Esox lucius) is an economically important cold-water fish species in northern China. It exhibits pronounced sexual dimorphism, yet the molecular mechanism underlying its sex differentiation remains unclear, which hinders the development of aquaculture. Whole-transcriptome sequencing is a powerful approach for screening sex-related genes; however, no such study has been reported for this species to date. In this study, gonadal tissues from three female and three male E. lucius were collected for whole-transcriptome sequencing. A total of 14,941 differentially expressed messengerRNAs, 119 differentially expressed microRNAs, 229 differentially expressed circularRNAs, and 2055 differentially expressed long non-codingRNAs were identified. Functional enrichment analysis revealed that the differentially expressed genes were significantly enriched in pathways closely associated with sex differentiation, such as steroid hormone biosynthesis and oocyte meiosis. Several key sex-biased genes were identified, including female-biased genes (FANCL, DDX5, SRSF5B) and male-biased genes (STAR, FDX1B, ITGA2B). Furthermore, a competing endogenous RNA (ceRNA) regulatory network involving dre-miR-107b was constructed, which may represent a candidate for further investigation into sex differentiation in E. lucius. This study provides the first comprehensive whole-transcriptome dataset of female and male gonads in E. lucius, identifies key sex-biased genes and core pathways involved in its sex differentiation, and thereby identifies the dre-miR-107b-centered ceRNA network and key sex-biased genes (FANCL, DDX5, SRSF5B, STAR, FDX1B, ITGA2B) as core molecular players in sex differentiation of this species.