Genome Wide Scan to Identify Potential Genomic Regions Associated With Milk Protein and Minerals in Vrindavani Cattle

In this study, genome-wide association study (GWAS) was conducted for identifying significantly associated genomic regions/SNPs with milk protein and minerals in the 96 taurine-indicine crossbred (Vrindavani) cows using 50K SNP Chip. After quality control, a total of 41,427 SNPs were retained and were further analyzed using a single-SNP additive linear model. Lactation stage, parity, test day milk yield and proportion of exotic inheritance were included as fixed effects in GWAS model. Across all traits, 13 genome-wide significant (p < 1.20 x 10(−06)) and 49 suggestive significant (p < 2.41 x 10(−05)) SNPs were identified which were located on 18 different autosomes. The strongest association for protein percentage, calcium (Ca), phosphorus (P), copper (Cu), zinc (Zn), and iron (Fe) were found on BTA 18, 7, 2, 3, 14, and 2, respectively. No significant SNP was detected for manganese (Mn). Several significant SNPs identified were within or close proximity to CDH13, BHLHE40, EDIL3, HAPLN1, INHBB, USP24, ZFAT, and IKZF2 gene, respectively. Enrichment analysis of the identified candidate genes elucidated biological processes, cellular components, and molecular functions involved in metal ion binding, ion transportation, transmembrane protein, and signaling pathways. This study provided a groundwork to characterize the molecular mechanism for the phenotypic variation in milk protein percentage and minerals in crossbred cattle. Further work is required on a larger sample size with fine mapping of identified QTL to validate potential candidate regions.