Neutrophil trapping and nexocytosis, mast cell-mediated processes for inflammatory signal relay.

Neutrophils are sentinel immune cells with essential roles for antimicrobial defense. Most of our knowledge on neutrophil tissue navigation derived from wounding and infection models, whereas allergic conditions remained largely neglected. Here, we analyzed allergen-challenged mouse tissues and discovered that degranulating mast cells (MCs) trap living neutrophils inside them. MCs release the attractant leukotriene B4 to re-route neutrophils toward them, thus exploiting a chemotactic system that neutrophils normally use for intercellular communication. After MC intracellular trap (MIT) formation, neutrophils die, but their undigested material remains inside MC vacuoles over days. MCs benefit from MIT formation, increasing their functional and metabolic fitness. Additionally, they are more pro-inflammatory and can exocytose active neutrophilic compounds with a time delay (nexocytosis), eliciting a type 1 interferon response in surrounding macrophages. Together, our study highlights neutrophil trapping and nexocytosis as MC-mediated processes, which may relay neutrophilic features over the course of chronic allergic inflammation.

Effective mechanical potential of cell-cell interaction in tissues harboring cavity and in cell sheet toward morphogenesis.

Measuring mechanical forces of cell-cell interactions is important for studying morphogenesis in multicellular organisms. We previously reported an image-based statistical method for inferring effective mechanical potentials of pairwise cell-cell interactions by fitting cell tracking data with a theoretical model. However, whether this method is applicable to tissues with non-cellular components such as cavities remains elusive. Here we evaluated the applicability of the method to cavity-harboring tissues. Using synthetic data generated by simulations, we found that the effect of expanding cavities was added to the pregiven potentials used in the simulations, resulting in the inferred effective potentials having an additional repulsive component derived from the expanding cavities. Interestingly, simulations by using the effective potentials reproduced the cavity-harboring structures. Then, we applied our method to the mouse blastocysts, and found that the inferred effective potentials can reproduce the cavity-harboring structures. Pairwise potentials with additional repulsive components were also detected in two-dimensional cell sheets, by which curved sheets including tubes and cups were simulated. We conclude that our inference method is applicable to tissues harboring cavities and cell sheets, and the resultant effective potentials are useful to simulate the morphologies.

Ischemia Reperfusion Responses in Human Lung Transplants at the Single Cell Resolution.

Ischemia reperfusion is an unavoidable step of organ transplantation. Development of therapeutics for lung injury during transplantation has proved challenging; understanding lung injury from human data at the single cell resolution is required to accelerate the development of therapeutics. Donor lung biopsies from six human lung transplant cases were collected at the end of cold preservation and 2-hour reperfusion and underwent single cell RNA sequencing. Donor and recipient origin of cells from the reperfusion timepoint were deconvolved. Gene expression profiles were (1) compared between each donor cell type between timepoints and (2) compared between donor and recipient cells. Inflammatory responses from donor lung macrophages were found after reperfusion with upregulation of multiple cytokines and chemokines, especially IL-1ß and IL-1α. Significant inflammatory responses were found in alveolar epithelial cells (featured by CXCL8) and lung endothelial cells (featured by IL-6 upregulation). Different inflammatory responses were noted between donor and recipient monocytes and CD8+ T cells. The inflammatory signals and differences between donor and recipient cells observed provide insight into the cellular and molecular mechanisms of ischemia reperfusion induced lung injury. Further investigations may lead to the development of novel targeted therapeutics.

Cell-cell interactions and communication dynamics in lung fibrosis.

Cell-cell interactions are essential components of coordinated cell function in lung homeostasis. Lung diseases involve altered cell-cell interactions and communication between different cell types, as well as between subsets of cells of the same type. The identification and understanding of intercellular signaling in lung fibrosis offer insights into the molecular mechanisms underlying these interactions and their implications in the development and progression of lung fibrosis. A comprehensive cell atlas of the human lung, established with the facilitation of single-cell RNA transcriptomic analysis, has enabled the inference of intercellular communications using ligand-receptor databases. In this review, we provide a comprehensive overview of the modified cell-cell communications in lung fibrosis. We highlight the intricate interactions among the major cell types within the lung and their contributions to fibrogenesis. The insights presented in this review will contribute to a better understanding of the molecular mechanisms underlying lung fibrosis and may guide future research efforts in developing targeted therapies for this debilitating disease.

