_1757587072 20250911103752000 Rovedar Daryoushbabazadeh@gmail.com Rovedar Small Animal Advances Small Anim. Adv. 2821-2363 09 01 2025 4 3 Umar Aziz https://orcid.org/0009-0009-6266-4340 Abdul Rehman https://orcid.org/0009-0009-3221-2984 Muhammad Hanzalah Yousaf https://orcid.org/0009-0005-7241-5701 Fasih Ur Rehman https://orcid.org/0000-0002-7750-3453 Muhammad Mushahid https://orcid.org/0000-0001-8670-3011 Nauman Khan https://orcid.org/0009-0008-1289-6859 Jiayuan Li Xugan Wang Hanbing Yan Xiaopeng An Introduction: Goat milk production is a vital economic trait, driven by rising global demand due to its digestibility, nutritional benefits, and hypoallergenic properties. To explore the genetic basis of milk yield, a genome-wide association study (GWAS) was conducted using whole-genome sequencing (WGS) data from two dairy goat breeds, Guanzhong (China) and Beetal (Pakistan). The present study aimed to identify genomic variants linked to milk yield in Guanzhong (China) and Beetal (Pakistan) goat breeds by performing an in silico GWAS using available WGS data.Materials and methods: Raw sequencing reads from both breeds were retrieved from public repositories and processed through an established bioinformatics pipeline. A GWAS was performed using a linear mixed model with GCTA and GEMMA, accounting for population structure and polygenic background. After quality control and alignment to the ARS1 goat reference genome using BWA, single-nucleotide polymorphisms (SNPs) were identified, and variants were filtered using SAMtools/BCFtools by applying thresholds of minor allele frequency (> 5%) and genotype call rate (> 90%). Population structure was assessed through principal component analysis and a genomic kinship matrix, both conducted using GCTA software, to control for stratification. The Manhattan plot revealed several genome-wide significant peaks, including loci near LALBA, PRLR, and SPP1, which are associated with lactation traits in dairy goats.Results: The GWAS revealed significant SNPs near LALBA on chromosome 19 (p = 1 × 10-¹⁰) and PRLR on chromosome X (p = 3.2 × 10-⁹) strongly associated with milk yield in Guanzhong and Beetal goats. In Guanzhong goats, SNPs near ANPEP, ADRA1A, and PRKG1 exhibited significant allele frequency differences, while in Beetal goats, SNPs near IGFBP3 and LEPR were linked to lactation traits. These loci provided robust genomic markers for enhancing the dairy goat breeding program.Conclusion: The present study demonstrated the feasibility of WGS-based GWAS in goats and identified candidate loci such as SPP1, ERBB4, and LALBA, previously linked to lactation, that may serve as genomic markers in future selection programs. 09 01 2025 11 19 https://creativecommons.org/licenses/by/4.0 10.58803/saa.v4i3.36 https://saa.rovedar.com/index.php/SAA/article/view/36 https://saa.rovedar.com/index.php/SAA/article/download/36/62 https://saa.rovedar.com/index.php/SAA/article/download/36/62 10.1111/age.13380 10.3389/fvets.2024.1409282 10.1038/s41597-025-04880-6 10.1038/s41598-023-35973-0 10.3389/fgene.2022.1044017 10.3168/jds.2022-22223 10.1017/S0022029924000013 10.1038/s41598-023-49125-x 10.3390/foods13121913 10.3390/ani14081246 10.1080/09712119.2023.2273945 10.3168/jds.2023-23284 10.1186/s40168-023-01651-6 10.1007/s11250-023-03606-3 10.1038/s41598-023-28644-7 10.1080/09712119.2023.2300624 10.1080/1828051X.2022.2160665 10.1186/s12891-023-06136-z