Skin Network

Skin Network - Overview

Skin undergoes constant self-renewal, and its functional decline is a visible consequence of aging. Understanding human skin aging requires in-depth knowledge of the molecular and functional properties of various skin cell types. We performed single-cell RNA sequencing of human eyelid skin from healthy individuals across different ages and identified eleven canonical cell types, as well as six subpopulations of basal cells. Further analysis revealed progressive accumulation of photoaging-related changes and increased chronic inflammation with age. Transcriptional factors involved in the developmental process underwent early-onset decline during aging. Furthermore, inhibition of key transcription factors HES1 in fibroblasts and KLF6 in keratinocytes not only compromised cell proliferation, but also increased inflammation and cellular senescence during aging. Lastly, we found that genetic activation of HES1 or pharmacological treatment with quercetin alleviated cellular senescence of dermal fibroblasts. These findings provide a single-cell molecular framework of human skin aging, providing a rich resource for developing therapeutic strategies against aging-related skin disorders.

A Single-Cell Transcriptomic Atlas of Human Skin Aging.

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for approximately 1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems. Overall design: Applying single cell chromatin accessibility assays to diverse human tissue types from four donors.

A single-cell atlas of chromatin accessibility in the human genome

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.

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

The skin serves as a physical barrier and an immunological interface that protects the body from the external environment1-3. Aberrant activation of immune cells can induce common skin autoimmune diseases such as vitiligo, which are often characterized by bilateral symmetric lesions in certain anatomic regions of the body4-6. Understanding what orchestrates the activities of cutaneous immune cells at an organ level is necessary for the treatment of autoimmune diseases. Here we identify subsets of dermal fibroblasts that are responsible for driving patterned autoimmune activity, by using a robust mouse model of vitiligo that is based on the activation of endogenous auto-reactive CD8+ T cells that target epidermal melanocytes. Using a combination of single-cell analysis of skin samples from patients with vitiligo, cell-type-specific genetic knockouts and engraftment experiments, we find that among multiple interferon-γ (IFNγ)-responsive cell types in vitiligo-affected skin, dermal fibroblasts are uniquely required to recruit and activate CD8+ cytotoxic T cells through secreted chemokines. Anatomically distinct human dermal fibroblasts exhibit intrinsic differences in the expression of chemokines in response to IFNγ. In mouse models of vitiligo, regional IFNγ-resistant fibroblasts determine the autoimmune pattern of depigmentation in the skin. Our study identifies anatomically distinct fibroblasts with permissive or repressive IFNγ responses as the key determinant of body-level patterns of lesions in vitiligo, and highlights mesenchymal subpopulations as therapeutic targets for treating autoimmune diseases.

Anatomically distinct fibroblast subsets determine skin autoimmune patterns.

Inflammatory conditions represent the largest class of chronic skin disease, but the molecular dysregulation underlying many individual cases remains unclear. Single-cell RNA sequencing (scRNA-seq) has increased precision in dissecting the complex mixture of immune and stromal cell perturbations in inflammatory skin disease states. We single-cell-profiled CD45+ immune cell transcriptomes from skin samples of 31 patients (7 atopic dermatitis, 8 psoriasis vulgaris, 2 lichen planus (LP), 1 bullous pemphigoid (BP), 6 clinical/histopathologically indeterminate rashes, and 7 healthy controls). Our data revealed active proliferative expansion of the Treg and Trm components and universal T cell exhaustion in human rashes, with a relative attenuation of antigen-presenting cells. Skin-resident memory T cells showed the greatest transcriptional dysregulation in both atopic dermatitis and psoriasis, whereas atopic dermatitis also demonstrated recurrent abnormalities in ILC and CD8+ cytotoxic lymphocytes. Transcript signatures differentiating these rash types included genes previously implicated in T helper cell (TH2)/TH17 diatheses, segregated in unbiased functional networks, and accurately identified disease class in untrained validation data sets. These gene signatures were able to classify clinicopathologically ambiguous rashes with diagnoses consistent with therapeutic response. Thus, we have defined major classes of human inflammatory skin disease at the molecular level and described a quantitative method to classify indeterminate instances of pathologic inflammation. To make this approach accessible to the scientific community, we created a proof-of-principle web interface (RashX), where scientists and clinicians can visualize their patient-level rash scRNA-seq-derived data in the context of our TH2/TH17 transcriptional framework.

Classification of human chronic inflammatory skin disease based on single-cell immune profiling

4 independent pools of iPS cells were differentiated and used to generate kidney organoids. Each pool was used to generate each of the 3 types of conformations: Manual organoids, bioprinted organoids composed of a single deposition of cells (R0) and bioprinted organoids composed a thin line of cells (R40). Manual and two types of bioprinted kidney organoids were analysed by single cell RNA-seq.

Comparison of manual and two types of bioprinted kidney organoids by single cell RNA-seq.

Single-cell RNA sequencing was used to determine the cell type composition of all major human organs and construct a basic scheme for the human cell landscape (HCL).

Construction of A Human Cell Landscape by Single-cell mRNA-seq

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells/multipotent progenitors (HSC/MPPs), yet remains poorly defined in humans. Using single cell transcriptome profiling of ~133,000 fetal liver and ~65,000 fetal skin and kidney cells, we identify the repertoire of blood and immune cells in first and early second trimesters of development. From this data, we infer differentiation trajectories from HSC/MPPs, and evaluate the impact of tissue microenvironment on blood and immune cell development. We predict coupling of mast cell differentiation with erythro-megakaryopoiesis and identify physiological erythropoiesis in fetal skin. We demonstrate a shift in fetal liver haematopoietic composition during gestation away from being erythroid-predominant, accompanied by a parallel change in HSC/MPP differentiation potential, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a valuable reference for understanding and harnessing the therapeutic potential of HSC/MPPs.

Decoding human fetal liver haematopoiesis

The epidermis and its appendage, the hair follicle, represent an elegant developmental system in which cells are replenished with regularity because of controlled proliferation, lineage specification, and terminal differentiation. Although transcriptome data exists for human epidermal and dermal cells, the hair follicle remains poorly characterized. Through single-cell resolution profiling of the epidermis and anagen hair follicle, we characterized the anatomical, transcriptional, functional, and pathological profiles of distinct epidermal, hair follicle, and hair follicle-associated cell subpopulations including melanocytes, endothelial cells, and immune cells. We additionally traced the differentiation trajectory of interfollicular and matrix cell progenitors and explored the association of specific cell subpopulations to known molecular signatures of common skin conditions. These data simultaneously corroborate prior murine and human studies while offering new insights into epidermal and hair follicle differentiation and pathogenesis.

Defining Transcriptional Signatures of Human Hair Follicle Cell States

We generated a comprehensive atlas of human skin in early prenatal life, in adulthood, and during inflammatory skin disease by profiling the transcriptomes of more than 500,000 single cells. We analyzed human embryonic skin between 7 and 10 post-conception weeks, healthy adult skin surplus from mammoplasty surgery, and skin biopsies from patients affected by atopic dermatitis (AD) and psoriasis, two common inflammatory skin diseases.

Developmental cell programs are co-opted in inflammatory skin disease.

Palmoplantar skin is structurally and functionally unique, but the transcriptional programs driving this specialization are unclear. Here, we use bulk and single-cell RNA sequencing of human palm, sole, and hip skin to describe the distinguishing characteristics of palmoplantar and non-palmoplantar skin while also uncovering differences between palmar and plantar sites. Our approach reveals an altered immune environment in palmoplantar skin, with downregulation of diverse immunological processes and decreased immune cell populations. Further, we identify specific fibroblast populations that appear to orchestrate key differences in cell-cell communication in palm, sole, and hip. Dedicated keratinocyte analysis highlights major differences in basal cell fraction among the three sites and demonstrates the existence of two spinous keratinocyte populations constituting parallel, site-selective epidermal differentiation trajectories. In summary, this deep characterization of highly adapted palmoplantar skin contributes key insights into the fundamental biology of human skin and provides a valuable data resource for further investigation.

Differential cell composition and split epidermal differentiation in human palm, sole, and hip skin

