A prenatal skin atlas reveals immune regulation of human skin morphogenesis.

Human prenatal skin is populated by innate immune cells, including macrophages, but whether they act solely in immunity or have additional functions in morphogenesis is unclear. Here we assembled a comprehensive multi-omics reference atlas of prenatal human skin (7-17 post-conception weeks), combining single-cell and spatial transcriptomics data, to characterize the microanatomical tissue niches of the skin. This atlas revealed that crosstalk between non-immune and immune cells underpins the formation of hair follicles, is implicated in scarless wound healing and is crucial for skin angiogenesis. We systematically compared a hair-bearing skin organoid (SkO) model derived from human embryonic stem cells and induced pluripotent stem cells to prenatal and adult skin1. The SkO model closely recapitulated in vivo skin epidermal and dermal cell types during hair follicle development and expression of genes implicated in the pathogenesis of genetic hair and skin disorders. However, the SkO model lacked immune cells and had markedly reduced endothelial cell heterogeneity and quantity. Our in vivo prenatal skin cell atlas indicated that macrophages and macrophage-derived growth factors have a role in driving endothelial development. Indeed, vascular network remodelling was enhanced following transfer of autologous macrophages derived from induced pluripotent stem cells into SkO cultures. Innate immune cells are therefore key players in skin morphogenesis beyond their conventional role in immunity, a function they achieve through crosstalk with non-immune cells.

Neutral evolution of snoRNA Host Gene long non-coding RNA affects cell fate control.

A fundamental challenge in molecular biology is to understand how evolving genomes can acquire new functions. Actively transcribed, non-coding parts of the genome provide a potential platform for the development of new functional sequences, but their biological and evolutionary roles remain largely unexplored. Here, we show that a set of neutrally evolving long non-coding RNAs (lncRNAs) whose introns encode small nucleolar RNAs (snoRNA Host Genes, SNHGs) are highly expressed in skin and dysregulated in inflammatory conditions. Using SNHG7 and human epidermal keratinocytes as a model, we describe a mechanism by which these lncRNAs can increase self-renewal and inhibit differentiation. The activity of SNHG7 lncRNA has been recently acquired in the primate lineage and depends on a short sequence required for microRNA binding. Taken together, our results highlight the importance of understanding the role of fast-evolving transcripts in normal and diseased epithelia, and show how poorly conserved, actively transcribed non-coding sequences can participate in the evolution of genomic functionality.

Recommendations for accelerating open preprint peer review to improve the culture of science.

Peer review is an important part of the scientific process, but traditional peer review at journals is coming under increased scrutiny for its inefficiency and lack of transparency. As preprints become more widely used and accepted, they raise the possibility of rethinking the peer-review process. Preprints are enabling new forms of peer review that have the potential to be more thorough, inclusive, and collegial than traditional journal peer review, and to thus fundamentally shift the culture of peer review toward constructive collaboration. In this Consensus View, we make a call to action to stakeholders in the community to accelerate the growing momentum of preprint sharing and provide recommendations to empower researchers to provide open and constructive peer review for preprints.

Multiscale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma.

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this, we have created a multiscale 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.

Genome-wide meta-analysis implicates variation affecting mast cell biology in urticaria.

BACKGROUND: Urticaria is characterized by inappropriate mast cell degranulation leading to the development of wheals and/or angioedema. Twin and family studies indicate that there is a substantial heritable component to urticaria risk. OBJECTIVE: Our aim was to identify genomic loci at which common genetic variation influences urticaria susceptibility. METHODS: Genome-wide association studies of urticaria (including all subtypes) from 3 European cohorts (UK Biobank, FinnGen, and the Trøndelag Health Study [HUNT]) were combined through statistical meta-analysis (14,306 urticaria cases and 650,664 controls). Cases were identified via electronic health care records from primary and/or secondary care. To identify putative causal variants and genes, statistical fine-mapping, colocalization, and interrogation of publicly available single-cell transcriptome sequencing resources were performed. RESULTS: Genome-wide significant associations (P < 5 × 10-8) were identified at 6 independent loci. These included 2 previously reported association signals at 1q44 and the human leucocyte antigen region on chromosome 6. Genes with expected or established roles in mast cell biology were associated with the 4 other genome-wide association signals (GCSAML, FCER1A, TPSAB1, and CBLB). Colocalization of association signals consistent with the presence of shared causal variants was observed between urticaria susceptibility and increased expression of GCSAML (posterior probability of colocalization [PPcoloc] = 0.89) and FCER1A (PPcoloc = 0.91) in skin. CONCLUSION: Common genetic variation influencing the risk of developing urticaria was identified at 6 genomic loci. The relationship between genes with roles in mast cell biology and several association signals implicates genetic variability of specific components of mast cell function in the development of urticaria.

