Apoptosis recognition receptors regulate skin tissue repair in mice

Apoptosis and clearance of apoptotic cells via efferocytosis are evolutionarily conserved processes that drive tissue repair. However, the mechanisms by which recognition and clearance of apoptotic cells regulate repair are not fully understood. Here, we use single-cell RNA sequencing to provide a map of the cellular dynamics during early inflammation in mouse skin wounds. We find that apoptotic pathways and efferocytosis receptors are elevated in fibroblasts and immune cells, including resident Lyve1+ macrophages, during inflammation. Interestingly, human diabetic foot wounds upregulate mRNAs for efferocytosis pathway genes and display altered efferocytosis signaling via the receptor Axl and its ligand Gas6. During early inflammation in mouse wounds, we detect upregulation of Axl in dendritic cells and fibroblasts via TLR3-independent mechanisms. Inhibition studies in vivo in mice reveal that Axl signaling is required for wound repair but is dispensable for efferocytosis. By contrast, inhibition of another efferocytosis receptor, Timd4, in mouse wounds decreases efferocytosis and abrogates wound repair. These data highlight the distinct mechanisms by which apoptotic cell detection coordinates tissue repair and provides potential therapeutic targets for chronic wounds in diabetic patients.

Lineage commitment of dermal fibroblast progenitors is controlled by Kdm6b-mediated chromatin demethylation

Dermal Fibroblast Progenitors (DFPs) differentiate into distinct fibroblast lineages during skin development. However, the epigenetic mechanisms that regulate DFP differentiation are not known. Our objective was to use multimodal single-cell approaches, epigenetic assays, and allografting techniques to define a DFP state and the mechanism that governs its differentiation potential. Our initial results indicated that the overall transcription profile of DFPs is repressed by H3K27me3 and has inaccessible chromatin at lineage-specific genes. Surprisingly, the repressive chromatin profile of DFPs renders them unable to reform the skin in allograft assays despite their multipotent potential. We hypothesized that chromatin derepression was modulated by the H3K27me3 demethylase, Kdm6b/Jmjd3. Dermal fibroblast-specific deletion of Kdm6b/Jmjd3 in mice resulted in adipocyte compartment ablation and inhibition of mature dermal papilla functions, confirmed by additional single-cell RNA-seq, ChIP-seq, and allografting assays. We conclude that DFPs are functionally derepressed during murine skin development by Kdm6b/Jmjd3. Our studies therefore reveal a multimodal understanding of how DFPs differentiate into distinct fibroblast lineages and provide a novel publicly available multiomics search tool.

Lineage Commitment of Dermal Fibroblast Progenitors is Mediated by Chromatin De-repression

Dermal Fibroblast Progenitors (DFPs) differentiate into distinct fibroblast lineages during skin development. However, the mechanisms that regulate lineage commitment of naive dermal progenitors to form niches around the hair follicle, dermis, and hypodermis, are unknown. In our study, we used multimodal single-cell approaches, epigenetic assays, and allografting techniques to define a DFP state and the mechanisms that govern its differentiation potential. Our results indicate that the overall chromatin profile of DFPs is repressed by H3K27me3 and has inaccessible chromatin at lineage specific genes. Surprisingly, the repressed chromatin profile of DFPs renders them unable to reform skin in allograft assays despite their multipotent potential. Distinct fibroblast lineages, such as the dermal papilla and adipocytes contained specific chromatin profiles that were de-repressed during late embryogenesis by the H3K27-me3 demethylase, Kdm6b/Jmjd3. Tissue-specific deletion of Kdm6b/Jmjd3 resulted in ablating the adipocyte compartment and inhibiting mature dermal papilla functions in single-cell-RNA-seq, ChIPseq, and allografting assays. Altogether our studies reveal a mechanistic multimodal understanding of how DFPs differentiate into distinct fibroblast lineages, and we provide a novel multiomic search-tool within skinregeneration.org.

Parallel single cell multi-omics analysis of neonatal skin reveals transitional fibroblast states that restricts differentiation into distinct fates