Cellular heterogeneity can emerge from the expression of only one parental allele, however it has remained unknown to what degree patterns of random monoallelic expression of autosomal genes (aRME) are mitotically inherited (clonal) or stochastic (dynamic) in somatic cells. Here, we resolved this by applying allele-sensitive single-cell RNA-seq on primary mouse fibroblasts and in vivo human T-cells to simultaneously investigate clonal and dynamic aRME. Dynamic aRME affected a considerable portion of the transcriptome, with levels dependent on the cell’s transcriptional activity, but clonal aRME was surprisingly scarce (<1% of genes) and affected mainly lowly expressed genes. Consequently, the overwhelming portion of aRME occurs transiently and is scattered throughout somatic populations rather than, as previously hypothesized, confined to spatially restricted patches of clonally related cells.

Dissecting the clonal nature of allelic expression in somatic cells by single-cell RNA-seq

Direct comparisons of human and non-human primate brain tissue have the potential to reveal molecular pathways underlying remarkable specializations of the human brain. However, chimpanzee tissue is largely inaccessible during neocortical neurogenesis when differences in brain size first appear. To identify human-specific features of cortical development, we leveraged recent innovations that permit generating pluripotent stem cell-derived cerebral organoids from chimpanzee. First, we systematically evaluated the fidelity of organoid models to primary human and macaque cortex, finding organoid models preserve gene regulatory networks related to cell types and developmental processes but exhibit increased metabolic stress. Second, we identified 261 genes differentially expressed in human compared to chimpanzee organoids and macaque cortex. Many of these genes overlap with human-specific segmental duplications and a subset suggest increased PI3K/AKT/mTOR activation in human outer radial glia. Together, our findings establish a platform for systematic analysis of molecular changes contributing to human brain development and evolution. Overall design: Single cell mRNA sequencing of iPS-derived neural and glial progenitor cells using the Fluidigm C1 system This series includes re-analysis of publicly available data in accessions: phs000989.v3, GSE99951, GSE86207, GSE75140. Sample metadata and accession IDs for the re-analyzed samples are included in the file "GSE124299_metadata_on_processed_samples.xlsx" available on the foot of this record. The following samples have no raw data due to data loss: GSM3569728, GSM3569738, GSM3571601, GSM3571606, GSM3571615, GSM3571621, GSM3571625, and GSM3571631

Establishing Cerebral Organoids as Models of Human-Specific Brain Evolution

Reprogramming, CRISPR/Cas9 gene-editing, organoid technologies and single-cell RNA sequencing technology were combined to study the nephron lineage in a human context.

Fate-mapping within human iPSC-derived kidney organoids reveals conserved mammalian nephron progenitor lineage relationships.

Unbiased whole-genome methylome was studied in the wound-edge (WE) tissue of chronic wound patients. Methylation status of proximal promoter (1Kb) was calculated using MethylCap-Seq and PrEMeR-CG data analyses. A total of 4689 differentially methylated regions (DMRs) were identified in chronic WE compared to unwounded (UW) human skin. Hypermethylation was more frequently observed (3661 DMRs) in the chronic WE compared to hypomethylation (1028 DMRs). Twenty-six hypermethylated DMRs were involved in epithelial to mesenchymal transition (EMT). Bisulfite sequencing validated hypermethylation of a predicted specific upstream regulator TP53. Whole genome RNA sequencing analysis was performed to qualify findings from methylome analysis. Hierarchical clustering analyses identified a large set of genes, the expression of which were significantly downregulated in chronic WE compared to UW skin. Analysis of the downregulated genes identified the TP53 signaling pathway, including P63, P73, FOXO3, SIRT1 and HDAC1, as being significantly silenced. Direct comparison of hypermethylation and downregulated gene expression identified four genes, ADAM17, NOTCH, TWIST1 and SMURF1, that were common to both data sets and functionally represented the EMT pathway. Single cell RNA sequencing studies identified that these effects on gene expression were limited to the keratinocyte cell compartment. Overall design: In this work, scRNA-seq analyses of human chronic WE characterized its heterogenous cellular composition as compared to unwounded skin.

Genome-Wide DNA Hypermethylation in the Wound-Edge of Chronic Wound Patients Opposes Closure by Impairing Epithelial to Mesenchymal Transition

Cerebral organoids – three-dimensional cultures of human cerebral tissue derived from pluripotent stem cells – have emerged as models of human cortical development. However, the extent to which in vitro organoid systems recapitulate neural progenitor cell proliferation and neuronal differentiation programs observed in vivo remains unclear. Here we use single-cell RNA sequencing (scRNA-seq) to dissect and compare cell composition and progenitor-to-neuron lineage relationships in human cerebral organoids and fetal neocortex. Covariation network analysis using the fetal neocortex data reveals known and novel interactions among genes central to neural progenitor proliferation and neuronal differentiation. In the organoid, we detect diverse progenitors and differentiated cell types of neuronal and mesenchymal lineages, and identify cells that derived from regions resembling the fetal neocortex. We find that these organoid cortical cells use gene expression programs remarkably similar to those of the fetal tissue in order to organize into cerebral cortex-like regions. Our comparison of in vivo and in vitro cortical single cell transcriptomes illuminates the genetic features underlying human cortical development that can be studied in organoid cultures.

Human cerebral organoids recapitulate gene expression programs of fetal neocortex development.

In humans, epidermal melanocytes are responsible for skin pigmentation, defence against ultraviolet radiation and the deadliest common skin cancer, melanoma. Although there is substantial overlap in melanocyte development pathways between different model organisms, species-dependent differences are frequent and the conservation of these processes in human skin remains unresolved. Here, we used a single-cell enrichment and RNA-sequencing pipeline to study human epidermal melanocytes directly from the skin, capturing transcriptomes across different anatomical sites, developmental age, sexes and multiple skin tones. We uncovered subpopulations of melanocytes that exhibit anatomical site-specific enrichment that occurs during gestation and persists through adulthood. The transcriptional signature of the volar-enriched subpopulation is retained in acral melanomas. Furthermore, we identified human melanocyte differentiation transcriptional programs that are distinct from gene signatures generated from model systems. Finally, we used these programs to define patterns of dedifferentiation that are predictive of melanoma prognosis and response to immune checkpoint inhibitor therapy.

Human melanocyte development and melanoma dedifferentiation at single-cell resolution

Brain organoids represent a powerful tool for studying human neurological diseases, particularly those that affect brain growth and structure. However, many diseases manifest with clear evidence of physiological and network abnormality in the absence of anatomical changes, raising the question of whether organoids possess sufficient neural network complexity to model these conditions. Here, we explore the network-level functions of brain organoids using calcium sensor imaging and extracellular recording approaches that together reveal the existence of complex network dynamics reminiscent of intact brain preparations. We demonstrate highly abnormal and epileptiform-like activity in organoids derived from induced pluripotent stem cells from individuals with Rett syndrome, accompanied by transcriptomic differences revealed by single-cell analyses. We also rescue key physiological activities with an unconventional neuroregulatory drug, pifithrin-α. Together, these findings provide an essential foundation for the utilization of brain organoids to study intact and disordered human brain network formation and illustrate their utility in therapeutic discovery.

Identification of neural oscillations and epileptiform changes in human brain organoids.

Bud tip progenitor cells give rise to all murine lung epithelial lineages and have been described in the developing human lung; however, the mechanisms controlling human bud tip differentiation into specific lineages are unclear. Here, we used homogeneous human bud tip organoid cultures and identified SMAD signaling as a key regulator of the bud tip-to-airway transition. SMAD induction led to the differentiation of airway-like organoids possessing functional basal cells capable of clonal expansion and multilineage differentiation. To benchmark in vitro-derived organoids, we developed a single-cell mRNA sequencing atlas of the human lung from 11.5 to 21 weeks of development, which revealed high degrees of similarity between the in vitro-derived and in vivo airway. Together, this work sheds light on human airway differentiation in vitro and provides a single-cell atlas of the developing human lung.

In Vitro and In Vivo Development of the Human Airway at Single-Cell Resolution

Induced pluripotent stem cell (iPSC)-derived organoids provide models to study human organ development. Single-cell transcriptomics enable highly resolved descriptions of cell states within these systems; however, approaches are needed to directly measure lineage relationships. Here we establish iTracer, a lineage recorder that combines reporter barcodes with inducible CRISPR–Cas9 scarring and is compatible with single-cell and spatial transcriptomics. We apply iTracer to explore clonality and lineage dynamics during cerebral organoid development and identify a time window of fate restriction as well as variation in neurogenic dynamics between progenitor neuron families.

