Single-cell transcriptome profiling reveals vascular endothelial cell heterogeneity in human skin

Abnormal dermal microvascular endothelial cells in psoriatic excessive angiogenesis

Psoriasis is characterized by excessive angiogenesis, with increased distortion and dilation of the dermal blood vessels. These vascular alterations are ascribed, at least in part, to the changes in dermal microvascular endothelial cell functions. However, despite the recognition of vascular normalization as an emerging strategy for the treatment of psoriasis, in-depth studies of human dermal microvascular endothelial cells (HDMECs) have been missing. The difficulty of isolation and culture of HDMECs has impeded the study of endothelial dysfunction in psoriasis. Researchers have done a great deal of work to study the abnormal characteristics of keratinocytes, fibroblasts, and leukocytes in psoriatic skin tissue. Recently, with successful isolation of HDMECs from psoriasis, great progress has been made in the elucidation of the pathogenic role of these cells in psoriasis. It is of great therapeutic significance to study the molecular mechanism of HDMECs in psoriasis. We review here the abnormalities of HDMECs in psoriasis.

Integrative transcriptome analysis reveals the molecular events underlying impaired T-cell responses in EGFR-mutant lung cancer

EGFR mutations are critical oncogenic drivers in lung adenocarcinoma (LUAD). However, the mechanisms by which they impact the tumor microenvironment (TME) and tumor immunity are unclear. Furthermore, the reasons underlying the poor response of EGFR-mutant (EGFR-MU) LUADs to immunotherapy with PD-1/PD-L1 inhibitors are unknown. Utilizing single-cell RNA (sc-RNA) and bulk RNA sequencing datasets, we conducted high-dimensional weighted gene coexpression network analysis to identify key genes and immune-related pathways contributing to the immunosuppressive TME. EGFR-MU cancer cells downregulated MHC class I genes to evade CD8+ cytotoxic T cells, expressed substantial levels of MHC class II molecules, and engaged with CD4+ regulatory T cells (Tregs). EGFR-MU tumors may recruit Tregs primarily through the CCL17/CCL22/CCR4 axis, leading to a Treg-enriched TME. High levels of MHC class II-positive cancer-associated fibroblasts and tumor endothelial cells were found within EGFR-MU tumors. Owing to the absence of costimulatory factors, they may inhibit rather than activate the tumor antigen-specific CD4+ T-cell response, contributing further to immune suppression. Multiplex immunohistochemistry analyses in a LUAD cohort confirmed increased expression of MHC class II molecules in cancer cells and fibroblasts in EGFR-MU tumors. Our research elucidates the highly immunosuppressive TME in EGFR-MU LUAD and suggests potential targets for effective immunotherapy.

Endothelial cell diversity: the many facets of the crystal

Endothelial cells (ECs) form the inner lining of blood vessels and play crucial roles in angiogenesis. While it has been known for a long time that there are considerable differences among ECs from lymphatic and blood vessels, as well as among arteries, veins and capillaries, the full repertoire of endothelial diversity is only beginning to be elucidated. It has become apparent that the role of ECs is not just limited to their exchange functions. Indeed, a multitude of organ-specific functions, including release of growth factors, regulation of immune functions, have been linked to ECs. Recent years have seen a surge into the identification of spatiotemporal molecular and functional heterogeneity of ECs, supported by technologies such as single-cell RNA sequencing (scRNA-seq), lineage tracing and intersectional genetics. Together, these techniques have spurred the generation of epigenomic, transcriptomic and proteomic signatures of ECs. It is now clear that ECs across organs and in different vascular beds, but even within the same vessel, have unique molecular identities and employ specialized molecular mechanisms to fulfil highly specialized needs. Here, we focus on the molecular heterogeneity of the endothelium in different organs and pathological conditions.

Generation and characterisation of scalable and stable human pluripotent stem cell-derived microvascular-like endothelial cells for cardiac applications

Coronary microvascular disease (CMD) and its progression towards major adverse coronary events pose a significant health challenge. Accurate in vitro investigation of CMD requires a robust cell model that faithfully represents the cells within the cardiac microvasculature. Human pluripotent stem cell-derived endothelial cells (hPSC-ECs) offer great potential; however, they are traditionally derived via differentiation protocols that are not readily scalable and are not specified towards the microvasculature. Here, we report the development and comprehensive characterisation of a scalable 3D protocol enabling the generation of phenotypically stable cardiac hPSC-microvascular-like ECs (hPSC-CMVECs) and cardiac pericyte-like cells. These were derived by growing vascular organoids within 3D stirred tank bioreactors and subjecting the emerging 3D hPSC-ECs to high-concentration VEGF-A treatment (3DV). Not only did this promote phenotypic stability of the 3DV hPSC-ECs; single cell-RNA sequencing (scRNA-seq) revealed the pronounced expression of cardiac endothelial- and microvascular-associated genes. Further, the generated mural cells attained from the vascular organoid exhibited markers characteristic of cardiac pericytes. Thus, we present a suitable cell model for investigating the cardiac microvasculature as well as the endothelial-dependent and -independent mechanisms of CMD. Moreover, owing to their phenotypic stability, cardiac specificity, and high angiogenic potential, the cells described within would also be well suited for cardiac tissue engineering applications.

Novel insights into the heterogeneity of FOXP3 + Treg cells in drug-induced allergic reactions through single-cell transcriptomics

This study uncovers the novel heterogeneity of FOXP3 + regulatory T (Treg) cells and their pivotal role in modulating immune responses during drug-induced allergic reactions, employing cutting-edge single-cell transcriptomics. We established a mouse model for drug-induced allergic reactions and utilized single-cell RNA sequencing (scRNA-seq) to analyze the transcriptomic landscapes of FOXP3 + Treg cells isolated from affected tissues. The study involved both in vitro and in vivo approaches to evaluate the impact of FOXP3 expression levels on the immunoregulatory functions of Treg cells during allergic responses. Techniques included flow cytometry, cluster analysis, principal component analysis (PCA), CCK8 and CSFE assays for cell proliferation, LDH release assays for toxicity, ELISA for cytokine profiling, and CRISPR/Cas9 technology for gene editing. Our findings revealed significant transcriptomic heterogeneity among FOXP3 + Treg cells in the context of drug-induced allergic reactions, with distinct subpopulations exhibiting unique gene expression profiles. This heterogeneity suggests specialized roles in immune regulation. We observed a decrease in the proliferative capacity and cytokine secretion of FOXP3 + Treg cells following allergic stimulation, alongside an increase in reaction toxicity. Manipulating FOXP3 expression levels directly influenced these outcomes, where FOXP3 deletion exacerbated allergic responses, whereas its overexpression mitigated them. Notably, in vivo experiments demonstrated that FOXP3 overexpression significantly reduced the severity of allergic skin reactions in mice. Our study presents novel insights into the heterogeneity and crucial immunoregulatory role of FOXP3 + Treg cells during drug-induced allergic reactions. Overexpression of FOXP3 emerges as a potential therapeutic strategy to alleviate such allergic responses. These findings contribute significantly to our understanding of immune regulation and the development of targeted treatments for drug-induced allergies.