Gradient boosting approaches can outperform logistic regression for risk prediction in cutaneous allergy.

BACKGROUND: Contact allergy is a major clinical and public health challenge. It is important to identify individuals who are at risk and perform patch testing to identify relevant allergens. Predicting clinical risk on the basis of input parameters is common in clinical medicine and traditionally has been achieved with linear models. OBJECTIVES: We hypothesized that the risk of a clinically relevant positive patch test could be predicted according to clinical and demographic parameters. METHODS: We compared the predictive accuracy of logistic regression with more sophisticated machine learning approaches such as gradient boosting, in the prediction of patch testing results. RESULTS: We found that both logistic regression and more sophisticated machine learning approaches were able to predict the risk of positive patch tests. For certain predictions, including the overall risk of a clinically relevant positive patch test, gradient boosting approaches can outperform logistic regression. CONCLUSIONS: These findings suggest that complex nonlinear interactions between input variables are relevant in risk prediction. While a risk prediction model cannot replace the judgment of an experienced clinician, quantifying the risk of a clinically relevant positive patch test result has the potential to assist in decision making and to inform discussions with patients.

Angiotensin-Converting Enzyme 2 Expression Is Detectable in Keratinocytes, Cutaneous Appendages, and Blood Vessels by Multiplex RNA In Situ Hybridization.

OBJECTIVE: The angiotensin-converting enzyme 2 (ACE2) receptor mediates uptake of SARS-CoV-2, the virus responsible for COVID-19. Previous work analyzing publicly available bulk RNA-sequencing data sets has shown the expression of ACE2 in human keratinocytes. This finding is potentially relevant for the etiology of COVID-19-associated rashes and might also suggest a possible entry mechanism for the SARS-CoV-2 virus. In this study, the authors examined the spatial localization of ACE2 mRNA in vivo. METHODS AND RESULTS: The authors analyzed several publicly available single-cell RNA-sequencing data sets. They determined spatial localization of ACE2 mRNA using multiplex RNA in situ hybridization in human skin. CONCLUSIONS: Both analyses supported ACE2 expression in keratinocytes and skin vasculature, which could reflect a potential cutaneous entry point for SARS-CoV-2, particularly in damaged or broken skin. Moreover, ACE2 expression in vascular endothelial cells may support direct, virally mediated mechanisms in the etiology of the chilblain-like acral eruption that is seen in patients with COVID-19.

Mapping epidermal and dermal cellular senescence in human skin aging.

Single-cell RNA sequencing and spatial transcriptomics enable unprecedented insight into cellular and molecular pathways implicated in human skin aging and regeneration. Senescent cells are individual cells that are irreversibly cell cycle arrested and can accumulate across the human lifespan due to cell-intrinsic and -extrinsic stressors. With an atlas of single-cell RNA-sequencing and spatial transcriptomics, epidermal and dermal senescence and its effects were investigated, with a focus on melanocytes and fibroblasts. Photoaging due to ultraviolet light exposure was associated with higher burdens of senescent cells, a sign of biological aging, compared to chronological aging. A skin-specific cellular senescence gene set, termed SenSkin™, was curated and confirmed to be elevated in the context of photoaging, chronological aging, and non-replicating CDKN1A+ (p21) cells. In the epidermis, senescent melanocytes were associated with elevated melanin synthesis, suggesting haphazard pigmentation, while in the dermis, senescent reticular dermal fibroblasts were associated with decreased collagen and elastic fiber synthesis. Spatial analysis revealed the tendency for senescent cells to cluster, particularly in photoaged skin. This work proposes a strategy for characterizing age-related skin dysfunction through the lens of cellular senescence and suggests a role for senescent epidermal cells (i.e., melanocytes) and senescent dermal cells (i.e., reticular dermal fibroblasts) in age-related skin sequelae.

