Cell-cell Communication Networks Research Articles

Abstract 1587: Multimodal spatial transcriptomics uncover distinct tumor microenvironment states and cell-cell communication networks in molecular pancreatic cancer subtypes

Abstract The tumor microenvironment (TME) presents a complex ecosystem comprising a diversity of immune and stromal cell populations that influence tumor progression, drug delivery and therapy outcome. High levels of inter- and intra-tumor heterogeneity in TME state are driven by molecular tumor phenotypes. Pancreatic ductal adenocarcinoma (PDAC), one of the most lethal cancer types, displays a high genetic heterogeneity along with an immunosuppressive TME. However, it remains largely elusive how distinct TME states are mechanistically influenced, and which molecular processes dictate the emergence of different modes of immunosuppression within the TME of molecular PDAC subtypes.Here, we performed a systematic analysis of functional associations between the major molecular PDAC subtypes, the classical and mesenchymal subtype, and their TME states. We describe the TME composition and complex cell-cell communication networks within PDAC subtypes using a multimodal integration of spatial transcriptomics (ST) with complementary approaches, such as scRNA-seq, MS-based Secretomics and multiplexed histocytometry. We generated ST data sets (Visium) of a tumor cohort derived from a Kras-driven PDAC mouse model. To analyze TME composition as well as spatial niches and cell type communities, we computationally enhanced the spot-resolution of ST datasets, followed by cell type deconvolution and cell-cell communication analysis to identify subtype-specific TME communities and communication networks.To this end, we generated a large resource of PDAC mouse models which represent molecular subtypes of the disease and mimic the heterogeneity of TME states found in human PDAC cohorts. This analysis revealed that a set of subtype-specific secreted factors shape the immunosuppressive PDAC TME via direct and indirect communication networks with immunosuppressive myeloid and T cells in molecular PDAC subtypes. Multimodal ST analysis delineated spatial subtype-specific TME communities as well as spatial communication patterns. We functionally investigate the tumor cell intrinsic regulation of the identified subtype-specific secreted factors, providing potential therapeutic vulnerabilities for more effective combinatorial subtype-specific therapies including immunotherapeutic approaches. Citation Format: Stefanie Baerthel, Chiara Falcomatà, Vanessa Goelling, Daniele Lucarelli, Rushin Gindra, Constantin Schmitt, Jonathan Swietlik, Felix Meissner, Marc Schmidt-Supprian, Dieter Saur. Multimodal spatial transcriptomics uncover distinct tumor microenvironment states and cell-cell communication networks in molecular pancreatic cancer subtypes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1587.

Abstract 1152: Pan-cancer single-cell and subcellular characterization of lymphoid aggregates and tertiary lymphoid structures using high-plex spatial molecular imaging

