Revealing cell-cell communication pathways with their spatially coupled gene programs.

Inference of cell-cell communication (CCC) provides valuable information in understanding the mechanisms of many important life processes. With the rise of spatial transcriptomics in recent years, many methods have emerged to predict CCCs using spatial information of cells. However, most existing methods only describe CCCs based on ligand-receptor interactions, but lack the exploration of their upstream/downstream pathways. In this paper, we proposed a new method to infer CCCs, called Intercellular Gene Association Network (IGAN). Specifically, it is for the first time that we can estimate the gene associations/network between two specific single spatially adjacent cells. By using the IGAN method, we can not only infer CCCs in an accurate manner, but also explore the upstream/downstream pathways of ligands/receptors from the network perspective, which are actually exhibited as a new panoramic cell-interaction-pathway graph, and thus provide extensive information for the regulatory mechanisms behind CCCs. In addition, IGAN can measure the CCC activity at single cell/spot resolution, and help to discover the CCC spatial heterogeneity. Interestingly, we found that CCC patterns from IGAN are highly consistent with the spatial microenvironment patterns for each cell type, which further indicated the accuracy of our method. Analyses on several public datasets validated the advantages of IGAN.

Identification of ligand and receptor interactions in CKD and MASH through the integration of single cell and spatial transcriptomics.

BACKGROUND: Chronic Kidney Disease (CKD) and Metabolic dysfunction-associated steatohepatitis (MASH) are metabolic fibroinflammatory diseases. Combining single-cell (scRNAseq) and spatial transcriptomics (ST) could give unprecedented molecular disease understanding at single-cell resolution. A more comprehensive analysis of the cell-specific ligand-receptor (L-R) interactions could provide pivotal information about signaling pathways in CKD and MASH. To achieve this, we created an integrative analysis framework in CKD and MASH from two available human cohorts. RESULTS: The analytical framework identified L-R pairs involved in cellular crosstalk in CKD and MASH. Interactions between cell types identified using scRNAseq data were validated by checking the spatial co-presence using the ST data and the co-expression of the communicating targets. Multiple L-R protein pairs identified are known key players in CKD and MASH, while others are novel potential targets previously observed only in animal models. CONCLUSION: Our study highlights the importance of integrating different modalities of transcriptomic data for a better understanding of the molecular mechanisms. The combination of single-cell resolution from scRNAseq data, combined with tissue slide investigations and visualization of cell-cell interactions obtained through ST, paves the way for the identification of future potential therapeutic targets and developing effective therapies.

Integrating single-cell and spatial analysis reveals MUC1-mediated cellular crosstalk in mucinous colorectal adenocarcinoma.

BACKGROUND: Mucinous colorectal adenocarcinoma (MCA) is a distinct subtype of colorectal cancer (CRC) with the most aggressive pattern, but effective treatment of MCA remains a challenge due to its vague pathological characteristics. An in-depth understanding of transcriptional dynamics at the cellular level is critical for developing specialised MCA treatment strategies. METHODS: We integrated single-cell RNA sequencing and spatial transcriptomics data to systematically profile the MCA tumor microenvironment (TME), particularly the interactome of stromal and immune cells. In addition, a three-dimensional bioprinting technique, canonical ex vivo co-culture system, and immunofluorescence staining were further applied to validate the cellular communication networks within the TME. RESULTS: This study identified the crucial intercellular interactions that engaged in MCA pathogenesis. We found the increased infiltration of FGF7+/THBS1+ myofibroblasts in MCA tissues with decreased expression of genes associated with leukocyte-mediated immunity and T cell activation, suggesting a crucial role of these cells in regulating the immunosuppressive TME. In addition, MS4A4A+ macrophages that exhibit M2-phenotype were enriched in the tumoral niche and high expression of MS4A4A+ was associated with poor prognosis in the cohort data. The ligand-receptor-based intercellular communication analysis revealed the tight interaction of MUC1+ malignant cells and ZEB1+ endothelial cells, providing mechanistic information for MCA angiogenesis and molecular targets for subsequent translational applications. CONCLUSIONS: Our study provides novel insights into communications among tumour cells with stromal and immune cells that are significantly enriched in the TME during MCA progression, presenting potential prognostic biomarkers and therapeutic strategies for MCA. KEY POINTS: Tumour microenvironment profiling of MCA is developed. MUC1+ tumour cells interplay with FGF7+/THBS1+ myofibroblasts to promote MCA development. MS4A4A+ macrophages exhibit M2 phenotype in MCA. ZEB1+ endotheliocytes engage in EndMT process in MCA.

