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]

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis. [ABSTRACT FROM AUTHOR]