Open Access
How to translate text using browser tools
23 September 2022 Correlation analysis of four KRTAP gene polymorphisms and cashmere fiber diameters in two cashmere goat breeds
Cuiling Wua, Chongkai Qin, Xuefeng Fu, Bingru Zhao, Yujiang Wu, Junmin He, Jingyi Mao, Jing Liu, Xixia Huang, Kechuan Tian
Author Affiliations +

Fiber diameter, a quantitative trait, is controlled by minor effect polygenes. Keratin-associated proteins (KRTAPs) are an important part of hair, and their rich polymorphisms facilitate the mining of cashmere trait molecular markers. In this study, Jiangnan and Tibetan cashmere goats were taken as the research object; multiplex PCR and exome sequencing technology were used to identify the exon regional polymorphisms of cashmere goats KRTAP15-1, KRTAP13.1, KRTAP27-1, and KRTAP24-1. The effects of mutation sites on the fiber diameter of cashmere were analyzed by least square method. The results showed that there were 28 mutation sites in the four KRTAP genes in Jiangnan cashmere goats and Tibetan cashmere goat populations. Among them, the KRTAP13.1, KRTAP27-1, and KRTAP24-1 gene polymorphisms were found to be significantly related to the fiber diameter of Jiangnan cashmere goats. The exploration of molecular markers in this study will help to improve the fiber diameter of the down, while the identification of gene polymorphisms will provide original data for the utilization and protection of germplasm resources of cashmere goats.

Le diamètre de la fibre, une caractéristique quantitative, est contrôlé par des polygènes à effet mineur. Les protéines associées à la kératine (KRTAPs — «keratin-associated proteins») sont parties importantes des poils, et leurs polymorphismes riches facilitent le minage des marqueurs moléculaires des caractéristiques du cachemire. Dans cette étude, les chèvres à cachemire de Jiangnan et du Tibet ont été utilisées comme objets de recherche. La PCR multiplexe et la technologie de séquençage d’exons ont été utilisées afin de déterminer les polymorphismes régionaux des exons des chèvres à cachemire KRTAP15-1, KRTAP13.1, KRTAP27-1, et KRTAP24-1. Les effets des sites de mutation sur le diamètre de la fibre du cachemire ont été analysés par la méthode des moindres carrés. Les résultats ont montré qu’il y avait 28 sites de mutation dans les quatre gènes KRTAP dans les populations de chèvres à cachemire de Jiangnan et du Tibet. Parmi ceux-ci, les polymorphismes de gènes KRTAP13.1, KRTAP27-1, et KRTAP24-1 se sont avérés reliés de façon significative au diamètre de la fibre des chèvres à cachemire de Jiangnan. L’exploration des marqueurs moléculaires dans cette étude aidera à améliorer le diamètre de la fibre du duvet, tandis que l’identification des polymorphismes des gènes offrira des données originales pour l’utilisation et la protection des ressources du patrimoine génétique des chèvres à cachemire. [Traduit par la Rédaction]


Jiangnan cashmere goat is a white cashmere goat breed cultivated in China, and its main area of production is Aksu, Xinjiang. Xinjiang Cashmere Goat Core Breeding Base (Akesu Comprehensive Experimental Station) is located in arid desert and semidesert grassland at an altitude of 2000 m. The annual average temperature is 6.5°C, the rainfall is 200 mm, 90% of the grassland is desert grassland, and the vegetation coverage rate is about 10%–30% (Qin et al. 2016). Tibetan cashmere goats are a local goat breed that has been retained after long-term natural and artificial selection and adapted to the special environment of high-altitude areas. They are also known as Tibetan Kashmir cashmere goats (Ci 2012; Fu 2021).

In recent years, with the increase in demand in the mutton market and the fall in grassroots breeding work, there has been a loss of original fine cashmere traits in both Jiangnan cashmere and Tibetan cashmere goats, seriously affecting the quality of cashmere and the sales price of cashmere (Tonin et al. 2002). Cashmere goats are the main source of income for local herdsmen. The breeding and popularization of high-quality cashmere goats has become an important means for these people to rise above poverty and become rich based on local resources. It is thus necessary for local people to breed ultrafine cashmere goats that are suitable for the extreme climates of Xinjiang and Tibet. Therefore, it is particularly important to find an effective and convenient breeding method which can accelerate the breeding process of cashmere goats in a short period of time. The classical breeding method mainly employs crossing, but this method has a slower genetic progress and a longer breeding cycle. With the rapid development of molecular biology theory and technology, molecular genetics and genetic engineering methods combined with traditional marker-assisted selection methods can provide better solutions to breeding problems (Niu and Hong-Bin 2008; Yang et al. 2017). The search for molecular markers related to cashmere traits is our primary focus (Hailemariam and Yadeta 2020).

