Bioinformatics analysis of candidate genes for milk production traits in water buffalo (Bubalus bubalis)

Water buffalo (Bubalus bubalis) is an important livestock species in developing countries due to its contribution to meat, milk production, and a certain form of labor. However, the genetic potential of buffalo milk production traits has not been fully exploited. To date, 516 candidate genes associated with milk production traits of buffalo have been identified. The present study aimed to explore the possible molecular mechanisms underlying milk production traits of this species through functional genomics analysis of these candidate genes by using different bioinformatics tools. Gene ontology (GO) analysis indicated that these candidate genes were associated with complex biological processes, such as cell proliferation and mitotic nuclear division. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis demonstrated that these candidate genes were enriched in multiple signaling pathways, such as AMPK, ErbB, Toll-like receptor, and Jak-STAT. In addition, one function module consisting of 57 nodes and 139 edges were identified from the protein-protein interaction (PPI) network. GO analysis showed that the 57 candidate genes in this function module were enriched in three main biological processes, including homeostasis, metabolism, and cell response. These three distinct biological processes are well known for regulating mammary gland activities, which explained clearly the mechanism underlying milk production traits. This study provides a novel perspective for better understanding of the biological processes linked with milk production traits. This knowledge is conducive to the improvement of milk yield and composition of this species.