Identification of N6-Methyladenosine-Related Factors and the Prediction of the Regulatory Mechanism of Hair Follicle Development in Rex and Hycole Rabbits

Abstract

Hair follicle development directly affects the development of the rabbit fur industry. The growth and development of a hair follicle is modified and regulated by many genes and mechanisms. M6A is an important RNA modification. However, there are few studies on the effects of the regulation of m6A on hair follicle growth and development. In this study, hematoxylin–eosin (HE) staining was used to explore the difference in hair follicle development between Rex rabbits and Hycole rabbits, and we performed m6A sequencing to identify the key genes with m6A modification in hair follicle growth. The results showed that the hair length, coarse hair percentage, primary hair follicle ratio, and skin thickness of Hycole rabbits were significantly higher than those of Rex rabbits. However, the proportion of secondary hair follicles in Hycole rabbits was significantly lower than that in Rex rabbits. In addition, we found five differential methylases, 20 differential genes, and 24 differential signaling pathways related to hair growth and development. The results of the Sankey diagram showed that 12 genes were related to 13 signal pathways. Finally, we found that five methylases regulated the development of hair follicles through differential genes/signal pathways. These findings laid a molecular foundation for the function of m6A modification in hair development.