Application of single-cell RNA sequencing on human skin: Technical evolution and challenges

Single-cell transcriptomics reveals aberrant skin-resident cell populations and identifies fibroblasts as a determinant in rosacea

Rosacea is a chronic inflammatory skin disorder, whose underlying cellular and molecular mechanisms remain obscure. Here, we generate a single-cell atlas of facial skin from female rosacea patients and healthy individuals. Among keratinocytes, a subpopulation characterized by IFNγ-mediated barrier function damage is found to be unique to rosacea lesions. Blocking IFNγ signaling alleviates rosacea-like phenotypes and skin barrier damage in mice. The papulopustular rosacea is featured by expansion of pro-inflammatory fibroblasts, Schwann, endothelial and macrophage/dendritic cells. The frequencies of type 1/17 and skin-resident memory T cells are increased, and vascular mural cells are characterized by activation of inflammatory pathways and impaired muscle contraction function in rosacea. Most importantly, fibroblasts are identified as the leading cell type producing pro-inflammatory and vasodilative signals in rosacea. Depletion of fibroblasts or knockdown of PTGDS, a gene specifically upregulated in fibroblasts, blocks rosacea development in mice. Our study provides a comprehensive understanding of the aberrant alterations of skin-resident cell populations and identifies fibroblasts as a key determinant in rosacea development.

Microfluidic device for the high-throughput and selective encapsulation of single target cells

Droplet-based microfluidic technologies for encapsulating single cells have rapidly evolved into powerful tools for single-cell analysis. In conventional passive single-cell encapsulation techniques, because cells arrive randomly at the droplet generation section, to encapsulate only a single cell with high precision, the average number of cells per droplet has to be decreased by reducing the average frequency at which cells arrive relative to the droplet generation rate. Therefore, the encapsulation efficiency for a given droplet generation rate is very low. Additionally, cell sorting operations are required prior to the encapsulation of target cells for specific cell type analysis. To address these challenges, we developed a cell encapsulation technology with a cell sorting function using a microfluidic chip. The microfluidic chip is equipped with an optical detection section to detect the optical information of cells and a sorting section to encapsulate cells into droplets by controlling a piezo element, enabling active encapsulation of only the single target cells. For a particle population including both targeted and non-targeted particles arriving at an average frequency of up to 6000 particles per s, with an average number of particles per droplet of 0.45, our device maintained a high purity above 97.9% for the single-target-particle droplets and achieved an outstanding throughput, encapsulating up to 2900 single target particles per s. The proposed encapsulation technology surpasses the encapsulation efficiency of conventional techniques, provides high efficiency and flexibility for single-cell research, and shows excellent potential for various applications in single-cell analysis.

Optimization of Tissue Digestion Methods for Characterization of Photoaged Skin by Single Cell RNA Sequencing Reveals Preferential Enrichment of T Cell Subsets

Healthy human skin tissue is often used as a control for comparison to diseased skin in patients with skin pathologies, including skin cancers or other inflammatory conditions such as atopic dermatitis or psoriasis. Although non-affected skin from these patients is a more appropriate choice for comparison, there is a paucity of studies examining such tissue. This lack is exacerbated by the difficulty of processing skin tissue for experimental analysis. In addition, choosing a processing protocol for skin tissue which preserves cell viability and identity while sufficiently dissociating cells for single-cell analysis is not a trivial task. Here, we compare three digestion methods for human skin tissue, evaluating the cell yield and viability for each protocol. We find that the use of a sequential dissociation method with multiple enzymatic digestion steps produces the highest cell viability. Using single-cell sequencing, we show this method results in a relative increase in the proportion of non-antigen-presenting mast cells and CD8 T cells as well as a relative decrease in the proportion of antigen-presenting mast cells and KYNU+ CD4 T cells. Overall, our findings support the use of this sequential digestion method on freshly processed human skin samples for optimal cell yield and viability.

The molecular landscape of neurological disorders: insights from single-cell RNA sequencing in neurology and neurosurgery

