Genome-wide CNV analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis)

Abstract

Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variant (CNVs) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 copy number variant regions (CNVRs), with 1,993 shared CNVRs being found within the studied buffalo types. Analyzing CNVRs highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVRs that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci (eQTL) analysis revealed differential expression of CNVR-driven genes (DECGs) associated with milk production traits. Notably, known milk production-related genes were among these DECGs, validating their relevance. Last, a genome-wide association study (GWAS) identified 3 CNVRs significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.