Rajawat, Divya, Ghildiyal, Kanika, Sonejita Nayak, Sonali, Sharma, Anurodh, Parida, Subhashree, Kumar, Shive, Ghosh, A.K., Singh, Umesh, Sivalingam, Jayakumar, Bhushan, Bharat, Dutt, Triveni, and Panigrahi, Manjit
• The study aimed to elucidate genetic diversity and selective sweeps in 82 Sahiwal cattle using ddRAD data. • A declining trend in the effective population size (74 in most recent generation) observed in the population which necessitates a breeding plan to maintain a sufficiently large Ne. The genomic inbreeding coefficients ranged from −0.0289 to 0.0725. • The LD estimates revealed a minimum of 114,800 SNP markers and at least 14,350 informative markers necessary for GWAS analysis and genomic selection studies, respectively. • We found 146 regions undergoing selective sweeps using five distinct statistical tests and these regions, located in non-overlapping. • We identified various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK) , as well as production traits (NRXN1, PLCE1, GHR). The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (F HOM , F UNI , F ROH , and F GROM), ranged from −0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075–35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population. [ABSTRACT FROM AUTHOR]