Advancing immunotherapy using biomaterials to control tissue, cellular, and molecular level immune signaling in skin

Immunotherapies have been transformative in many areas, including cancer treatments, allergies, and autoimmune diseases. However, significant challenges persist in extending the reach of these technologies to new indications and patients. Some of the major hurdles include narrow applicability to patient groups, transient efficacy, high cost burdens, poor immunogenicity, and side effects or off-target toxicity that results from lack of disease-specificity and inefficient delivery. Thus, there is a significant need for strategies that control immune responses generated by immunotherapies while targeting infection, cancer, allergy, and autoimmunity. Being the outermost barrier of the body and the first line of host defense, the skin presents a unique immunological interface to achieve these goals. The skin contains a high concentration of specialized immune cells, such as antigen-presenting cells and tissue-resident memory T cells. These cells feature diverse and potent combinations of immune receptors, providing access to cellular and molecular level control to modulate immune responses. Thus, skin provides accessible tissue, cellular, and molecular level controls that can be harnessed to improve immunotherapies. Biomaterial platforms - microneedles, nano- and micro-particles, scaffolds, and other technologies - are uniquely capable of modulating the specialized immunological niche in skin by targeting these distinct biological levels of control. This review highlights recent pre-clinical and clinical advances in biomaterial-based approaches to target and modulate immune signaling in the skin at the tissue, cellular, and molecular levels for immunotherapeutic applications. We begin by discussing skin cytoarchitecture and resident immune cells to establish the biological rationale for skin-targeting immunotherapies. This is followed by a critical presentation of biomaterial-based pre-clinical and clinical studies aimed at controlling the immune response in the skin for immunotherapy and therapeutic vaccine applications in cancer, allergy, and autoimmunity.

Single-cell transcriptome profiling reveals cell type-specific variation and development in HLA expression of human skin

Skin, the largest organ of body, is a highly immunogenic tissue with a diverse collection of immune cells. Highly polymorphic human leukocyte antigen (HLA) molecules have a central role in coordinating immune responses as recognition molecules. Nevertheless, HLA gene expression patterns among diverse cell types within a specific organ, like the skin, have yet to be thoroughly investigated, with stromal cells attracting much less attention than immune cells. To illustrate HLA expression profiles across different cell types in the skin, we performed single-cell RNA sequencing (scRNA-seq) analyses on skin datasets, covering adult and fetal skin, and hair follicles as the skin appendages. We revealed the variation in HLA expression between different skin populations by examining normal adult skin datasets. Moreover, we evaluated the potential immunogenicity of multiple skin populations based on the expression of classical HLA class I genes, which were well represented in all cell types. Furthermore, we generated scRNA-seq data of developing skin from fetuses of 15 post conception weeks (PCW), 17 PCW, and 22 PCW, delineating the dynamic expression of HLA genes with cell type-dependent variation among various cell types during development. Notably, the pseudotime trajectory analysis unraveled the significant variance in HLA genes during the evolution of vascular endothelial cells. Moreover, we uncovered the immune-privileged properties of hair follicles at single-cell resolution. Our study presents a comprehensive single-cell transcriptomic landscape of HLA genes in the skin, which provides new insights into variation in HLA molecules and offers a clue for allogeneic skin transplantation.

A high-resolution view of the heterogeneous aging endothelium

Vascular endothelial cell (EC) aging has a strong impact on tissue perfusion and overall cardiovascular health. While studies confined to the investigation of aging-associated vascular readouts in one or a few tissues have already drastically expanded our understanding of EC aging, single-cell omics and other high-resolution profiling technologies have started to illuminate the intricate molecular changes underlying endothelial aging across diverse tissues and vascular beds at scale. In this review, we provide an overview of recent insights into the heterogeneous adaptations of the aging vascular endothelium. We address critical questions regarding tissue-specific and universal responses of the endothelium to the aging process, EC turnover dynamics throughout lifespan, and the differential susceptibility of ECs to acquiring aging-associated traits. In doing so, we underscore the transformative potential of single-cell approaches in advancing our comprehension of endothelial aging, essential to foster the development of future innovative therapeutic strategies for aging-associated vascular conditions.

Revealing the pathogenesis of keloids based on the status: Active vs inactive

Recently, the pathomechanisms of keloids have been extensively researched using transcriptomic analysis, but most studies did not consider the activity of keloids. We aimed to profile the transcriptomics of keloids according to their clinical activity and location within the keloid lesion, compared with normal and mature scars. Tissue samples were collected (keloid based on its activity (active and inactive), mature scar from keloid patients and normal scar (NS) from non-keloid patients). To reduce possible bias, all keloids assessed in this study had no treatment history and their location was limited to the upper chest or back. Multiomics assessment was performed by using single-cell RNA sequencing and multiplex immunofluorescence. Increased mesenchymal fibroblasts (FBs) was the main feature in keloid patients. Noticeably, the proportion of pro-inflammatory FBs was significantly increased in active keloids compared to inactive ones. To explore the nature of proinflammatory FBs, trajectory analysis was conducted and CCN family associated with mechanical stretch exhibited higher expression in active keloids. For vascular endothelial cells (VECs), the proportion of tip and immature cells increased in keloids compared to NS, especially at the periphery of active keloids. Also, keloid VECs highly expressed genes with characteristics of mesenchymal activation compared to NS, especially those from the active keloid center. Multiomics analysis demonstrated the distinct expression profile of active keloids. Clinically, these findings may provide the future appropriate directions for development of treatment modalities of keloids. Prevention of keloids could be possible by the suppression of mesenchymal activation between FBs and VECs and modulation of proinflammatory FBs may be the key to the control of active keloids.

Characterization of Vascular Niche in Systemic Sclerosis by Spatial Proteomics

BACKGROUND

