Abstract
Present study aimed to explore genomic basis of adaptation of Indian native cattle and to predict the impact of key SNPs on amino acid changes that affect protein function. Four native cattle breeds belonging to contrasting landscape and climatic conditions were genotyped using Illumina 777 K BovineHD BeadChip: Siri & Ladakhi from cold hilly areas, and Kankrej and Hallikar from hot arid and semi-arid regions, respectively. The R.SamBada package in R was used to perform the genotype-environment association analysis. A total of 1,12,780 significant (q < 0.05), models with 30,350 unique SNPs were obtained. Significantly associated SNPs had impact on 4,435 genes and 141 pathways. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence.. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Identified adaptive alleles might be responsible for the local adaptation of these cattle breeds.