Lineage recording in human cerebral organoids.

The authors captured 301 single cells from 11 populations using microfluidics and analyzing single-cell transcriptomes across downsampled sequencing depths. They demonstrate that shallow single-cell mRNA sequencing (∼50,000 reads per cell) is sufficient for unbiased cell-type classification and biomarker identification in the developing cortex.

Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex.

Recent advances in single cell genomics technologies have facilitated studies on the developing immune system at unprecedented scale and resolution. However, these studies have focused on one or a few organs and were thus limited in understanding the developing immune system as a distributed network across tissues. Here, we profiled prenatal haematopoietic organs, lymphoid organs and non-lymphoid tissues using a combination of single-cell RNA sequencing, paired antigen-receptor sequencing and spatial transcriptomics to reconstruct the developing human immune system. Our analysis revealed the acquisition of immune effector transcriptome profiles in macrophages, mast cells and NK cells from the second trimester, and the transcriptomic changes accompanying the late-stage maturation of developing monocytes and T cells that extended from their organ of origin to peripheral tissues. We uncovered system-wide blood and immune cell development beyond the conventional primary haematopoietic organs. We further identified, extensively characterised and functionally validated the human prenatal B1 cells. Finally, we provide evidence for thymocyte-thymocyte selection origin for αβTCR- expressing unconventional T cells based on TCR gene usage and an in vitro artificial thymic organoid culture model. Our comprehensive atlas of the developing human immune system provides both valuable data resources and biological insights that will facilitate cell engineering, regenerative medicine and disease understanding. 

One-Sentence Summary 
By performing a comprehensive single-cell RNA sequencing atlas of human developing immune system together with antigen-receptor sequencing and spatial transcriptomics, we explored the cross-gestation and cross-organ variability in immune cells, discovered system-wide blood and immune cell development, identified, characterised and functionally validated the properties of human prenatal B1 cells and the origin of unconventional T cells.

Mapping the developing human immune system across organs

The cancer microenvironment is a complex ecosystem characterized by dynamic interactions between diverse cell types, including malignant, immune and stromal cells. Here, we performed single-cell RNA sequencing on CD45+ and CD45- cells isolated from tumour and lymph nodes during a mouse model of melanoma. The transcriptional profiles of these individual cells taken at different time points coupled with assembled T cell receptor sequences, allowed us to identify distinct immune subpopulations and delineate their developmental trajectory. Our study provides insights into the complex interplay among cells within the tumour microenvironment and presents a valuable resource for future translational applications.

Melanoma infiltration of stromal and immune cells

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this we have created a multi-scale spatial atlas of healthy skin and basal cell carcinoma (BCC), incorporating in vivo optical coherence tomography, single cell RNA sequencing, spatial global transcriptional profiling and in situ sequencing. Computational spatial deconvolution and projection revealed the localisation of distinct cell populations to specific tissue contexts. Although cell populations were conserved between healthy anatomical sites and in BCC, mesenchymal cell populations including fibroblasts and pericytes retained signatures of developmental origin. Spatial profiling and in silico lineage tracing support a hair follicle origin for BCC and demonstrate that cancer-associated fibroblasts are an expansion of a POSTN + subpopulation associated with hair follicles in healthy skin. RGS5+ pericytes are also expanded in BCC suggesting a role in vascular remodelling. We propose that the identity of mesenchymal cell populations is regulated by signals emanating from adjacent structures and that these signals are repurposed to promote the expansion of skin cancer stroma. The resource we have created is publicly available in an interactive format for the research community. Significance statement Single cells RNA sequencing has revolutionised cell biology, enabling high resolution analysis of cell types and states within human tissues. Here, we report a comprehensive spatial atlas of adult human skin across different anatomical sites and basal cell carcinoma (BCC) - the most common form of skin cancer - encompassing in vivo optical coherence tomography, single cell RNA sequencing, global spatial transcriptomic profiling and in situ sequencing. In combination these modalities have allowed us to assemble a comprehensive nuclear-resolution atlas of cellular identity in health and disease.

Multi-scale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma

Skin and lung fibrosis in systemic sclerosis (SSc) is driven by myofibroblasts, alpha-smooth muscle actin expressing cells. The number of myofibroblasts in SSc skin correlates with the modified Rodnan skin score, the most widely used clinical measure of skin disease severity. Murine fibrosis models indicate that myofibroblasts can arise from a variety of different cell types, but their origin in SSc skin has remained uncertain. Utilizing single cell RNA-sequencing, we define different dermal fibroblast populations and transcriptome changes, comparing SSc to healthy dermal fibroblasts. Here, we show that SSc dermal myofibroblasts arise in two steps from an SFRP2hi/DPP4-expressing progenitor fibroblast population. In the first step, SSc fibroblasts show globally upregulated expression of transcriptome markers, such as PRSS23 and THBS1. A subset of these cells shows markers indicating that they are proliferating. Only a fraction of SFRP2hi SSc fibroblasts differentiate into myofibroblasts, as shown by expression of additional markers, SFRP4 and FNDC1. Bioinformatics analysis of the SSc fibroblast transcriptomes implicated upstream transcription factors, including FOSL2, RUNX1, STAT1, FOXP1, IRF7 and CREB3L1, as well as SMAD3, driving SSc myofibroblast differentiation.

Myofibroblast transcriptome indicates SFRP2hi fibroblast progenitors in systemic sclerosis skin

22q11.2 deletion syndrome (22q11DS) is a highly penetrant and common genetic cause of neuropsychiatric disease. Here we generated induced pluripotent stem cells from 15 individuals with 22q11DS and 15 control individuals and differentiated them into three-dimensional (3D) cerebral cortical organoids. Transcriptional profiling across 100 days showed high reliability of differentiation and revealed changes in neuronal excitability-related genes. Using electrophysiology and live imaging, we identified defects in spontaneous neuronal activity and calcium signaling in both organoid- and 2D-derived cortical neurons. The calcium deficit was related to resting membrane potential changes that led to abnormal inactivation of voltage-gated calcium channels. Heterozygous loss of DGCR8 recapitulated the excitability and calcium phenotypes and its overexpression rescued these defects. Moreover, the 22q11DS calcium abnormality could also be restored by application of antipsychotics. Taken together, our study illustrates how stem cell derived models can be used to uncover and rescue cellular phenotypes associated with genetic forms of neuropsychiatric disease.

Neuronal defects in a human cellular model of 22q11.2 deletion syndrome.

The purpose of this project is to assess the relevance of pluripotent stem cell-derived cerebral and liver organoids to recapitulate the variation in cell-type specific gene expression programs between individuals. Towards this aim, we will generate reference atlases of the developing cortex and liver from multiple individuals, derive iPSC lines from these same individuals, and determine if inter-individual gene expression variation is recapitulated in cerebral and liver organoids from the same individual from which we have reference maps. In parallel we will assess the genetic contribution to variablity between organoids from different iPSCs of multiple human individuals that are available in existing iPSC resources (e.g. HipSci).

Organoid single-cell genomic atlas uncovers human-specific features of brain development

An efficient pooling strategy was used to differentiate 215 iPS cell lines towards a midbrain neural fate, including dopaminergic neurons. Over 1 million cells were sampled across three differentiation timepoints and profiled using single cell RNA sequencing.

Population-scale single-cell RNA-seq profiling across dopaminergic neuron differentiation.

Here we demonstrate single-cell RNA-5′end-sequencing (sc-end5-seq) methods can detect transcribed CREs (tCREs), enabling simultaneous quantification of gene expression and enhancer activities in a single assay at no extra cost. We showed enhancer RNAs can be detected using sc-end5-seq methods with either random or oligo(dT) priming.

Profiling of transcribed cis-regulatory elements in single cells.

Transcriptomic profiling of skin wounds in a time-course manner. Full thickness punch biopsies were performed on the flank of healthy donors to create skin wounds. This was followed by larger biopsies, to encompass the original wounds, taken between one to four days later. The skin biopsies and blood samples were subjected to single-cell RNA sequencing analysis.

Re-evaluation of human BDCA-2+ DC during acute sterile skin inflammation.