Systemic sclerosis (SSc) is a connective tissue disease that can serve as a model to study vascular changes in response to inflammation, autoimmunity, and fibrotic remodeling. Although microvascular changes are the earliest histopathologic manifestation of SSc, the vascular pathophysiology remains poorly understood.

METHODS

We applied spatial proteomic approaches to deconvolute the heterogeneity of vascular cells at the single-cell level in situ and characterize cellular alterations of the vascular niches of patients with SSc. Skin biopsies of patients with SSc and control individuals were analyzed by imaging mass cytometry, yielding a total of 90 755 cells including 2987 endothelial cells and 4096 immune cells.

RESULTS

We identified 7 different subpopulations of blood vascular endothelial cells (VECs), 2 subpopulations of lymphatic endothelial cells, and 3 subpopulations of pericytes. A novel population of CD34+;αSMA+ (α-smooth muscle actin);CD31+ VECs was more common in SSc, whereas endothelial precursor cells were decreased. Co-detection by indexing and tyramide signal amplification confirmed these findings. The microenvironment of CD34+;αSMA+;CD31+ VECs was enriched for immune cells and myofibroblasts, and CD34+;αSMA+;CD31+ VECs expressed markers of endothelial-to-mesenchymal transition. The density of CD34+;αSMA+;CD31+ VECs was associated with clinical progression of fibrosis in SSc.

CONCLUSIONS

Using spatial proteomics, we unraveled the heterogeneity of vascular cells in control individuals and patients with SSc. We identified CD34+;αSMA+;CD31+ VECs as a novel endothelial cell population that is increased in patients with SSc, expresses markers for endothelial-to-mesenchymal transition, and is located in close proximity to immune cells and myofibroblasts. CD34+;αSMA+;CD31+ VEC counts were associated with clinical outcomes of progressive fibrotic remodeling, thus providing a novel cellular correlate for the crosstalk of vasculopathy and fibrosis.

Endothelial Dysfunction in Psoriasis: An Integrative Review

Vascular endothelial cells (ECs), the inner layer of blood vessels, were previously considered to be a passive lining that facilitates cellular and molecular exchange. However, recent studies have revealed that ECs can respond to various stimuli and actively regulate vascular function and skin inflammation. Specific subtypes of ECs are known to have significant roles in a diverse range of physiological and pathological processes in the skin. This review suggests that EC dysfunction is both causal and consequential in the pathogenesis of psoriasis. Further investigations into dysregulated pathways in EC dysfunction may provide new insights for the treatment of psoriasis.

The BulkECexplorer compiles endothelial bulk transcriptomes to predict functional versus leaky transcription

Transcriptomic data can be mined to understand the molecular activity of cell types. Yet, functional genes may remain undetected in RNA sequencing (RNA-seq) experiments for technical reasons, such as insufficient read depth or gene dropout. Conversely, RNA-seq experiments may detect lowly expressed mRNAs thought to be biologically irrelevant products of leaky transcription. To represent a cell type's functional transcriptome more accurately, we propose compiling many bulk RNA-seq datasets into a compendium and applying established classification models to predict whether detected transcripts are likely products of active or leaky transcription. Here, we present the BulkECexplorer (bulk RNA-seq endothelial cell explorer) compendium of 240 bulk RNA-seq datasets from five vascular endothelial cell subtypes. This resource reports transcript counts for genes of interest and predicts whether detected transcripts are likely the products of active or leaky gene expression. Beyond its usefulness for vascular biology research, this resource provides a blueprint for developing analogous tools for other cell types.

The intricate cellular ecosystem of human peripheral veins as revealed by single-cell transcriptomic analysis

The venous system has been historically understudied despite its critical roles in blood distribution, heart function, and systemic immunity. This study dissects the microanatomy of upper arm veins at the single cell level, and how it relates to wall structure, remodeling processes, and inflammatory responses to injury. We applied single-cell RNA sequencing to 4 non-diseased human veins (3 basilic, 1 cephalic) obtained from organ donors, followed by bioinformatic and histological analyses. Unsupervised clustering of 20,006 cells revealed a complex ecosystem of endothelial cell (EC) types, smooth muscle cell (SMCs) and pericytes, various types of fibroblasts, and immune cell populations. The venous endothelium showed significant upregulation of cell adhesion genes, with arteriovenous zonation EC phenotypes highlighting the heterogeneity of vasa vasorum (VV) microvessels. Venous SMCs had atypical contractile phenotypes and showed widespread localization in the intima and media. MYH11+DESlo SMCs were transcriptionally associated with negative regulation of contraction and pro-inflammatory gene expression. MYH11+DEShi SMCs showed significant upregulation of extracellular matrix genes and pro-migratory mediators. Venous fibroblasts ranging from secretory to myofibroblastic phenotypes were 4X more abundant than SMCs and widely distributed throughout the wall. Fibroblast-derived angiopoietin-like factors were identified as versatile signaling hubs to regulate angiogenesis and SMC proliferation. An abundant monocyte/macrophage population was detected and confirmed by histology, including pro-inflammatory and homeostatic phenotypes, with cell counts positively correlated with age. Ligand-receptor interactome networks identified the venous endothelium in the main lumen and the VV as a niche for monocyte recruitment and infiltration. This study underscores the transcriptional uniqueness of venous cells and their relevance for vascular inflammation and remodeling processes. Findings from this study may be relevant for molecular investigations of upper arm veins used for vascular access creation, where single-cell analyses of cell composition and phenotypes are currently lacking.

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

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

The mycosis fungoides cutaneous microenvironment shapes dysfunctional cell trafficking, antitumor immunity, matrix interactions, and angiogenesis

Malignant T lymphocyte proliferation in mycosis fungoides (MF) is largely restricted to the skin, implying that malignant cells are dependent on their specific cutaneous tumor microenvironment (TME), including interactions with non-malignant immune and stromal cells, cytokines, and other immunomodulatory factors. To explore these interactions, we performed a comprehensive transcriptome analysis of the TME in advanced-stage MF skin tumors by single-cell RNA sequencing. Our analysis identified cell-type compositions, cellular functions, and cell-to-cell interactions in the MF TME that were distinct from those from healthy skin and benign dermatoses. While patterns of gene expression were common among patient samples, high transcriptional diversity was also observed in immune and stromal cells, with dynamic interactions and crosstalk between these cells and malignant T lymphocytes. This heterogeneity mapped to processes such as cell trafficking, matrix interactions, angiogenesis, immune functions, and metabolism that affect cancer cell growth, migration, and invasion, as well as antitumor immunity. By comprehensively characterizing the transcriptomes of immune and stromal cells within the cutaneous microenvironment of individual MF tumors, we have identified patterns of dysfunction common to all tumors that represent a resource for identifying candidates with therapeutic potential as well as patient-specific heterogeneity that has important implications for personalized disease management.

