Investigation on the growth factor regulatory network of dermal fibroblasts in mouse full-thickness skin defect wounds based on single-cell RNA sequencing

Objective To explore the heterogeneity and growth factor regulatory network of dermal fibroblasts (dFbs) in mouse full-thickness skin defect wounds based on single-cell RNA sequencing.

Methods The experimental research methods were adopted. The normal skin tissue from 5 healthy 8-week-old male C57BL/6 mice (the same mouse age, sex, and strain below) was harvested, and the wound tissue of another 5 mice with full-thickness skin defect on the back was harvested on post injury day (PID) 7. The cell suspension was obtained by digesting the tissue with collagenase D and DNase Ⅰ, sequencing library was constructed using 10x Genomics platform, and single-cell RNA sequencing was performed by Illumina Novaseq6000 sequencer. The gene expression matrices of cells in the two kinds of tissue were obtained by analysis of Seurat 3.0 program of software R4.1.1, and two-dimensional tSNE plots classified by cell group, cell source, and gene labeling of major cells in skin were used for visual display. According to the existing literature and the CellMarker database searching, the expression of marker genes in the gene expression matrices of cells in the two kinds of tissue was analyzed, and each cell group was numbered and defined. The gene expression matrices and cell clustering information were introduced into CellChat 1.1.3 program of software R4.1.1 to analyze the intercellular communication in the two kinds of tissue and the intercellular communication involving vascular endothelial growth factor (VEGF), platelet-derived growth factor (PDGF), epidermal growth factor (EGF), and fibroblast growth factor (FGF) signal pathways in the wound tissue, the relative contribution of each pair of FGF subtypes and FGF receptor (FGFR) subtypes (hereinafter referred to as FGF ligand receptor pairs) to FGF signal network in the two kinds of tissue, and the intercellular communication in the signal pathway of FGF ligand receptor pairs with the top 2 relative contributions in the two kinds of tissue. The normal skin tissue from one healthy mouse was harvested, and the wound tissue of one mouse with full-thickness skin defect on the back was harvested on PID 7. The multiple immunofluorescence staining was performed to detect the expression and distribution of FGF7 protein and its co-localized expression with dipeptidyl peptidase 4 (DPP4), stem cell antigen 1 (SCA1), smooth muscle actin (SMA), and PDGF receptor α (PDGFRα) protein.

Results Both the normal skin tissue of healthy mice and the wound tissue of full-thickness skin defected mice on PID 7 contained 25 cell groups, but the numbers of cells in each cell group between the two kinds of tissue were different. Genes PDGFRα, platelet endothelial cell adhesion molecule 1, lymphatic endothelial hyaluronic acid receptor 1, receptor protein tyrosine phosphatase C, keratin 10, and keratin 79 all had distinct distributions on two-dimensional tSNE plots, indicating specific cell groups respectively. The 25 cell groups were numbered by C0-C24 and divided into 9 dFb subgroups and 16 non-dFb groups. dFb subgroups included C0 as interstitial progenitor cells, C5 as adipose precursor cells, and C13 as contractile muscle cells related fibroblasts, etc. Non-dFb group included C3 as neutrophils, C8 as T cells, and C18 as erythrocytes, etc. Compared with that of the normal skin tissue of healthy mice, the intercellular communication in the wound tissue of full-thickness skin defected mice on PID 7 was more and denser, and the top 3 cell groups in intercellular communication intensity were dFb subgroups C0, C1, and C2, of which all communicated with other cell groups in the wound tissue. In the wound tissue of full-thickness skin defected mice on PID 7, VEGF signals were mainly sent by the dFb subgroup C0 and received by vascular related cell groups C19 and C21, PDGF signals were mainly sent by peripheral cells C14 and received by multiple dFb subgroups, EGF signals were mainly sent by keratinocyte subgroups C9 and C11 and received by the dFb subgroup C0, and the main sender and receiver of FGF signals were the dFb subgroup C6. In the relative contribution rank of FGF ligand receptor pairs to FGF signal network in the normal skin tissue of healthy mice and the wound tissue of full-thickness skin defected mice on PID 7, FGF7-FGFR1 was the top 1, and FGF7-FGFR2 or FGF10-FGFR1 was in the second place, respectively; compared with those in the normal skin tissue, there was more intercellular communication in FGF7-FGFR1 signal pathway, while the intercellular communication in FGF7-FGFR2 and FGF10-FGFR1 signal pathways decreased slightly or did not change significantly in the wound tissue; the intercellular communication in FGF7-FGFR1 signal pathway in the wound tissue was stronger than that in FGF7-FGFR2 or FGF10-FGFR1 signal pathway; in the two kinds of tissue, FGF7 signal was mainly sent by dFb subgroups C0, C1, and C2, and received by dFb subgroups C6 and C7. Compared with that in the normal skin tissue of healthy mouse, the expression of FGF7 protein was higher in the wound tissue of full-thickness skin defected mouse on PID 7; in the normal skin tissue, FGF7 protein was mainly expressed in the skin interstitium and also expressed in the white adipose tissue near the dermis layer; in the two kinds of tissue, FGF7 protein was co-localized with DPP4 and SCA1 proteins and expressed in the skin interstitium, co-localized with PDGFRα protein and expressed in dFbs, but was not co-localized with SMA protein, with more co-localized expression of FGF7 in the wound tissue than that in the normal skin tissue.

Conclusions In the process of wound healing of mouse full-thickness skin defect wound, dFbs are highly heterogeneous, act as potential major secretory or receiving cell populations of a variety of growth factors, and have a close and complex relationship with the growth factor signal pathways. FGF7-FGFR1 signal pathway is the main FGF signal pathway in the process of wound healing, which targets and regulates multiple dFb subgroups.