Bioinformatics analysis and verification of m6A related genes based on the construction of keloid diagnostic model

Keloids are formed due to abnormal hyperplasia of the skin connective tissue. We explored the relationship between N6-methyladenosine (m6A)-related genes and keloids. The transcriptomic datasets (GSE44270 and GSE185309) of keloid and normal skin tissues samples were obtained from the Gene Expression Omnibus database. We constructed the m6A landscape and verified the corresponding genes using immunohistochemistry. We extracted hub genes for unsupervised clustering analysis using protein-protein interaction (PPI) network; gene ontology enrichment analysis was performed to determine the biological processes or functions affected by the differentially expressed genes (DEGs). We performed immune infiltration analysis to determine the relationship between keloids and the immune microenvironment using single-sample gene set enrichment analysis and CIBERSORT. Differential expression of several m6A genes was observed between the two groups; insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) was significantly upregulated in keloid patients. PPI analysis elucidated six genes with significant differences between the two keloid sample groups. Enrichment analysis revealed that the DEGs were mainly enriched in cell division, proliferation, and metabolism. Moreover, significant differences in immunity-related pathways were observed. Therefore, the results of this study will provide a reference for the elucidation of the pathogenesis and therapeutic targets of keloids.