Characterization of biological pathways associated with semen traits in the Thai multibreed dairy population

The objective of this research was to characterize biological pathways associated with semen volume (VOL), number of sperm (NS), and sperm motility (MOT) of dairy bulls in the Thai multibreed dairy population. Phenotypes for VOL (n = 13,535), NS (n = 12,773), and MOT (n = 12,660) came from 131 bulls of the Dairy Farming Promotion Organization of Thailand. Genotypic data consisted of 76,519 imputed and actual single nucleotide polymorphisms (SNP) from 72 animals. The SNP variances for VOL, NS, and MOT were estimated using a three-trait genomic-polygenic repeatability model. Fixed effects were contemporary group, ejaculate order, age of bull, ambient temperature, and heterosis. Random effects were animal additive genetic, permanent environmental, and residual. Individual SNP explaining at least 0.001% of the total genetic variance for each trait were selected to identify associated genes in the NCBI database (UMD Bos taurus 3.1 assembly) using the R package Map2NCBI. A set of 1,999 NCBI genes associated with all three semen traits was utilized for the pathway analysis conducted with the ClueGO plugin of Cytoscape using information from the Kyoto Encyclopedia of Genes and Genomes database. The pathway analysis revealed seven significant biological pathways involving 127 genes that explained 1.04% of the genetic variance for VOL, NS, and MOT. These genes were known to affect cell structure, motility, migration, proliferation, differentiation, survival, apoptosis, signal transduction, oxytocin release, calcium channel, neural development, and immune system functions related to sperm morphology and physiology during spermatogenesis.