Single-cell analysis of psoriasis resolution demonstrates an inflammatory fibroblast state targeted by IL-23 blockade.

Biologic therapies targeting the IL-23/IL-17 axis have transformed the treatment of psoriasis. However, the early mechanisms of action of these drugs remain poorly understood. Here, we perform longitudinal single-cell RNA-sequencing in affected individuals receiving IL-23 inhibitor therapy. By profiling skin at baseline, day 3 and day 14 of treatment, we demonstrate that IL-23 blockade causes marked gene expression shifts, with fibroblast and myeloid populations displaying the most extensive changes at day 3. We also identify a transient WNT5A+/IL24+ fibroblast state, which is only detectable in lesional skin. In-silico and in-vitro studies indicate that signals stemming from these WNT5A+/IL24+ fibroblasts upregulate multiple inflammatory genes in keratinocytes. Importantly, the abundance of WNT5A+/IL24+ fibroblasts is significantly reduced after treatment. This observation is validated in-silico, by deconvolution of multiple transcriptomic datasets, and experimentally, by RNA in-situ hybridization. These findings demonstrate that the evolution of inflammatory fibroblast states is a key feature of resolving psoriasis skin.

Single-cell RNA sequencing of human epidermis identifies Lunatic fringe as a novel regulator of the stem cell compartment.

Single-cell RNA sequencing (scRNA-seq) of human skin provides a tool for validating observations from in vitro experimental models. By analyzing a published dataset of healthy adult epidermis, we confirm that the basal epidermal layer is heterogeneous, and three subpopulations of non-dividing cells can be distinguished. We show that Delta-like ligand 1 (DLL1) is expressed in a subset of basal cells previously identified as stem cells in cultured human keratinocytes and map the distribution of other Notch ligands and receptors to specific epidermal cell compartments. Although DLL1 is expressed at low levels, it is expressed in the same cell state as the Notch regulator, Lunatic -fringe (LFNG, O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase). Overexpression of LFNG amplifies the effects of DLL1 in cultured keratinocytes, increasing proliferation and colony-forming ability. We conclude that using scRNA-seq resources from healthy human skin not only validates previous experimental data but allows formulation of testable new hypotheses.

Myc-dependent dedifferentiation of Gata6+ epidermal cells resembles reversal of terminal differentiation.

Dedifferentiation is the process by which terminally differentiated cells acquire the properties of stem cells. During mouse skin wound healing, the differentiated Gata6-lineage positive cells of the sebaceous duct are able to dedifferentiate. Here we have integrated lineage tracing and single-cell mRNA sequencing to uncover the underlying mechanism. Gata6-lineage positive and negative epidermal stem cells in wounds are transcriptionally indistinguishable. Furthermore, in contrast to reprogramming of induced pluripotent stem cells, the same genes are expressed in the epidermal dedifferentiation and differentiation trajectories, indicating that dedifferentiation does not involve adoption of a new cell state. We demonstrate that dedifferentiation is not only induced by wounding, but also by retinoic acid treatment or mechanical expansion of the epidermis. In all three cases, dedifferentiation is dependent on the master transcription factor c-Myc. Mechanotransduction and actin-cytoskeleton remodelling are key features of dedifferentiation. Our study elucidates the molecular basis of epidermal dedifferentiation, which may be generally applicable to adult tissues.

Fibroblast Heterogeneity in Healthy and Wounded Skin.