In the hair cortex, hair keratins are embedded in an interfilamentous matrix consisting of keratin-associated proteins (KRTAPs), which are important for the formation of hair shaft as a result of disulphide bonds between cysteine residues. The genes encoding KRTAP are usually rich in polymorphism, which may be related to the cashmere traits of cashmere goats. At present, the polymorphisms of KRTAP15-1 (Wang et al. 2017b), KRTAP13.1 (Wu et al. 2018), KRTAP27-1 (Zhao et al. 2020), and KRTAP24-1 (Zhou et al. 2012) have been confirmed to affect fiber traits in some sheep and goat breeds, including Xinji fine-wool sheep, Merino × Southdown-cross sheep, Longdong cashmere goats, and Xinjiang local goat. However, it is unclear whether these genes can affect the fiber diameter of Jiangnan and Tibetan cashmere goats. Therefore, this study analyzed the effects of four gene polymorphism (KRTAP15-1, KRTAP13.1, KRTAP27-1, and KRTAP24-1) on cashmere fiber diameter of Jiangnan and Tibetan cashmere goats. The aim of this study was to obtain the molecular markers of cashmere fineness of cashmere goats, which would lay the foundation for the protection, development, and utilization of cashmere genetic resources, and provide theoretical basis for cashmere molecular breeding.

Materials and methods

Animal care

All animal experiments were strictly performed according to the guidelines established by the Animal Care and Use Committee of Xinjiang Academy of Animal Science (Approval number 2019009). Sample collection was carried out under license in accordance with the Guidelines for Care and Use of Laboratory Animals of China.

Experimental animals

A total of 353 two-year-old female Jiangnan cashmere goats were selected from the breeding center of Wenshu County, Aksu Prefecture, Xinjiang Province. They were from the Aerken group (n = 144), Samusak group (n = 79), Tuniazi group (n = 54), and Yiming group (n = 79). A total of 299 two-year-old female Tibetan cashmere goats were selected from the original breed of cashmere goats in Ngari Prefecture (Ritu County, n = 132; Gaize County, n = 13) and Nagqu Prefecture (Nima County, n = 57). Before the experiment, all cashmere goats are healthy and raised by grazing.

Sample collection

Cashmere samples were collected from Jiangnan and Tibetan cashmere goats at the 10 cm posterior margin of the scapula above the left midline of the body. Cashmere is naturally dried after washing according to its conventional washing process. A fiber diameter optical analyzer (OFDA2000) was used to determine the cashmere mean fiber diameter (MFD), fiber diameter standard deviation (FDSD), and coefficient of variation of fiber diameter (CVFD) under the conditions of a constant temperature of 20 ± 2 °C and a humidity of 65% ± 4%.

Corresponding to the cashmere sample, we collected 5 mL of experimental goat blood in an anticoagulation tube and stored it in a refrigerator at ?20 °C. A blood genomic DNA extraction kit (TIANGENG, USA) was used to extract DNA from cashmere goats. The quality and concentration of DNA were determined using 1.0% agarose gel electrophoresis and the Qubit 2.0 (Thermo, USA).

Multiple PCR and exome sequencing

According to the sequence published in the NCBI goat chromosome 1 accession number (NC_030808.1), the exon regions of KRTAP15-1, KRTAP13.1, KRTAP27-1, and KRTAP24-1 were selected. A Primer pool containing four exon regions of target genes was designed using Primer 5.0 and synthesized by Shanghai Sangong. The list of primers is shown in  Table S1 (cjas-2021-0112_suppla.docx). Then, the target SNP sequences of 353 Jiangnna cashmere goats and 299 Tibetan cashmere goats were amplified by two-step PCR and Illumina sequencing library was prepared. The two-step PCR system is shown in  Tables S2 (cjas-2021-0112_suppla.docx) and  S4 (cjas-2021-0112_suppla.docx), and the reaction procedure is shown in  Tables S3 (cjas-2021-0112_suppla.docx) and  S5 (cjas-2021-0112_suppla.docx). The final PCR product was purified and recovered using AMPure XP magnetic beads. Each PCR product was equally mixed and sequenced using a Hiseq XTEN sequencer (Illumina, San Diego, CA, US).

Sequencing data analysis and validation