Single-cell ribonucleic acid sequencing (scRNA-seq) has emerged as a transformative technology in neurological and neurosurgical research, revolutionising our comprehension of complex neurological disorders. In brain tumours, scRNA-seq has provided valuable insights into cancer heterogeneity, the tumour microenvironment, treatment resistance, and invasion patterns. It has also elucidated the brain tri-lineage cancer hierarchy and addressed limitations of current models. Neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis have been molecularly subtyped, dysregulated pathways have been identified, and potential therapeutic targets have been revealed using scRNA-seq. In epilepsy, scRNA-seq has explored the cellular and molecular heterogeneity underlying the condition, uncovering unique glial subpopulations and dysregulation of the immune system. ScRNA-seq has characterised distinct cellular constituents and responses to spinal cord injury in spinal cord diseases, as well as provided molecular signatures of various cell types and identified interactions involved in vascular remodelling. Furthermore, scRNA-seq has shed light on the molecular complexities of cerebrovascular diseases, such as stroke, providing insights into specific genes, cell-specific expression patterns, and potential therapeutic interventions. This review highlights the potential of scRNA-seq in guiding precision medicine approaches, identifying clinical biomarkers, and facilitating therapeutic discovery. However, challenges related to data analysis, standardisation, sample acquisition, scalability, and cost-effectiveness need to be addressed. Despite these challenges, scRNA-seq has the potential to transform clinical practice in neurological and neurosurgical research by providing personalised insights and improving patient outcomes.

Single-nuclei RNA-seq reveals skin cell responses to Aeromonas hydrophila infection in Chinese longsnout catfish Leiocassis longirostris

As the primary natural barrier that protects against adverse environmental conditions, the skin plays a crucial role in the innate immune response of fish, particularly in relation to bacterial infections. However, due to the diverse functionality and intricate anatomical and cellular composition of the skin, deciphering the immune response of the host is a challenging task. In this study, single nuclei RNA-sequencing (snRNA-seq) was performed on skin biopsies obtained from Chinese longsnout catfish (Leiocassis longirostris), comparing Aeromonas hydrophila-infected subjects to healthy control subjects. A total of 19,581 single nuclei cells were sequenced using 10x Genomics (10,400 in the control group and 9,181 in the treated group). Based on expressed unique transcriptional profiles, 33 cell clusters were identified and classified into 12 cell types including keratinocyte (KC), fibroblast (FB), endothelial cells (EC), secretory cells (SC), immune cells, smooth muscle cells (SMC), and other cells such as pericyte (PC), brush cell (BC), red blood cell (RBC), neuroendocrine cell (NDC), neuron cells (NC), and melanocyte (MC). Among these, three clusters of KCs, namely, KC1, KC2, and KC5 exhibited significant expansion after A. hydrophila infection. Analysis of pathway enrichment revealed that KC1 was primarily involved in environmental signal transduction, KC2 was primarily involved in endocrine function, and KC5 was primarily involved in metabolism. Finally, our findings suggest that neutrophils may play a crucial role in combating A. hydrophila infections. In summary, this study not only provides the first detailed comprehensive map of all cell types present in the skin of teleost fish but also sheds light on the immune response mechanism of the skin following A. hydrophila infection in Chinese longsnout catfish.

New insight into the role of fibroblasts in the epithelial immune microenvironment in the single-cell era

The skin is exposed to environmental challenges and contains heterogeneous cell populations such as epithelial cells, stromal cells, and skin-resident immune cells. As the most abundant type of stromal cells, fibroblasts have been historically considered silent observers in the immune responses of the cutaneous epithelial immune microenvironment (EIME), with little research conducted on their heterogeneity and immune-related functions. Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) have overcome the limitations of bulk RNA sequencing and help recognize the functional and spatial heterogeneity of fibroblasts, as well as their crosstalk with other types of cells in the cutaneous EIME. Recently, emerging single-cell sequencing data have demonstrated that fibroblasts notably participate in the immune responses of the EIME and impact the initiation and progression of inflammatory skin diseases. Here, we summarize the latest advances in the role of fibroblasts in the cutaneous EIME of inflammatory skin diseases and discuss the distinct functions and molecular mechanisms of activated fibroblasts in fibrotic skin diseases and non-fibrotic inflammatory skin diseases. This review help unveil the multiple roles of fibroblasts in the cutaneous EIME and offer new promising therapeutic strategies for the management of inflammatory skin diseases by targeting fibroblasts or the fibroblast-centered EIME.