Single-cell RNA-Seq reveals the heterogeneity of fibroblasts within the tympanojugular paraganglioma microenvironment.

Tympanojugular paragangliomas (TJP) originate from the parasympathetic ganglia in the lateral base of the skull. Although the cellular composition and oncogenic mechanisms of paragangliomas have been evaluated, a comprehensive transcriptomic atlas specific to TJP remains to be established to facilitate further investigations. In this study, single-cell RNA sequencing and whole-exome sequencing were conducted on six surgically excised TJP samples to determine their cellular composition and intratumoral heterogeneity. Fibroblasts were sub-classified into two distinct groups: myofibroblasts and fibroblasts associated with bone remodeling. Additionally, an elaborate regulatory and cell-cell communication network was determined, highlighting the multifaceted role of fibroblasts, which varies depending on expression transitions. The Kit receptor (KIT) signaling pathway mediated interactions between fibroblasts and mast cells, whereas robust connections with endothelial and Schwann cell-like cells were facilitated through the platelet-derived growth factor signaling pathway. These findings establish a foundation for studying the mechanisms underlying protumor angiogenesis and the specific contributions of fibroblasts within the TJP microenvironment. IL6 signaling pathway of fibroblasts interacting with macrophages and endothelial cells may be involved in tumor regrowth. These results enhance our understanding of fibroblast functionality and provide a resource for future therapeutic targeting of TJP.

Elucidating bacterial coaggregation through a physicochemical and imaging surface characterization.

Bacterial coaggregation is a highly specific type of cell-cell interaction, well-documented among oral bacteria, and involves specific characteristics of the cell surface of the coaggregating strains. However, the understanding of the mechanisms promoting coaggregation in aquatic systems remains limited. This gap is critical to address, given the broad implications of coaggregation for multispecies biofilm formation, water quality, the performance of engineered systems, and diverse biotechnological applications. Therefore, this study aims to comprehensively characterize the cell surface of the coaggregating strain Delftia acidovorans 005P, isolated from drinking water, alongside a non-coaggregating strain, D. acidovorans 009P. By analyzing two strains of the same species, we aim to identify the factors contributing to the coaggregation ability of strain 005P. To achieve this, we employed a combination of physicochemical characterization, Fourier-transform infrared spectroscopy (FTIR), and advancing imaging techniques [transmission electron microscopy and cryo-electron tomography (cryo-ET)]. The coaggregating strain (005P) exhibited higher surface hydrophobicity, negative surface charge, and cell surface and co-adhesion energies than the non-coaggregating strain (009P). The chemical characterization of bacterial surfaces through FTIR revealed subtle differences, particularly in spectral regions linked to carbohydrates and phosphodiesters/amide III of proteins (860-930 cm-1 and 1212-1240 cm-1, respectively). Cryo-ET highlighted significant differences in pili structures between the strains, such as variations in length, frequency, and arrangement. The pili in the 005P strain, identified as pili-like adhesins, serve as key mediators of coaggregation. By integrating physicochemical analyses and high-resolution imaging techniques, this study conclusively links the coaggregation ability of D. acidovorans 005P to its unique pili characteristics, emphasizing their crucial role in microbial coaggregation in aquatic environments.

A 3D micro-printed single cell micro-niche with asymmetric niche signals - An in vitro model for asymmetric cell division study.

Intricate microenvironment signals orchestrate to affect cell behavior and fate during tissue morphogenesis. However, the underlying mechanisms on how specific local niche signals influence cell behavior and fate are not fully understood, owing to the lack of in vitro platform able to precisely, quantitatively, spatially, and independently manipulate individual niche signals. Here, microarrays of protein-based 3D single cell micro-niche (3D-SCµN), with precisely engineered biophysical and biochemical niche signals, are micro-printed by a multiphoton microfabrication and micropatterning technology. Mouse embryonic stem cell (mESC) is used as the model cell to study how local niche signals affect stem cell behavior and fate. By precisely engineering the internal microstructures of the 3D SCµNs, we demonstrate that the cell division direction can be controlled by the biophysical niche signals, in a cell shape-independent manner. After confining the cell division direction to a dominating axis, single mESCs are exposed to asymmetric biochemical niche signals, specifically, cell-cell adhesion molecule on one side and extracellular matrix on the other side. We demonstrate that, symmetry-breaking (asymmetric) niche signals successfully trigger cell polarity formation and bias the orientation of asymmetric cell division, the mitosis process resulting in two daughter cells with differential fates, in mESCs.