Raw reads were filtered according to two steps: (i) removing any adaptor sequence the reads contained using cutadapt (v 1.2.1); (ii) removing low-quality bases from reads 3? to 5? (Q < 20) using PRINSEQ-lite (v 0.20.3). Additionally, the remaining clean data were mapped to the reference genome by BWA (version 0.7.13-r1126) with default parameters. Samtools (Version: 0.1.18) was used to calculate each genotype of the target site. Annovar was used to detect genetic variants.

The heterozygous individuals with four mutation sites selected for each gene were verified using first-generation sequencing technology. The first-generation sequencing results were assembled and corrected using the SeqMan program of the DNASTAR software, and the peak maps were compared with the BioEdit software.

Statistical analysis

The Popgene software was used to calculate the minor allele frequency (MAF) and Hardy–Weinberg equilibrium (HWE) of SNPs. The Linkage Format function of the Haploview software was used to analyze the linkage disequilibrium of SNPs. The GLM model in the SAS 9.2 software was used to analyze the influence of different SNPs genotypes on the cashmere fiber diameter. The results are expressed in the form of least squares mean (±standard error), and the linear model is


In the formula, Yick denotes the individual phenotypic value of cashmere goats, μ the population mean, Gi the genotype SNP effect, Fc the group effect, and eick denotes the random error.


Descriptive statistics of fiber diameter traits

A descriptive statistical analysis of the MFD, FDSD, and CVFD of Jiangnan cashmere goats and Tibet cashmere goats was performed. The basic statistics are shown in Table 1. The MFD, FDSD, and CVFD of Jiangnan cashmere goats were 15.69 µm, 3.26 µm, and 20.86%, respectively. The MFD, FDSD, and CVFD of Tibetan cashmere goats were 15.10 µm, 3.25 µm, and 21.53%, respectively. It is not difficult to see that, compared with Jiangnan cashmere goats, Tibetan cashmere goats have a finer MFD and a higher CVFD.

Table 1.

Descriptive statistics on fiber diameter traits of two Chinese cashmere goat breeds.


Quality control of sequencing data

The genomic DNA of cashmere goats was detected using 1% agarose gel electrophoresis, and the DNA bands were bright, as shown in  Fig. S1 (cjas-2021-0112_suppla.docx). A nucleic acid protein detector was used to detect DNA, and the OD ratio at 260–280 nm was between 1.8 and 2.1, indicating that the quality and purity of DNA meet the requirements for subsequent library construction. Quality control was performed on the data after the HiSeq XTen sequencer was sequenced. The average coverage ratio (Coverage) of each fragment and the target area sequence comparison was 96.09%; the average coverage depth (Mean_depth) was 4838.75; and the sequencing error percentage (Error_ratio) was 4.52%. It can be seen that the quality of the sequencing data is relatively high, which can satisfy the follow-up experiment. In addition, our sequencing results were submitted to the NCBI public database (PRJNA738549).

Classification results annotation and verification

Combined with multiplex PCR technology and exome sequencing, a total of 28 mutation sites were obtained in the four genes of Jiangnan cashmere goats and Tibet cashmere goats, including 19 missense mutations and 9 synonymous mutations (Table 2). KRTAP15-1, KRTAP13.1, KRTAP27, and KRTAP24 genes had 5, 10, 7, and 6 mutation sites, respectively. It is worth noting that the SNP12 and SNP18 mutations were missing in the Jiangnan cashmere goat population. There was no SNP17 mutation in the Tibetan cashmere goat population. Heterozygous individuals with mutation sites were randomly selected for first-generation sequencing; the results are shown in Fig. 1. It can be seen from the figure that the high-throughput typing results are consistent with the first-generation sequencing results, indicating that the exome sequencing results are reliable.

Fig. 1.

Validation of high-throughput sequencing results. (A) SNP1, SNP3, SNP4, and SNP5 of the KRTAP15-1 gene. (B) SNP6, SNP8, SNP14, and SNP15 of the KRTAP13.1 gene. (C) SNP19, SNP20, SNP21, and SNP22 of the KRTAP27 gene. (D) SNP23, SNP205, SNP27, and SNP28 of the KRTAP24-1 gene.


Table 2.

Nucleotide variant information of candidate genes.


Statistical analysis of SNP

