Construction of transcriptome atlas of white yak hair follicle during anagen and catagen using single-cell RNA sequencing

Abstract Background As the direct organ of villus, hair follicles have obvious seasonal cycles. The hair follicle cycle is orchestrated by multiple cell types that together direct cell renewal and differentiation. But the regulation property of hair follicle cells from anagen to catagen in yak is still unknown. Results In this study, single-cell RNA sequencing was performed on 24,124 single cells of the scapular skin from white yak. Based on tSNE cluster analysis, the cell types of IFE-DC, epidermal cell lines, fibroblasts, keratinocytes, IRS, DS, INFU, and other cells in yak hair follicles during anagen and catagen were successfully identified, and the gene expression profiles were described. The GO enrichment analysis indicated the different cells characteristic genes to be mainly enriched in the epidermal development, epithelial cell differentiation and wound healing pathways. The pseudotime trajectory analysis described the differentiation trajectory of the epidermal lineage and dermal lineage of the hair follicle during anagen and catagen. Moreover, the dynamic changes of the genes like LHX2, KRT25, and KRT71 were found to be highly expressed in HS and IRS, but not in the IFE-DC, INFU, and keratinocyte during differentiation. Conclusions Our results analyzed the time-varying process of gene expression in the dermal cell lineage and epidermal cell lineage of hair follicles during anagen and catagen during fate differentiation was expounded at the single cell level, revealing the law of fate specialization of different types of cells. In addition, based on the enrichment analysis, the transcriptional regulatory factors involved in the different cell fates were also revealed. These results will help to enhance our understanding of yak hair follicle cycle and promote the development and utilization of yak villus.