Single-cell analysis of human prepuce reveals dynamic changes in gene regulation and cellular communications

BACKGROUND

The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells.

RESULTS

We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults.

CONCLUSIONS

Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

Assessment of Treatment-Relevant Immune Biomarkers in Psoriasis and Atopic Dermatitis: Toward Personalized Medicine in Dermatology

Immunologically targeted therapies have revolutionized the treatment of inflammatory dermatoses, including atopic dermatitis and psoriasis. Although immunologic biomarkers hold great promise for personalized classification of skin disease and tailored therapy selection, there are no approved or widely used approaches for this in dermatology. This review summarizes the translational immunologic approaches to measuring treatment-relevant biomarkers in inflammatory skin conditions. Tape strip profiling, microneedle-based biomarker patches, molecular profiling from epidermal curettage, RNA in situ hybridization tissue staining, and single-cell RNA sequencing have been described. We discuss the advantages and limitations of each and open questions for the future of personalized medicine in inflammatory skin disease.

Pellino-1 promotes intrinsic activation of skin-resident IL-17A-producing T cells in psoriasis

BACKGROUND

Psoriasis is a chronically relapsing inflammatory skin disease primarily perpetuated by skin-resident IL-17-producing T (T17) cells. Pellino-1 (Peli1) belongs to a member of E3 ubiquitin ligase mediating immune receptor signaling cascades, including nuclear factor kappa-light-chain enhancer of activated B cells (NF-κB) pathway.

OBJECTIVE

We explored the potential role of Peli1 in psoriatic inflammation in the context of skin-resident T17 cells.

METHODS

We performed single-cell RNA sequencing of relapsing and resolved psoriatic lesions with analysis for validation data set of psoriasis. Mice with systemic and conditional depletion of Peli1 were generated to evaluate the role of Peli1 in imiquimod-induced psoriasiform dermatitis. Pharmacologic inhibition of Peli1 in human CD4⁺ T cells and ex vivo human skin cultures was also examined to evaluate its potential therapeutic implications.

RESULTS

Single-cell RNA sequencing analysis revealed distinct T-cell subsets in relapsing psoriasis exhibiting highly enriched gene signatures for (1) tissue-resident T cells, (2) T17 cells, and (3) NF-κB signaling pathway including PELI1. Peli1-deficient mice were profoundly protected from psoriasiform dermatitis, with reduced IL-17A production and NF-κB activation in γδ T17 cells. Mice with conditional depletion of Peli1 treated with FTY720 revealed that Peli1 was intrinsically required for the skin-resident T17 cell immune responses. Notably, pharmacologic inhibition of Peli1 significantly ameliorated murine psoriasiform dermatitis and IL-17A production from the stimulated human CD4⁺ T cells and ex vivo skin explants modeling psoriasis.

CONCLUSION

Targeting Peli1 would be a promising therapeutic strategy for psoriasis by limiting skin-resident T17 cell immune responses.

Applications of single-cell RNA sequencing in atopic dermatitis and psoriasis

Single-cell RNA sequencing (scRNA-seq) is a novel technology that characterizes molecular heterogeneity at the single-cell level. With the development of more automated, sensitive, and cost-effective single-cell isolation methods, the sensitivity and efficiency of scRNA-seq have improved. Technological advances in single-cell analysis provide a deeper understanding of the biological diversity of cells present in tissues, including inflamed skin. New subsets of cells have been discovered among common inflammatory skin diseases, such as atopic dermatitis (AD) and psoriasis. ScRNA-seq technology has also been used to analyze immune cell distribution and cell-cell communication, shedding new light on the complex interplay of components involved in disease responses. Moreover, scRNA-seq may be a promising tool in precision medicine because of its ability to define cell subsets with potential treatment targets and to characterize cell-specific responses to drugs or other stimuli. In this review, we briefly summarize the progress in the development of scRNA-seq technologies and discuss the latest scRNA-seq-related findings and future trends in AD and psoriasis. We also discuss the limitations and technical problems associated with current scRNA-seq technology.

