Genome-wide analysis of dysregulated RNA-binding proteins and alternative splicing genes in keloid

Introduction: The pathogenesis of keloids remains unclear. Methods: In this study, we analyzed RNA-Seq data (GSE113619) of the local skin tissue of 8 keloid-prone individuals (KPI) and 6 healthy controls (HC) before and 42 days after trauma from the gene expression omnibus (GEO) database. The differential alternative splicing (AS) events associated with trauma healing between KPIs and HCs were identifified, and their functional differences were analyzed by gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathways. The co-expression relationship of differentially alternative splicing genes and differentially expressed RNA binding proteins (RBPs) was established subsequently. Results: A total of 674 differential AS events between the KD42 and the KD0 and 378 differential AS events between the HD42 and the HD0 were discovered. Notably, most of the differential genes related to keloids are enriched in actin, microtubule cells, and cortical actin cytoskeletal tissue pathway. We observed a signifificant association between AS genes (EPB41, TPM1, NF2, PARD3) and trauma healing in KPIs and HCs. We also found that the differential expression of healthy controls-specifific trauma healing-related RBPs (TKT, FDPS, SAMHD1) may affect the response of HCs to trauma healing by regulating the AS of downstream trauma healing-related genes such as DCN and DST. In contrast, KPIs also has specifific differential expression of trauma healing related RBPs (S100A9, HspB1, LIMA1, FBL), which may affect the healing response of KPIs to trauma by regulating the AS of downstream trauma healing-related genes such as FN1 and TPM1. Discussion: Our results were innovative in revealing early wound healing-related genes (EPB41, TPM1, NF2, PARD3) in KPI from the perspective of AS regulated by RBPs.