From Mendel laws to whole genetic association study to decipher the swine mulefoot phenotype

The swine mulefoot (SM) is a rare condition characterized by a non-cloven hoof due to the partial or total fusion of the phalanges. No comprehensive study has been conducted to identify associated markers with this phenotype until now. We aimed to characterize the association between SNP and the mulefoot phenotype using a Genome-Wide Association Study (GWAS). An experimental population was produced using a half-sib mating where the male had the mulefoot phenotype and the females (n = 6) had cloven hoofs. The cross resulted in 27 (47%) animals with the mulefoot characteristic and 30 (53%) normal animals, indicating the possible dominant gene action. Animals were further genotyped using the Illumina PorcineSNP50k BeadChip, and SNPs were tested for associations. Twenty-nine SNPs located on the SSC15, SSC4, and SSCX were associated with the mulefoot phenotype (p-value <5 × 10^-5). Six markers were found in the intronic regions of VWC2L, CATIP, PDK3, PCYT1B, and POLA1 genes. The marker rs81277626, on SSC15:116,886,110 bp, is located in the Von Willebrand Factor C Domain (VWC2L), a possible functional candidate gene. The VWC2L is part of a biological process involved with the bone morphogenetic protein (BMP) signaling pathway, previously associated with syndactyly in other species. In conclusion, the identified markers suggest the involvement of the VWC2L gene in the SM phenotype in this population.