To assess the transcriptomic profile of disease-specific cell populations, fibroblasts from patients with primary open-angle glaucoma (POAG) were reprogrammed into induced pluripotent stem cells (iPSCs) before being differentiated into retinal organoids and compared to those from healthy individuals. We performed single-cell RNA-sequencing of a total of 247,520 cells and identified cluster-specific molecular signatures. Comparing the gene expression profile between cases and controls, we identified novel genetic associations for this blinding disease. Expression quantitative trait mapping identified a total of 4,443 significant loci across all cell types, 496 of which are specific to the retinal ganglion cell subpopulations, which ultimately degenerate in POAG. Transcriptome-wide association analysis identified genes at loci previously associated with POAG, and analysis, conditional on disease status, implicated 97 statistically significant retinal ganglion cell-specific expression quantitative trait loci. This work highlights the power of large-scale iPSC studies to uncover context-specific profiles for a genetically complex disease.

Retinal ganglion cell-specific genetic regulation in primary open angle glaucoma

Myeloid cells are central to homeostasis and immunity. Characterising in vitro myelopoiesis protocols is imperative for their use in research, immunotherapies, and understanding human myelopoiesis. The authors generated a >470K cells molecular map of human induced pluripotent stem cells (iPSC) differentiation into macrophages. Integration with in vivo single-cell atlases showed in vitro differentiation recapitulates features of yolk sac hematopoiesis, before definitive hematopoietic stem cells (HSC) emerge. The diversity of myeloid cells generated, including mast cells and monocytes, suggests that HSC-independent hematopoiesis can produce multiple myeloid lineages.

Robust temporal map of human in vitro myelopoiesis using single-cell genomics.

Lupus nephritis is a leading cause of mortality among systemic lupus erythematosus (SLE) patients, and its heterogeneous nature poses a significant challenge to the development of effective diagnostics and treatments. Single cell RNA sequencing (scRNA-seq) offers a potential solution to dissect the heterogeneity of the disease and enables the study of similar cell types distant from the site of renal injury to identify novel biomarkers. We applied scRNA-seq to human renal and skin biopsy tissues and demonstrated that scRNA-seq can be performed on samples obtained during routine care. Chronicity index, IgG deposition, and quantity of proteinuria correlated with a transcriptomic-based score composed of IFN-inducible genes in renal tubular cells. Furthermore, analysis of cumulative expression profiles of single cell keratinocytes dissociated from nonlesional, non–sun-exposed skin of patients with lupus nephritis also revealed upregulation of IFN-inducible genes compared with keratinocytes isolated from healthy controls. This indicates the possible use of scRNA-seq analysis of skin biopsies as a biomarker of renal disease. These data support the potential utility of scRNA-seq to provide new insights into the pathogenesis of lupus nephritis and pave the way for exploiting a readily accessible tissue to reflect injury in the kidney.

Single cell RNA sequencing to dissect the molecular heterogeneity in lupus nephritis

These files represent single cell RNA-Seq data generated on a 10x Chromium genomics platform from three biological replicates from the embryonic day (E)18.5 developing mouse kidney and three biological replicates of iPSC-derived human kidney organoids differentiated according to our published protocol (Takasato et al., Nature Protocols 2016). When aggregated, the mouse data represents >6000 cells that passed our QC, containing most major cell types known to exist in the developing mouse kidney. The aggregated human organoid data contains of >7000 cells that passed our QC and contains populations representing endothelial cells, podocytes, stroma, nephron, and off-target populations with similarity to neurons.

Single cell RNA-Seq of E18.5 developing mouse kidney and human kidney organoids

Here, we studied transcriptome profiles of human nail units from polydactly samples. Single cell RNAseq with 11,541 cells from 4 extra digits revealed nail-specific mesenchymal and epithelial cell populations Overall design: We profiled transcriptomes of nail unit from 4 fresh human polydactyly specimens using single cell RNAseq (scRNAseq). Extra-digits from 6-12 months old patients were delivered to the dermatologic clinic immediately after the digit removal. All samples containing nail mesenchyme and epithelium were minced and dissociated within 3 hours of digit removal.

Single cell RNAseq of human nail unit defines RSPO4 onychofibroblasts and SPINK6 nail epithelium

Endometriosis is characterized by growth of endometrial-like tissue outside of the uterus affecting many women during their reproductive age, causing years of pelvic pain and potential infertility. Its pathophysiology remains largely unknown, limiting early diagnosis and treatment. We characterized peritoneal and ovarian lesions at single-cell transcriptome resolution and compared to matched eutopic endometrium, control endometrium, and organoids derived from these tissues, generating data on over 122,000 cells across 14 individuals. We spatially localized many of the cell types using imaging mass cytometry. We identify a perivascular mural cell unique to the peritoneal lesions with dual roles in angiogenesis promotion and immune cell trafficking. We define an immunotolerant peritoneal niche, fundamental differences in eutopic endometrium and between lesions microenvironments, and an unreported progenitor-like epithelial cell subpopulation. Altogether, this study provides a holistic view of the endometriosis microenvironment representing a comprehensive cell atlas of the disease in hormonally treated patients, essential information for advancing therapeutics and diagnostics. Overall design: Biopsies were collected from a total of 27 patients, 19 with confirmed endometriosis and 8 patients control. Biopsies were subjected to single cell experiments, bulk sequencing and cell hashing on patient-derived organoid. Single cells experiement were conducted on 31 biopsies consisting of 3 control and 9 eutopic endometrium, 8 ectopic peritoneal, 6 ectopic peritoneal adjacent, 4 ectopic ovary of endometriosis patient. A total of 24 biopsies were subjected to bulk sequencing, consist of 5 control and 7 eutopic endometrium, 6 ectopic peritoneal, 6 ectopic ovary of endometriosis patient. Single-cell cell hashing experiment were conducted on patient-derived organoid from 2 control and 4 eutopic endometrium, 4 ectopic peritoneum, 1 ectopic peritoneal adjacent of endometriosis patient. All single cell experiments were performed on 10x Genomics platform.

Single cell analysis of endometriosis reveals a coordinated transcriptional program driving immunotolerance and angiogenesis across eutopic and ectopic tissues

Psoriasis is the most common skin disease in adults. Current experimental and clinical evidences suggested the infiltrating immune cells could target local skin cells and thus induce psoriatic phenotype. However, recent studies indicated the existence of a potential feedback signaling loop from local resident skin cells to infiltrating immune cells. Here, we deconstructed the full-thickness human skins of both healthy donors and patients with psoriasis vulgaris at single cell transcriptional level, and further built a neural-network classifier to evaluate the evolutional conservation of skin cell types between mouse and human. Last, we systematically evaluated the intrinsic and intercellular molecular alterations of each cell type between healthy and psoriatic skin. Cross-checking with psoriasis susceptibility gene loci, cell-type based differential expression, and ligand-receptor communication revealed that the resident psoriatic skin cells including mesenchymal and epidermis cell types, which specifically harbored the target genes of psoriasis susceptibility loci, intensively evoked the expression of major histocompatibility complex (MHC) genes, upregulated interferon (INF), tumor necrosis factor (TNF) signalling and increased cytokine gene expression for primarily aiming the neighboring dendritic cells in psoriasis. The comprehensive exploration and pathological observation of psoriasis patient biopsies proposed an uncovered immunoregulatory axis from skin local resident cells to immune cells, thus provided a novel insight for psoriasis treatment. In addition, we published a user-friendly website to exhibit the transcriptional change of each cell type between healthy and psoriatic human skin.

Single cell transcriptional zonation of human psoriasis skin identifies an alternative immunoregulatory axis conducted by skin resident cells

How stem cells give rise to epidermis is unclear despite the crucial role the epidermis plays in barrier and appendage formation. Here we use single cell-RNA sequencing to interrogate basal stem cell heterogeneity of human interfollicular epidermis and find four spatially distinct stem cell populations at the top and bottom of rete ridges and transitional positions between the basal and suprabasal epidermal layers. Cell-cell communication modeling suggests that basal cell populations serve as crucial signaling hubs to maintain epidermal communication. Combining pseudotime, RNA velocity, and cellular entropy analyses point to a hierarchical differentiation lineage supporting multi-stem cell interfollicular epidermal homeostasis models and suggest that transitional basal stem cells are stable states essential for proper stratification. Finally, alterations in differentially expressed transitional basal stem cell genes result in severe thinning of human skin equivalents, validating their essential role in epidermal homeostasis and reinforcing the critical nature of basal stem cell heterogeneity.