A new integrative analysis of histopathology and single cell RNA-seq reveals the CCL5 mediated T and NK cell interaction with vascular cells in idiopathic pulmonary arterial hypertension.

BACKGROUND: Inflammation and dysregulated immunity play vital roles in idiopathic pulmonary arterial hypertension (IPAH), while the mechanisms that initiate and promote these processes are unclear. METHODS: Transcriptomic data of lung tissues from IPAH patients and controls were obtained from the Gene Expression Omnibus database. Weighted gene co-expression network analysis (WGCNA), differential expression analysis, protein-protein interaction (PPI) and functional enrichment analysis were combined with a hemodynamically-related histopathological score to identify inflammation-associated hub genes in IPAH. The monocrotaline-induced rat model of pulmonary hypertension was utilized to confirm the expression pattern of these hub genes. Single-cell RNA-sequencing (scRNA-seq) data were used to identify the hub gene-expressing cell types and their intercellular interactions. RESULTS: Through an extensive bioinformatics analysis, CXCL9, CCL5, GZMA and GZMK were identified as hub genes that distinguished IPAH patients from controls. Among these genes, pulmonary expression levels of Cxcl9, Ccl5 and Gzma were elevated in monocrotaline-exposed rats. Further investigation revealed that only CCL5 and GZMA were highly expressed in T and NK cells, where CCL5 mediated T and NK cell interaction with endothelial cells, smooth muscle cells, and fibroblasts through multiple receptors. CONCLUSIONS: Our study identified a new inflammatory pathway in IPAH, where T and NK cells drove heightened inflammation predominantly via the upregulation of CCL5, providing groundwork for the development of targeted therapeutics.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

SEGCECO: Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication.

Recent advances in single-cell RNA sequencing technology have eased analyses of signaling networks of cells. Recently, cell-cell interaction has been studied based on various link prediction approaches on graph-structured data. These approaches have assumptions about the likelihood of node interaction, thus showing high performance for only some specific networks. Subgraph-based methods have solved this problem and outperformed other approaches by extracting local subgraphs from a given network. In this work, we present a novel method, called Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication (SEGCECO), which uses an attributed graph convolutional neural network to predict cell-cell communication from single-cell RNA-seq data. SEGCECO captures the latent and explicit attributes of undirected, attributed graphs constructed from the gene expression profile of individual cells. High-dimensional and sparse single-cell RNA-seq data make converting the data into a graphical format a daunting task. We successfully overcome this limitation by applying SoptSC, a similarity-based optimization method in which the cell-cell communication network is built using a cell-cell similarity matrix which is learned from gene expression data. We performed experiments on six datasets extracted from the human and mouse pancreas tissue. Our comparative analysis shows that SEGCECO outperforms latent feature-based approaches, and the state-of-the-art method for link prediction, WLNM, with 0.99 ROC and 99% prediction accuracy. The datasets can be found at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84133 and the code is publicly available at Github https://github.com/sheenahora/SEGCECO and Code Ocean https://codeocean.com/capsule/8244724/tree.

Dynamic YAP expression in the non-parenchymal liver cell compartment controls heterologous cell communication.

