Genome-wide association study identifies candidate genes affecting body conformation traits of Zhongwei goat

Identifying markers or genes crucial for growth traits in Zhongwei goats is pivotal for breeding. Pinpointing genetic factors linked to body size gain enhances breeding efficiency and economic value. In this study, we used the MGISEQ-T7 platform to re-sequence 240 Zhongwei goats (133 male; 107 female) belonging to 5 metrics of growth traits at different growth stages (40 days and 6 months, here in after referred to as 40d and 6 m), namely, Body Weight (BW), Body Height (BH), Body Length (BL), Chest Circumference (CC), Tube Circumference (TC) were examined. A Genome-wide association study (GWAS) was conducted to identify candidate genes associated with the five indicators of body conformation traits, thereby establishing a foundation for subsequent investigations into the biological functions of these genes. A total of 19.89 Tb of raw data was generated with an average sequencing depth of about 20×. After quality control, 15,958,716 SNPs were available for the analysis. A total of 342 genome-wide significant SNPs were obtained. Among them, in the two physiological stages of 40d and 6 m, 147 and 32 SNPs were significantly associated with BW; 1 and 4 SNPs were significantly associated with BH; 19 and 6 SNPs were significantly associated with BL; 33 and 64 SNPs were significantly associated with CC, 34 and 2 SNPs were significantly associated with TC. These SNPs were annotated to 425 candidate genes. Finally, A total of 39 candidate genes are closely related to biological processes such as skeletal muscle development, skeletal formation, carcass quality, and embryonic development, where ADIPOQ, CCDD39, PTPRT, ZNF215, VRTN, ABCD4, DLST, ADAMTS2, ROBO1, AKAP13, AQPI, SOX2, and AHSG were identified as an important component of the genetic framework that may control somatic conformational traits in Zhongwei goats. which warrants further investigation and review. We verified the polymorphism of 11 SNPs by KASP, and found that Chr13_g.11,700,438 A > G, Chr15_g.37,120,328 A > G, Chr6_g.7,209,383 C > T, Chr20_g.51277932T > A, Chr19_g.17,078,199 A > G, and Chr1_g.79,943,276 C > T were significantly genotyped in verified populations (P