Single cell transcriptomics of human epidermis identifies basal stem cell transition states.

Gastric cancer heterogeneity represents a barrier to disease management. We generated a comprehensive single-cell atlas of gastric cancer (>200,000 cells) comprising 48 samples from 31 patients across clinical stages and histologic subtypes. We identified 34 distinct cell-lineage states including novel rare cell populations. Many lineage states exhibited distinct cancer-associated expression profiles, individually contributing to a combined tumor-wide molecular collage. We observed increased plasma cell proportions in diffuse-type tumors associated with epithelial-resident KLF2 and stage-wise accrual of cancer-associated fibroblast subpopulations marked by high INHBA and FAP coexpression. Single-cell comparisons between patient-derived organoids (PDO) and primary tumors highlighted inter- and intralineage similarities and differences, demarcating molecular boundaries of PDOs as experimental models. We complemented these findings by spatial transcriptomics, orthogonal validation in independent bulk RNA-sequencing cohorts, and functional demonstration using in vitro and in vivo models. Our results provide a high-resolution molecular resource of intra- and interpatient lineage states across distinct gastric cancer subtypes.SignificanceWe profiled gastric malignancies at single-cell resolution and identified increased plasma cell proportions as a novel feature of diffuse-type tumors. We also uncovered distinct cancer-associated fibroblast subtypes with INHBA-FAP-high cell populations as predictors of poor clinical prognosis. Our findings highlight potential origins of deregulated cell states in the gastric tumor ecosystem. This article is highlighted in the In This Issue feature, p. 587.

Single-Cell Atlas of Lineage States, Tumor Microenvironment, and Subtype-Specific Expression Programs in Gastric Cancer

PurposeThe heterogeneity of tumor cells presents a major challenge to cancer diagnosis and therapy. Cutaneous T-cell lymphomas (CTCL) are a group of T lymphocyte malignancies that primarily affect skin. Lack of highly specific markers for malignant lymphocytes prevents early diagnosis, while only limited treatment options are available for patients with advanced stage CTCL. Droplet-based single-cell transcriptome analysis of CTCL skin biopsies opens avenues for dissecting patient-specific T lymphocyte heterogeneity, providing a basis for identifying specific markers for diagnosis and cure of CTCL.Experimental designSingle-cell RNA-sequencing was performed by Droplet-based sequencing (10X Genomics), focusing on 14,056 CD3+ lymphocytes (448 cells from normal and 13,608 cells from CTCL skin samples) from skin biopsies of 5 patients with advanced-stage CTCL and 4 healthy donors. Protein expression of identified genes was validated in advanced stage CTCL skin tumors by immunohistochemistry and confocal immunofluorescence microscopy.ResultsOur analysis revealed a large inter- and intratumor gene expression heterogeneity in the T lymphocyte subset, as well as a common gene expression signature in highly proliferating lymphocytes that was validated in multiple advanced-stage skin tumors. In addition, we established the immunologic state of reactive lymphocytes and found heterogeneity in effector and exhaustion programs across patient samples.ConclusionsSingle-cell analysis of CTCL skin tumor samples reveals patient-specific landscapes of malignant and reactive lymphocytes within the local microenvironment of each tumor, giving an unprecedented view of lymphocyte heterogeneity and identifying tumor-specific molecular signatures, with important implications for diagnosis and personalized disease treatment.

Single-Cell Lymphocyte Heterogeneity in Advanced Cutaneous T-cell Lymphoma Skin Tumors

We induced differentiation from human induced pluripotent stem cell (hiPSC)-derived alveolar type2 (AT2) cells into hiPSC-derived 2D alveoalr type 1 (2D-iAT1) cells. Single cell RNA sequencing revealed that 2D-iAT1 cells were composed of three clusters. Overall design: After isolation of hiPSC-derived AT2 cells, hiPSC-dervied AT2 cells were cultured in two-demensional condition for 7 days, and then single cell libraries were generated using Fluidigm C1 system. Single "2D_EpCAM+" cells were captured and separated on the each chamber of medium-sized (17-25 µm) C1 integrated fluidic circuit (IFC) with the C1 system (Fluidgm) Then, first-strand cDNA was synthesized and amplified on the independent chamber using SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing.

Single-Cell RNA Sequencing Reveals Human iPS Cell-derived Alveolar Type 1 Cells in Alveolar Organoids

Background In early-stage mycosis fungoides (MF), the most common primary cutaneous T-cell lymphoma, limited skin involvement with patches and plaques is associated with a favorable prognosis. Nevertheless, approximately 20–30% of cases progress to tumors or erythroderma, resulting in poor outcome. At present, factors contributing to this switch from indolent to aggressive disease are only insufficiently understood. Methods In patients with advanced-stage MF, we compared patches with longstanding history to newly developed plaques and tumors by using single-cell RNA sequencing, and compared results with early-stage MF as well as nonlesional MF and healthy control skin. Results Despite considerable inter-individual variability, lesion progression was uniformly associated with downregulation of the tissue residency markers CXCR4 and CD69, the heat shock protein HSPA1A, the tumor suppressors and immunoregulatory mediators ZFP36 and TXNIP, and the interleukin 7 receptor (IL7R) within the malignant clone, but not in benign T cells. This phenomenon was not only found in conventional TCR-αβ MF, but also in a case of TCR-γδ MF, suggesting a common mechanism across MF subtypes. Conversely, malignant cells in clinically unaffected skin from MF patients showed upregulation of these markers. Conclusions Our data reveal a specific panel of biomarkers that might be used for monitoring MF disease progression. Altered expression of these genes may underlie the switch in clinical phenotype observed in advanced-stage MF.

Single-cell RNA sequencing reveals markers of disease progression in primary cutaneous T-cell lymphoma

Human kidney organoids hold promise for studying development, disease modelling and drug screening. However, the utility of stem cell-derived kidney tissues will depend on how faithfully these replicate normal fetal development at the level of cellular identity and complexity. Here, we present an integrated analysis of single cell datasets from human kidney organoids and human fetal kidney to assess similarities and differences between the component cell types.

Single-cell analysis reveals congruence between kidney organoids and human fetal kidney

BackgroundAlthough ample knowledge exists about phenotype and function of cutaneous T lymphocytes, much less is known about the lymphocytic components of the skin's innate immune system.ObjectiveTo better understand the biologic role of cutaneous innate lymphoid cells (ILCs), we investigated their phenotypic and molecular features under physiologic (normal human skin [NHS]) and pathologic (lesional skin of patients with atopic dermatitis [AD]) conditions.MethodsSkin punch biopsies and reduction sheets as well as blood specimens were obtained from either patients with AD or healthy individuals. Cell and/or tissue samples were analyzed by flow cytometry, immunohistochemistry, and single-cell RNA sequencing or subjected to in vitro/ex vivo culture.ResultsNotwithstanding substantial quantitative differences between NHS and AD skin, we found that the vast majority of cutaneous ILCs belong to the CRTH2+ subset and reside in the upper skin layers. Single-cell RNA sequencing of cutaneous ILC-enriched cell samples confirmed the predominance of biologically heterogeneous group 2 ILCs and, for the first time, demonstrated considerable ILC lineage infidelity (coexpression of genes typical of either type 2 [GATA3 and IL13] or type 3/17 [RORC, IL22, and IL26] immunity within individual cells) in lesional AD skin, and to a much lesser extent, in NHS. Similar events were demonstrated in ILCs from skin explant cultures and in vitro expanded ILCs from the peripheral blood.ConclusionThese findings support the concept that instead of being a stable entity with well-defined components, the skin immune system consists of a network of highly flexible cellular players that are capable of adjusting their function to the needs and challenges of the environment.

Single-cell analysis reveals innate lymphoid cell lineage infidelity in atopic dermatitis

