Transcriptomic analysis of the HPG axis-related tissues reveals potential candidate genes and regulatory pathways associated with testicular size in Hu sheep.

Testicular size is an excellent proxy for selecting high-fertility rams. The hypothalamus-pituitary-gonadal (HPG) axis plays an important role in regulating reproductive capacity in vertebrates, while key genes and regulatory pathways within the HPG axis associated with testicular size remain largely unknown in sheep. This study comprehensively compared the transcriptomic profiles in the hypothalamus, pituitary and testis of rams after sexual maturity between the large-testis group (LTG, testicular weight = 454.29 ± 54.24 g) and the small-testis group (STG, testicular weight = 77.29 ± 10.76 g). In total, 914, 795 and 10518 differentially expressed genes (DEGs) were identified in the hypothalamus, pituitary and testis between LTG and STG, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DEGs were mainly involved in the biological processes of reproduction, biological regulation, and development process. Notably, the neuroactive ligand-receptor interaction and cAMP signaling pathways, commonly enriched by the DEGs in the hypothalamus and pituitary between two groups, were considered as two key signal pathways regulating testicular development through the HPGs axis. Weighted gene co-expression network analysis (WGCNA) identified two modules that were significantly associated with testicular size, and 97 key genes were selected with high module membership (MM) and gene significance (GS) in these two modules. Finally, a protein-protein interaction (PPI) network was constructed, and ten genes with the highest degree were represented as hub genes, including FOS, NPY, SST, F2, AGT, NTS, OXT, EDN1, VIP and TAC1. Taken together, these results provide new insights into the molecular mechanism underlying the HPG axis regulating testicular size of Hu sheep.