Genome-wide analysis reveals genetic diversity, linkage disequilibrium, and selection for milk traits in Chinese buffalo breeds

Abstract

Water buffalo holds the tremendous potential of milk and meat that widespread throughout central and southern China. However, characterization of the population genetics of Chinese buffalo is poorly understood. Using Axiom(R) buffalo genotyping array, we performed the genetic diversity, linkage disequilibrium (LD) pattern and signature of selection in the 176 Chinese buffaloes from thirteen breeds. A total of 35,547 SNPs passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between the swamp and river types. Ten Chinese indigenous breeds clustered into the swamp group, Murrah and Nili-Ravi breeds were the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosities compared to the river group. LD decay distance was much shorter in the swamp group compared with the river group with [Formula] value of approximately 50 Kb. Analysis of runs of homozygosity indicated that extensive remote and recent inbreeding activity was respectively found within swamp and river groups. Moreover, a total of 12 genomic regions under selection were detected between river and swamp groups. Further, 12 QTL regions were found associated with buffalo milk production traits. Some candidate genes within these QTLs were predicted to be involved in the cell structure and function, suggesting that these genes might play vital roles in the buffalo milk performance. Our data contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be utilized in buffalo breeding programs.\n\nAuthor SummaryIdentifying the causal genes or markers associated with important economic traits in livestock is critical to increasing the production level on the species. However, current understanding of the genetic basis for milk production traits in buffalo is limited. Here, we confirmed the divergent evolution, distinct population structure, and LD extent among Chinese buffalo breeds. We also identified 12 QTL regions associated with milk production traits in buffaloes using the selective sweeps and haplotype analysis. Further, a total of 7 genes involved in the cell structure and function were predicted within the identified QTLs. These findings suggested that these genes can serve as the candidate genes associated with buffalo milk production, which hold a vital role in the milk trait improvement of dairy buffalo industry.