Immunotherapies have been explored in treating solid tumors, albeit with disparate clinical effects in distinct cancer types. Systematic interrogation of immune cells in the tumor microenvironment (TME) is vital to the prediction of immunotherapy response and the development of innovative immunotherapeutics. To comprehensively characterize the immune microenvironment in advanced biliary tract cancer (BTC), we utilized single-cell RNA sequencing in unselected viable cells from 16 matched samples. We identified nineteen cell subsets from a total of 45,851 cells, in which exhausted CD8+ T cells, macrophages, and dendritic cells (DCs) in BTC were shown to augment and communicate within the TME. Transcriptional profiles coupled with T cell receptor (TCR) sequences revealed that exhausted CD8+ T cells retained clonal expansion and high proliferation in the TME, and some of them highly expressed the endoplasmic reticulum stress (ER) sensor XBP1, indicating its key role in remodeling TME. Functional assays demonstrated that XBP1 and common immune checkpoints (PD1, TIGIT) were significantly upregulated in CD8+ T cells cocultured within the TME of BTC cell lines (GBC-SD, HCCC-9810). When treating the coculture groups with the specific inhibitor of IRE1α-XBP1 (4μ8C), the downregulation of TIGIT was observed in the treatment group. Collectively, comprehensive transcriptome profiling provides deeper insights into the immune atlas in advanced BTC, which might be instrumental in exploring novel immunotherapy strategies. Overall design: Five patients diagnosed with cholangiocarcinoma in the Eastern Hepatobiliary Surgery Hospital were incorporated into the current study, including two intrahepatic cholangiocarcinoma (iCCA), two gallbladder carcinoma (GBC) and one dCCA patients. Of note, none of the patients have undergone any preoperative treatments (chemotherapy, radiation or anti-tumor medicines) prior to tumor resection. Tumor primary locus and metastasis as well as lymph node were consistently harvested for each participant, if available. Meanwhile, peripheral blood samples were collected prior to surgical procedures using anticoagulant tubes.

Single-cell atlas of diverse immune populations in the advanced biliary tract cancer microenvironment

The functions of innate lymphoid cells (ILCs) in immune system are increasingly appreciated, whereas the early development of ILCs in human remains elusive. In this study, we sorted humanhematopoietic stem progenitor cells, lymphoid progenitors, presumed ILC progenitor/precursors and mature ILCs in the fetal hematopoietic, lymphoid and non-lymphoid tissues, from 8 to 12 PCW,for single-cell RNA-sequencing, followed by computational analysis and functional validation. We delineated the early phase of ILC development, from hematopoietic stem progenitor cells to multipotent lymphoid progenitors and to ILC progenitors, which mainly occurred in fetal liver and intestine. We further unveiled interleukin-3 receptor alpha (IL-3RA) as the surface marker for the lymphoid progenitors with T cell, B cell and ILC potentials. Notably, we determined the heterogeneity and tissue distribution of each ILC subpopulation, revealing the shared proliferating characteristics of the precursors of each ILC subtype. Additionally, a novel unconventional ILC2 subpopulation (CRTH2-CCR9+ ILC2) was identified in fetal thymus. Taken together, our study illuminates the precise cellular and molecular features underlying the stepwise formation of human fetal ILC hierarchy with remarkable spatiotemporal heterogeneity. Overall design: Here, we provided an unprecedented single-cell transcriptome atlas of stepwise ILC ontogeny in early human fetus using cells from fetal hematopoietic, lymphoid and non-lymphoid tissues spanning 8 to 12 post-conception weeks (PCWs).

Single-cell deciphering of human fetal innate lymphoid cell development

The immune system of skin develops in stages in mice. However, the developmental dynamics of immune cells in human skin remains elusive. Here we perform transcriptome profiling of CD45+ hematopoietic cells in human fetal skin at an estimated gestational age of 10–17 weeks by single-cell RNA sequencing. A total of 13 immune cell types are identified. Skin macrophages show heterogeneity and dynamic over the course of skin development. A major shift in lymphoid cell developmental states occurs from the first to the second trimester that implied an in situ differentiation process. Gene expression analysis reveals a typical developmental program in immune cells in accordance with their functional maturation, possibly involving metabolic reprogramming. Finally, we identify transcription factors (TFs) that potentially regulate cellular transitions by comparing TFs and TF-target gene networks. These findings provide detailed insight into how the immune system of the human skin is established during development Overall design: Here we perform transcriptome profiling of CD45+ hematopoietic cells in human fetal skin at an estimated gestational age of 10–17 weeks by single-cell RNA sequencing

Single-cell transcriptome analysis reveals the dynamics of human immune cells during early fetal skin development

Background: As core units of organ tissues, cells of various types play their harmonious rhythms to maintain the homeostasis of the human body. It is essential to identify the characteristics of cells in human organs and their regulatory networks for understanding the biological mechanisms related to health and disease. However, a systematic and comprehensive single-cell transcriptional profile across multiple organs of a normal human adult is missing. Results: We perform single-cell transcriptomes of 84,363 cells derived from 15 tissue organs of one adult donor and generate an adult human cell atlas. The adult human cell atlas depicts 252 subtypes of cells, including major cell types such as T, B, myeloid, epithelial, and stromal cells, as well as novel COCH+ fibroblasts and FibSmo cells, each of which is distinguished by multiple marker genes and transcriptional profiles. These collectively contribute to the heterogeneity of major human organs. Moreover, T cell and B cell receptor repertoire comparisons and trajectory analyses reveal direct clonal sharing of T and B cells with various developmental states among different tissues. Furthermore, novel cell markers, transcription factors and ligand-receptor pairs are identified with potential functional regulations in maintaining the homeostasis of human cells among tissues. Conclusions: The adult human cell atlas reveals the inter- and intra-organ heterogeneity of cell characteristics and provides a useful resource in uncovering key events during the development of human diseases in the context of the heterogeneity of cells and organs. Overall design: We performed single-cell transcriptomes of 84,363 cells derived from 15 tissue organs of one adult donor and generated an adult human cell atlas (AHCA). Keywords: Human cell atlas, Single-cell RNA sequencing, TCR, BCR, Transcriptome

Single-cell transcriptome profiling of an adult human cell atlas of 15 major organs

Fibroblasts are the main dermal cell type and are essential for the architecture and function of human skin. Important differences have been described between fibroblasts localized in distinct dermal layers, and these cells are also known to perform varied functions. However, this phenomenon has not been analyzed comprehensively yet. Here we have used single-cell RNA sequencing to analyze >15,000 cells from a sun-protected area in young and old donors. Our results define four main fibroblast subpopulations that can be spatially localized and functionally distinguished. Importantly, intrinsic aging reduces this fibroblast ‘priming’, generates distinct expression patterns of skin aging-associated genes, and substantially reduces the interactions of dermal fibroblasts with other skin cell types. Our work thus provides comprehensive evidence for a functional specialization of human dermal fibroblasts and suggests that the age-related loss of fibroblast priming contributes to human skin aging. Overall design: Overall 5 samples were analyzed. Two of them are from young subjects (age=25) and serve as replicates, three of them are from old subjects (age >= 50) and serve as replicates.

Single-cell transcriptomes of the human skin reveal age-related loss of fibroblast priming

Neuroblastoma is a pediatric tumor of the developing sympathetic nervous system. However, the cellular origin of neuroblastoma has yet to be defined. Here we studied the single-cell transcriptomes of neuroblastomas and normal human developing adrenal glands at various stages of embryonic and fetal development. We defined normal differentiation trajectories from Schwann cell precursors over intermediate states to neuroblasts or chromaffin cells and showed that neuroblastomas transcriptionally resemble normal fetal adrenal neuroblasts. Importantly, neuroblastomas with varying clinical phenotypes matched different temporal states along normal neuroblast differentiation trajectories, with the degree of differentiation corresponding to clinical prognosis. Our work highlights the roles of oncogenic MYCN and loss of TFAP2B in blocking differentiation and may provide the basis for designing therapeutic interventions to overcome differentiation blocks.

Single-cell transcriptomic analyses provide insights into the developmental origins of neuroblastoma.