One of the keys to achieving skin regeneration lies within understanding the heterogeneity of neonatal fibroblasts, which support skin regeneration. However, the molecular underpinnings regulating the cellular states and fates of these cells are not fully understood. To investigate this, we performed a parallel multi-omics analysis by processing neonatal murine skin for single-cell ATAC-sequencing (scATAC-seq) and single-cell RNA-sequencing (scRNA-seq) separately. Our approach revealed that fibroblast clusters could be sorted into papillary and reticular lineages based on transcriptome profiling, as previously published. However, scATAC-seq analysis of neonatal fibroblast lineage markers, such as, Dpp4/CD26, Corin, and Dlk1 along with markers of myofibroblasts, revealed accessible chromatin in all fibroblast populations despite their lineage-specific transcriptome profiles. These results suggests that accessible chromatin does not always translate to gene expression and that many fibroblast lineage markers reflect a fibroblast state, which includes neonatal papillary, reticular, and myofibroblasts. This analysis also provides a possible explanation as to why these marker genes can be promiscuously expressed in different fibroblast populations under different conditions. Our scATAC-seq analysis also revealed that the functional lineage restriction between dermal papilla and adipocyte fates are regulated by distinct chromatin landscapes. Finally, we have developed a webtool for our multi-omics analysis: https://skinregeneration.org/scatacseq-and-scrnaseq-data-from-thompson-et-al-2021-2/.

Single-cell transcriptomic analysis of small and large wounds reveals the distinct spatial organization of regenerative fibroblasts

Wound-induced hair follicle neogenesis (WIHN) has been an important model to study hair follicle regeneration during wound repair. However, the cellular and molecular components of the dermis that make large wounds more regenerative are not fully understood. Here, we compare and contrast recently published scRNA-seq data of small scarring wounds to wounds that regenerate in hope to elucidate the role of fibroblasts lineages in WIHN. Our analysis revealed an over-representation of the newly identified upper wound fibroblasts in regenerative wound conditions, which express the retinoic acid binding protein Crabp1. This regenerative cell type shares a similar gene signature to the murine papillary fibroblast lineage, which are necessary to support hair follicle morphogenesis and homeostasis. RNA velocity analysis comparing scarring and regenerating wounds revealed the divergent trajectories towards upper and lower wound fibroblasts and that the upper populations were closely associated with the specialized dermal papilla. We also provide analyses and explanation reconciling the inconsistency between the histological lineage tracing and the scRNA-seq data from recent reports investigating large wounds. Finally, we performed a computational test to map the spatial location of upper wound fibroblasts in large wounds which revealed that upper peripheral fibroblasts might harbour equivalent regenerative competence as those in the centre. Overall, our scRNA-seq reanalysis combining multiple samples suggests that upper wound fibroblasts are required for hair follicle regeneration and that papillary fibroblasts may migrate from the wound periphery to the centre during wound re-epithelialization. Moreover, data from this publication are made available on our searchable web resource: https://skinregeneration.org/.

Lef1 expression in fibroblasts maintains developmental potential in adult skin to regenerate wounds

Scars are a serious health concern for burn victims and individuals with skin conditions associated with wound healing. Here, we identify regenerative factors in neonatal murine skin that transforms adult skin to regenerate instead of only repairing wounds with a scar, without perturbing development and homeostasis. Using scRNA-seq to probe unsorted cells from regenerating, scarring, homeostatic, and developing skin, we identified neonatal papillary fibroblasts that form a transient regenerative cell type that promotes healthy skin regeneration in young skin. These fibroblasts are defined by the expression of a canonical Wnt transcription factor Lef1 and using gain- and loss of function genetic mouse models, we demonstrate that Lef1 expression in fibroblasts primes the adult skin macroenvironment to enhance skin repair, including regeneration of hair follicles with arrector pili muscles in healed wounds. Finally, we share our genomic data in an interactive, searchable companion website (https://skinregeneration.org/). Together, these data and resources provide a platform to leverage the regenerative abilities of neonatal skin to develop clinically tractable solutions that promote the regeneration of adult tissue.