Deciphering the Transcriptomic Complexity of Yak Skin Across Different Ages and Body Sites.

Differences in skin and hair phenotypes between the scapular and ventral regions of yaks (Bos grunniens) are obvious and become more prominent with age. However, the genetic mechanism that causes differences in yak skin at different ages has not been reported. In this study, we investigated the transcriptomic profile of yak skin across different ages (0.5 years, 2.5 years, and 4.5 years) and body sites (scapular and ventral regions). Differential gene expression analysis was initially conducted to explore the transcriptomic differences in skin at different ages and different body sites. Subsequently, weighted gene co-expression network analysis (WGCNA) was employed to analyze the transcriptomic data comprehensively. The results showed that, among all comparison groups, the Y2.5_S vs. Y2.5_V group (regional comparison) exhibited the highest number of DEGs, with 491 genes (179 upregulated and 312 downregulated), followed by the Y2.5_V vs. Y0.5_V group (age comparison), which had 370 DEGs (103 upregulated and 267 downregulated). DEGs such as , , , and were found in all comparison groups of different ages, and multiple members of the HOX gene family including , , , , and were differentially expressed in comparison groups at different sites. Functional enrichment analysis showed that there were more obvious differences in immune function between different ages of skin and more obvious differences in endocrine function between different parts of skin. WGCNA revealed that genes related with immunity such as , , and were the core genes of the co-expression module associated with the scapula region, and multiple genes related to hair follicle development such as , , , , and were found to be the hub genes of the co-expression module associated with the ventral region. Overall, our study provides valuable insights into the transcriptomic complexity of yak skin across different ages and body sites. The differential gene expression patterns and co-expression network modules identified in this study lay the foundation for further research on skin biology and adaptation mechanisms in yaks.