CD31 defines a subpopulation of human adipose-derived regenerative cells with potent angiogenic effects

Cellular heterogeneity represents a major challenge for regenerative treatment using freshly isolated Adipose Derived Regenerative Cells (ADRCs). Emerging data suggest superior efficacy of ADRCs as compared to the ex vivo expanded and more homogeneous ADRCs (= ASCs) for indications involving (micro)vascular deficiency, however, it remains unknown which ADRC cell subtypes account for the improvement. Surprisingly, we found regarding erectile dysfunction (ED) that the number of injected CD31+  ADRCs correlated positively with erectile function 12 months after one bolus of autologous ADRCs. Comprehensive in vitro and ex vivo analyses confirmed superior pro-angiogenic and paracrine effects of human CD31+ enriched ADRCs compared to the corresponding CD31- and parent ADRCs. When CD31+, CD31- and ADRCs were co-cultured in aortic ring- and corpus cavernous tube formation assays, the CD31+  ADRCs induced significantly higher tube development. This effect was corroborated using conditioned medium (CM), while quantitative mass spectrometric analysis suggested that this is likely explained by secretory pro-angiogenic proteins including DKK3, ANGPT2, ANAX2 and VIM, all enriched in CD31+  ADRC CM. Single-cell RNA sequencing showed that transcripts of the upregulated and secreted proteins were present in 9 endothelial ADRC subsets including endothelial progenitor cells in the heterogenous non-cultured ADRCs. Our data suggest that the vascular benefit of using ADRCs in regenerative medicine is dictated by CD31+  ADRCs.

Single-cell RNA-seq and bulk-seq identify RAB17 as a potential regulator of angiogenesis by human dermal microvascular endothelial cells in diabetic foot ulcers

Background

Angiogenesis is crucial in diabetic wound healing and is often impaired in diabetic foot ulcers (DFUs). Human dermal microvascular endothelial cells (HDMECs) are vital components in dermal angiogenesis; however, their functional and transcriptomic characteristics in DFU patients are not well understood. This study aimed to comprehensively analyse HDMECs from DFU patients and healthy controls and find the potential regulator of angiogenesis in DFUs.

Methods

HDMECs were isolated from skin specimens of DFU patients and healthy controls via magnetic-activated cell sorting. The proliferation, migration and tube-formation abilities of the cells were then compared between the experimental groups. Both bulk RNA sequencing (bulk-seq) and single-cell RNA-seq (scRNA-seq) were used to identify RAB17 as a potential marker of angiogenesis, which was further confirmed via weighted gene co-expression network analysis (WGCNA) and least absolute shrink and selection operator (LASSO) regression. The role of RAB17 in angiogenesis was examined through in vitro and in vivo experiments.

Results

The isolated HDMECs displayed typical markers of endothelial cells. HDMECs isolated from DFU patients showed considerably impaired tube formation, rather than proliferation or migration, compared to those from healthy controls. Gene set enrichment analysis (GSEA), fGSEA, and gene set variation analysis (GSVA) of bulk-seq and scRNA-seq indicated that angiogenesis was downregulated in DFU-HDMECs. LASSO regression identified two genes, RAB17 and CD200, as characteristic of DFU-HDMECs; additionally, the expression of RAB17 was found to be significantly reduced in DFU-HDMECs compared to that in the HDMECs of healthy controls. Overexpression of RAB17 was found to enhance angiogenesis, the expression of hypoxia inducible factor-1α and vascular endothelial growth factor A, and diabetic wound healing, partially through the mitogen-activated protein kinase/extracellular signal-regulated kinase signalling pathway.

Conclusions

Our findings suggest that the impaired angiogenic capacity in DFUs may be related to the dysregulated expression of RAB17 in HDMECs. The identification of RAB17 as a potential molecular target provides a potential avenue for the treatment of impaired angiogenesis in DFUs.

Mesenchymal Stem Cell Aggregation-Released Extracellular Vesicles Induce CD31+ EMCN+ Vessels in Skin Regeneration and Improve Diabetic Wound Healing

The blood vessel system is essential for skin homeostasis and regeneration. While the heterogeneity of vascular endothelial cells has been emergingly revealed, whether a regeneration-relevant vessel subtype exists in skin remains unknown. Herein, a specialized vasculature in skin featured by simultaneous CD31 and EMCN expression contributing to the regeneration process is identified, the decline of which functionally underlies the impaired angiogenesis of diabetic nonhealing wounds. Moreover, enlightened by the developmental process that mesenchymal condensation induces angiogenesis, it is demonstrated that mesenchymal stem/stromal cell aggregates (CAs) provide an efficacious therapy to enhance regrowth of CD31+ EMCN+ vessels in diabetic wounds, which is surprisingly suppressed by pharmacological inhibition of extracellular vesicle (EV) release. It is further shown that CAs promote secretion of angiogenic protein-enriched EVs by proteomic analysis, which directly exert high efficacy in boosting CD31+ EMCN+ vessels and treating nonhealing diabetic wounds. These results add to the current knowledge on skin vasculature and help establish feasible strategies to benefit wound healing under diabetic condition.

Endothelial-derived angiocrine factors as instructors of embryonic development

Blood vessels are well-known to play roles in organ development and repair, primarily owing to their fundamental function in delivering oxygen and nutrients to tissues to promote their growth and homeostasis. Endothelial cells however are not merely passive conduits for carrying blood. There is now evidence that endothelial cells of the vasculature actively regulate tissue-specific development, morphogenesis and organ function, as well as playing roles in disease and cancer. Angiocrine factors are growth factors, cytokines, signaling molecules or other regulators produced directly from endothelial cells to instruct a diverse range of signaling outcomes in the cellular microenvironment, and are critical mediators of the vascular control of organ function. The roles of angiocrine signaling are only beginning to be uncovered in diverse fields such as homeostasis, regeneration, organogenesis, stem-cell maintenance, cell differentiation and tumour growth. While in some cases the specific angiocrine factor involved in these processes has been identified, in many cases the molecular identity of the angiocrine factor(s) remain to be discovered, even though the importance of angiocrine signaling has been implicated. In this review, we will specifically focus on roles for endothelial-derived angiocrine signaling in instructing tissue morphogenesis and organogenesis during embryonic and perinatal development.