Fibroblasts are the main cell type in the dermis. They are responsible for the synthesis and deposition of structural proteins such as collagen and elastin, which are integrated into the extracellular matrix (ECM). Mouse and human studies using flow cytometry, cell culture, skin reconstitution, and lineage tracing experiments have shown the existence of different subpopulations of fibroblasts, including papillary fibroblasts, reticular fibroblasts, and fibroblasts comprising the dermal papilla at the base of the hair follicle. In recent years, the technological advances in single-cell sequencing have allowed researchers to study the repertoire of cells present in full-thickness skin including the dermis. Multiple groups have confirmed that distinct fibroblast populations can be identified in mouse and human dermis on the basis of differences in the transcriptional profile. Here, we discuss the current state of knowledge regarding dermal fibroblast heterogeneity in healthy mouse and human skin, highlighting the similarities and differences between mouse and human fibroblast subpopulations. We also discuss how fibroblast heterogeneity may provide insights into physiological wound healing and its dysfunction in pathological states such as hypertrophic and keloid scars.

Ex Vivo COL7A1 Correction for Recessive Dystrophic Epidermolysis Bullosa Using CRISPR/Cas9 and Homology-Directed Repair.

Recessive dystrophic epidermolysis bullosa is a rare and severe genetic skin disease resulting in blistering of the skin and mucosa. Recessive dystrophic epidermolysis bullosa (RDEB) is caused by a wide variety of mutations in COL7A1-encoding type VII collagen, which is essential for dermal-epidermal adhesion. Here we demonstrate the feasibility of ex vivo COL7A1 editing in primary RDEB cells and in grafted 3D skin equivalents through CRISPR/Cas9-mediated homology-directed repair. We designed five guide RNAs to correct a RDEB causative null mutation in exon 2 (c.189delG; p.Leu64Trpfs*40). Among the site-specific guide RNAs tested, one showed significant cleavage activity in primary RDEB keratinocytes and in fibroblasts when delivered as integration-deficient lentivirus. Genetic correction was detected in transduced keratinocytes and fibroblasts by allele-specific highly sensitive TaqMan-droplet digital PCR (ddPCR), resulting in 11% and 15.7% of corrected COL7A1 mRNA expression, respectively, without antibiotic selection. Grafting of genetically corrected 3D skin equivalents onto nude mice showed up to 26% re-expression and normal localization of type VII collagen as well as anchoring fibril formation at the dermal-epidermal junction. Our study provides evidence that precise genome editing in primary RDEB cells is a relevant strategy to genetically correct COL7A1 mutations for the development of future ex vivo clinical applications.

Gene-Corrected Fibroblast Therapy for Recessive Dystrophic Epidermolysis Bullosa using a Self-Inactivating COL7A1 Retroviral Vector.

Patients with recessive dystrophic epidermolysis bullosa (RDEB) lack type VII collagen and therefore have severely impaired dermal-epidermal stability causing recurrent skin and mucosal blistering. There is currently no specific approved treatment for RDEB. We present preclinical data showing that intradermal injections of genetically corrected patient-derived RDEB fibroblasts using a Good Manufacturing Practices grade self-inactivating COL7A1 retroviral vector reverse the disease phenotype in a xenograft model in nude mice. We obtained 50% transduction efficiency in primary human RDEB fibroblasts with an average low copy number (range = 1-2) of integrated provirus. Transduced fibroblasts showed strong type VII collagen re-expression, improved adhesion properties, normal proliferative capabilities, and viability in vitro. We show that a single intradermal injection of 3 × 10(6) genetically corrected RDEB fibroblasts beneath RDEB skin equivalents grafted onto mice allows type VII collagen deposition, anchoring fibril formation at the dermal-epidermal junction, and improved dermal-epidermal adherence 2 months after treatment, supporting functional correction in vivo. Gene-corrected fibroblasts previously showed no tumorigenicity. These data show the efficacy and safety of gene-corrected fibroblast therapy using a self-inactivating vector that has now been good manufacturing grade-certified and pave the way for clinical translation to treat nonhealing wounds in RDEB patients.