T cell and airway smooth muscle interaction: a key driver of asthmatic airway inflammation and remodelling.

Crosstalk between T cells and airway smooth muscle (ASM) may play a role in modulating asthmatic airway inflammation and remodelling. Infiltrating T cells have been observed within the ASM bundles of asthmatics, and a wide range of direct and indirect interactions between T cells and ASM have been demonstrated using various in vitro and in vivo model systems. Contact-dependent mechanisms such as ligation and activation of cellular adhesion and costimulatory molecules, as well as the formation of lymphocyte-derived membrane conduits, facilitate the adhesion, bidirectional communication and transfer of materials between T and ASM cells. T cell-derived cytokines, particularly of the Th1, Th2 and Th17 subsets, modulate the secretome, proliferation and contractility of ASM cells. This review summarizes the mechanisms governing T cell-ASM crosstalk in the context of asthma. Understanding the underlying mechanistic basis is important for directing future research and developing therapeutic interventions targeted towards this complex interaction.

Invasion of the stigma by oomycete pathogenic hyphae or pollen tubes : striking similarities and differences.

Both filamentous pathogens' hyphae and pollen tube penetrate the host's outer layer and involve growth within the host tissues. Early epidermal responses are decisive for the outcome of these two-cell interaction processes. We identified a single cell type, the papilla of Arabidospis thaliana's stigma, as a tool to conduct a comprehensive comparative analysis on how an epidermal cell responds to the invasion of an unwanted pathogen or a welcomed pollen tube. We showed that Phytophtora parasitica, a root oomycete, effectively breaches the stigmatic cell wall and develops as a biotroph within the papilla cytoplasm. These invasive features resemble the behaviour exhibited by the pathogen within its natural host cells, but diverge from the manner in which the pollen tube progresses, being engulfed within the papilla cell wall. Quantitative analysis revealed that both invaders trigger reorganisation of the stigmatic endomembrane system and the actin cytoskeleton. While some remodelling processes are shared between the two interactions, others appear more specific towards the respective invader. These findings underscore the remarkable ability of an epidermal cell to differentiate between two types of invaders, thereby enabling it to trigger the most suitable response during the onset of invasion.

Mapping Single-Cell Transcriptomes of Endometrium Reveals Potential Biomarkers in Endometrial Cancer.

Background: The heterogeneity and dynamic changes of endometrial cells have a significant impact on health as they determine the normal function of the endometrium during the menstrual cycle. Dysfunction of the endometrium can lead to the occurrence of various gynecological diseases. Therefore, deconvolution of immune microenvironment that drives transcriptional programs throughout the menstrual cycle is key to understand regulatory biology of endometrium. Methods: Herein, we comprehensively analyzed single-cell transcriptome of 59,397 cells across ten human endometrium samples and revealed the dynamic cellular heterogeneity throughout the menstrual cycle. Results: We identified two perivascular cell subtypes, four epithelial subtypes and four fibroblast cell types in endometrium. Moreover, we inferred the cell type-specific transcription factor (TF) activities and linked critical TFs to transcriptional output of diverse immune cell types, highlighting the importance of transcriptional regulation in endometrium. Dynamic interactions between various types of cells in endometrium contribute to a range of biological pathways regulating differentiation of secretory. Integration of the molecular biomarkers identified in endometrium and bulk transcriptome of 535 endometrial cancers (EC), we revealed five RNA-based molecular subtypes of EC with highly intratumoral heterogeneity and different clinical manifestations. Mechanism analysis uncovered clinically relevant pathways for pathogenesis of EC. Conclusion: In summary, our results revealed the dynamic immune microenvironment of endometrium and provided novel insights into future development of RNA-based treatments for endometriosis and endometrial carcinoma.

A global view of altered ligand-receptor interactions in bone marrow aging based on single-cell sequencing.

