Comparative RNA-Seq Analysis of Differentially Expressed Genes in the Testis and Ovary of Mudskipper, Boleophthalmus pectinirostris

Sex determination and differentiation in teleosts are governed by complex genetic regulatory networks that include evolutionarily conserved mechanisms. In this study, we investigated Boleophthalmus pectinirostris, a Gobiidae species lacking heterogametic sex chromosomes, using comparative gonadal transcriptome analysis to identify sex differentially expressed genes (DEGs). RNA sequencing of ovarian and testicular tissues identified 17,214 DEGs, including 14,302 upregulated in males and 2912 upregulated in females. These DEGs were primarily associated with male (e.g., dmrt1, amh, amhr2) or female (e.g., bmp15, gdf9, rspo1) sex determination and differentiation, steroidogenesis (e.g., hsd17b1, hsd3b1, cyp17a1), and meiosis (e.g., cyp26b1, aldh1a2, piwil2). Functional enrichment analysis revealed that male upregulated DEGs were involved in spermatogenesis pathways such as calcium signaling, while female upregulated DEGs were associated with oogenesis pathways including oocyte meiosis and progesterone-mediated oocyte maturation. Conserved regulators, notably dmrt1 and amh, were predicted to act as key hubs in protein–protein interaction networks, being primarily associated with reproductive processes and sex differentiation in B. pectinirostris. The amh gene produces two alternatively spliced isoforms that differ by a partial deletion in the second exon, both expressed in ovaries and testes. Collectively, this study provides the first comprehensive molecular framework of sex determination and differentiation in Gobiidae species, offering critical insights into the regulatory mechanisms of B. pectinirostris reproductive development.