Fig 1: 3′ RNA-seq analysis identified keloid specific differentially expressed genes. (a) Scheme of dissections of the TAGLN-positive area, and workflow of the analysis performed in this study. (b) 3′ RNA-seq analysis for immature scars (n = 2) and keloids (n = 3). Volcano plots showing the differentially 112 upregulated genes and 108 downregulated genes in keloids (fold change > 2, FDR P-value ≤ 0.1). (c) Venn diagram showing the eight shared differentially expressed genes among the aforementioned 220 genes identified though the current 3′ RNA-seq analysis and 414 genes detected in the GSE113619 dataset (fold change > 2, FDR P-value ≤ 0.1). (d) Relative mRNA expressions of the eight differentially expressed genes in immature scars (n = 10) and keloids (n = 12) compared with GAPDH. Error bars show the mean ± SD. P-value was determined using the two-tailed t-test. *P < 0.05, **P < 0.01, ***P < 0.001. ELN elastin, FDR false discovery rate, GAPDH glyceraldehyde-3-phosphate dehydrogenase, KANK4 KN motif and ankyrin repeat domains 4, KRT6A keratin 6A, OPCML opioid binding protein/cell adhesion molecule like, PTPRD protein tyrosine phosphatase receptor type D, qPCR quantitative polymerase chain reaction, RNA-seq RNA sequencing, S100A7 S100 calcium binding protein A7, TNXB tenascin XB, S100A8 S100 calcium binding protein A8, SD standard deviation, TAGLN transgelin.