Integrated Genome-Wide Association Studies and Selection Signature Analysis Reveal the Genetic Architecture of the Wattle Trait in Goats.

Abstract

Wattles are finger-like appendages on the ventral neck of goats, serving as a distinctive morphological marker for breed identification that serves potential implications for production performance. However, their genetic basis remains incompletely characterized. Here, we integrated a genome-wide association study (GWAS) and selection signature analysis to identify candidate genes and genomic regions associated with the wattle trait in goats. Using a linear mixed model, GWAS on 463 goats (23 wattled and 440 non-wattled) identified 385 quantitative trait loci (QTLs) at a 5% false discovery rate, yielding 346 candidate genes. The most significant association signal was observed on chromosome 10 (72.61-73.48 Mb), where the lead SNP (rs636481767) is located within a region containing GJD2, GREM1, and FMN1, showing strong linkage disequilibrium (r2 > 0.6) with surrounding loci. Subsequent selection signature analysis (23 wattled and 23 non-wattled) identified 83 genomic regions harboring 119 candidate genes. The strongest signals were detected at MFSD14B on chromosome 8 (FST = 0.154, log2π-ratio = 2.611) and PDLIM7 on chromosome 7 (FST = 0.144, log2π-ratio = 0.806). KEGG pathway enrichment analysis revealed that GWAS-associated genes were involved in glycosylation and immune responses, whereas selection-signature genes were enriched in DNA repair and the Hippo, Notch, and Wnt pathways. Furthermore, cross-species PheWAS revealed that human FMN1 is associated with dermatological, skeletal, and metabolic phenotypes, while porcine FMN1 is associated with backfat thickness and loin muscle depth. Overall, this study provides molecular markers of potential value for goat breeding and pinpoints key candidate genes for future functional validation of wattle development.