Exploring Genomic Variability in the Mediterranean Buffalo Breed: A Step Towards Custom SNP Array

Genomic improvement and sustainable breeding of Mediterranean buffalo are hampered by the lack of breed-specific genomic resources compared with bovine. To address this gap, we aimed to identify a comprehensive set of high-confidence single-nucleotide polymorphisms (SNPs) in Mediterranean buffalo and evaluate their informativeness across other buffalo populations. A total of 58 whole-genome sequencing samples, from three different sets, were merged through the integration of short- and long-read sequencing technologies. Variants, both unique and shared among datasets, were identified using a combination of bioinformatic tools to increase the reliability of the dataset. From these, over 11 million high-confidence biallelic SNPs were identified in the Mediterranean breed. The same SNPs were also tested in additional populations (other rivers and swamps) to evaluate their polymorphism and missingness. These results provide a robust genomic resource for Mediterranean buffalo, overcoming the limitations of bovine-derived genotyping tools available today. The identified SNPs dataset lays the basis for the development of a breedspecific SNP array, providing a dense and informative set of markers that could support cost-effective SNP chip development compared with existing arrays and sequencing technologies. This resource will facilitate more accurate genomic selection, precision breeding, and the conservation of genetic diversity in the Mediterranean buffalo population.