INTRODUCTION: The Hippo pathway and its transcriptional effectors yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) are targets for cancer therapy. It is important to determine if the activation of one factor compensates for the inhibition of the other. Moreover, it is unknown if YAP/TAZ-directed perturbation affects cell-cell communication of non-malignant liver cells. MATERIALS AND METHODS: To investigate liver-specific phenotypes caused by YAP and TAZ inactivation, we generated mice with hepatocyte (HC) and biliary epithelial cell (BEC)-specific deletions for both factors (YAPKO, TAZKO and double knock-out (DKO)). Immunohistochemistry, single-cell sequencing, and proteomics were used to analyze liver tissues and serum. RESULTS: The loss of BECs, liver fibrosis, and necrosis characterized livers from YAPKO and DKO mice. This phenotype was weakened in DKO tissues compared to specimens from YAPKO animals. After depletion of YAP in HCs and BECs, YAP expression was induced in non-parenchymal cells (NPCs) in a cholestasis-independent manner. YAP positivity was detected in subgroups of Kupffer cells (KCs) and endothelial cells (ECs). The secretion of pro-inflammatory chemokines and cytokines such as C-X-C motif chemokine ligand 11 (CXCL11), fms-related receptor tyrosine kinase 3 ligand (FLT3L), and soluble intercellular adhesion molecule-1 (ICAM1) was increased in the serum of YAPKO animals. YAP activation in NPCs could contribute to inflammation via TEA domain transcription factor (TEAD)-dependent transcriptional regulation of secreted factors. CONCLUSION: YAP inactivation in HCs and BECs causes liver damage, and concomitant TAZ deletion does not enhance but reduces this phenotype. Additionally, we present a new mechanism by which YAP contributes to cell-cell communication originating from NPCs.

Pancreatic cancer environment: from patient-derived models to single-cell omics.

Pancreatic cancer (PC) is a highly malignant cancer characterized by poor prognosis, high heterogeneity, and intricate heterocellular systems. Selecting an appropriate experimental model for studying its progression and treatment is crucial. Patient-derived models provide a more accurate representation of tumor heterogeneity and complexity compared to cell line-derived models. This review initially presents relevant patient-derived models, including patient-derived xenografts (PDXs), patient-derived organoids (PDOs), and patient-derived explants (PDEs), which are essential for studying cell communication and pancreatic cancer progression. We have emphasized the utilization of these models in comprehending intricate intercellular communication, drug responsiveness, mechanisms underlying tumor growth, expediting drug discovery, and enabling personalized medical approaches. Additionally, we have comprehensively summarized single-cell analyses of these models to enhance comprehension of intercellular communication among tumor cells, drug response mechanisms, and individual patient sensitivities.

The diversification of methods for studying cell-cell interactions and communication.

No cell lives in a vacuum, and the molecular interactions between cells define most phenotypes. Transcriptomics provides rich information to infer cell-cell interactions and communication, thus accelerating the discovery of the roles of cells within their communities. Such research relies heavily on algorithms that infer which cells are interacting and the ligands and receptors involved. Specific pressures on different research niches are driving the evolution of next-generation computational tools, enabling new conceptual opportunities and technological advances. More sophisticated algorithms now account for the heterogeneity and spatial organization of cells, multiple ligand types and intracellular signalling events, and enable the use of larger and more complex datasets, including single-cell and spatial transcriptomics. Similarly, new high-throughput experimental methods are increasing the number and resolution of interactions that can be analysed simultaneously. Here, we explore recent progress in cell-cell interaction research and highlight the diversification of the next generation of tools, which have yielded a rich ecosystem of tools for different applications and are enabling invaluable discoveries.

Comparative analysis of cell-cell communication at single-cell resolution.

Inference of cell-cell communication from single-cell RNA sequencing data is a powerful technique to uncover intercellular communication pathways, yet existing methods perform this analysis at the level of the cell type or cluster, discarding single-cell-level information. Here we present Scriabin, a flexible and scalable framework for comparative analysis of cell-cell communication at single-cell resolution that is performed without cell aggregation or downsampling. We use multiple published atlas-scale datasets, genetic perturbation screens and direct experimental validation to show that Scriabin accurately recovers expected cell-cell communication edges and identifies communication networks that can be obscured by agglomerative methods. Additionally, we use spatial transcriptomic data to show that Scriabin can uncover spatial features of interaction from dissociated data alone. Finally, we demonstrate applications to longitudinal datasets to follow communication pathways operating between timepoints. Our approach represents a broadly applicable strategy to reveal the full structure of niche-phenotype relationships in health and disease.

Identification of cytokine-induced cell communications by pan-cancer meta-analysis.