Background
Atopic dermatitis (AD) is the most common chronic inflammatory skin disease, but its complex pathogenesis is only insufficiently understood, resulting in still limited treatment options.ObjectiveWe sought to characterize AD on both transcriptomic and proteomic levels in humans.Methods
We used skin suction blistering, a painless and nonscarring procedure that can simultaneously sample skin cells and interstitial fluid. We then compared results with conventional biopsies.Results
Suction blistering captured epidermal and most immune cells equally well as biopsies, except for mast cells and nonmigratory CD163+ macrophages that were only present in biopsy isolates. Using single-cell RNA sequencing, we found comparable transcriptional profiles of key inflammatory pathways between blister and biopsy AD, but suction blistering was superior in cell-specific resolution for high-abundance transcripts (KRT1/KRT10, KRT16/KRT6A, S100A8/S100A9), which showed some background signals in biopsy isolates. Compared with healthy controls, we found characteristic upregulation of AD-typical cytokines such as IL13 and IL22 in Th2 and Th22 cells, respectively, but we also discovered these mediators in proliferating T cells and natural killer T cells, that also expressed the antimicrobial cytokine IL26. Overall, not T cells, but myeloid cells were most strongly enriched in AD, and we found dendritic cell (CLEC7A, amphiregulin/AREG, EREG) and macrophage products (CCL13) among the top upregulated proteins in AD blister fluid proteomic analyses.ConclusionThese data show that by using cutting-edge technology, suction blistering offers several advantages over conventional biopsies, including better transcriptomic resolution of skin cells, combined with proteomic information from interstitial fluid, unraveling novel inflammatory players that shape the cellular and proteomic microenvironment of AD.

Single-cell transcriptomics combined with interstitial fluid proteomics defines cell type-specific immune regulation in atopic dermatitis.

Cutaneous lupus is commonly present in patients with systemic lupus erythematosus (SLE). T cells have been strongly suspected to contribute to the pathology of cutaneous lupus; however, our understanding of the relevant T cell phenotypes and functions remains incomplete. Here, we present a detailed single-cell RNA-Seq profile of T and NK cell populations present within lesional and nonlesional skin biopsies of patients with cutaneous lupus. T cells across clusters from lesional and nonlesional skin biopsies expressed elevated levels of IFN-simulated genes (ISGs). Compared with T cells from control skin, however, T cells from cutaneous lupus lesions did not show elevated expression profiles of activation, cytotoxicity, or exhaustion. Integrated analyses indicated that skin lymphocytes appeared less activated and lacked the expanded cytotoxic populations prominent in lupus nephritis kidney T/NK cells. Comparison of skin T cells from lupus and systemic sclerosis skin biopsies further revealed an elevated ISG signature specific to cells from lupus biopsies. Overall, these data represent the first detailed transcriptomic analysis to our knowledge of the T and NK cells in cutaneous lupus at the single-cell level and have enabled a cross-tissue comparison that highlights stark differences in composition and activation of T/NK cells in distinct tissues in lupus.

Single-cell transcriptomics reveals distinct effector profiles of infiltrating T cells in lupus skin and kidney

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Inflammatory skin diseases including atopic dermatitis (AD) and psoriasis (PSO) are underpinned by dendritic cell (DC)-mediated T cell responses. Currently, the heterogeneous human cutaneous DC population is incompletely characterized, and its contribution to these diseases remains unclear. Here, we performed index-sorted single-cell flow cytometry and RNA sequencing of lesional and nonlesional AD and PSO skin to identify macrophages and all DC subsets, including the newly described mature LAMP3+BIRC3+ DCs enriched in immunoregulatory molecules (mregDC) and CD14+ DC3. By integrating our indexed data with published skin datasets, we generated a myeloid cell universe of DC and macrophage subsets in healthy and diseased skin. Importantly, we found that CD14+ DC3s increased in PSO lesional skin and co-produced IL1B and IL23A, which are pathological in PSO. Our study comprehensively describes the molecular characteristics of macrophages and DC subsets in AD and PSO at single-cell resolution, and identifies CD14+ DC3s as potential promoters of inflammation in PSO.

Skin single cell universe identifies IL-1B/IL23A co-producing CD14+ type 3 dendritic cells in psoriasis

Fibroblasts synthesize the extracellular matrix of connective tissue and play an essential role in maintaining tissue integrity. We have previously shown that mouse skin connective tissue, the dermis, is comprised of functionally distinct fibroblast lineages. However, the extent of fibroblast heterogeneity in human skin is unknown. Here, using a combination of spatial transcriptional profiling of human and mouse dermis and single cell transcriptional profiling of human dermal fibroblasts, we show that there are at least four distinct fibroblast populations in adult human skin. We define markers permitting prospective isolation of these cells and show that although marker expression is rapidly lost in culture, different fibroblast subpopulations retain distinct functionality in terms of Wnt signalling, T cell communication and the ability to support human epidermal reconstitution in organotypic culture. Furthermore, while some fibroblast subpopulations are spatially segregated, others are not. These findings have profound implications for normal wound healing and diseases characterized by excessive fibrosis, and suggest that ex vivo expansion or in vivo ablation of specific fibroblast subpopulations may have therapeutic applications.

Spatial and single-cell transcriptional profiling identifies functionally distinct human dermal fibroblast subpopulations

We have created a compendium of single cell transcriptome data from the model
organism Mus musculus comprising more than 100,000 cells from 20 organs and
tissues. These data represent a new resource for cell biology, revealing gene
expression in poorly characterized cell populations and allowing for direct and
controlled comparison of gene expression in cell types shared between tissues, such
as T-lymphocytes and endothelial cells from distinct anatomical locations. Two
distinct technical approaches were used for most tissues: one approach, microfluidic
droplet-based 3’-end counting, enabled the survey of thousands of cells at relatively
low coverage, while the other, FACS-based full length transcript analysis, enabled
characterization of cell types with high sensitivity and coverage. The cumulative
data provide the foundation for an atlas of transcriptomic cell biology.

Tabula Muris: Transcriptomic characterization of 20 organs and tissues from Mus musculus at single cell resolution

Cardiac fibrosis is causally linked to heart failure pathogenesis and adverse clinical outcomes. However, the precise fibroblast populations that drive fibrosis in the human heart and the mechanisms that govern their emergence remain incompletely defined. Here, we performed Cellular Indexing of Transcriptomes and Epitomes by sequencing (CITE-seq) in 22 explanted human hearts from healthy donors, acute myocardial infarction (MI), and chronic ischemic and non-ischemic cardiomyopathy patients. We identified a fibroblast trajectory marked by fibroblast activator protein (FAP) and periostin (POSTN) expression that was independent of myofibroblasts, peaked early after MI, remained elevated in chronic heart failure, and displayed a transcriptional signature consistent with fibrotic activity. We assessed the applicability of cardiac fibrosis models and demonstrated that mouse MI, angiotensin II/phenylephrine infusion, and pressure overload models were superior compared to cultured human heart and dermal fibroblasts in recapitulating cardiac fibroblast diversity including pathogenic cell states. Ligand-receptor analysis and spatial transcriptomics predicted interactions between macrophages, T cells, and fibroblasts within spatially defined niches. CCR2 + monocyte and macrophage states were the dominant source of ligands targeting fibroblasts. Inhibition of IL-1β signaling to cardiac fibroblasts was sufficient to suppress fibrosis, emergence, and maturation of FAP + POSTN + fibroblasts. Herein, we identify a human fibroblast trajectory marked by FAP and POSTN expression that is associated with cardiac fibrosis and identify macrophage-fibroblast crosstalk mediated by IL-1β signaling as a key regulator of pathologic fibroblast differentiation and fibrosis.

Targeting the Immune-Fibrosis Axis in Myocardial Infarction and Heart Failure

In recent years there has been tremendous progress towards deep molecular characterization of cell types using single cell transcriptome sequencing. Here we report a single cell transcriptomic atlas comprising nearly 500,000 cells from 24 different human tissues and organs. In several instances multiple organs were analyzed from the same donor. Analyzing organs from the same individual controls for genetic background, age, environment, and epigenetic effects, and enables a detailed comparison of cell types that are shared between tissues. This resource provides a rich molecular characterization of more than 400 cell types, their distribution across tissues, and detailed information about tissue specific variation in gene expression. We have used the fact that multiple tissues came from the same donor to study the clonal distribution of T cells between tissues, to understand the tissue specific mutation rate in B cells, and to analyze the cell cycle state and proliferative potential of shared cell types across tissues. Finally, we have also used this data to characterize cell type specific RNA splicing and how such splicing varies across tissues within an individual.

The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans

