Using WGCNA (weighted gene co-expression network analysis) to identify the hub genes of skin hair follicle development in fetus stage of Inner Mongolia cashmere goat

Objective Mature hair follicles represent an important stage of hair follicle development, which determines the stability of hair follicle structure and its ability to enter the hair cycle. Here, we used weighted gene co-expression network analysis (WGCNA) to identify hub genes of mature skin and hair follicles in Inner Mongolian cashmere goats. Methods We used transcriptome sequencing data for the skin of Inner Mongolian cashmere goats from fetal days 45–135 days, and divided the co expressed genes into different modules by WGCNA. Characteristic values were used to screen out modules that were highly expressed in mature skin follicles. Module hub genes were then selected based on the correlation coefficients between the gene and module eigenvalue, gene connectivity, and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The results were confirmed by quantitative polymerase chain reaction (qPCR). Results Ten modules were successfully defined, of which one, with a total of 3166 genes, was selected as a specific module through sample and gene expression pattern analyses. A total of 584 candidate hub genes in the module were screened by the correlation coefficients between the genes and module eigenvalue and gene connectivity. Finally, GO/KEGG functional enrichment analyses detected WNT10A as a key gene in the development and maturation of skin hair follicles in fetal Inner Mongolian cashmere goats. qPCR showed that the expression trends of 13 genes from seven fetal skin samples were consistent with the sequencing results, indicating that the sequencing results were reliable.n