Endothelial Cells and Mitochondria: Two Key Players in Liver Transplantation

Building the inner layer of our blood vessels, the endothelium forms an important line communicating with deeper parenchymal cells in our organs. Previously considered passive, endothelial cells are increasingly recognized as key players in intercellular crosstalk, vascular homeostasis, and blood fluidity. Comparable to other cells, their metabolic function strongly depends on mitochondrial health, and the response to flow changes observed in endothelial cells is linked to their mitochondrial metabolism. Despite the direct impact of new dynamic preservation concepts in organ transplantation, the impact of different perfusion conditions on sinusoidal endothelial cells is not yet explored well enough. This article therefore describes the key role of liver sinusoidal endothelial cells (LSECs) together with their mitochondrial function in the context of liver transplantation. The currently available ex situ machine perfusion strategies are described with their effect on LSEC health. Specific perfusion conditions, including perfusion pressure, duration, and perfusate oxygenation are critically discussed considering the metabolic function and integrity of liver endothelial cells and their mitochondria.

A single-cell transcriptome atlas of pig skin characterizes anatomical positional heterogeneity

Different anatomical locations of the body skin show differences in their gene expression patterns depending on different origins, and the inherent heterogeneous information can be maintained in adults. However, highly resolvable cellular specialization is less well characterized in different anatomical regions of the skin. Pig is regarded as an excellent model animal for human skin research in view of its similar physiology to human. In this study, single-cell RNA sequencing was performed on pig skin tissues from six different anatomical regions of Chenghua (CH) pigs, with a superior skin thickness trait, and the back site of large white (LW) pigs. We obtained 233,715 cells, representing seven cell types, among which we primarily characterized the heterogeneity of the top three cell types, including smooth muscle cells (SMCs), endothelial cells (ECs), and fibroblasts (FBs). Then, we further identified several subtypes of SMCs, ECs, and FBs, and discovered the expression patterns of site-specific genes involved in some important pathways such as the immune response and extracellular matrix (ECM) synthesis in different anatomical regions. By comparing differentially expressed genes of skin FBs among different anatomical regions, we considered TNN, COL11A1, and INHBA as candidate genes for facilitating ECM accumulation. These findings of heterogeneity in the main skin cell types from different anatomical sites will contribute to a better understanding of inherent skin information and place the potential focus on skin generation, transmission, and transplantation, paving the foundation for human skin priming.

The First Transcriptomic Atlas of the Adult Lacrimal Gland Reveals Epithelial Complexity and Identifies Novel Progenitor Cells in Mice

The lacrimal gland (LG) secretes aqueous tears. Previous studies have provided insights into the cell lineage relationships during tissue morphogenesis. However, little is known about the cell types composing the adult LG and their progenitors. Using scRNAseq, we established the first comprehensive cell atlas of the adult mouse LG to investigate the cell hierarchy, its secretory repertoire, and the sex differences. Our analysis uncovered the complexity of the stromal landscape. Epithelium subclustering revealed myoepithelial cells, acinar subsets, and two novel acinar subpopulations: Tfrchi and Car6hi cells. The ductal compartment contained Wfdc2+ multilayered ducts and an Ltf+ cluster formed by luminal and intercalated duct cells. Kit+ progenitors were identified as: Krt14+ basal ductal cells, Aldh1a1+ cells of Ltf+ ducts, and Sox10+ cells of the Car6hi acinar and Ltf+ epithelial clusters. Lineage tracing experiments revealed that the Sox10+ adult populations contribute to the myoepithelial, acinar, and ductal lineages. Using scRNAseq data, we found that the postnatally developing LG epithelium harbored key features of putative adult progenitors. Finally, we showed that acinar cells produce most of the sex-biased lipocalins and secretoglobins detected in mouse tears. Our study provides a wealth of new data on LG maintenance and identifies the cellular origin of sex-biased tear components.

Single-cell transcriptomics reveals that tumor-infiltrating natural killer cells are activated by localized ablative immunotherapy and share anti-tumor signatures induced by immune checkpoint inhibitors

Rationale

Natural killer (NK) cells provide protective anti-cancer immunity. However, the cancer therapy induced activation gene signatures and pathways in NK cells remain unclear.

Methods

We applied a novel localized ablative immunotherapy (LAIT) by synergizing photothermal therapy (PTT) with intra-tumor delivering of the immunostimulant N-dihydrogalactochitosan (GC), to treat breast cancer using a mammary tumor virus-polyoma middle tumor-antigen (MMTV-PyMT) mouse model. We performed single-cell RNA sequencing (scRNAseq) analysis to unveil the cellular heterogeneity and compare the transcriptional alterations induced by PTT, GC, and LAIT in NK cells within the tumor microenvironment (TME).

Results

ScRNAseq showed that NK subtypes, including cycling, activated, interferon-stimulated, and cytotoxic NK cells. Trajectory analysis revealed a route toward activation and cytotoxicity following pseudotime progression. Both GC and LAIT elevated gene expression associated with NK cell activation, cytolytic effectors, activating receptors, IFN pathway components, and cytokines/chemokines in NK subtypes. Single-cell transcriptomics analysis using immune checkpoint inhibitor (ICI)-treated animal and human samples revealed that ICI-induced NK activation and cytotoxicity across several cancer types. Furthermore, ICI-induced NK gene signatures were also induced by LAIT treatment. We also discovered that several types of cancer patients had significantly longer overall survival when they had higher expression of genes in NK cells that were also specifically upregulated by LAIT.

Conclusion

Our findings show for the first time that LAIT activates cytotoxicity in NK cells and the upregulated genes positively correlate with beneficial clinical outcomes for cancer patients. More importantly, our results further establish the correlation between the effects of LAIT and ICI on NK cells, hence expanding our understanding of mechanism of LAIT in remodeling TME and shedding light on the potentials of NK cell activation and anti-tumor cytotoxic functions in clinical applications.

Are circulating endothelial cells the next target for transcriptome-level pathway analysis in ARDS?

