Screening and identification of hub genes of scar physique via weighted gene co-expression network analysis

Abstract

Scar physique refers to the abnormal repair of skin injury in some people, which may easily lead to keloid or hypertrophic scar. However, the mechanism of scar physique is still unclear. GSE108110 was obtained from the gene expression omnibus database. Differently expression genes (DEGs) between normal skin tissue of non-scar physique individuals and normal skin tissue of scar physique individuals were screened by R package “limma”. Weighted gene co-expression network analysis was performed to find highly relevant gene modules. Functional annotation of DEGs was made. Protein-protein interaction network was constructed, and the identification and analysis of hub DEGs were performed, including identification of hub DEGs associated with scar diseases, MiRNA of hub DEGs prediction, and functional annotation of miRNA. A total of 1389 up-regulate DEGs and 1672 down-regulate DEGs were screened. weighted gene co-expression network analysis analysis showed that the dendrogram and heatmap were used to quantify module similarity by correlation. The associations between clinic traits and the modules were identified based on the correlation between module and scar physique. Eight common hub genes were obtained. The comparative toxicogenomics database shows common hub genes associated with scar tissue. Gene ontology and Kyoto encyclopedia of genes and genomes analysis were significantly enriched in “fibroblast growth factor receptor signaling pathway”, “epidermal growth factor receptor signaling pathway”, “G1/S transition of mitotic cell cycle”, protein polyubiquitination”, and others. The 8 hub genes might be involved in the development of scarring and used as early diagnosis, prevention and treatment of scar physique.