Altered cell-cell communication is a hallmark of aging, but its impact on bone marrow aging remains poorly understood. Based on a common and effective pipeline and single-cell transcriptome sequencing, we detected 384,124 interactions including 2575 ligand-receptor pairs and 16 non-adherent bone marrow cell types in old and young mouse and identified a total of 5560 significantly different interactions, which were then verified by flow cytometry and quantitative real-time PCR. These differential ligand-receptor interactions exhibited enrichment for the senescence-associated secretory phenotypes. Further validation demonstrated supplementing specific extracellular ligands could modify the senescent signs of hematopoietic stem cells derived from old mouse. Our work provides an effective procedure to detect the ligand-receptor interactions based on single-cell sequencing, which contributes to understand mechanisms and provides a potential strategy for intervention of bone marrow aging.

Endothelial and macrophage interactions in the angiogenic niche.

The interactions between vascular cells, especially endothelial cells, and macrophages play a pivotal role in maintaining the subtle balance of vascular biology, which is crucial for angiogenesis in both healthy and diseased states. These cells are central to ensuring a harmonious balance between tissue repair and preventing excessive angiogenic activity, which could lead to pathological conditions. Recent advances in sophisticated genetic engineering vivo models and novel sequencing approaches, such as single-cell RNA-sequencing, in immunobiology have significantly enhanced our understanding of the gene expression and behavior of macrophages. These insights offer new perspectives on the role macrophages play not only in development but also across various health conditions. In this review, we explore the complex interactions between multiple types of macrophages and endothelium, focusing on their impact on new blood vessel formation. By understanding these intricate interactions, we aim to provide insights into new methods for managing angiogenesis in various diseases, thereby offering hope for the development of novel therapeutic approaches.

Minimal cellular automaton model with heterogeneous cell sizes predicts epithelial colony growth.

Regulation of cell proliferation is a crucial aspect of tissue development and homeostasis and plays a major role in morphogenesis, wound healing, and tumor invasion. A phenomenon of such regulation is contact inhibition, which describes the dramatic slowing of proliferation, cell migration and individual cell growth when multiple cells are in contact with each other. While many physiological, molecular and genetic factors are known, the mechanism of contact inhibition is still not fully understood. In particular, the relevance of cellular signaling due to interfacial contact for contact inhibition is still debated. Cellular automata (CA) have been employed in the past as numerically efficient mathematical models to study the dynamics of cell ensembles, but they are not suitable to explore the origins of contact inhibition as such agent-based models assume fixed cell sizes. We develop a minimal, data-driven model to simulate the dynamics of planar cell cultures by extending a probabilistic CA to incorporate size changes of individual cells during growth and cell division. We successfully apply this model to previous in-vitro experiments on contact inhibition in epithelial tissue: After a systematic calibration of the model parameters to measurements of single-cell dynamics, our CA model quantitatively reproduces independent measurements of emergent, culture-wide features, like colony size, cell density and collective cell migration. In particular, the dynamics of the CA model also exhibit the transition from a low-density confluent regime to a stationary postconfluent regime with a rapid decrease in cell size and motion. This implies that the volume exclusion principle, a mechanical constraint which is the only inter-cellular interaction incorporated in the model, paired with a size-dependent proliferation rate is sufficient to generate the observed contact inhibition. We discuss how our approach enables the introduction of effective bio-mechanical interactions in a CA framework for future studies.

Tissue module discovery in single-cell-resolution spatial transcriptomics data via cell-cell interaction-aware cell embedding.

Computational methods are desired for single-cell-resolution spatial transcriptomics (ST) data analysis to uncover spatial organization principles for how individual cells exert tissue-specific functions. Here, we present ST data analysis via interaction-aware cell embedding (SPACE), a deep-learning method for cell-type identification and tissue module discovery from single-cell-resolution ST data by learning a cell representation that captures its gene expression profile and interactions with its spatial neighbors. SPACE identified spatially informed cell subtypes defined by their special spatial distribution patterns and distinct proximal-interacting cell types. SPACE also automatically discovered "cell communities"-tissue modules with discernible boundaries and a uniform spatial distribution of constituent cell types. For each cell community, SPACE outputs a characteristic proximal cell-cell interaction network associated with physiological processes, which can be used to refine ligand-receptor-based intercellular signaling analyses. We envision that SPACE can be used in large-scale ST projects to understand how proximal cell-cell interactions contribute to emergent biological functions within cell communities. A record of this paper's transparent peer review process is included in the supplemental information.

Exploring coaggregation mechanisms involved in biofilm formation in drinking water through a proteomic-based approach.