Acute respiratory distress syndrome (ARDS) has had no mortality-improving pharmacological intervention despite 50 years of high-caliber research due to its heterogeneity (Huppert LA, Matthay MA, Ware LB. Semin Respir Crit Care Med 40: 31-39, 2019). For the field to advance, better definitions for ARDS subgroups that more uniformly respond to therapies are needed (Bos LDJ, Scicluna BP, Ong DSY, Cremer O, van der Poll T, Schultz MJ. Am J Respir Crit Care Med 200: 42-50, 2019; Dickson RP, Schultz MJ, T van der P, Schouten LR, Falkowski NR, Luth JE, Sjoding MW, Brown CA, Chanderraj R, Huffnagle GB, Bos LDJ, Biomarker Analysis in Septic ICU Patients (BASIC) Consortium. Am J Respir Crit Care Med 201: 555-563, 2020; Sinha P, Calfee CS. Am J Respir Crit Care Med 200: 4-6, 2019; Calfee CS, Delucchi K, Parsons PE, Thompson BT, Ware LB, Matthay MA, NHLBI ARDS Network. Lancet Respir Med 2: 611-620, 2014; Hendrickson CM, Matthay MA. Pulm Circ 8: 1-12, 2018). A plethora of high-quality clinical research has uncovered the next generation of soluble biomarkers that provide the predictive enrichment necessary for trial recruitment; however, plasma-soluble markers do not specify the damaged organ of origin nor do they provide insight into disease mechanisms. In this perspective, we make the case for querying the transcriptome of circulating endothelial cells (CECs), which when shed from vessels after inflammatory insult, become heralds of site-specific inflammatory damage. We review the application of CEC quantification to multiple disease phenotypes (including myocardial infarction, vasculitides, cancer, and ARDS), in each case supporting the association of CEC number with disease severity. We also argue for the utility of single-cell RNA transcriptomics to the understanding of cell-specific contributions to disease pathophysiology and its potential to uncover novel insight on signals contributing to CEC shedding in ARDS.

Single-cell transcriptomic analysis of tumor heterogeneity and intercellular networks in human urothelial carcinoma

BACKGROUND

Heterogeneity of tumor cells and the tumor microenvironment (TME) is significantly associated with clinical outcomes and treatment responses in patients with urothelial carcinoma (UC). Comprehensive profiling of the cellular diversity and interactions between malignant cells and TME may clarify the mechanisms underlying UC progression and guide the development of novel therapies. This study aimed to extend our understanding of intra-tumoral heterogeneity and the immunosuppressive TME in UC and provide basic support for the development of novel UC therapies.

METHODS

Seven patients with UC were included who underwent curative surgery at our hospital between July 2020 and October 2020. We performed single-cell RNA sequencing (scRNA-seq) analysis in seven tumors with six matched adjacent normal tissues and integrated the results with two public scRNA-seq datasets. The functional properties and intercellular interactions between single cells were characterized, and the results were validated using multiplex immunofluorescence staining, flow cytometry, and bulk transcriptomic datasets. All statistical analyses were performed using the R package with two-sided tests. Wilcoxon-rank test, log-rank test, one-way analysis of variance test, and Pearson correlation analysis were used properly.

RESULTS

Unsupervised t-distributed stochastic neighbor embedding clustering analysis identified ten main cellular subclusters in urothelial tissues. Of them, seven urothelial subtypes were noted, and malignant urothelial cells were characterized with enhanced cellular proliferation and reduced immunogenicity. CD8 + T cell subclusters exhibited enhanced cellular cytotoxicity activities along with increased exhaustion signature in UC tissues, and the recruitment of CD4 + T regulatory cells was also increased in tumor tissues. Regarding myeloid cells, coordinated reprogramming of infiltrated neutrophils, M2-type polarized macrophages, and LAMP3 + dendritic cells contribute to immunosuppressive TME in UC tissues. Tumor tissues demonstrated enhanced angiogenesis mediated by KDR + endothelial cells and RGS5 + /ACTA2 + pericytes. Through deconvolution analysis, we identified multiple cellular subtypes may influence the programmed death-ligand 1 (PD-L1) immunotherapy response in patients with UC.

CONCLUSION

Our scRNA-seq analysis clarified intra-tumoral heterogeneity and delineated the pro-tumoral and immunosuppressive microenvironment in UC tissues, which may provide novel therapeutic targets.

Significance of anti-transcobalamin receptor antibodies in cutaneous arteritis revealed by proteome-wide autoantibody screening

Cutaneous arteritis (CA) is a single-organ vasculitis that exclusively affects the small to medium-sized arteries of the skin. Diagnosis depends on a histological investigation with skin biopsy, which could be burdensome for both patients and clinicians. Moreover, the pathogenesis of CA remains unstudied, and treatment has not yet been established. Herein, we applied our proteome-wide autoantibody screening method to explore autoantibodies in the serum of CA patients. As a result, anti-transcobalamin receptor (TCblR) antibodies (Abs) were specifically detected in 24% of CA patients. Patients with positive anti-TCblR Abs were spared from peripheral neuropathy compared to those with negative anti-TCblR Abs, showing characteristics as CA confined to the skin. In addition, we revealed that anti-TCblR Abs trigger the autocrine loop of interleukin-6 mediated by tripartite motif-containing protein 21 in human endothelial cells and induce periarterial inflammation in murine skin. Furthermore, we demonstrated that methylcobalamin, a ligand of TCblR, ameliorates inflammation caused by anti-TCblR Abs both in vitro and in vivo. Collectively, our investigation unveils the pathologic significance of anti-TCblR Abs in CA and their potential as a diagnostic marker and a pathophysiology-oriented therapeutic target.

CD146 expression profile in human skin and pre-vascularized dermo-epidermal skin substitutes in vivo

BACKGROUND

CD146 is a cell adhesion molecule whose expression profile in human skin has not yet been elucidated. Here, we characterize CD146 expression pattern in human skin, in particular in blood endothelial cells (BECs) and lymphatic endothelial cells (LECs), which constitute human dermal microvascular endothelial cells (HDMECs), as well as in perivascular cells.

RESULTS

We demonstrated that CD146 is a specific marker of BECs, but not of LECs. Moreover, we found CD146 expression also in human pericytes surrounding blood capillaries in human skin. In addition, we demonstrated that CD146 expression is up-regulated by the TNFα-IL-1β/NF-kB axis in both BECs and pericytes. Finally, we engineered 3D collagen hydrogels composed of HDMECs, CD146⁺ pericytes, and fibroblasts which developed, in vitro and in vivo, a complete microvasculature network composed of blood and lymphatic capillaries with pericytes investing blood capillaries.

CONCLUSIONS

Overall, our results proved that CD146 is a specific marker of BECs and pericytes, but not LECs in human skin. Further, the combination of CD146⁺ pericytes with HDMECs in skin substitutes allowed to bioengineer a comprehensive 3D in vitro and in vivo model of the human dermal microvasculature.

