Identification of genomic regions associated with differences in fleece type in Huacaya and Suri alpacas (Vicugna pacos)

Abstract

The difference in fleece type is the distinguishing trait between the two types of alpacas (Vicugna pacos), Huacaya and Suri. The Suri fleece type has been found to be inherited dominantly over the Huacaya type, resulting in offspring with the Suri phenotype. The aim of our study was to map genomic regions associated with the two different fleece types. In this study, 91 alpacas (54 Huacayas and 37 Suris) from Germany and Switzerland were genotyped using the 76k alpaca SNP array. Only 59k chromosome‐localised markers map to the alpaca reference assembly VicPac3.1, and after quality control 49 866 SNPs, were retained for population structure assessment and to conduct a genome‐wide association study. Both principal component and neighbour‐joining tree analysis showed that the two fleece‐type cohorts overlapped rather than forming two distinct clusters. Genome‐wide significantly associated markers were observed in the scaffold region of chromosome 16 (NW_021964192.1), which contains a cluster of keratin genes. A haplotype predominantly found in Suri alpacas has been identified which supports dominant inheritance. Variant filtering of nine whole‐genome sequenced alpacas from both fleece types in the critical interval of 0.4 Mb did not reveal perfect segregation of either fleece type for specific variants. To our knowledge, this is the first study to use the recently developed species‐specific SNP array to identify genomic regions associated with differences in fleece type in alpacas. There are still some limitations, such as the preliminary status of the reference assembly and the incomplete annotation of the alpaca genome.