The results of the genotyping and MAF and HWE analyses of 28 mutation sites in two cashmere goat populations are shown in Table 3. The MAFs of SNP11, SNP12, SNP16, SNP17, SNP18, SNP24, and SNP26 in the Jiangnan cashmere goat population were all less than 0.03. The MAFs of SNP11, SNP17, SNP18, SNP24, and SNP26 in the Tibetan cashmere goat population were all less than 0.03. SNP12 and SNP18 were not classified in Jiangnan cashmere goats, while SNP17 was not classified in Tibetan cashmere goats. Excluding SNP7, the remaining 27 SNPs complied with the Hardy–Weinberg equilibrium (P > 0.05).

Table 3.

The basic characters of 28 SNPs located in four genes.


Linkage disequilibrium analysis

We used the Haploview software to analyze the linkage disequilibrium of 28 SNP loci. When the MAF > 0.03 and HWE > 0.05, the linkage disequilibrium analysis of 21 SNPs was as shown in Fig. 2. The results indicate that we detected the existence of five blocks. SNP1, SNP2, SNP3, and SNP4 constitute block 1; SNP5, SNP6, SNP8, and SNP10 constitute block 2; SNP13 and SNP14 constitute block 3; SNP19, SNP20, SNP21, and SNP22 constitute block 4; and SNP23, SNP25, SNP27, and SNP28 constitute block 5.

Fig. 2.

Blocks found in the linkage disequilibrium (LD) analysis of SNPs. The block with a black border is a completely linked haplotype block. The color of each block in the picture changes from blue to red, indicating that the degree of linkage is becoming higher and higher. Each square number indicates the D? value.


The Tag SNP can be predicted using the linkage disequilibrium relationship between SNPs. The predicted Tag SNP can cover all SNP sites. When the R2 value >0.8, the 12 Tag SNPs are as listed in Table 4. The two Tag SNPs predicted by the KRTAP15-1 gene were SNP2 and SNP5. The KRTAP13.1 gene has four Tag SNPs, which are SNP6, SNP9, SNP10, and SNP14. The KRTAP27-1 gene has five Tag SNPs—namely, SNP16, SNP19, SNP20, SNP21, and SNP22. The KRTAP24-1 gene had only one Tag SNP, which was SNP28.

Table 4.

Summary of predictive tag SNP information.


SNP effect analysis

The correlation analysis between Tag SNP and the diameter of cashmere fiber using the SAS 9.2 software is shown in Table 5. In Jiangnan cashmere goats, Tag SNP6, Tag SNP19, and Tag SNP22 significantly affected the MFD (P < 0.01); Tag SNP9, Tag SNP14, and Tag SNP21 significantly affected the MFD (P < 0.05); and Tag SNP21 significantly affected the FDSD (P < 0.05). Tag SNP2 and Tag SNP19 significantly affected the CVFD (P < 0.05). In Tibetan cashmere goats, none of the mutation sites had any significant effects on MFD, FDSD, and CVFD.

Table 5.

Correlation analysis between mutation sites and cashmere fiber diameter.



The price of cashmere is influenced by the color and fineness of cashmere. Age, feeding management, and environment can affect the fiber diameter of a cashmere goat (Zhou et al. 2003). However, in addition to these nongenetic factors, genes are the fundamental causes of the fineness of cashmere. Recently, major cashmere goat breeds in China, such as the Liaoning cashmere goat (Yu et al. 2014; Zheng et al. 2019, 2020b), the Inner Mongolia cashmere goat (Zheng et al. 2020a), and the Tibetan cashmere goat (Fu et al. 2020), have been studied in terms of their fiber diameter at the transcriptomic level. In terms of genome, we prefer the use of genome-wide association analysis (Li et al. 2017; Zheng et al. 2020) and candidate gene polymorphism analysis (Zhao et al. 2008; Wang et al. 2017a, 2018; Zhang et al. 2019) to explore genes and molecular markers associated with cashmere fiber diameter.

In the hair cortex, hair keratin intermediate filaments are embedded in an interfilamentous matrix, consisting of hair KRTAP, which is essential for the formation of a rigid and resistant hair shaft through extensive disulfide bond cross-linking with the abundant cysteine residues of hair keratins (Strasser et al. 2015). KRTAP13, KRTAP15, KRTAP24, and KRTAP27 are all high-sulfur keratin-associated proteins (HS-KRTAP; <30 mol% cysteine) (Gong et al. 2012, 2016). Studies have shown that a change in amino acids may cause the loss of phosphorylation sites in the process of the post-translational modification of proteins, as well as leading to changes in the net charge of proteins (Gong et al. 2011).