Multimodal human thymic profiling reveals trajectories and cellular milieu for T agonist selection

To prevent autoimmunity, thymocytes expressing self-reactive T cell receptors (TCRs) are negatively selected, however, divergence into tolerogenic, agonist selected lineages represent an alternative fate. As thymocyte development, selection, and lineage choices are dependent on spatial context and cell-to-cell interactions, we have performed Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) and spatial transcriptomics on paediatric human thymu​​s. Thymocytes expressing markers of strong TCR signalling diverged from the conventional developmental trajectory prior to CD4⁺ or CD8⁺ lineage commitment, while markers of different agonist selected T cell populations (CD8αα(I), CD8αα(II), T_(agonist), T_reg(diff), and T_reg) exhibited variable timing of induction. Expression profiles of chemokines and co-stimulatory molecules, together with spatial localisation, supported that dendritic cells, B cells, and stromal cells contribute to agonist selection, with different subsets influencing thymocytes at specific developmental stages within distinct spatial niches. Understanding factors influencing agonist T cells is needed to benefit from their immunoregulatory effects in clinical use.

Mitochondrial ATP Production is Required for Endothelial Cell Control of Vascular Tone

Arteries and veins are lined by nonproliferating endothelial cells that play a critical role in regulating blood flow. Endothelial cells also regulate tissue perfusion, metabolite exchange, and thrombosis. It is thought that endothelial cells rely on ATP generated via glycolysis, rather than mitochondrial oxidative phosphorylation, to fuel each of these energy-demanding processes. However, endothelial metabolism has mainly been studied in the context of proliferative cells, and little is known about energy production in endothelial cells within the fully formed vascular wall. Using intact arteries isolated from rats and mice, we show that inhibiting mitochondrial respiration disrupts endothelial control of vascular tone. Basal, mechanically activated, and agonist-evoked calcium activity in intact artery endothelial cells are each prevented by inhibiting mitochondrial ATP synthesis. Agonist-evoked calcium activity was also inhibited by blocking the transport of pyruvate, the master fuel for mitochondrial energy production, through the mitochondrial pyruvate carrier. The role for mitochondria in endothelial cell energy production is independent of species, sex, or vascular bed. These data show that a mitochondrial ATP supply is necessary for calcium-dependent, nitric oxide-mediated endothelial control of vascular tone, and identifies the critical role of endothelial mitochondrial energy production in fueling perfused blood vessel function.

Comprehensive bioinformatics analysis reveals biomarkers of DNA methylation-related genes in varicose veins

Background: Patients with Varicose veins (VV) show no obvious symptoms in the early stages, and it is a common and frequent clinical condition. DNA methylation plays a key role in VV by regulating gene expression. However, the molecular mechanism underlying methylation regulation in VV remains unclear. Methods: The mRNA and methylation data of VV and normal samples were obtained from the Gene Expression Omnibus (GEO) database. Methylation-Regulated Genes (MRGs) between VV and normal samples were crossed with VV-associated genes (VVGs) obtained by weighted gene co-expression network analysis (WGCNA) to obtain VV-associated MRGs (VV-MRGs). Their ability to predict disease was assessed using receiver operating characteristic (ROC) curves. Biomarkers were then screened using a random forest model (RF), support vector machine model (SVM), and generalized linear model (GLM). Next, gene set enrichment analysis (GSEA) was performed to explore the functions of biomarkers. Furthermore, we also predicted their drug targets, and constructed a competing endogenous RNAs (ceRNA) network and a drug target network. Finally, we verified their mRNA expression using quantitative real-time polymerase chain reaction (qRT-PCR). Results: Total three VV-MRGs, namely Wnt1-inducible signaling pathway protein 2 (WISP2), Cysteine-rich intestinal protein 1 (CRIP1), and Odd-skipped related 1 (OSR1) were identified by VVGs and MRGs overlapping. The area under the curves (AUCs) of the ROC curves for these three VV-MRGs were greater than 0.8. RF was confirmed as the optimal diagnostic model, and WISP2, CRIP1, and OSR1 were regarded as biomarkers. GSEA showed that WISP2, CRIP1, and OSR1 were associated with oxidative phosphorylation, extracellular matrix (ECM), and respiratory system functions. Furthermore, we found that lncRNA MIR17HG can regulate OSR1 by binding to hsa-miR-21-5p and that PAX2 might treat VV by targeting OSR1. Finally, qRT-PCR results showed that the mRNA expression of the three genes was consistent with the results of the datasets. Conclusion: This study identified WISP2, CRIP1, and OSR1 as biomarkers of VV through comprehensive bioinformatics analysis, and preliminary explored the DNA methylation-related molecular mechanism in VV, which might be important for VV diagnosis and exploration of potential molecular mechanisms.

Metabolic Reprogramming in Tumor Endothelial Cells

The dynamic crosstalk between the different components of the tumor microenvironment is critical to determine cancer progression, metastatic dissemination, tumor immunity, and therapeutic responses. Angiogenesis is critical for tumor growth, and abnormal blood vessels contribute to hypoxia and acidosis in the tumor microenvironment. In this hostile environment, cancer and stromal cells have the ability to alter their metabolism in order to support the high energetic demands and favor rapid tumor proliferation. Recent advances have shown that tumor endothelial cell metabolism is reprogrammed, and that targeting endothelial metabolic pathways impacts developmental and pathological vessel sprouting. Therefore, the use of metabolic antiangiogenic therapies to normalize the blood vasculature, in combination with immunotherapies, offers a clinical niche to treat cancer.

Single-cell transcriptome atlas of the human corpus cavernosum

The corpus cavernosum is the most important structure for penile erection, and its dysfunction causes many physiological and psychological problems. However, its cellular heterogeneity and signalling networks at the molecular level are poorly understood because of limited access to samples. Here, we profile 64,993 human cavernosal single-cell transcriptomes from three males with normal erection and five organic erectile dysfunction patients. Cell communication analysis reveals that cavernosal fibroblasts are central to the paracrine signalling network and regulate microenvironmental homeostasis. Combining with immunohistochemical staining, we reveal the cellular heterogeneity and describe a detailed spatial distribution map for each fibroblast, smooth muscle and endothelial subcluster in the corpus cavernosum. Furthermore, comparative analysis and related functional experiments identify candidate regulatory signalling pathways in the pathological process. Our study provides an insight into the human corpus cavernosum microenvironment and a reference for potential erectile dysfunction therapies.