To delineate the ontogeny of the mammalian endoderm, we generated 112,217 single-cell transcriptomes representing all endoderm populations within the mouse embryo until midgestation. By employing graph-based approaches, we modelled differentiating cells for spatio-temporal characterization of developmental trajectories and defined the transcriptional architecture that accompanies the emergence of the first (primitive or extra-embryonic) endodermal population and its sister pluripotent (embryonic) epiblast lineage. We uncovered a relationship between descendants of these two lineages, whereby epiblast cells differentiate into endoderm at two distinct time-points, before and during gastrulation. Trajectories of endoderm cells were mapped as they acquired embryonic versus extra-embryonic fates, and as they spatially converged within the nascent gut endoderm; revealing them to be globally similar but retaining aspects of their lineage history. We observed the regionalized identity of cells along the anterior-posterior axis of the emergent gut tube, reflecting their embryonic or extra-embryonic origin, and their coordinate patterning into organ-specific territories.

The emergent landscape of the mouse gut endoderm at single-cell resolution

Our current knowledge of the different immune cells in neuroblastoma is based on in vitro and in vivo studies mainly focusing on a single cell type. Importantly, different studies have conveyed conflicting results. Moreover, a comprehensive immune cell overview at the single-cell level is still missing and understanding the complete immune cell composition of neuroblastoma will be crucial for the development of novel immunotherapeutics against the disease. In this study, we performed single-cell RNA-sequencing on nineteen human neuroblastoma samples coupled with multiplex immunohistochemistry and survival analysis using additional datasets to provide a comprehensive cellular and molecular immune cell landscape of human neuroblastoma. Further, we contrasted our data with single-cell RNA-sequencing data from normal fetal adrenal gland to characterize cell-state changes from normal tissue to cancerous neuroblastoma. Our analysis revealed 27 immune cell subtypes including distinct subpopulations of myeloid, NK, B and T cells not identified in neuroblastoma before. Several immune cell subtypes demonstrated a survival benefit such as inflammatory monocytes, tumor associated macrophages, various T cell populations, and Active NK cells. Furthermore, in contrast to adult cancers and previous neuroblastoma studies, we demonstrated an increase in inflammatory monocyte cell-state when contrasting normal and tumor tissue, while we do not observe differences in cytotoxicity and exhaustion score for cytotoxic T cells, nor in Treg activity. Finally, we performed a systemic receptor-ligand interaction analysis between tumor, stroma and immune cells, where we showed the neuroblastoma tumor microenvironment is highly complex and strongly correlated to survival. In addition, we highlighted several interactions that we suggest to be tested in future studies as a therapeutic option in human neuroblastoma. Taken together, our study significantly adds to the in depth understanding of the immune cell landscape, the complexity of the tumor microenvironment and it provides a resource for the development of novel immunotherapeutics for neuroblastoma.

The immune cell atlas of human neuroblastoma

Rhabdomyosarcoma (RMS) is the most soft tissue sarcoma of childhood and histologically mimics developing skeletal muscle. To better understand RMS, we performed single-cell RNA-seq of one embryonal subtype RMS (ERMS) and one alveolar subtype RMS (ARMS). We expanded our observations by performing single-nucleus RNA-seq of 18 patient tumors (12 ERMS, 6 ARMS). All 18 of these tumors were used to generate orthotopic patient-derived xenografts, and we performed single-cell RNA-seq of all 18 O-PDXs to compare heterogeneity between the original tumor and the patient-derived model. 15 of the O-PDXs underwent lentiviral barcoding, and lineage could be tracked using a separate dial-out PCR step during the single-cell RNA-seq workflow. We performed single-cell RNA-seq of an ex vivo organoid, showing that we can preserve heterogeneity in this model. We expanded our study to the epigenetic level by using single-cell ATAC-sequencing of 7 O-PDXs. Finally, we longitudinally tracked shifts in heterogeneity in an O-PDX model treated with chemotherapy. Overall design: We performed single-cell RNA-seq from 2 patient RMS tumors, 18 O-PDXs, and 1 organoid. Additionally, we performed single-nucleus RNA-seq from 18 patient RMS tumors and single-cell ATAC-seq from 7 O-PDXs. We also provide dialout PCR data from 15 RMS O-PDX that were used for lentiviral barcode experiments to track clonal lineages within O-PDXs.

The myogenesis program drives clonal selection and drug resistance in rhabdomyosarcoma

The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and under different conditions. Human induced pluripotent stem cells (hiPSCs) are uniquely suited to study these context-dependent effects, but to do so requires cell lines from hundreds or potentially thousands of individuals. Village cultures, where multiple hiPSC lines are cultured and differentiated together in a single dish, provide an elegant solution for scaling hiPSC experiments to the necessary sample sizes required for population-scale studies. Here, we show the utility of village models, demonstrating how cells can be assigned back to a donor line using single cell sequencing, and addressing whether line-specific signaling alters the transcriptional profiles of companion lines in a village culture. We generated single cell RNA sequence data from hiPSC lines cultured independently (uni-culture) and in villages at three independent sites. We show that the transcriptional profiles of hiPSC lines are highly consistent between uni- and village cultures for both fresh (0.46 < R < 0.88) and cryopreserved samples (0.46 < R < 0.62). Using a mixed linear model framework, we estimate that the proportion of transcriptional variation across cells is predominantly due to donor effects, with minimal evidence of variation due to culturing in a village system. We demonstrate that the genetic, epigenetic or hiPSC line-specific effects on gene expression are consistent whether the lines are uni- or village-cultured (0.82 < R < 0.94). Finally, we identify the consistency in the landscape of cell states between uni- and village-culture systems. Collectively, we demonstrate that village methods can be effectively used to detect hiPSC line-specific effects including sensitive dynamics of cell states.

Village in a dish: a model system for population-scale hiPSC studies

Regenerative medicine relies on basic research to find safe and useful outcomes that are only practical when cost-effective. The human eyeball requires the retinal pigment epithelium (RPE) for support and maintenance that interfaces the neural retina and the choroid at large. Nearly 200 million people suffer from age-related macular degeneration (AMD), a blinding multifactor genetic disease among other retinal pathologies related to RPE degradation. Recently, autologous pluripotent stem cell-derived RPE cells were prohibitively expensive due to time, therefore we developed a new simplified cell reprogramming system. We stably induced RPE-like cells (iRPE) from human fibroblasts by conditional overexpression of broad plasticity and lineage-specific pioneering transcription factors and treated some with nicotinamide and chetomin (iRPENC). iRPE cells showed features of modern RPE benchmarks and significant in-vivo integration in transplanted chimeric hosts. Herein, we provide the single-cell RNA (scRNA) sequencing profiling of the foreskin fibroblasts, iRPE, iRPENC, and model iPSC.RPE cells. Overall design: Cell culture samples were prepared as approximately 3,000 single cells and then sampled for scRNA-seq 10x Genomics 5' (v1.0 Chemistry) VDJ kit library construction with standard protocol for 5' Gene Expressoin libraries that skips Steps 4&5. Libraries were sequenced with paired-end reads on Illumina HiSeq X at one library per lane and then mapped to a custom reference human genome (hg38) using Cell Ranger (Version 3.1). related study: https://doi.org/10.1101/2020.07.27.215103

scRNA-seq experiments on human foreskin fibroblasts, foreskin fibroblasts induced directly to retina pigment epithelium, and induced pluripotent stem cell derived retina pigment epithelium.

T cells in human skin play an important role in the immune defense against pathogens and tumors. T cells are present already in fetal skin, where little is known about their cellular phenotype and biological function. Using single-cell analyses, we identified a naive T cell population expressing αβ and γδ T cell receptors (TCRs) that was enriched in fetal skin and intestine but not detected in other fetal organs and peripheral blood. TCR sequencing data revealed that double-positive (DP) αβγδ T cells displayed little overlap of CDR3 sequences with single-positive αβ T cells. Gene signatures, cytokine profiles and in silico receptor-ligand interaction studies indicate their contribution to early skin development. DP αβγδ T cells were phosphoantigen responsive, suggesting their participation in the protection of the fetus against pathogens in intrauterine infections. Together, our analyses unveil a unique cutaneous T cell type within the native skin microenvironment and point to fundamental differences in the immune surveillance between fetal and adult human skin.

Skin Network

Network Description

Skin Network atlases are still under development.

Network Coordinators