Cancer immune responses are complex cellular processes in which cytokine-receptor interactions play central roles in cancer development and response to therapy; dysregulated cytokine-receptor communication may lead to pathological processes, including cancer, autoimmune diseases, and cytokine storm; however, our knowledge regarding cytokine-mediated cell-cell communication (CCI) in different cancers remains limited. The present study presents a single-cell and pan-cancer-level transcriptomics integration of 41,900 cells across 25 cancer types. We developed a single-cell method to actively express 62 cytokine-receptor pairs to reveal stable cytokine-mediated cell communications involving 84 cytokines and receptors. The correlation between the sample-based CCI profile and the interactome analysis indicates multiple cytokine-receptor modules including TGFB1, IL16ST, IL15, and the PDGF family. Some isolated cytokine interactions, such as FN1-IL17RC, displayed diverse functions within over ten single-cell transcriptomics datasets. Further functional enrichment analysis revealed that the constructed cytokine-receptor interaction map is associated with the positive regulation of multiple immune response pathways. Using public TCGA pan-cancer mutational data, co-mutational analysis of the cytokines and receptors provided significant co-occurrence features, implying the existence of cooperative mechanisms. Analysis of 10,967 samples from 32 TCGA cancer types revealed that the 84 cytokine and receptor genes are significantly associated with clinical survival time. Interestingly, the tumor samples with mutations in any of the 84 cytokines and receptors have a substantially higher mutational burden, offering insights into antitumor immune regulation and response. Clinical cancer stage information revealed that tumor samples with mutations in any of the 84 cytokines and receptors stratify into earlier tumor stages, with unique cellular compositions and clinical outcomes. This study provides a comprehensive cytokine-receptor atlas of the cellular architecture in multiple cancers at the single-cell level.

Advances in Intercellular Communication Mediated by Exosomal ncRNAs in Cardiovascular Disease.

Cardiovascular diseases are a leading cause of worldwide mortality, and exosomes have recently gained attention as key mediators of intercellular communication in these diseases. Exosomes are double-layered lipid vesicles that can carry biomolecules such as miRNAs, lncRNAs, and circRNAs, and the content of exosomes is dependent on the cell they originated from. They can be involved in the pathophysiological processes of cardiovascular diseases and hold potential as diagnostic and monitoring tools. Exosomes mediate intercellular communication, stimulate or inhibit the activity of target cells, and affect myocardial hypertrophy, injury and infarction, ventricular remodeling, angiogenesis, and atherosclerosis. Exosomes can be released from various types of cells, including endothelial cells, smooth muscle cells, cardiomyocytes, fibroblasts, platelets, adipocytes, immune cells, and stem cells. In this review, we highlight the communication between different cell-derived exosomes and cardiovascular cells, with a focus on the roles of RNAs. This provides new insights for further exploring targeted therapies in the clinical management of cardiovascular diseases.

Exploration of the molecular mechanism of intercellular communication in paediatric neuroblastoma by single-cell sequencing.

Neuroblastoma (NB) is an embryonic tumour that originates in the sympathetic nervous system and occurs most often in infants and children under 2 years of age. Moreover, it is the most common extracranial solid tumour in children. Increasing studies suggest that intercellular communication within the tumour microenvironment is closely related to tumour development. This study aimed to construct a prognosis-related intercellular communication-associated genes model by single-cell sequencing and transcriptome sequencing to predict the prognosis of patients with NB for precise management. Single-cell data from patients with NB were downloaded from the gene expression omnibus database for comprehensive analysis. Furthermore, prognosis-related genes were screened in the TARGET database based on epithelial cell marker genes through a combination of Cox regression and Lasso regression analyses, using GSE62564 and GSE85047 for external validation. The patients' risk scores were calculated, followed by immune infiltration analysis, drug sensitivity analysis, and enrichment analysis of risk scores, which were conducted for the prognostic model. I used the Lasso regression feature selection algorithm to screen characteristic genes in NB and developed a 21-gene prognostic model. The risk scores were highly correlated with multiple immune cells and common anti-tumour drugs. Furthermore, the risk score was identified as an independent prognostic factor for NB. In this study, I constructed and validated a prognostic signature based on epithelial marker genes, which may provide useful information on the development and prognosis of NB.