The corpus cavernosum is the most important structure for penile erection, and its dysfunction causes physiological and psychological problems. Here the authors perform single-cell RNA-sequencing on corpus cavernosum samples from males with normal erection and erectile dysfunction patients, providing insights into this pathology.

Claudin5 protects the peripheral endothelial barrier in an organ and vessel type-specific manner

Dysfunctional and leaky blood vessels resulting from disruption of the endothelial cell (EC) barrier accompanies numerous diseases. The EC barrier is established through endothelial cell tight and adherens junctions. However, the expression pattern and precise contribution of different junctional proteins to the EC barrier is poorly understood. Here, we focus on organs with continuous endothelium to identify structural and functional in vivo characteristics of the EC barrier. Assembly of multiple single-cell RNAseq datasets into a single integrated database revealed the variability and commonalities of EC barrier patterning. Across tissues, Claudin5 exhibited diminishing expression along the arteriovenous axis, correlating with EC barrier integrity. Functional analysis identified tissue-specific differences in leakage properties and response to the leakage agonist histamine. Loss of Claudin5 enhanced histamine-induced leakage in an organotypic and vessel type-specific manner in an inducible, EC-specific, knock-out mouse. Mechanistically, Claudin5 loss left junction ultrastructure unaffected but altered its composition, with concomitant loss of zonula occludens-1 and upregulation of VE-Cadherin expression. These findings uncover the organ-specific organisation of the EC barrier and distinct importance of Claudin5 in different vascular beds, providing insights to modify EC barrier stability in a targeted, organ-specific manner.

Novel Blood Vascular Endothelial Subtype-Specific Markers in Human Skin Unearthed by Single-Cell Transcriptomic Profiling

Ample evidence pinpoints the phenotypic diversity of blood vessels (BVs) and site-specific functions of their lining endothelial cells (ECs). We harnessed single-cell RNA sequencing (scRNA-seq) to dissect the molecular heterogeneity of blood vascular endothelial cells (BECs) in healthy adult human skin and identified six different subpopulations, signifying arterioles, post-arterial capillaries, pre-venular capillaries, post-capillary venules, venules and collecting venules. Individual BEC subtypes exhibited distinctive transcriptomic landscapes associated with diverse biological pathways. These functionally distinct dermal BV segments were characterized by their unique compositions of conventional and novel markers (e.g., arteriole marker GJA5; arteriole capillary markers ASS1 and S100A4; pre-venular capillary markers SOX17 and PLAUR; venular markers EGR2 and LRG1), many of which have been implicated in vascular remodeling upon inflammatory responses. Immunofluorescence staining of human skin sections and whole-mount skin blocks confirmed the discrete expression of these markers along the blood vascular tree in situ, further corroborating BEC heterogeneity in human skin. Overall, our study molecularly refines individual BV compartments, whilst the identification of novel subtype-specific signatures provides more insights for future studies dissecting the responses of distinct vessel segments under pathological conditions.

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

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

Antigen presentation and interferon signatures in B cells driven by localized ablative cancer immunotherapy correlate with extended survival

Rationale: B cells have emerged as key regulators in protective cancer immunity. However, the activation pathways induced in B cells during effective immunotherapy are not well understood. Methods: We used a novel localized ablative immunotherapy (LAIT), combining photothermal therapy (PTT) with intra-tumor delivery of the immunostimulant N-dihydrogalactochitosan (GC), to treat mice bearing mouse mammary tumor virus-polyoma middle tumor-antigen (MMTV-PyMT). We used single-cell RNA sequencing to compare the transcriptional changes induced by PTT, GC and PTT+GC in B cells within the tumor microenvironment (TME). Results: LAIT significantly increased survival in the tumor-bearing mice, compared to the treatment by PTT and GC alone. We found that PTT, GC and PTT+GC increased the proportion of tumor-infiltrating B cells and induced gene expression signatures associated with B cell activation. Both GC and PTT+GC elevated gene expression associated with antigen presentation, whereas GC elevated transcripts that regulate B cell activation and GTPase function and PTT+GC induced interferon response genes. Trajectory analysis, where B cells were organized according to pseudotime progression, revealed that both GC and PTT+GC induced the differentiation of B cells from a resting state towards an effector phenotype. The analyses confirmed upregulated interferon signatures in the differentiated tumor-infiltrating B cells following treatment by PTT+GC but not by GC. We also observed that breast cancer patients had significantly longer survival time if they had elevated expression of genes in B cells that were induced by PTT+GC therapy in the mouse tumors. Conclusion: Our findings show that the combination of local ablation and local application of immunostimulant initiates the activation of interferon signatures and antigen-presentation in B cells which is associated with positive clinical outcomes for breast cancer. These findings broaden our understanding of LAIT's regulatory roles in remodeling TME and shed light on the potentials of B cell activation in clinical applications.

Single-Cell RNA Sequencing of Human Corpus Cavernosum Reveals Cellular Heterogeneity Landscapes in Erectile Dysfunction

PURPOSE

To assess the diverse cell populations of human corpus cavernosum in patients with severe erectile dysfunction (ED) at the single-cell level.

METHODS

Penile tissues collected from three patients were subjected to single-cell RNA sequencing using the BD Rhapsody™ platform. Common bioinformatics tools were used to analyze cellular heterogeneity and gene expression profiles from generated raw data, including the packages Seurat, Monocle, and CellPhoneDB.

RESULTS

Disease-related heterogeneity of cell types was determined in the cavernous tissue such as endothelial cells (ECs), smooth muscle cells, fibroblasts, and immune cells. Reclustering analysis of ECs identified an arteriole ECs subcluster and another one with gene signatures of fibroblasts. The proportion of fibroblasts was higher than the other cell populations and had the most significant cellular heterogeneity, in which a distinct subcluster co-expressed endothelial markers. The transition trajectory of differentiation from smooth muscle cells into fibroblasts was depicted using the pseudotime analysis, suggesting that the expansion of corpus cavernosum is possibly compromised as a result of fibrosis. Cell-cell communications among ECs, smooth muscle cells, fibroblasts, and macrophages were robust, which indicated that inflammation may also have a crucial role in the development of ED.

CONCLUSIONS

Our study has demonstrated a comprehensive single-cell atlas of cellular components in human corpus cavernosum of ED, providing in-depth insights into the pathogenesis. Future research is warranted to explore disease-specific alterations for individualized treatment of ED.