Single-Cell and Bulk Transcriptome Data Integration Reveals Dysfunctional Cell Types and Aberrantly Expressed Genes in Hypertrophic Scar

Abstract

Hypertrophic scar (HS) is a common skin disorder characterized by excessive extracellular matrix (ECM) deposition. However, it is still unclear how the cellular composition, cell-cell communications, and crucial transcriptionally regulatory network were changed in HS. In the present study, we found that FB-1, which was identified a major type of fibroblast and had the characteristics of myofibroblast, was significantly expanded in HS by integrative analysis of the single-cell and bulk RNA sequencing (RNA-seq) data. Moreover, the proportion of KC-2, which might be a differentiated type of keratinocyte (KC), was reduced in HS. To decipher the intercellular signaling, we conducted the cell-cell communication analysis between the cell types, and found the autocrine signaling of HB-1 through COL1A1/2-CD44 and CD99-CD99 and the intercellular contacts between FB-1/FB-5 and KC-2 through COL1A1/COL1A2/COL6A1/COL6A2-SDC4. Almost all the ligands and receptors involved in the autocrine signaling of HB-1 were upregulated in HS by both scRNA-seq and bulk RNA-seq data. In contrast, the receptor of KC-2, SDC4, which could bind to multiple ligands, was downregulated in HS, suggesting that the reduced proportion of KC-2 and apoptotic phenotype of KC-2 might be associated with the downregulation of SDC4. Furthermore, we also investigated the transcriptionally regulatory network involved in HS formation. The integrative analysis of the scRNA-seq and bulk RNA-seq data identified CREB3L1 and TWIST2 as the critical TFs involved in the myofibroblast of HS. In summary, the integrative analysis of the single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data greatly improved our understanding of the biological characteristics during the HS formation.