Screening the key genes of hair follicle growth cycle in Inner Mongolian Cashmere goat based on RNA sequencing.

Inner Mongolian Cashmere goat is an excellent local breed selected for the dual-purpose of cashmere and meat. There are three lines of Inner Mongolian Cashmere goat: Erlangshan, Alashan and Aerbasi. Cashmere is a kind of precious textile raw material with a high price. Cashmere is derived from secondary hair follicle (SHF), while hair is derived from primary hair follicle (PHF). The growth cycle of SHF of cashmere goat is 1 year, and it can be divided into three different stages: anagen, catagen and telogen. In this study, we tried to find some important influence factors of SHF growth cycle in skin tissue from Inner Mongolian Cashmere goats by RNA sequencing (RNA-Seq). Three female Aerbasi Inner Mongolian Cashmere goats (2 years old) were used as experimental samples in this study. Skin samples were collected in September (anagen), December (catagen) and March (telogen) at dorsal side from cashmere goats. Results showed that over 511 396 044 raw reads and 487 729 890 clean reads were obtained from sequence data. In total, 51 different expression genes (DEGs) including 29 downregulated genes and 22 upregulated genes were enriched in anagen-catagen comparing group. The 443 DEGs contained 117 downregulated genes and 326 upregulated genes that were enriched in catagen-telogen comparing group. In telogen-anagen comparing group, 779 DEGs were enriched including 582 downregulated genes and 197 upregulated genes. The result of gene ontology (GO) annotation showed that DEGs are in different growth cycle periods, and enriched GO items are mostly related to the transformation of cell and protein. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment result indicated that metabolic process has a great impact on SHF growth cycle. Based on the results of a comprehensive analysis of differentially expressed genes, GO enrichment and KEGG enrichment, we found that