AIM: Coaggregation, a highly specific cell-cell interaction mechanism, plays a pivotal role in multispecies biofilm formation. While it has been mostly studied in oral environments, its occurrence in aquatic systems is also acknowledged. Considering biofilm formation's economic and health-related implications in engineered water systems, it is crucial to understand its mechanisms. Here, we hypothesized that traceable differences at the proteome level might determine coaggregation ability. METHODS AND RESULTS: Two strains of Delftia acidovorans, isolated from drinking water were studied. First, in vitro motility assays indicated more swarming and twitching motility for the coaggregating strain (C+) than non-coaggregating strain (C-). By transmission electronic microscopy, we confirmed the presence of flagella for both strains. By proteomics, we detected a significantly higher expression of type IV pilus twitching motility proteins in C+, in line with the motility assays. Moreover, flagellum ring proteins were more abundant in C+, while those involved in the formation of the flagellar hook (FlE and FilG) were only detected in C-. All the results combined suggested structural and conformational differences between stains in their cell appendages. CONCLUSION: This study presents an alternative approach for identifying protein biomarkers to detect coaggregation abilities in uncharacterized strains.

Intratumor transforming growth factor-ß signaling with extracellular matrix-related gene regulation marks chemotherapy-resistant gastric cancer.

Drug-tolerant persister (DTP) cells remain following chemotherapy and can cause cancer relapse. However, it is unclear when acquired resistance to chemotherapy emerges. Here, we compared the gene expression profiles of gastric cancer patient-derived cells (GC PDCs) and their respective xenograft tumors with different sensitivities to 5-fluorouracil (5-FU) by using immunodeficient female BALB/c-nu mice. RNA sequencing analysis of 5-FU-treated PDCs demonstrated that DNA replication/cell cycle-related genes were transiently induced in the earlier phase of DTP cell emergence, while extracellular matrix (ECM)-related genes were sustainably upregulated during long-term cell survival in 5-FU-resistant residual tumors. NicheNet analysis, which uncovers cell-cell signal interactions, indicated the transforming growth factor-ß (TGF-ß) pathway as the upstream regulator in response to 5-FU treatment. This induced ECM-related gene expression in the 5-FU-resistant tumor model. In the 5-FU-resistant residual tumors, there was a marked upregulation of cancer cell-derived TGF-ß1 expression and increased phosphorylation of SMAD3, a downstream regulator of the TGF-ß receptor. By contrast, these responses were not observed in a 5-FU-sensitive tumor model. We further found that TGF-ß-related upregulation of ECM genes was preferentially observed in non-responders to chemotherapy with 5-FU and/or oxaliplatin among 22 patient-derived xenograft tumors. These observations suggest that chemotherapy-induced activation of the TGF-ß1/SMAD3/ECM-related gene axis is a potential biomarker for the emergence of drug resistance in GCs.

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.

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.

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.

Fibroblast regulates angiogenesis in assembled oral cancer organoid: A possible role of NNMT.

OBJECTIVE: Tumour angiogenesis is affected by various cell types in the tumour microenvironment (TME), including cancer cells and cancer-associated fibroblasts (CAFs). Here, an assembled organoid model was generated to investigate the mechanism by which the TME regulates angiogenesis in oral squamous cell carcinoma (OSCC). MATERIALS AND METHODS: Secretion of vascular endothelial growth factor-A (VEGFA) was analysed to compare the proangiogenic properties of OSCC cells and corresponding CAFs. Cell aggregates consisting of endothelial cells (ECs), CAFs and cancer cells were generated to construct assembled organoids. Nicotinamide N-methyltransferase (NNMT) was pharmacologically or genetically inhibited to block the activation of CAFs. ATAC-seq was employed to test the transcriptional network of fibroblasts overexpressing NNMT. RESULTS: Compared with cancer cells, CAFs secreted more VEGFA. Coculture with CAFs more effectively promoted the sprouting of ECs. Blockade of CAF activation via inhibition of NNMT drastically reduced the expression of CD31 in the assembled organoids. Overexpression of NNMT enhanced the transcription of genes related to angiogenesis in fibroblasts. Specifically, NNMT orchestrated the enrichment of the transcription factor JUNB at the promoter of VEGFA. CONCLUSIONS: We clarify that stromal NNMT enables the steady reproduction of angiogenesis in assembled oral cancer organoids, providing a novel target for exploiting antiangiogenic therapy.