Transcriptome profiling reveals transcriptional and alternative splicing regulation in the early embryonic development of hair follicles in the cashmere goat

The undercoat fiber of the cashmere goat, from the secondary hair follicle (HF), possesses commercial value. However, very few studies have focused on the molecular details of primary and secondary HF initiation and development in goat embryos. In this study, skin samples at embryonic day 45, 55, and 65 (E45, E55, and E65) were collected and prepared for RNA sequencing (RNA-seq). We found that the HF probably initiated from E55 to E65 by analyzing the functional pathways of differentially expressed genes (DEGs). Most key genes in canonical signaling pathways, including WNT, TGF-beta, FGF, Hedgehog, NOTCH, and other factors showed clear expression changes from E55 to E65. We, for the first time, explored alternative splicing (AS) alterations, which showed distinct patterns among these three stages. Functional pathways of AS-regulated genes showed connections to HF development. By comparing the published RNA-seq samples from the E60, E120, and newborn (NB) stages, we found the majority of WNT/beta-catenin signaling genes were important in the initiation of HF development, while other factors including FOXN1, GATA3, and DLX3 may have a consistent influence on HF development. Our investigation supported the time points of embryonic HF initiation and identified genes that have potential functions of embryonic HF initiation and development. We further explored the potential regulatory roles of AS in HF initiation, which extended our knowledge about the molecular mechanisms of HF development.