Abstract:BackgroundKeloid is a benign proliferative disease characterized by excessive deposition of extracellular matrix collagen during skin wound healing. The mechanisms of keloid formation have not been fully elucidated, and the current treatment methods are not effective for all keloid patients. Therefore, there is an urgent need to find more effective therapies, and our research focused on identifying characteristic molecular signatures of keloid to explore potential therapeutic targets.MethodsGene expression profiles of keloid and control group samples were retrieved from the GEO database. Taking the GSE113619 dataset as the training set, the dataset collected skin tissues from non-lesion sites of healthy and keloid-prone individuals, denoted as Day0. The second sampling was performed 42 days later at the original sampling site of control and keloid groups, denoted as Day42.The 'limma' package and Venn diagram identified differentially expressed genes (DEGs) specific to keloid day42 versus day0 samples. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome pathway functional enrichment, and annotation of the characteristic genes were conducted on the Metascape website. Ingenuity canonical pathways, disease & function enrichment analysis and gene interaction network were performed and predicted in Ingenuity Pathway Analysis (IPA) software. Key module genes related to keloid were filtered out by Weighted Gene Co-expression Network Analysis (WGCNA). We utilized the Least Absolute Shrinkage and Selection Operator (LASSO) algorithm to screen the characteristic genes in keloid by the ‘glmnet’ package. The area under the curve (AUC) of receiver operating characteristic (ROC) was utilized to determine the effectiveness of potential signatures in discriminating keloid samples from normal samples and performed by using the ‘pROC’ package. The enrich scores of 24 immune cells in each sample were calculated by the single-sample gene set enrichment analysis (ssGSEA) algorithm, and then the Gene Set Variation Analysis (GSVA) was performed. Finally, RNA from 4 normal and 6 keloid samples was extracted, and RT-qPCR and Western Blot validated the expression of characteristic genes.ResultsA total of 640 DEGs specific to keloid day42 versus day0 samples were detected. 69 key module genes were uncovered and implicated in ‘NCAM signaling for neurite out-growth’, ‘oncogenic MAPK signaling’, ‘transmission across chemical synapses’ pathways, and the mitotic cell cycle-related processes. Five characteristic genes (MTUS1, UNC5C, CEP57, NAA35, and HOXD3) of keloid were identified by LASSO, and among which UNC5C and HOXD3 were validated by ROC plot in external dataset, RT-qPCR and Western Blot in validation samples. The result of ssGSEA indicated that the infiltration of neutrophils showed a relatively higher abundance and natural killer cells with relatively low enrichment in the keloid group compared to the control group. UNC5C was correlated with more immune cells compared with other characteristic genes.ConclusionIn this study, characteristic genes associated with keloid were identified by bioinformatic approaches and verified in clinical validation samples, providing potential targets for the diagnosis and treatment of keloid.

Comprehensive Analysis of Circular RNA Expression in Cerna Networks and Identification of the Effects of Hsa_circ_0006867 in Keloid Dermal Fibroblasts

Comprehensive Analysis of Differentially Expressed Circular RNAs in Keloid Dermal Tissues

Gene Profiling of Circular RNAs in Keloid-prone Individuals and ceRNA Network Construction During Wound Healing

Circular RNAs: Emerging Players in the Pathogenesis of Keloid

Circular RNA CircCOL5A1 Sponges the MiR-7-5p/Epac1 Axis to Promote the Progression of Keloids Through Regulating PI3K/Akt Signaling Pathway

Detection and Analysis of Long Noncoding RNA Expression Profiles Related to Epithelial-Mesenchymal Transition in Keloids

Circular RNA Microarray Expression Profile and Potential Function of Circ0005875 in Clear Cell Renal Cell Carcinoma.

Circular RNA Hsa_circ_0084443 is Upregulated in Diabetic Foot Ulcer and Modulates Keratinocyte Migration and Proliferation.

MicroRNAs Associated with Keloids Identified by Microarray Analysis and In Vitro Experiments

Comparative Study of Microrna Profiling in Keloid Fibroblast and Annotation of Differential Expressed Micrornas.

The biomarkers associated with epithelial-mesenchymal transition in human keloids

Integrated bioinformatics analysis of core regulatory elements involved in keloid formation

Identification and validation of HOXD3 and UNC5C as molecular signatures in keloid based on weighted gene co-expression network analysis

Keloid Patient Plasma-Derived Exosomal Hsa_circ_0020792 Promotes Normal Skin Fibroblasts Proliferation, Migration, and Fibrogenesis Via Modulating Mir-193A-5p and Activating TGF-β1/Smad2/3 Signaling

High‑throughput sequencing reveals differentially expressed lncRNAs and circRNAs, and their associated functional network, in human hypertrophic scars.

Silencing circular RNAPTPN12 promoted the growth of keloid fibroblasts by activating Wnt signaling pathway via targeting microRNA-21-5p

Identification of potential therapeutic target SPP1 and related RNA regulatory pathway in keloid based on bioinformatics analysis

Identification of Crucial Noncoding RNAs and Mrnas in Hypertrophic Scars Via RNA Sequencing.

Combined Analyses of GWAS, RNA-Sequence and Hi-C—A Novel Approach to Dissect Keloid Genetic Mechanism

The study on circRNA profiling uncovers the regulatory function of the hsa_circ_0059665/miR-602 pathway in breast cancer

Long Non-Coding RNA Expression Profiling in the Lesional Tissue and Derived Fibroblasts of Keloid