Abstract Tertiary lymphoid structures (TLSs), clusters of immune cells that form in non-lymphoid tissues, are frequently identified in a wide range of solid tumors. They function as immunological “hubs”, driving anti-tumor responses. Previous studies have demonstrated that the presence of TLSs is often associated with improved responses to immunotherapy and superior clinical outcomes. Despite these insights, our understanding of their phenotypic diversity, spatial organization, intercellular communication, and interactions with cancer cells and the tumor microenvironment remains incomplete. Existing methods lack the resolution, depth, and scale required for a thorough understanding of these complex structures. We have obtained in-house and public spatial multi-omics data generated on various tissues across 9 different cancer types using the cutting-edge single-cell spatial transcriptomics and high-plex spatial molecular imaging platforms. We have developed innovative computational approaches for the detailed characterization of TLSs at single-cell and subcellular resolution and carried out an extensive analysis of lymphoid aggregates (LAs) and TLSs in tissue sections across these different solid cancer types. We defined the phenotypic state of each individual cell within TLSs including B/plasma cells, follicular helper/regulatory/stressed/exhausted T cells, follicular dendritic cells, fibroblasts and high endothelial venules. Then, we performed unbiased phenotyping and classification of TLSs based on their cellular states, compositions and organization patterns across various types of tissues. Additionally, we performed spatial neighborhood analysis, characterized the neighboring cells of each TLS and analyzed their cell-cell communication networks. Our approach allowed an unprecedented examination of the spatial phenotypes of TLSs in large, complex tissues at super resolution, in ways that were not previously described. We comprehensively characterized TLSs across different solid tumor types. Our analysis revealed a high degree of spatial and phenotypic heterogeneity in LAs and TLSs. We created a comprehensive spatial atlas of LAs and TLSs across solid tumor types, profiled TLSs at various maturation states with distinct cellular composition and organization patterns, and discovered new TLS states. Notably, we observed distinct communication networks associated with TLS phenotypic states. Our innovative methods and analytical approaches offer a more comprehensive way to study LAs and TLSs, enhancing detection capabilities and providing deeper insights into their spatial interaction and functional phenotypes. This research holds the potential to significantly advance our understanding of TLSs and their roles in immunobiology and immunotherapy responses across a broad range of cancers and diseases. Citation Format: Yunhe Liu, Yuanyuan Zhang, Guangchun Han, Guangsheng Pei, Kyung Serk Cho, Enyu Dai, Yanshuo Chu, Jiahui Jiang, Ansam Sinjab, Luisa M. Soto, Alejandra Serrano, Jianjun Gao, Humam Kadara, Catherine Fridman, Wolf Herve Fridman, Linghua Wang. Pan-cancer single-cell and subcellular characterization of lymphoid aggregates and tertiary lymphoid structures using high-plex spatial molecular imaging [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1152.

ICELLNET v2: a versatile method for cell-cell communication analysis from human transcriptomic data.

Several methods have been developed in the past years to infer cell-cell communication networks from transcriptomic data based on ligand and receptor expression. Among them, ICELLNET is one of the few approaches to consider the multiple subunits of ligands and receptors complexes to infer and quantify cell communication. In here, we present a major update of ICELLNET. As compared to its original implementation, we (i) drastically expanded the ICELLNET ligand-receptor database from 380 to 1669 biologically curated interactions, (ii) integrated important families of communication molecules involved in immune crosstalk, cell adhesion, and Wnt pathway, (iii) optimized ICELLNET framework for single-cell RNA sequencing data analyses, (iv) provided new visualizations of cell-cell communication results to facilitate prioritization and biological interpretation. This update will broaden the use of ICELLNET by the scientific community in different biological fields. ICELLNET package is implemented in R. Source code, documentation and tutorials are available on GitHub (https://github.com/soumelis-lab/ICELLNET).

Single-cell transcriptomics identifies senescence-associated secretory phenotype (SASP) features of testicular aging in human.

The male reproductive system experiences degradation with age, predominantly impacting the testes. Testicular aging can result in failure to produce physiological testosterone levels, normal sperm concentrations, or both. However, we cannot predict the onset of testicular aging in advance. Using single-cell RNA sequencing (scRNA-seq) from Gene Expression Omnibus (GEO) database, we conducted cell-cell communication network of human testis between older and young group, indicating Leydig cells' potential role in spermatogenesis microenvironment of aging testis. And we depicted the senescence-Associated Secretory Phenotype (SASP) features of aging testis by identifying differentially expressed senescence-associated secretory phenotype (SASP)-related genes between two group. Notably, IGFBP7 mainly expressed in Leydig cells of those differentially expressed SASP-related genes in aging testis. Furthermore, IGFBP7 protein located in the interstitial compartment of older mice confirmed by immunofluorescence and highly expressed in both human seminal plasma and mouse testis in the older group confirmed through Western blot. Together, our findings suggest that IGFBP7 may be a new biomarker of testicular aging.

Cell-cell communication network-based interpretable machine learning predicts cancer patient response to immune checkpoint inhibitors.

Immune checkpoint inhibitors (ICIs) have revolutionized cancer treatment. However, only some patients respond to ICIs, and current biomarkers for ICI efficacy have limited performance. Here, we devised an interpretable machine learning (ML) model trained using patient-specific cell-cell communication networks (CCNs) decoded from the patient's bulk tumor transcriptome. The model could (i) predict ICI efficacy for patients across four cancer types (median AUROC: 0.79) and (ii) identify key communication pathways with crucial players responsible for patient response or resistance to ICIs by analyzing more than 700 ICI-treated patient samples from 11 cohorts. The model prioritized chemotaxis communication of immune-related cells and growth factor communication of structural cells as the key biological processes underlying response and resistance to ICIs, respectively. We confirmed the key communication pathways and players at the single-cell level in patients with melanoma. Our network-based ML approach can be used to expand ICIs' clinical benefits in cancer patients.

Characterization of CD4+ and CD8+ T cells responses in the mixed lymphocyte reaction by flow cytometry and single cell RNA sequencing.

The Mixed Lymphocyte Reaction (MLR) consists in the allogeneic co-culture of monocytes derived dendritic cells (MoDCs) with T cells from another donor. This in vitro assay is largely used for the assessment of immunotherapy compounds. Nevertheless, the phenotypic changes associated with lymphocyte responsiveness under MLR have never been thoroughly evaluated. Here, we used multiplex cytokine and chemokine assays, multiparametric flow cytometry and single cell RNA sequencing to deeply characterize T cells activation and function in the context of CD4+- and CD8+-specific MLR kinetics. We showed that CD4+ and CD8+ T cells in MLR share common classical markers of response such as polyfunctionality, increased proliferation and CD25 expression but differ in their kinetics and amplitude of activation as well as their patterns of cytokines secretion and immune checkpoints expression. The analysis of immunoreactive Ki-67+CD25+ T cells identified PBK, LRR1 and MYO1G as new potential markers of MLR response. Using cell-cell communication network inference and pathway analysis on single cell RNA sequencing data, we also highlighted key components of the immunological synapse occurring between T cells and the stimulatory MoDCs together with downstream signaling pathways involved in CD4+ and CD8+ T cells activation. These results provide a deep understanding of the kinetics of the MLR assay for CD4+ or CD8+ T cells and may allow to better characterize compounds impacting MLR and eventually identify new strategies for immunotherapy in cancer.

Bioengineered hydrogels enhance ex vivo preservation of patient-derived tumor explants for drug evaluation

Ex vivo patient-derived tumor slices (PDTS) are currently limited by short-term viability in culture. Here, we show how bioengineered hydrogels enable the identification of key matrix parameters that significantly enhance PDTS viability compared to conventional culture systems. As demonstrated using single-cell RNA sequencing and high-dimensional flow cytometry, hydrogel-embedded PDTS tightly preserved cancer, cancer-associated fibroblast, and various immune cell populations and subpopulations in the corresponding original tumor. Cell-cell communication networks within the tumor microenvironment, including immune checkpoint ligand-receptor interactions, were also maintained. Remarkably, our results from a co-clinical trial suggest hydrogel-embedded PDTS may predict sensitivity to immune checkpoint inhibitors (ICIs) in head and neck cancer patients. Further, we show how these longer term-cultured tumor explants uniquely enable the sampling and detection of temporal evolution in molecular readouts when treated with ICIs. By preserving the compositional heterogeneity and complexity of patient tumors, hydrogel-embedded PDTS provide a valuable tool to facilitate experiments targeting the tumor microenvironment.

Elucidating T cell dynamics and molecular mechanisms in syngeneic and allogeneic islet transplantation through single-cell RNA sequencing.

Islet transplantation is a promising therapy for diabetes treatment. However, the molecular underpinnings governing the immune response, particularly T-cell dynamics in syngeneic and allogeneic transplant settings, remain poorly understood. Understanding these T cell dynamics is crucial for enhancing graft acceptance and managing diabetes treatment more effectively. This study aimed to elucidate the molecular mechanisms, gene expression differences, biological pathway alterations, and intercellular communication patterns among T-cell subpopulations after syngeneic and allogeneic islet transplantation. Using single-cell RNA sequencing, we analyzed cellular heterogeneity and gene expression profiles using the Seurat package for quality control and dimensionality reduction through t-SNE. Differentially expressed genes (DEGs) were analyzed among different T cell subtypes. GSEA was conducted utilizing the HALLMARK gene sets from MSigDB, while CellChat was used to infer and visualize cell-cell communication networks. Our findings revealed genetic variations within T-cell subpopulations between syngeneic and allogeneic islet transplants. We identified significant DEGs across these conditions, highlighting molecular discrepancies that may underpin rejection or other immune responses. GSEA indicated activation of the interferon-alpha response in memory T cells and suppression in CD4+ helper and γδ T cells, whereas TNFα signaling via NFκB was particularly active in regulatory T cells, γδ T cells, proliferating T cells, and activated CD8+ T cells. CellChat analysis revealed complex communication patterns within T-cell subsets, notably between proliferating T cells and activated CD8+ T cells. In conclusion, our study provides a comprehensive molecular landscape of T-cell diversity in islet transplantation. The insights into specific gene upregulation in xenotransplants suggest potential targets for improving graft tolerance. The differential pathway activation across T-cell subsets underscores their distinct roles in immune responses posttransplantation.

An advanced comprehensive muti-cell-type-specific model for predicting anti-PD-1 therapeutic effect in melanoma.

Rationale: Immune checkpoint inhibitors targeting the programmed cell death (PD)-1/PD-L1 pathway have promise in patients with advanced melanoma. However, drug resistance usually results in limited patient benefits. Recent single-cell RNA sequencing studies have elucidated that MM patients display distinctive transcriptional features of tumor cells, immune cells and interstitial cells, including loss of antigen presentation function of tumor cells, exhaustion of CD8+T and extracellular matrix secreted by fibroblasts to prevents immune infiltration, which leads to a poor response to immune checkpoint inhibitors (ICIs). However, cell subgroups beneficial to anti-tumor immunity and the model developed by them remain to be further identified. Methods: In this clinical study of neoadjuvant therapy with anti-PD-1 in advanced melanoma, tumor tissues were collected before and after treatment for single-nucleus sequencing, and the results were verified using multicolor immunofluorescence staining and public datasets. Results: This study describes four cell subgroups which are closely associated with the effectiveness of anti-PD-1 treatment. It also describes a cell-cell communication network, in which the interaction of the four cell subgroups contributes to anti-tumor immunity. Furthermore, we discuss a newly developed predictive model based on these four subgroups that holds significant potential for assessing the efficacy of anti-PD-1 treatment. Conclusions: These findings elucidate the primary mechanism of anti-PD-1 resistance and offer guidance for clinical drug administration for melanoma.

Dictionary of immune responses to cytokines at single-cell resolution

Cytokines mediate cell–cell communication in the immune system and represent important therapeutic targets1–3. A myriad of studies have highlighted their central role in immune function4–13, yet we lack a global view of the cellular responses of each immune cell type to each cytokine. To address this gap, we created the Immune Dictionary, a compendium of single-cell transcriptomic profiles of more than 17 immune cell types in response to each of 86 cytokines (>1,400 cytokine–cell type combinations) in mouse lymph nodes in vivo. A cytokine-centric view of the dictionary revealed that most cytokines induce highly cell-type-specific responses. For example, the inflammatory cytokine interleukin-1β induces distinct gene programmes in almost every cell type. A cell-type-centric view of the dictionary identified more than 66 cytokine-driven cellular polarization states across immune cell types, including previously uncharacterized states such as an interleukin-18-induced polyfunctional natural killer cell state. Based on this dictionary, we developed companion software, Immune Response Enrichment Analysis, for assessing cytokine activities and immune cell polarization from gene expression data, and applied it to reveal cytokine networks in tumours following immune checkpoint blockade therapy. Our dictionary generates new hypotheses for cytokine functions, illuminates pleiotropic effects of cytokines, expands our knowledge of activation states of each immune cell type, and provides a framework to deduce the roles of specific cytokines and cell–cell communication networks in any immune response.

Early Intervention and Favorable Biologic Effects of Personalized Neoantigen Vaccines on the Tumor Immune Microenvironment in Smoldering Waldenstrom Macroglobulinemia

Lymphoplasmacytic lymphoma (LPL) is an incurable low-grade B-cell lymphoma, characterized by the presence of clonal tumor cells infiltrating the bone marrow. Despite a cumulative risk of progression, there is no approved therapy for patients in the asymptomatic phase. We conducted a first-in-human clinical trial of a novel individualized therapeutic DNA idiotype neoantigen vaccine as an early intervention in nine patients with untreated asymptomatic LPL. The vaccine platform encoded the autologous LPL-derived Ig single chain variable fragment fused to human chemokine CCL20 (macrophage inflammatory protein-3, MIP-3a), which was designed to trigger T-cell immunity by targeting delivery of the expressed fusion protein to antigen presenting cells. The vaccine was well tolerated with no dose limiting toxicities. One patient achieved a minor response, and all remaining patients experienced stable disease, with median time to disease progression of 61+ months, which compares favorably to median time to progression of 3.9 years in published series. Direct interrogation of the tumor microenvironment by single-cell transcriptome analysis revealed explanations for lack of more robust objective clinical responses. Comparing paired pre- and post- vaccine bone marrow samples, we observed a striking dichotomous pattern of significantly reduced numbers of clonal tumor cells, identified by their unique BCR sequences, and downregulation of genes involved in signaling pathways critical for B-cell survival among the mature B cells subpopulation post-vaccine in the majority of patients, but no change in clonally related plasma cell-like clusters of any patient (Figure 1). Downregulation of HLA molecule expression suggested intrinsic resistance by tumor plasma cell subpopulations and cell-cell interaction analyses predicted paradoxical upregulation of IGF signaling post vaccine by plasma cell, but not mature B-cell subpopulations, suggesting a potential mechanism of acquired resistance. Vaccine therapy induced dynamic changes in bone marrow T-cells, including upregulation of signaling pathways involved in T-cell activation, expansion of T-cell clonotypes (Figure 2), increased T-cell clonal diversity, and functional tumor idiotype-specific cytokine production, with little change in co-inhibitory pathways or Treg. Vaccine therapy also globally altered cell-cell communication networks across various bone marrow cell types. CellChat analysis inferred significant downregulation of signaling pathways post-vaccine that likely directly promote growth of LPL cells, such as APRIL and IL-6 which are known to promote B- or plasma cell survival. Other pathways were reduced post-vaccine, such as RESISTIN, which has a known role in supporting proliferation of solid cancers. A role for RESISTIN in supporting LPL has not been previously inferred, but it has been reported to induce multidrug resistance in multiple myeloma. Bioinformatic analysis also identified a predominant role for myeloid cells in the tumor microenvironment as a source of vaccine-induced, downregulated pro-tumoral signaling to LPL cells. The pathways affected were primarily associated with monocytic, rather than granulocytic or dendritic cell subpopulations, particularly non-classical CD14 -CD16 + monocytes, and to a lesser extent classical CD14 +CD16 - monocytes . Recent reports indicate that an increased pro-inflammatory myeloid signature is an early step in the development of LPL and in monoclonal gammopathy of undetermined significance. Taken together, these results suggest that this prototype neoantigen vaccine favorably perturbed the tumor immune microenvironment, resulting in reduction of clonal tumor mature B-cell, but not plasma cell subpopulations. Future functional studies of the pathways affected are needed to confirm mechanisms of resistance elucidated and to design combination strategies to circumvent them. Such strategies could include adding IFN gamma or epigenetic drugs, designed to increase HLA molecule expression on plasma cell-like LPL subpopulations and combining neoantigen vaccines with agents that specifically target plasma cells or pathways known to promote their growth, such as IGF-1 receptor inhibitors. Finally, our data suggest that combinations of these vaccines with myeloid cell checkpoint blockade may be worthwhile.

Galectin-1 Fuels Monocyte-Driven Hyperinflammation and Represents a Novel Therapeutic Target in Myeloproliferative Neoplasms

JDP and STO co-corresponding authors Previous studies by our group and others have demonstrated that monocytes play a key role in driving hyperinflammation in myeloproliferative neoplasms (MPNs), and that aberrant inflammatory cytokine signaling contributes to disease progression and poor prognosis in MPNs. To interrogate relationships between cellular sources and targets of inflammatory cytokines, we performed single cell RNA-seq (scRNA-seq) in PBMCs and progenitor cells from MPN patients and healthy controls. Monocytes from MPN patients displayed significant enrichment of inflammation-relevant genes as expected. Cell-cell communication networks inferred from expression of ligands and receptors predicted monocytes as a pivotal mediator of cell interactions. Further analysis revealed galectin signaling as one of the most robust input/output pathways for monocytes and significant enrichment of LGALS1 (galectin-1) expression in MPN monocytes. Galectins are a class of proteins that bind specifically to β-galactoside carbohydrates, such as N-linked or O-linked glycosylated proteins. Abnormal expression of galectins in various cancers has been found to mediate broad biological functions, including cell proliferation, apoptosis, adhesion, and inflammation. However, the mechanisms by which galectins contribute to MPN pathogenesis remain incompletely understood. To corroborate our sRNA-seq findings, we performed flow cytometry analysis which demonstrated elevated expression of galectin-1 in MPN patient CD14+ monocytes. Increased plasma levels of galectin-1 were also identified in MPN patients compared to healthy individuals. Higher galectin-1 protein levels in both Kit+ progenitor cells and CD11b+ myeloid cells from Jak2 V617F knock-in mice were also observed when compared to wild type mice. Additionally, galectin-1 expression was induced by MPL W515L and JAK2 V617F and inhibited by ruxolitinib in Ba/F3 cells. Thus, our results show enrichment of galectin-1 in both MPN patients and mouse models, with evidence of direct contribution from specific MPN driver mutations. Our scRNA-seq dataset analysis showed significant correlations of galectin-1 expression with inflammatory and oxidative phosphorylation (OXPHOS) genes in MPN patient monocytes. To investigate the effects of galectin-1 in monocyte activation, we incubated CD14+ monocytes from MPN patients with recombinant galectin-1 (rGal-1) and observed markedly stimulated secretion of inflammatory cytokines, such as IL-1α, IL-1β, IL-6, IL-8 and TNF-α. Furthermore, OTX008, a galectin-1 inhibitor, inhibited the transcription and secretion of inflammatory cytokines in MPN patient CD14+ monocytes and monocytic cell lines. Notably, both genetic and pharmacologic inhibition of galectin-1 reduced cellular ATP level and oxygen consumption rate, suggesting a metabolic reprogramming in MPNs. Although galectin-1 mediates broad biological functions, the molecular mechanisms by which galectin-1 regulates signaling pathways and transcription in MPN cells had not been delineated. Our scRNA-seq heterogeneity analysis suggested that genes in the PI3K-AKT-mTOR pathway were differentially expressed between galectin-1 high- and low-expressing monocytes. We therefore incubated both monocytic cell lines and monocytes from MPN patients with rGal-1 and observed stimulation of the PI3K-AKT-mTOR pathway, represented by increased levels of phosphorylated mTOR, AKT and S6. Consistently, genetic and pharmacologic inhibition of galectin-1 inhibited activation of the PI3K-AKT-mTOR pathway. In summary, our results demonstrate upregulation of galectin-1 in MPN monocytes, potentially due to activation of JAK2-driven signaling. We further demonstrate that galectin-1 fuels inflammation potentially via metabolic reprogramming of monocytes and activation of PI3K-AKT-mTOR signaling. Our data also suggest galectin-1 as a putative therapeutic target in MPNs.

Single-nucleus RNA sequencing reveals heterogenous microenvironments and specific drug response between cervical squamous cell carcinoma and adenocarcinoma

Cervical squamous cell carcinoma (CSCC) and adenocarcinoma (CAde) are two major pathological types of cervical cancer (CC), but their high-resolution heterogeneity of tumor and immune microenvironment remains elusive. Here, we performed single-nucleus RNA sequencing (snRNA-seq) from five CSCC and three CAde samples, and systematically outlined their specific transcriptome atlas. We found CD8+ T cells in CSCC were more cytotoxic but lower exhausted compared to those in CAde, and phagocytic MRC1+ macrophages were specifically enriched in CSCC. Interestingly, we discovered that pro-tumoral cancer-associated myofibroblasts (myoCAFs) and cancer-associated vascular-fibroblasts (vCAFs) were more abundant in CSCC, and further verified their pro-metastatic roles invitro. Furthermore, we also identified some specific chemotherapy drugs for CSCC (Dasatinib and Doramapimod) and CAde (Pyrimethamine and Lapatinib) by revealing their heterogeneity in transcriptomic profiles of malignant epithelial cells, and further verified their specific sensitivity in cell lines and constructed CC-derived organoids. Cell-cell communication networks revealed that the pathways of NRG1-ERBB2, and FN1-ITAG3 were specific for CAde and CSCC, respectively, which may partly explain the specificities of identified chemotherapy drugs. Our study described the immune heterogeneity and specific cellular interactions between CSCC and CAde, which could provide insights for uncovering pathogenesis and designing personalized treatment. National Key R&D Program of China (2021YFC2701201), National Natural Science Foundation of China (82072895, 82141106, 82103134, 81903114).

Single-cell RNA-seq public data reveal the gene regulatory network landscape of respiratory epithelial and peripheral immune cells in COVID-19 patients.

Infection with SARS-CoV-2 leads to coronavirus disease 2019 (COVID-19), which can result in acute respiratory distress syndrome and multiple organ failure. However, its comprehensive influence on pathological immune responses in the respiratory epithelium and peripheral immune cells is not yet fully understood. In this study, we analyzed multiple public scRNA-seq datasets of nasopharyngeal swabs and peripheral blood to investigate the gene regulatory networks (GRNs) of healthy individuals and COVID-19 patients with mild/moderate and severe disease, respectively. Cell-cell communication networks among cell types were also inferred. Finally, validations were conducted using bulk RNA-seq and proteome data. Similar and dissimilar regulons were identified within or between epithelial and immune cells during COVID-19 severity progression. The relative transcription factors (TFs) and their targets were used to construct GRNs among different infection sites and conditions. Between respiratory epithelial and peripheral immune cells, different TFs tended to be used to regulate the activity of a cell between healthy individuals and COVID-19 patients, although they had some TFs in common. For example, XBP1, FOS, STAT1, and STAT2 were activated in both the epithelial and immune cells of virus-infected individuals. In contrast, severe COVID-19 cases exhibited activation of CEBPD in peripheral immune cells, while CEBPB was exclusively activated in respiratory epithelial cells. Moreover, in patients with severe COVID-19, although some inflammatory genes, such as S100A8/A9, were found to be upregulated in both respiratory epithelial and peripheral immune cells, their relative regulators can differ in terms of cell types. The cell-cell communication analysis suggested that epidermal growth factor receptor signaling among epithelia contributes to mild/moderate disease, and chemokine signaling among immune cells contributes to severe disease. This study identified cell type- and condition-specific regulons in a wide range of cell types from the initial infection site to the peripheral blood, and clarified the diverse mechanisms of maladaptive responses to SARS-CoV-2 infection.

Single-Cell RNA Sequencing Analysis of the Early Postnatal Mouse Lens Epithelium.

The lens epithelium maintains the overall health of the organ. We used single-cell RNA sequencing (scRNA-seq) technology to assess transcriptional heterogeneity between cells in the postnatal day 2 (P2) epithelium and identify distinct epithelial cell subtypes. Analysis of these data was used to better understand lens growth, differentiation, and homeostasis on P2. scRNA-seq on P2 mouse lenses was performed using the 10x Genomics Chromium Single Cell 3' Kit (v3.1) and short-read Illumina sequencing. Sequence alignment and preprocessing of data were conducted using 10x Genomics Cell Ranger software. Seurat was employed for preprocessing, quality control, dimensionality reduction, and cell clustering, and Monocle was utilized for trajectory analysis to understand the developmental progression of the lens cells. CellChat and GO analyses were used to explore cell-cell communication networks and signaling interactions. Lens epithelial cells (LECs) were divided into seven subclusters, classified by specific gene markers. The expression of crystallin, cell-cycle, and metabolic genes was not uniform, indicating distinct functional roles of LECs. Trajectory analysis predicted a bifurcation of differentiating and cycling cells from an Igfbp5+ progenitor pool. We also identified heterogeneity in signaling molecules and pathways, suggesting that cycling and progenitor subclusters have prominent roles in coordinating crosstalk. scRNA-seq corroborated many known markers of epithelial differentiation and proliferation while providing further insight into the pathways and genes directing these processes. Interestingly, we demonstrated that the developing epithelium can be divided into distinct subpopulations. These clusters reflect the transcriptionally diverse roles of the epithelium in proliferation, signaling, and maintenance.

CellCommuNet: an atlas of cell-cell communication networks from single-cell RNA sequencing of human and mouse tissues in normal and disease states.

Cell-cell communication, as a basic feature of multicellular organisms, is crucial for maintaining the biological functions and microenvironmental homeostasis of cells, organs, and whole organisms. Alterations in cell-cell communication contribute to many diseases, including cancers. Single-cell RNA sequencing (scRNA-seq) provides a powerful method for studying cell-cell communication by enabling the analysis of ligand-receptor interactions. Here, we introduce CellCommuNet (http://www.inbirg.com/cellcommunet/), a comprehensive data resource for exploring cell-cell communication networks in scRNA-seq data from human and mouse tissues in normal and disease states. CellCommuNet currently includes 376 single datasets from multiple sources, and 118 comparison datasets between disease and normal samples originating from the same study. CellCommuNet provides information on the strength of communication between cells and related signalling pathways and facilitates the exploration of differences in cell-cell communication between healthy and disease states. Users can also search for specific signalling pathways, ligand-receptor pairs, and cell types of interest. CellCommuNet provides interactive graphics illustrating cell-cell communication in different states, enabling differential analysis of communication strength between disease and control samples. This comprehensive database aims to be a valuable resource for biologists studying cell-cell communication networks.

Single-Cell Transcriptomics Reveals Cellular Heterogeneity and Complex Cell-Cell Communication Networks in the Mouse Cornea.

To generate a single-cell RNA-sequencing (scRNA-seq) map and construct cell-cell communication networks of mouse corneas. C57BL/6 mouse corneas were dissociated to single cells and subjected to scRNA-seq. Cell populations were clustered and annotated for bioinformatic analysis using the R package "Seurat." Differential expression patterns were validated and spatially mapped with whole-mount immunofluorescence staining. Global intercellular signaling networks were constructed using CellChat. Unbiased clustering of scRNA-seq transcriptomes of 14,732 cells from 40 corneas revealed 17 cell clusters of six major cell types: nine epithelial cell, three keratocyte, two corneal endothelial cell, and one each of immune cell, vascular endothelial cell, and fibroblast clusters. The nine epithelial cell subtypes included quiescent limbal stem cells, transit-amplifying cells, and differentiated cells from corneas and two minor conjunctival epithelial clusters. CellChat analysis provided an atlas of the complex intercellular signaling communications among all cell types. We constructed a complete single-cell transcriptomic map and the complex signaling cross-talk among all cell types of the cornea, which can be used as a foundation atlas for further research on the cornea. This study also deepens the understanding of the cellular heterogeneity and heterotypic cell-cell interaction within corneas.

Single-cell RNA sequencing of retina revealed novel transcriptional landscape in high myopia and underlying cell-type-specific mechanisms.

High myopia is a leading cause of blindness worldwide with increasing prevalence. Retina percepts visual information and triggers myopia development, but the underlying etiology is not fully understood because of cellular heterogeneity. In this study, single-cell RNA sequencing analysis was performed on retinas of mouse highly myopic and control eyes to dissect the involvement of each cell type during high myopia progression. For highly myopic photoreceptors, Hk2 inhibition underlying metabolic remodeling from aerobic glycolysis toward oxidative phosphorylation and excessive oxidative stress was identified. Importantly, a novel Apoe + rod subpopulation was specifically identified in highly myopic retina. In retinal neurons of highly myopic eyes, neurodegeneration was generally discovered, and the imbalanced ON/OFF signaling driven by cone-bipolar cells and the downregulated dopamine receptors in amacrine cells were among the most predominant findings, indicating the aberrant light processing in highly myopic eyes. Besides, microglia exhibited elevated expression of cytokines and TGF-β receptors, suggesting enhanced responses to inflammation and the growth-promoting states involved in high myopia progression. Furthermore, cell-cell communication network revealed attenuated neuronal interactions and increased glial/vascular interactions in highly myopic retinas. In conclusion, this study outlines the transcriptional landscape of highly myopic retina, providing novel insights into high myopia development and prevention.

Spatially distinct otic mesenchyme cells show molecular and functional heterogeneity patterns before hearing onset

The cochlea consists of diverse cellular populations working in harmony to convert mechanical stimuli into electrical signals for the perception of sound. Otic mesenchyme cells (OMCs), often considered a homogeneous cell type, are essential for normal cochlear development and hearing. Despite being the most numerous cell type in the developing cochlea, OMCs are poorly understood. OMCs are known to differentiate into spatially and functionally distinct cell types, including fibrocytes of the lateral wall and spiral limbus, modiolar osteoblasts, and specialized tympanic border cells of the basilar membrane. Here, we show that OMCs are transcriptionally and functionally heterogeneous and can be divided into four distinct populations that spatially correspond to OMC-derived cochlear structures. We also show that this heterogeneity and complexity of OMCs commences during early phases of cochlear development. Finally, we describe the cell-cell communication network of the developing cochlea, inferring a major role for OMC in outgoing signaling.

Analysis of matrisome expression patterns in murine and human dorsal root ganglia.

The extracellular matrix (ECM) is a dynamic structure of molecules that can be divided into six different categories and are collectively called the matrisome. The ECM plays pivotal roles in physiological processes in many tissues, including the nervous system. Intriguingly, alterations in ECM molecules/pathways are associated with painful human conditions and murine pain models. Nevertheless, mechanistic insight into the interplay of normal or defective ECM and pain is largely lacking. The goal of this study was to integrate bulk, single-cell, and spatial RNA sequencing (RNAseq) datasets to investigate the expression and cellular origin of matrisome genes in male and female murine and human dorsal root ganglia (DRG). Bulk RNAseq showed that about 65% of all matrisome genes were expressed in both murine and human DRG, with proportionally more core matrisome genes (glycoproteins, collagens, and proteoglycans) expressed compared to matrisome-associated genes (ECM-affiliated genes, ECM regulators, and secreted factors). Single cell RNAseq on male murine DRG revealed the cellular origin of matrisome expression. Core matrisome genes, especially collagens, were expressed by fibroblasts whereas matrisome-associated genes were primarily expressed by neurons. Cell-cell communication network analysis with CellChat software predicted an important role for collagen signaling pathways in connecting vascular cell types and nociceptors in murine tissue, which we confirmed by analysis of spatial transcriptomic data from human DRG. RNAscope in situ hybridization and immunohistochemistry demonstrated expression of collagens in fibroblasts surrounding nociceptors in male and female human DRG. Finally, comparing human neuropathic pain samples with non-pain samples also showed differential expression of matrisome genes produced by both fibroblasts and by nociceptors. This study supports the idea that the DRG matrisome may contribute to neuronal signaling in both mouse and human, and that dysregulation of matrisome genes is associated with neuropathic pain.