We identified a total of 28 mutant sites in both Jiangnan and Tibetan cashmere goats. The SNP12 and SNP18 mutations were absent in the Jiangnan cashmere goat population. At present, the SNP18 mutation site has not been reported in other goat breeds. The MAF of SNP12 and SNP18 in the Tibetan cashmere goat population was low. Similarly, the SNP17 mutation was found to be missing in the Tibetan cashmere goat population, and the MAF was also low in Jiangnan cashmere goats. On the one hand, this suggests that we need to expand our sample size in future research; on the other hand, this also reflects the diversity of species.

There are five mutation sites (SNP1-SNP5) in the KRTAP15-1 gene of both Jiangnan and Tibetan cashmere goats, among which SNP1 is a synonymous mutation and SNP2-SNP5 is a missense mutation. These five mutations were also found in Tan sheep and Hu sheep (Wang et al. 2017b). Among them, SNP1–SNP4 are strongly linked (D? = 1). The analysis of Tag SNP2 shows that it is significantly correlated with the CVFD of the Jiangnan cashmere goat population (P > 0.05), and the haplotype “GTGA” is the dominant genotype. Among 396 Merino × Southdown-cross sheep, the KTAP15-1 polymorphism was found to be significantly correlated with wool yield and fiber diameter standard deviation (P < 0.05) (Li et al. 2018). This suggests that KRTAP15-1 may have an effect on the villus traits of sheep and goats.

Recently, KRTAP 13.1 in both sheep and goats has received a lot of research. Sun et al. (2014) found that three synonymous mutations in the KRTAP13.1 gene were significantly related to wool length and wool yield on Xinji fine wool sheep (P < 0.05). Yu et al. (2014) used Solexa technology to find the differential expression of the KRTAP13-1 gene in the skin tissues of the fine cashmere group and the coarse cashmere group in Liaoning cashmere goats. This suggests that the expression level of KRTAP13-1 may affect the fiber diameter of cashmere goats. Li et al. (2013) found five mutation sites in the KRTAP13.1 gene in Hexi mountain cashmere goats, Huanxian cashmere goats, Liaoning cashmere goats, and Inner Mongolia cashmere goats, including the SNP11 mutation site in our research results. This has no significant effect on the fiber diameter of Liaoning cashmere goats (P > 0.05) but a significant effect on the level of cashmere production and body weight after fleecing (P < 0.05). Fang et al. (2010) used PCR-RFLP technology to study the correlation between the KRTAP13.1 gene and the cashmere traits of Jiangnan goats, and the results showed that there was no significant difference in cashmere fineness between the different genotypes (P > 0.05). Similarly, Shanaz et al. (2020) used PCR-RFLP technology and found that the KRTAP13.1 gene polymorphism had no significant effect on the cashmere traits of Changthangi goats (P > 0.05). The above two loci are not present in our research results, which indicates the rich polymorphism of the KRTAP13.1 gene. Wu et al. (2018) used the PCR-SSCP method to find that the polymorphism of the KRTAP13.1 gene was significantly correlated with the cashmere yield and hair length of Tibetan cashmere goats.

We found 10 mutation sites (SNP6-SNP15) in the KRTAP13.1 gene, among which SNP8-SNP10 was synonymous and strongly linked (D? = 1), while SNP6, SNP7, and SNP11–SNP15 were missense mutations. These mutation sites were not annotated by the ensemble database (ENSCHIT00000017674.1), including SNP7, SNP8, and SNP11-15. It is worth noting that SNP7 has no heterozygous genotype in Jiangnan cashmere goat and Tibetan cashmere goat populations; the only two genotypes are CC and TT. And MAF of SNP11 in Jiangnan and Tibetan cashmere goats was less than 0.03; there was no correlation analysis between SNP11 and cashmere fiber diameter.

KRTAP15-1 and KRTAP24-1, like most KRTAP genes, do not contain introns. However, KRTAP13.1 and KRTAP27-1 have an intron. KRTAP27-1 is mainly expressed in skin tissues and is slightly or not expressed in longissimus dorsi muscle, heart, kidney, liver, lung, and spleen tissues (Zhao et al. 2020). This suggests that KRTAP27-1 may play a unique role in the skin. We found seven mutation sites on KRTAP27-1, of which SNP20 and SNP22 were synonymous mutations, while SNP16–SNP19 and SNP21 were missense mutations. Interestingly, we found a novel mutation site SNP18 in Tibetan cashmere goats. This suggests that the genetic resources of Tibetan cashmere goats have abundant genetic resources. Zhao et al. (2020) found SNP21 and SNP22 sites on Longdong cashmere goats, among which the SNP21 site significantly affected the fiber diameter of Longdong cashmere goats, which was consistent with the research results found for Jiangnan cashmere goats. This indicates that the KTAP27-1 polymorphism has a significant effect on the cashmere fiber diameter of cashmere goats.