Revisiting the significance of keratin expression in complex epithelia

A large group of keratin genes (n=54 in the human genome) code for intermediate filament (IF)-forming proteins and show differential regulation in epithelial cells and tissues. Keratin expression can be highly informative about the type of epithelial tissue, differentiation status of constituent cells and biological context (e.g. normal versus diseased settings). The foundational principles underlying the use of keratin expression to gain insight about epithelial cells and tissues primarily originated in pioneering studies conducted in the 1980s. The recent emergence of single cell transcriptomics provides an opportunity to revisit these principles and gain new insight into epithelial biology. Re-analysis of single-cell RNAseq data collected from human and mouse skin has confirmed long-held views regarding the quantitative importance and pairwise regulation of specific keratin genes in keratinocytes of surface epithelia. Furthermore, such analyses confirm and extend the notion that changes in keratin gene expression occur gradually as progenitor keratinocytes commit to and undergo differentiation, and challenge the prevailing assumption that specific keratin combinations reflect a mitotic versus a post-mitotic differentiating state. Our findings provide a blueprint for similar analyses in other tissues, and warrant a more nuanced approach in the use of keratin genes as biomarkers in epithelia.

How "Neuronal" Are Human Skin Mast Cells?

Mast cells are evolutionarily old cells and the principal effectors in allergic responses and inflammation. They are seeded from the yolk sac during embryogenesis or are derived from hematopoietic progenitors and are therefore related to other leukocyte subsets, even though they form a separate clade in the hematopoietic system. Herein, we systematically bundle information from several recent high-throughput endeavors, especially those comparing MCs with other cell types, and combine such information with knowledge on the genes’ functions to reveal groups of neuronal markers specifically expressed by MCs. We focus on recent advances made regarding human tissue MCs, but also refer to studies in mice. In broad terms, genes hyper-expressed in MCs, but largely inactive in other myelocytes, can be classified into subcategories such as traffic/lysosomes (MLPH and RAB27B), the dopamine system (MAOB, DRD2, SLC6A3, and SLC18A2), Ca²⁺-related entities (CALB2), adhesion molecules (L1CAM and NTM) and, as an overall principle, the transcription factors and modulators of transcriptional activity (LMO4, PBX1, MEIS2, and EHMT2). Their function in MCs is generally unknown but may tentatively be deduced by comparison with other systems. MCs share functions with the nervous system, as they express typical neurotransmitters (histamine and serotonin) and a degranulation machinery that shares features with the neuronal apparatus at the synapse. Therefore, selective overlaps are plausible, and they further highlight the uniqueness of MCs within the myeloid system, as well as when compared with basophils. Apart from investigating their functional implications in MCs, a key question is whether their expression in the lineage is due to the specific reactivation of genes normally silenced in leukocytes or whether the genes are not switched off during mastocytic development from early progenitors.

Single-Cell Transcriptional Analysis Deciphers the Inflammatory Response of Skin-Resident Stromal Cells

The skin is the outermost barrier of the body. It has developed a sophisticated system against the ever-changing environment. The application of single-cell technologies has revolutionized dermatology research and unraveled the changes and interactions across skin resident cells in the healthy and inflamed skin. Single-cell technologies have revealed the critical roles of stromal cells in an inflammatory response and explained a series of plausible previous findings concerning skin immunity. Here, we summarized the functional diversity of skin stromal cells defined by single-cell analyses and how these cells orchestrated events leading to inflammatory diseases, including atopic dermatitis, psoriasis, vitiligo, and systemic lupus erythematosus.

An Optimized Tissue Dissociation Protocol for Single-Cell RNA Sequencing Analysis of Fresh and Cultured Human Skin Biopsies

