Bioinformatic identification of Single Nucleotide Polymorphisms (SNPs) in keratin-associated protein genes in alpacas (Vicugna pacos)

Abstract

This study aims to identify single nucleotide polymorphisms (SNPs) within KRTAP genes in alpacas (Vicugna pacos), which play a fundamental role in defining their physico-mechanical properties and potentially the quality of alpaca fiber, the primary product of their breeding. Thirty-four KRTAP genes, as annotated in the reference genome VicPac3.1, were investigated. Utilizing the reference genome, along with nine additional genomes and reads from 300 reduced representation DNA libraries of alpacas, SNPs were identified. Minor allele frequency (MAF) and genotyping rates were computed using PLINK software, while Illumina Scores were determined for each SNP using Illumina Design Studio software. Markers meeting the criteria of MAF ≥ 0.05, genotyping rate >45%, and Illumina Score ≥ 0.6 per SNP were selected. A total of 67 SNPs were identified within intronic, exonic, and untranslated regions of KRTAP genes. Among these, 35 SNPs were incorporated into the 76K Alpaca SNP microarray, with 32 SNPs subsequently validated in a population of 936 alpacas. In conclusion, our findings delineate SNPs within KRTAPs that hold potential utility in genome-wide association studies, thereby facilitating the integration of modern breeding technologies into alpaca breeding programs.