Starvation induces changes in abundance and small RNA cargo of extracellular vesicles released from Plasmodium falciparum infected red blood cells.

The lethal malaria parasite Plasmodium falciparum needs to constantly respond and adapt to changes within the human host in order to survive and transmit. One such change is composed of nutritional limitation, which is augmented with increased parasite loads and intimately linked to severe disease development. Extracellular vesicles released from infected red blood cells have been proposed as important mediators of disease pathogenesis and intercellular communication but whether important for the parasite response to nutritional availability is unknown. Therefore, we investigated the abundance and small RNA cargo of extracellular vesicles released upon short-term nutritional starvation of P. falciparum in vitro cultures. We show that primarily ring-stage parasite cultures respond to glucose and amino acid deprivation with an increased release of extracellular vesicles. Small RNA sequencing of these extracellular vesicles further revealed human miRNAs and parasitic tRNA fragments as the main constituent biotypes. Short-term starvations led to alterations in the transcriptomic profile, most notably in terms of the over-represented biotypes. These data suggest a potential role for extracellular vesicles released from P. falciparum infected red blood cells in the response to nutritional perturbations, their potential as prognostic biomarkers and point towards an evolutionary conserved role among protozoan parasites.

Decoding functional cell-cell communication events by multi-view graph learning on spatial transcriptomics.

Cell-cell communication events (CEs) are mediated by multiple ligand-receptor (LR) pairs. Usually only a particular subset of CEs directly works for a specific downstream response in a particular microenvironment. We name them as functional communication events (FCEs) of the target responses. Decoding FCE-target gene relations is: important for understanding the mechanisms of many biological processes, but has been intractable due to the mixing of multiple factors and the lack of direct observations. We developed a method HoloNet for decoding FCEs using spatial transcriptomic data by integrating LR pairs, cell-type spatial distribution and downstream gene expression into a deep learning model. We modeled CEs as a multi-view network, developed an attention-based graph learning method to train the model for generating target gene expression with the CE networks, and decoded the FCEs for specific downstream genes by interpreting trained models. We applied HoloNet on three Visium datasets of breast cancer and liver cancer. The results detangled the multiple factors of FCEs by revealing how LR signals and cell types affect specific biological processes, and specified FCE-induced effects in each single cell. We conducted simulation experiments and showed that HoloNet is more reliable on LR prioritization in comparison with existing methods. HoloNet is a powerful tool to illustrate cell-cell communication landscapes and reveal vital FCEs that shape cellular phenotypes. HoloNet is available as a Python package at https://github.com/lhc17/HoloNet.

ESICCC as a systematic computational framework for evaluation, selection, and integration of cell-cell communication inference methods.

Cell-cell communication (CCC) is critical for determining cell fates and functions in multicellular organisms. With the advent of single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), an increasing number of CCC inference methods have been developed. Nevertheless, a thorough comparison of their performances is yet to be conducted. To fill this gap, we developed a systematic benchmark framework called ESICCC to evaluate 18 ligand-receptor (LR) inference methods and five ligand/receptor-target inference methods using a total of 116 data sets, including 15 ST data sets, 15 sets of cell line perturbation data, two sets of cell type-specific expression/proteomics data, and 84 sets of sampled or unsampled scRNA-seq data. We evaluated and compared the agreement, accuracy, robustness, and usability of these methods. Regarding accuracy evaluation, RNAMagnet, CellChat, and scSeqComm emerge as the three best-performing methods for intercellular ligand-receptor inference based on scRNA-seq data, whereas stMLnet and HoloNet are the best methods for predicting ligand/receptor-target regulation using ST data. To facilitate the practical applications, we provide a decision-tree-style guideline for users to easily choose best tools for their specific research concerns in CCC inference, and develop an ensemble pipeline CCCbank that enables versatile combinations of methods and databases. Moreover, our comparative results also uncover several critical influential factors for CCC inference, such as prior interaction information, ligand-receptor scoring algorithm, intracellular signaling complexity, and spatial relationship, which may be considered in the future studies to advance the development of new methodologies.