We present an optimized dissociation protocol for preparing high-quality skin cell suspensions for in-depth single-cell RNA-sequencing (scRNA-seq) analysis of fresh and cultured human skin. Our protocol enabled the isolation of a consistently high number of highly viable skin cells from small freshly dissociated punch skin biopsies, which we use for scRNA-seq studies. We recapitulated not only the main cell populations of existing single-cell skin atlases, but also identified rare cell populations, such as mast cells. Furthermore, we effectively isolated highly viable single cells from ex vivo cultured skin biopsy fragments and generated a global single-cell map of the explanted human skin. The quality metrics of the generated scRNA-seq datasets were comparable between freshly dissociated and cultured skin. Overall, by enabling efficient cell isolation and comprehensive cell mapping, our skin dissociation-scRNA-seq workflow can greatly facilitate scRNA-seq discoveries across diverse human skin pathologies and ex vivo skin explant experimentations.

Single-cell transcriptomics in human skin research: available technologies, technical considerations, and disease applications

Single‐cell technologies have revolutionized research in the last decade, including for skin biology. Single‐cell RNA sequencing has emerged as a powerful tool allowing the dissection of human disease pathophysiology at unprecedented resolution by assessing cell‐to‐cell variation, facilitating identification of rare cell populations and elucidating cellular heterogeneity. In dermatology, this technology has been widely applied to inflammatory skin disorders, fibrotic skin diseases, wound healing complications and cutaneous neoplasms. Here, we discuss the available technologies and technical considerations of single‐cell RNA sequencing and describe its applications to a broad spectrum of dermatological diseases.

Single-Cell Transcriptome Reveals the Metabolic and Clinical Features of a Highly Malignant Cell Subpopulation in Pancreatic Ductal Adenocarcinoma

Background: Pancreatic ductal adenocarcinoma (PDAC) is a malignant tumor with a high mortality rate. PDAC exhibits significant heterogeneity as well as alterations in metabolic pathways that are associated with its malignant progression. In this study, we explored the metabolic and clinical features of a highly malignant subgroup of PDAC based on single-cell transcriptome technology.

Methods: A highly malignant cell subpopulation was identified at single-cell resolution based on the expression of malignant genes. The metabolic landscape of different cell types was analyzed based on metabolic pathway gene sets. In vitro experiments to verify the biological functions of the marker genes were performed. PDAC patient subgroups with highly malignant cell subpopulations were distinguished according to five glycolytic marker genes. Five glycolytic highly malignant-related gene signatures were used to construct the glycolytic highly malignant-related genes signature (GHS) scores.

Results: This study identified a highly malignant tumor cell subpopulation from the single-cell RNA sequencing (scRNA-seq) data. The analysis of the metabolic pathway revealed that highly malignant cells had an abnormally active metabolism, and enhanced glycolysis was a major metabolic feature. Five glycolytic marker genes that accounted for the highly malignant cell subpopulations were identified, namely, EN O 1, LDHA, PKM, PGK1, and PGM1. An in vitro cell experiment showed that proliferation rates of PANC-1 and CFPAC-1 cell lines decreased after knockdown of these five genes. Patients with metabolic profiles of highly malignant cell subpopulations exhibit clinical features of higher mortality, higher mutational burden, and immune deserts. The GHS score evaluated using the five marker genes was an independent prognostic factor for patients with PDAC.

Conclusion: We revealed a subpopulation of highly malignant cells in PDAC with enhanced glycolysis as the main metabolic feature. We obtained five glycolytic marker gene signatures, which could be used to identify PDAC patient subgroups with highly malignant cell subpopulations, and proposed a GHS prognostic score.

Isolating and cryopreserving pig skin cells for single-cell RNA sequencing study

The pig skin architecture and physiology are similar to those of humans. Thus, the pig model is very valuable for studying skin biology and testing therapeutics. The single-cell RNA sequencing (scRNA-seq) technology allows quantitatively analyzing cell types, compositions, states, signaling, and receptor-ligand interactome at single-cell resolution and at high throughput. scRNA-seq has been used to study mouse and human skins. However, studying pig skin with scRNA-seq is still rare. A critical step for successful scRNA-seq is to obtain high-quality single cells from the pig skin tissue. Here we report a robust method for isolating and cryopreserving pig skin single cells for scRNA-seq. We showed that pig skin could be efficiently dissociated into single cells with high cell viability using the Miltenyi Human Whole Skin Dissociation kit and the Miltenyi gentleMACS Dissociator. Furthermore, the obtained single cells could be cryopreserved using 90% FBS + 10% DMSO without causing additional cell death, cell aggregation, or changes in gene expression profiles. Using the developed protocol, we were able to identify all the major skin cell types. The protocol and results from this study are valuable for the skin research scientific community.

