Parallel single cell multi-omics analysis of neonatal skin reveals transitional fibroblast states that restrict 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 adipocytes 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/.

Trunk skin from 4 wild-type postnatal day 0 mouse littermates were used to generate single-cell suspensions of skin cells. These single-cell suspensions were pooled and parallel scRNA-seq and scATAC-seq experiments were conducted on the same pool of cells. The resulting sequencing data was aligned to the mm10 genome.