Search

• This is the first systematic study of the keratin family in yak. 53 keratins are identified in yak. • The results show that yak keratins have the same conserved domain and sequence as other species keratins. • The study demonstrates the expression pattern of keratins genes throughout the hair cycle and identifies five different expression patterns. • Nine hub keratin genes are identified by WGCNA method in the whole hair cycle. As the largest subgroup of intermediate filament proteins, keratins are divided into two types of subfamily. Currently, the molecular mechanism of keratins in several animals has been reported but is limited in yak. Here, 53 different kinds of keratins were identified in the yak genome, including 23 type I and 30 type II keratins. Bioinformatics analysis in this study revealed that multiple phosphorylation sites were identified among all the family members. And the subcellular localization of these proteins was predicted to be in the nucleus, cytoskeleton, and cytoplasm. All keratin family proteins were unstable and the scores of instability coefficient were higher than 40. Phylogenetic analysis showed that high consistency results of the sequence conservation and grouping were found in the genomes of yak, sheep, cattle, mouse, rat, and human. Based on the expression patterns obtained from the transcriptome data, keratin genes (KRTs) were grouped into five clusters, and results also showed that KRTs were highly activated in skin tissues during the hair cycle in yak. Among the five clusters, Cluster II contained the most KRTs , which was the main expression pattern of the yak hair follicle cycle, followed by Cluster III. These results indicated the transition period from telogen to anagen and catagen to telogen were highly dynamic in yak. Gene expression correlation analysis showed that KRTs exhibited a strong correlation (mainly positive correlation) throughout the hair follicle development cycle. And the identification of hub KRTs in specific modules related to hair follicle development in this study was performed using the Weight Gene Co-Expression Network Analysis (WGCNA). Specific modules that include KRTs were darkgreen (KRT40), darkgrey (KRT5), turquoise (KRT1 , KRT2 , KRT10), bisque4 (KRT4), thistle2 (KRT9 , KRT39), and yellowgreen (KRT24). The interaction network showed that these genes were found to be related to the regulation of cell cycle, melanogenesis, hair follicle development, keratinocyte proliferation. Our study provides theoretical support for the study of the evolutionary relationship and molecular mechanism of keratin family in B. grunnien. [ABSTRACT FROM AUTHOR]