The role of non-coding RNAs in extracellular vesicles in breast cancer and their diagnostic implications.

Breast Cancer (BC) is the most common form of cancer worldwide, responsible for 25% of cancers in women. Whilst treatment is effective and often curative in early BC, metastatic disease is incurable, highlighting the need for early detection. Currently, early detection relies on invasive procedures, however recent studies have shown extracellular vesicles (EVs) obtained from liquid biopsies may have clinical utility. EVs transport diverse bioactive cargos throughout the body, play major roles in intercellular communication and, importantly, mirror their cell of origin. In cancer cells, EVs alter the behaviour of the tumour microenvironment (TME), forming a bridge of communication between cancerous and non-cancerous cells to alter all aspects of cancer progression, including the formation of a pre-metastatic niche. Through gene regulatory frameworks, non-coding RNAs (ncRNAs) modulate vital molecular and cellular processes and can act as both tumour suppressors and oncogenic drivers in various cancer types. EVs transport and protect ncRNAs, facilitating their use clinically as liquid biopsies for early BC detection. This review summarises current research surrounding ncRNAs and EVs within BC, focusing on their roles in cancer progression through bi-directional communication with the microenvironment and their diagnostic implications. The role of EV ncRNAs in breast cancer. A representation of the different EV ncRNAs involved in tumourigenic processes in breast cancer. Pro-tumourigenic ncRNAs displayed in green and ncRNAs which inhibit oncogenic processes are shown in red.

TimeTalk uses single-cell RNA-seq datasets to decipher cell-cell communication during early embryo development.

Early embryonic development is a dynamic process that relies on proper cell-cell communication to form a correctly patterned embryo. Early embryo development-related ligand-receptor pairs (eLRs) have been shown to guide cell fate decisions and morphogenesis. However, the scope of eLRs and their influence on early embryo development remain elusive. Here, we developed a computational framework named TimeTalk from integrated public time-course mouse scRNA-seq datasets to decipher the secret of eLRs. Extensive validations and analyses were performed to ensure the involvement of identified eLRs in early embryo development. Process analysis identified that eLRs could be divided into six temporal windows corresponding to sequential events in the early embryo development process. With the interpolation strategy, TimeTalk is powerful in revealing paracrine settings and studying cell-cell communication during early embryo development. Furthermore, by using TimeTalk in the blastocyst and blastoid models, we found that the blastoid models share the core communication pathways with the epiblast and primitive endoderm lineages in the blastocysts. This result suggests that TimeTalk has transferability to other bio-dynamic processes. We also curated eLRs recognized by TimeTalk, which may provide valuable clues for understanding early embryo development and relevant disorders.

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

BACKGROUND: The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells. RESULTS: We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults. CONCLUSIONS: Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

miRNA-126a plays important role in myoblast and endothelial cell interaction.

Muscle satellite cells (SCs) are stem cells and the main players in skeletal muscle reconstruction. Since satellite cells are located near or in direct contact with blood vessels their niche is formed, inter alia, by endothelial cells. The cross-talk between satellite cells and endothelial cells determines quiescence or proliferation of these cells. However, little is known about the role of miRNA in these interactions. In the present study we identified miRNA that were up-regulated in SC-derived myoblasts treated with stromal derived factor-1 (SDF-1) and/or down-regulated in cells in which the expression of CXCR4 or CXCR7, that is, SDF-1 receptors, was silenced. SDF-1 is one of the important regulators of cell migration, mobilization, skeletal muscle regeneration, and angiogenesis. We hypothesized that selected miRNAs affect SC-derived myoblast fate and interactions with endothelial cells. We showed that miR-126a-3p inhibited both, myoblast migration and fusion. Moreover, the levels of Cxcl12, encoding SDF-1 and Ackr3, encoding CXCR7, were reduced by miR-126a-3p mimic. Interestingly, the miR-126a-3p mimic significantly decreased the level of numerous factors involved in myogenesis and the miR-126a-5p mimic increased the level of Vefga. Importantly, the treatment of endothelial cells with medium conditioned by miR-126-5p mimic transfected SC-derived myoblasts promoted tubulogenesis.