Challenges and Opportunities for the Translation of Single-Cell RNA Sequencing Technologies to Dermatology

Skin is a complex and heterogeneous organ at the cellular level. This complexity is beginning to be understood through the application of single-cell genomics and computational tools. A large number of datasets that shed light on how the different human skin cell types interact in homeostasis-and what ceases to work in diverse dermatological diseases-have been generated and are publicly available. However, translation of these novel aspects to the clinic is lacking. This review aims to summarize the state-of-the-art of skin biology using single-cell technologies, with a special focus on skin pathologies and the translation of mechanistic findings to the clinic. The main implications of this review are to summarize the benefits and limitations of single-cell analysis and thus help translate the emerging insights from these novel techniques to the bedside.

Single-cell transcriptomic analysis reveals the critical molecular pattern of UV-induced cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is the second most common nonmelanoma skin cancer characterized by high invasiveness, heterogeneity, and mainly occurs in the ultraviolet (UV)-exposed regions of the skin, but its pathogenesis is still unclear. Here, we generated single-cell transcriptome profiles for 350 cells from six primary UV-induced cSCCs, together with matched adjacent skin samples, and three healthy control skin tissues by single-cell RNA-sequencing technology based on Smart-seq2 strategy. A series of bioinformatics analyses and in vitro experiments were used to decipher and validate the critical molecular pattern of cSCC. Results showed that cSCC cells and normal keratinocytes were significantly distinct in gene expression and chromosomal copy number variation. Furthermore, cSCC cells exhibited 18 hallmark pathways of cancer by gene set enrichment analysis. Differential expression analysis demonstrated that many members belonging to S100 gene family, SPRR gene family, and FABP5 were significantly upregulated in cSCC cells. Further experiments confirmed their upregulation and showed that S100A9 or FABP5 knockdown in cSCC cells inhibited their proliferation and migration through NF-κB pathway. Taken together, our data provide a valuable resource for deciphering the molecular pattern in UV-induced cSCC at a single-cell level and suggest that S100A9 and FABP5 may provide novel targets for therapeutic intervention of cSCC in the future.

Generation of Skin Organoids: Potential Opportunities and Challenges

Although several types of human skin substitutes are currently available, they usually do not include important skin appendages such as hair follicles and sweat glands, or various skin-related cells, such as dermal adipocytes and sensory neurons. This highlights the need to improve the in vitro human skin generation model for use as a tool for investigating skin diseases and as a source of cells or tissues for skin regeneration. Skin organoids are generated from stem cells and are expected to possess the complexity and function of natural skin. Here, we summarize the current literatures relating to the "niches" of the local skin stem cell microenvironment and the formation of skin organoids, and then discuss the opportunities and challenges associated with multifunctional skin organoids.

Human Dermal Fibroblast Subpopulations Are Conserved across Single-Cell RNA Sequencing Studies

On the basis of their differential location within the dermis and of discrete changes in gene and protein expression, two major fibroblast subtypes (papillary and reticular) have traditionally been distinguished. In the last 3 years, a number of research groups have begun to address transcriptomic heterogeneity of human skin cells at the single-cell level by determining mRNA levels of expressed genes through single-cell RNA sequencing technologies. However, the outcome of single-cell RNA sequencing studies is thus far confusing. Very little overlap was found in fibroblast subpopulations, which also varied in number and composition in each dataset. After a careful reappraisal of the transcriptomic data of 13,823 human adult dermal fibroblasts that have been sequenced to date, we show that fibroblasts may robustly be assigned to three major types (axes A‒C), which in turn are composed of 10 major subtypes (clusters), which we denominated A1‒A4, B1 and B2, and C1‒C4. These computationally determined axes and clusters represent the major fibroblast types and subtypes in adult healthy human skin across different datasets, accounting for 92.5% of the sequenced fibroblasts. They thus may provide the basis to improve our understanding of dermal homeostasis and cellular function at the transcriptomic level.