We found six mutation sites (SNP23-SNP28) in the 200 bp fragment within the exon of KRTAP24-1 gene, of which SNP23, SNP24, and SNP28 are synonymous mutations and SNP25–SNP27 are missense mutations. Zhou et al. (2012) found four missense mutations in New Zealand Romney-cross sheep. Sun et al. (2016) found that two missense mutations in the KRTAP24-1 gene had a significant effect on the wool yield and hair length of Xinji fine wool sheep, but no significant effect on the fiber diameter. Comparing with our results, we found that these mutations in sheep do not exist in Jiangnan and Tibetan cashmere goats. Wang et al. (2019) found eight SNPs in the coding region of the KTAP24-1 gene of Longdong cashmere goats, which contained six mutation sites (SNP23-SNP28) that we found in Jiangnan and Tibetan cashmere goats. Similar to our results, SNP24 and SNP26 had a lower MAF in Longdong cashmere goats. Meanwhile, the combination of SNP23, SNP25, SNP27, and SNP28 in Longdong cashmere goats forms two unique band types using PCR-SSCP technology. This is consistent with our results of a strong linkage between SNP23, SNP25, SNP27, and SNP28.

According to linkage disequilibrium analysis, we selected a small number of SNPs from the set of SNP loci to represent the overall SNP. The purpose of this was to minimize the number of genotypes, thereby greatly reducing the number of SNPs used for further association studies, which will be helpful in accelerating the mining of functional genes (Wang et al. 2017c). Our results show that Tag SNP9, Tag SNP14, Tag SNP19, Tag SNP22, and Tag SNP28 significantly or extremely significantly affect the average fiber diameter of cashmere in the Jiangnan cashmere goat population. However, we did not find any mutation sites significantly related to cashmere fiber diameter in Tibetan cashmere goats. However, it is not difficult to see that the trend of cashmere fiber diameter of Tibetan cashmere goats and Jiangnan cashmere goats in different genotypes of significant mutation sites is roughly the same. This indicates that these Tag SNPs have a certain influence on the fiber diameter. On the other hand, this also explains the genetic differences between the two breeds.


The genetic polymorphisms related to fiber diameter in Jiangnan and Tibetan cashmere goats were tentatively explored in this study. These molecular markers can provide a theoretical scientific basis for the improvement of cashmere goat breeds. Our sequencing results were submitted to the NCBI public database (PRJNA738549). The accumulation of original genome data is helpful for research on the germplasm characteristics of local cashmere goats and the protection and utilization of resources.


This research was funded by Special Funds for Basic Scientific Research Operation of Public Welfare Scientific Research Institutes in Xinjiang Autonomous Region (2021), Xinjiang Autonomous Region Innovation Environment Construction Special Project (2020Q035, 2021D04008), Agriculture Research System of China (CARS-39), and Innovation Project of Shandong Academy of Agricultural Sciences (13200214443101).

Data availability

