It is estimated that 50% of men and 25% of women worldwide suffer from hair loss, and therefore it is of great significance to investigate the molecular pathways driving hair follicle de novo morphogenesis. However, due to high cellular heterogeneity and the asynchronous development of hair follicles, our current understanding of the molecular mechanisms involved in follicle development remains limited. Methods: Single-cell suspensions from the dorsal skin of E13.5 (induction stage), E16.5 (organogenesis) fetal mice, and newborn mice (cytodifferentiation stage, postnatal day 0, P0) were prepared for unbiased single-cell RNA sequencing. To delineate the single-cell transcriptional landscape during hair follicle de novo morphogenesis, we performed t-distributed Stochastic Neighbor Embedding (tSNE), pseudotime cell trajectory inference, and regulon enrichment analysis to dissect cellular heterogeneity and reveal the molecular pathways underlying major cell type cell fate decisions. To validate our analysis, we further performed immunohistochemistry analysis of the key molecules involved during hair follicle morphogenesis. Meanwhile, intercellular communication between different cell populations was inferred based on a priori knowledge of ligand-receptor pairs. Results: Based on tSNE analysis, we identified 14 cell clusters from skin tissue and delineated their cellular identity from specific gene expression profiles. By using pseudotime ordering analysis, we successfully constructed the epithelium/dermal cell lineage differentiation trajectory. For dermal cell lineage, our analysis here recapitulated the dynamic gene expression profiles during dermal condensate (DC) cell fate commitment and delineated the heterogeneity of the different dermal papilla (DP) cell populations during in utero hair follicle development. For the epithelium cell lineage, our analysis revealed the dynamic gene expression profiles of the underappreciated matrix, interfollicular epidermis (IFE), hair shaft and inner root sheath (IRS) cell populations. Furthermore, single-cell regulatory network inference and clustering analysis revealed key regulons during cell fate decisions. Finally, intercellular communication analysis demonstrated that strong intercellular communication was involved during early hair follicle development. Conclusions: Our findings here provide a molecular landscape during hair follicle epithelium/dermal cell lineage fate decisions, and recapitulate the sequential activation of core regulatory transcriptional factors (TFs) in different cell populations during hair follicle morphogenesis. More importantly, our study here represents a valuable resource for understanding the molecular pathways involved during hair follicle de novo morphogenesis, which will have implications for future hair loss treatments.

The role of melatonin in promoting the yield of Cashmere goat wool has been demonstrated for decades though there remains a lack of knowledge regarding melatonin mediated hair follicle growth. Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are widely transcribed in the genome and play ubiquitous roles in regulating biological processes. However, the role of lncRNAs in regulating melatonin mediated hair follicle growth remains unclear. In this study, we established an in vitro Cashmere goat secondary hair follicle culture system, and demonstrated that 500 ng/L melatonin exposure promoted hair follicle fiber growth. Based on long intergenic RNA sequencing, we demonstrated that melatonin promoted hair follicle elongation via regulating genes involved in focal adhesion and extracellular matrix receptor pathways and further cis predicting of lncRNAs targeted genes indicated that melatonin mediated lncRNAs mainly targeted vascular smooth muscle contraction and signaling pathways regulating the pluripotency of stem cells. We proposed that melatonin exposure not only perturbed key signals secreted from hair follicle stem cells to regulate hair follicle development, but also mediated lncRNAs mainly targeted to pathways involved in the microvascular system and extracellular matrix, which constitute the highly orchestrated microenvironment for hair follicle stem cell. Taken together, our findings here provide a profound view of lncRNAs in regulating Cashmere goat hair follicle circadian rhythms and broaden our knowledge on melatonin mediated hair follicle morphological changes.

Characterization of the morphological structure during hair follicle development has been well documented, while the current understanding of the molecular mechanisms involved in follicle development remain limited. Here, using unbiased single-cell RNA sequencing, we analyzed 15,086 single cell transcriptome profiles from E13.5 and E16.5 fetal mice, and newborn mouse (postnatal day 0, P0) dorsal skin cells. Based on t-distributed Stochastic Neighbor Embedding (tSNE) clustering, we identified 14 cell clusters from skin cells and delineated their cell identity gene expression profiles. Pseudotime ordering analysis successfully constructed epithelium/dermal cell lineage differentiation trajectory and revealed sequential activation of key regulons involved during cell fate decisions. Along with this, intercellular communication between different cell populations were inferred based on a priori knowledge of ligand-receptor pairs. Together, our findings here provide a molecular landscape during hair follicle epithelium/dermal cell lineage fate decisions, and more importantly, recapitulate sequential activation of core regulatory transcriptional factors for different cell populations during hair follicle morphogenesis.

Polymorphisms in miRNA genes could potentially alter various biological processes by influencing the processing and (or) target selection of miRNAs. The rs14120863 (C > G) mutation, which we characterized in a Gushi-Anka F2 resource population, resides in the precursor region of miR-1666. Association analysis with chicken carcass and growth traits showed that the SNP was significantly associated with carcass weight, evisceration weight, breast muscle weight, leg muscle weight, and body weight at 8 weeks of age, as well as some body size indexes including shank girth, chest breadth, breast bone length, and body slanting length, in the Gushi-Anka F2 resource population. Quantitative RT-PCR results showed that miR-1666 expression levels in muscle tissues differed within various genotypes. Experiment in DF1 cells further confirmed that the SNP in miR-1666 could significantly alter mature miRNA production. Subsequently, using dual-luciferase report assay, we verified that miR-1666 could perform its function through targeting of the CBFB gene. In conclusion, the SNP in the precursor of miR-1666 could significantly reduce mature miR-1666 production. It may further affect the function of miR-1666 through the target gene CBFB, hence it is associated with chicken growth traits.

Cashmere, also known as soft gold, is produced from the secondary hair follicles in Cashmere goats and it’s therefore of significance to investigate the molecular profiles during Cashmere goat hair follicle development. However, our current understanding of the machinery underlying Cashmere goat hair follicle remains largely unexplored and researches regarding hair follicle development mainly used the mouse as a research model. To provides comprehensively understanding on the cellular heterogeneity and cell lineage cell fate decisions, we performed single-cell RNA sequencing on 19,705 single cells from induction (embryonic day 60), organogenesis (embryonic day 90) and cytodifferentiation (embryonic day 120) stages of fetus Cashmere goat dorsal skin. Unsupervised clustering analysis identified 16 cell clusters and their corresponding cell types were also unprecedentedly characterized. Based on the lineage inference, we revealed detailed molecular landscape along the dermal and epidermal cell lineage developmental pathways. Notably, by cross-species comparasion of single cell data with murine model, we revelaed conserved programs during dermal condensate fate commitment and the heterochrony development of hair follicle development between mouse and Cashmere goat were also discussed here. Our work here delineate unparalleled molecular profiles of different cell populations during Cashmere goat hair follicle morphogenesis and provide a valuable resource for identifying biomarkers during Cashmere goat hair follicle development.