Comparative transcriptome analysis identifies crucial candidate genes and pathways in the hypothalamic-pituitary-gonadal axis during external genitalia development of male geese

BACKGROUND: All birds reproduce via internal fertilization, but only ~3% of male birds possess the external genitalia that allows for intromission. Waterfowl (e.g., duck and goose) are representatives of them, and the external genitalia development of male geese is directly related to mating ability. Notably, some male geese show abnormal external genitalia development during ontogenesis. However, until now little is known about the molecular mechanisms of the external genitalia development in goose. In the present study, comparative transcriptomic analyses were performed on the hypothalamus, pituitary gland, testis, and external genitalia isolated from the 245-day-old male Tianfu meat geese showing normal (NEGG, n = 3) and abnormal (AEGG, n = 3) external genitals in order to provide a better understanding of the mechanisms controlling the development of the external genitalia in aquatic bird species. RESULTS: There were 107, 284, 2192, and 1005 differentially expressed genes (DEGs) identified in the hypothalamus, pituitary gland, testis and external genitalia between NEGG and AEGG. Functional enrichment analysis indicated that the DEGs identified in the hypothalamus were mainly enriched in the ECM-receptor interaction pathway. The ECM-receptor interaction, focal adhesion, and neuroactive ligand-receptor interaction pathways were significantly enriched by the DEGs in the pituitary gland. In the testis, the DEGs were enriched in the neuroactive ligand-receptor interaction, cell cycle, oocyte meiosis, and purine metabolism. In the external genitalia, the DEGs were enriched in the metabolic, neuroactive ligand-receptor interaction, and WNT signaling pathways. Furthermore, through integrated analysis of protein-protein interaction (PPI) network and co-expression network, fifteen genes involved in the neuroactive ligand-receptor interaction and WNT signaling pathways were identified, including KNG1, LPAR2, LPAR3, NPY, PLCB1, AVPR1B, GHSR, GRM3, HTR5A, FSHB, FSHR, WNT11, WNT5A, WIF1, and WNT7B, which could play crucial roles in the development of goose external genitalia. CONCLUSIONS: This study is the first systematically comparing the hypothalamus, pituitary gland, testis, and external genitalia transcriptomes of male geese exhibiting normal and abnormal external genitals. Both bioinformatic analysis and validation experiments indicated that the neuroactive ligand-receptor interaction pathway could regulate the WNT signaling pathway through PLCB1 to control male goose external genitalia development. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12864-022-08374-2.