Abstract
Hypertrophic scar (HS) is a chronic inflammatory skin disorder characterized by excessive deposition of extracellular matrix, and the mechanisms underlying their formation remain poorly understood. We analysed scRNA-seq data from samples of normal skin and HS. Using the hdWGCNA method, key gene modules of fibroblasts in HS were identified. Non-negative matrix factorization was employed to perform subtype analysis of HS patients using these gene modules. Multiple machine learning algorithms were applied to screen and validate accurate gene signatures for identifying and predicting HS, and a convolutional neural network (CNN) based on deep learning was established and validated. Quantitative reverse transcription-polymerase chain reaction and western blotting were performed to measure mRNA and protein expression. Immunofluorescence was used for gene localization analysis, and biological features were assessed through CCK8 and wound healing assay. Single-cell sequencing revealed distinct subpopulations of fibroblasts in HS. HdWGCNA identified key gene characteristics of this population, and pseudotime analysis was conducted to investigate gene variation during fibroblast differentiation. By employing various machine learning algorithms, the gene range was narrowed down to three key genes. A CNN was trained using the expression of these key genes and immune cell infiltration, enabling diagnosis and prediction of HS. Functional experiments demonstrated that THBS2 is associated with fibroblast proliferation and migration in HS and affects the formation and development of HS through the TGFβ1/P-Smad2/3 pathway. Our study identifies unique fibroblast subpopulations closely associated with HS and provides biomarkers for the diagnosis and treatment of HS.