Data generated or analyzed during this study are available in the Sequence Read Archive (SRA) repository (PRJNA738549,

Author contributions

Data curation: QC and HJ

Formal analysis: ZB

Funding acquisition: FX, TK

Methodology: HX.

Resources: WC, QC, WY

Software: ZB

Supervision: TK

Validation: LJ, MJ

Visualization: HX

Writing—original draft: WC, FX

Writing—review and editing: WC

Supplementary material

Supplementary data are available with the article at

Supplementary Material 1 (DOCX / 132 KB).



Fu, X.F. 2021. Selection of regulatory factors related to cashmere fiber diameter trait from transcriptome and proteome profiles in Tibetan cashmere goats. Ph.D. Gansu Agricultural University, Lanzhou. Google Scholar


Gong, H., Zhou, H., Dyer, J.M., Plowman, J.E., and Hickford, J.G.H. 2011. Identification of the keratin-associated protein 13-3 (KAP13-3) gene in sheep. Open J. Genet. 1(3): 60–64. 13011Google Scholar


Gong, H., Zhou, H., Forrest, R.H.J., Li, S., Wang, J. Dyer, J.M., et al. 2016. Wool keratin-associated protein genes in sheep——a review. Genes, 7: 24. Google Scholar


Gong, H., Zhou, H., Mckenzie, G.W., Yu, Z., Clerens, S. Dyer, J.M., et al. 2012. An updated nomenclature for keratin-associated proteins (KAPs). Int J Biol Sci. 8: 258–264. 22298953. Google Scholar


Hailemariam, N., and Yadeta, W. 2020. Genetic polymorphism studies and their application on livestock breeding and health. Int. J. Genet. 10: 1–07. Google Scholar


Li, M.N., Liu, X., Wang, J.Q., Li, S.Q., and Luo, Y.Z. 2013. Genetic polymorphism of KAP13-1 gene of four cashmere goat populations in china. J. Northwest A & F Univ. (Nat. Sci. Ed.), 41: 7–13. jnwafu.2013.04.032Google Scholar


Li, W., Gong, H., Zhou, H., Wang, J., Liu, X. Li, S., et al. 2018. Variation in the ovine keratin-associated protein 15-1 gene affects wool yield. J. Agric. Sci. 156: 922–928. Google Scholar


Li, X., Su, R., Wan, W., Zhang, W., Jiang, H. Qiao, X., et al. 2017. Identification of selection signals by large-scale whole-genome resequencing of cashmere goats. Sci. Rep. 7: 15142. Google Scholar


Niu, L.L., and Hong-Bin, L.I. 2008. Molecular markers and their application in the study of domestic animal genetic diversity. Lett. Biotechnol. 02: 299–302 Google Scholar


Qin, C.K., Zhang, X.S., Zhang, F.Q., and Li, F.Z. 2016. Research on the growth and development rule of Nanjiang cashmere goat. Grass Feed. Livestock, 02: 12–16. Google Scholar


Shanaz, S., Ganai, A.S.T., Firdous Ahmad, S., Aalam Khan, A., Shabir, N. Shah, R., et al. 2020. Molecular characterization, polymorphism and association of Keratin-Associated Protein 13.1 (KAP13.1) gene in Changthangi goats. Biol Rhythm Res. 51: 915–924. 09291016.2018.1562122. Google Scholar


Strasser, B., Mlitz, V., Hermann, M., Tschachler, E., and Eckhart, L. 2015. Convergent evolution of cysteine-rich proteins in feathers and hair. BMC Evol. Biol. 15: 82. Google Scholar


Sun, F.L., Cao, Y., lu, C., Jin, H., and Zhang, L. 2014. Polymorphism of KAP13.1 gene and its effect on partial economic traits in Xiji fine wool sheep. Chin. J. Vet. Sci. 34: 2022–2026. Google Scholar


Sun, F.L., Jin, H.G., Lu, C., Cao, Y., Zhang, M.X., and Zhang, L.C. 2016. Correlation analysis of the polymorphisms of KAP24.1 gene and the wool traits in Xinji fine-wool sheep. Heilongjiang Anim. Sci. Vet. Med. 01: 70–73. Scholar


Tonin, C., Bianchetto, M., Vineis, C., and Festa Bianchet, M. 2002. Differentiating fine hairs from wild and domestic species: investigations of shatoosh, yangir, and cashmere fibers. Text. Res. J. 72: 701–705. Scholar


Wang, J., Che, L., Hickford, J.G.H., and Zhou, H. 2017a. Identification of the caprine keratin-associated protein 20-2 (KAP20-2) gene and its effect on cashmere traits. Genes, 8(11): 328. Scholar


Wang, J., Hao, Z., Zhou, H., Luo, Y., Hu, J. Liu, X., et al. 2018. A keratin-associated protein (KAP) gene that is associated with variation in cashmere goat fleece weight. Small Rumin. Res. 167: 104–109. j.smallrumres.2018.08.014. Google Scholar


Wang, J., Zhou, H., Luo, Y., Zhao, M., Gong, H. Hao, Z., et al. 2019. Variation in the caprine KAP24-1 gene affects cashmere fibre diameter. Animals, 9(1): 15. Scholar


Wang, J., Zhou, H., Zhu, J., Hu, J., Liu, X. Li, S., et al. 2017b. Identification of the ovine keratin-associated protein 15-1 gene (KRTAP15-1) and genetic variation in its coding sequence. Small Rumin. Res. 153: 131–136. Scholar


Wang, S., He, S., Yuan, F., and Zhu, X. 2017c. Tagging SNP-set selection with maximum information based on linkage disequilibrium structure in genome-wide association studies. Bioinformatics, 33: 2078–2081. Scholar


Ci, R.D.J. 2012. The study of cashmere performance in northwest Tibet. MA.Sc. Chinese Academy of Agricultural Sciences, Beijing. Google Scholar


Fang, Y., Liu, W.J., Zhang, F.Q., Shao, Y.G., and Yu, S.G. 2010. The polymorphism of a novel mutation of KAP13.1 gene and its associations with cashmere traits on Xinjiang local goat breed in China. Asian J. Anim. Vet. Adv. 5: 34–42. Scholar


Fu, X., Zhao, B., Tian, K., Wu, Y., Suo, L. Ba, G., et al. 2020. Integrated analysis of lncRNA and mRNA reveals novel insights into cashmere fineness in Tibetan cashmere goats. Peer J, 8: e10217. Google Scholar


Wu, Y.J., Suo, L.D., Ba, G., Ci, R.D.J, De, J., Song, T.Z., et al. 2018. Polymorphism of KAP13. 1 gene and relevance with production trait in Tibetan cashmere goat. Heilongjiang Anim. Sci. Vet. Med. 07: 79–82. Google Scholar


Yang, Y.-L., Zhou, R., and Li, K. 2017. Future livestock breeding: precision breeding based on multi-omics information and population personalization. J. Integr. Agric., 16: 2784–2791. Scholar


Yu, Y.S., Cao, Y., Piao, Q.L., Luo, X.T., and Jin, H.G. 2014. Screening of genes related to cashmere diameter in Liaoning cashmere goat. China Herbivore Sci. S1: 147–149. Scholar


Zhang, D., Song, T., Cheng, L., Du, X., Meng, F. Bu, G., et al. 2019. CDS cloning, exon polymorphism and expression analysis of KAP9.2 gene from Tibetan cashmere goat. Genomics Appl. Biol. 38(06): 2456–2465. Google Scholar


Zhao, M., Wang, X., Chen, H., Lan, X.Y., Guo, Y.K. Li, J.Y., et al. 2008. The PCR-SSCP and DNA sequencing methods detecting a large deletion mutation at KAP6.2 locus in the cashmere goat. Small Rumin. Res. 75(2–3): 243–246. pmid:32288216. Google Scholar


Zhao, M., Zhou, H., Luo, Y., Wang, J., Hu, J. Liu, X., et al. 2020. Variation in the caprine keratin-associated protein 27-1 gene is associated with cashmere fiber diameter. Genes, 11: 934. Google Scholar


Zheng, Y., Hui, T., Yue, C., Sun, J., Guo, D. Guo, S., et al. 2020a. Comprehensive analysis of circRNAs from cashmere goat skin by next generation RNA sequencing (RNA-seq). Sci. Rep. 10: 516. s41598-019-57404-9.pmid:31949277. Google Scholar


Zheng, Y., Sheng, S., Hui, T., Yue, C., Sun, J. Guo, D., et al. 2019. An integrated analysis of cashmere fineness lncRNAs in cashmere goats. Genes, 10: 266. Google Scholar


Zheng, Y., Wang, Y., Hui, T., Xu, Y., Zhang, Y. Wang, Z., et al. 2020b. 10×genomics single-cell sequencing reveals differential cell types in skin tissues of Liaoning cashmere goats and key genes related potentially to the fineness of cashmere fiber.Preprint version 1. 21203/ Google Scholar


Zhou, H., Gong, H., Yan, W., Luo, Y., and Hickford, J.G. 2012. Identification and sequence analysis of the Keratin-Associated Protein 24-1 (KAP24-1) gene homologue in sheep. Gene, 511: 62–65. 22995344. Google Scholar


Zhou, H.M., Allain, D., Li, J.Q., Zhang, W.G., and Yu, X.C. 2003. Effects of non-genetic factors on production traits of inner Mongolia cashmere goats in China. Small Rumin. Res. 47: 85–89. Scholar
Cuiling Wua, Chongkai Qin, Xuefeng Fu, Bingru Zhao, Yujiang Wu, Junmin He, Jingyi Mao, Jing Liu, Xixia Huang, and Kechuan Tian "Correlation analysis of four KRTAP gene polymorphisms and cashmere fiber diameters in two cashmere goat breeds," Canadian Journal of Animal Science 102(4), 561-570, (23 September 2022).
Received: 26 October 2021; Accepted: 9 May 2022; Published: 23 September 2022
cashmere fiber diameter
chèvre à cachemire de Jiangnan
chèvre à cachemire du Tibet
diamètre de la fibre de cachemire
Jiangnan cashmere goat
Back to Top