Abstract
BACKGROUND: Cashmere goats, as one of the important domesticated animal species, are known for their high-quality fiber. The growth of cashmere has seasonal variations caused by photoperiodic changes, but the molecular genetic mechanisms underlying this phenotype including the functional role of long non-coding RNAs (lncRNA) is still poorly understood. RESULTS: In this study, we analyzed the RNA-seq dataset of 39 Cashmere goat skin samples including all different growth stages and identified 1591 lncRNAs. These lncRNAs exhibited growth stage-specific expression patterns. Combining shortened light and hair follicle growth cycles, we found that 68% of differentially photo-responsive lncRNAs showed similar expression trends during transition phase I (early anagen to anagen phase). This suggests that the mechanism of light-controlled induction of hair follicles from early anagen to anagen is similar to that of transition phase I. According to weighted gene co-expression network analyses (WGCNA) analysis, it was found that two gene clusters and 10 hub lncRNAs participated in the transformation of hair follicle cycle, inducing hair follicles to enter the full growth phase in advance. These hub lncRNAs may regulate the development cycle of hair follicles through cis- or trans-regulation on clock genes, SLC superfamily genes, fibroblast growth factor genes. CONCLUSIONS: This study identified the key lncRNAs and target genes probably participating in the transformation of hair follicle cycle. This study will help further elucidate the role of lncRNAs in the hair follicle cycle and development.