SEARCH · PubMed Health

Results for “Cell Communication”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

Differential cell signaling testing for cell-cell communication inference from single-cell data by dominoSignal.

MOTIVATION: Algorithms for ligand-receptor network inference have emerged as commonly used tools to estimate cell-cell communication from reference single-cell data. Many studies employ these algorithms to compare signaling between conditions and lack methods to statistically identify signals that are significantly different. We previously developed the cell communication inference algorithm Domino, which considers ligand and receptor gene expression in association with downstream transcription factor activity scoring. We developed the dominoSignal software to innovate upon Domino and extend its functionality to test statistically differential cellular signaling. RESULTS: This new functionality includes the compilation of active signals as linkages from multiple subjects in a single-cell data set and testing condition-dependent signaling linkage. The software is applicable for analysis of single-cell data sets with multiple subjects as biological replicates as well as with bootstrapped replicates from data sets with few or pooled subjects. We use simulation studies to benchmark the number of subjects in compared groups and cells within an annotated cell type sufficient to accurately identify differential linkages. We demonstrate the application of the Differential Cell Signaling Test (DCST) in the dominoSignal software to investigate consequences of cancer cell phenotypes and immunotherapy on cell-cell communication in tumor microenvironments. These applications in cancer studies demonstrate the ability of differential cell signaling analysis to infer changes to cell communication networks from therapeutic or experimental perturbations, which is broadly applicable across biological systems. AVAILABILITY: dominoSignal is available through Bioconductor at https://www.bioconductor.org/packages/release/bioc/html/dominoSignal.html.

Cell Communication↗

PubMed2025

scACCorDiON: a clustering approach for explainable patient level cell-cell communication graph analysis.

MOTIVATION: Combining single-cell sequencing with ligand-receptor (LR) analysis paves the way for the characterization of cell communication events in complex tissues. In particular, directed weighted graphs naturally represent cell-cell communication events. However, current computational methods cannot yet analyze sample-specific cell-cell communication events, as measured in single-cell data produced in large patient cohorts. Cohort-based cell-cell communication analysis presents many challenges, such as the nonlinear nature of cell-cell communication and the high variability given by the patient-specific single-cell RNAseq datasets. RESULTS: Here, we present scACCorDiON (single-cell Analysis of Cell-Cell Communication in Disease clusters using Optimal transport in Directed Networks), an optimal transport algorithm exploring node distances on the Markov Chain as the ground metric between directed weighted graphs. Benchmarking indicates that scACCorDiON performs a better clustering of samples according to their disease status than competing methods that use undirected graphs. We provide a case study of pancreas adenocarcinoma, where scACCorDion detects a sub-cluster of disease samples associated with changes in the tumor microenvironment. Our study case corroborates that clusters provide a robust and explainable representation of cell-cell communication events and that the expression of detected LR pairs is predictive of pancreatic cancer survival. AVAILABILITY AND IMPLEMENTATION: The code of scACCorDiON is available at https://scaccordion.readthedocs.io/en/latest/. and https://doi.org/10.5281/zenodo.15267648. The survival analysis package can be found at https://github.com/CostaLab/scACCorDiON.su.

Humans↗

PubMed2026

Integrative Multi-Omics Analysis Identifies Thrombosis-Associated Molecular Features Linked to Germline Susceptibility and Immune Cell Communication in Gastric Cancer.

Emerging evidence indicates that coagulation-related molecular programs are associated with thrombosis, tumor progression, and molecular dysregulation in gastric cancer (GC). However, thrombosis-associated molecular features in GC and their potential links to inherited susceptibility remain insufficiently understood. Integrated analyses of transcriptomic data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were performed to identify thrombosis-associated genes and establish a machine learning-based prognostic signature. Genome-wide association study (GWAS), expression quantitative trait loci (eQTL), transcriptome-wide association study (TWAS), and Mendelian randomization (MR) analyses were conducted to investigate susceptibility-associated transcriptional programs in GC. Functional assays were used to evaluate candidate genes associated with malignant phenotypes. Single-cell RNA sequencing (scRNA-seq) and cell-cell communication analyses were further performed to characterize cell-type-specific expression patterns and potential intercellular interactions. A total of 22 differentially expressed thrombosis-associated genes were identified, and a prognostic signature comprising 14 genes was established. The signature stratified patients into high- and low-risk groups and showed prognostic performance in both the training and validation cohorts. Integrative GWAS, eQTL, and TWAS analyses identified susceptibility-associated transcriptional programs that were positively correlated with the thrombosis-associated risk score. Silencing ACTN2 and CRYAB significantly reduced GC cell migration and invasion. scRNA-seq analysis revealed relatively high CRYAB expression in neutrophils, and CellChat analysis suggested potential neutrophil-B cell interactions involving COLLAGEN-related signaling. This integrative multi-omics study identified a thrombosis-associated molecular signature linked to prognosis and germline susceptibility-associated transcriptional programs in GC. ACTN2 and CRYAB may represent candidate genes associated with GC cell migration and invasion, while single-cell analysis suggested potential immune-related communication features.

Humans↗

PubMed2025

Alcohol-induced KDM5B activation in hepatocytes drives pathogenic cell-cell communication, leading to loss of liver function.

BACKGROUND: Alcohol-associated liver disease (ALD) is a major cause of alcohol-associated mortality. Previously, we identified KDM5B as a sex-specific mediator of ALD development; however, the mechanism behind KDM5B-induced pathological changes is not established. METHODS: Kdm5b flox/flox female mice were fed a western diet and 20% alcohol in the drinking water for 8-16 weeks (WDA). To induce KO, mice received 2×1011 genome copies of AAV8-CMV-Cre, AAV8-TBG-Cre, or AAV8-control. To test the role of myeloid C/EBPβ, Cebpbfl/fl, or Cebpbfl/fl Lyz2-Cre mice were fed WDA for 16 weeks. RESULTS: We found that Kdm5b KO prevented alcohol-induced liver fibrosis and liver inflammation in female mice. These changes were in part mediated by hepatocyte-to-non-parenchymal cell communication changes. KDM5B in hepatocytes promoted pro-inflammatory and pro-fibrotic changes in liver macrophages, endothelial cells, and stellate cells. Moreover, KDM5B promoted alcohol-induced early increase in EpCAM-positive liver progenitors and loss of liver function at later time points of alcohol feeding. We found that loss of liver function was dependent on a hepatocyte-to-macrophage communication feedback loop. KDM5B in hepatocytes inhibited macrophage C/EBPβ expression, which in turn resulted in loss of the mature KCs phenotype and prevented the ability of KCs to support hepatocyte differentiation, ultimately leading to loss of liver synthetic function. CONCLUSIONS: KDM5B activation in hepatocytes drives pathogenic cell-cell communication, leading to alcohol-induced loss of liver function in ALD.

Animals↗

PubMed2024

MONOCYTES AND B CELLS MEDIATE ALTERATIONS IN THE GENETIC ASSOCIATION BETWEEN PLATELETS AND SEPSIS VIA CLEC SIGNALING PATHWAY.

Background: Sepsis is a life-threatening condition characterized by multiple organ dysfunction. Blood cells abnormalities play a significant role in the onset and progression of sepsis; however, the potential causal relationship between platelets and sepsis remains unclear, as does whether immune cells mediate the interaction between platelets and sepsis. This study aims to explore the potential causal relationship between platelets and sepsis and analyze the mediating effect of immune cells. In addition, cell-to-cell communication was analyzed to explore the interaction between blood cells and immune cells. Material and methods: In this study, genome-wide association study data were utilized to examine the association between blood cells and sepsis. Two-sample Mendelian randomization (MR) and reverse MR were performed to investigate the potential causal relationship between blood cells and sepsis, with a specific focus on the relationship between platelets and sepsis. Subsequently, two-step MR was employed to identify the immune cells that mediate the interaction between platelets and sepsis and to assess their potential mediating effects. Cellchat software was used to analyze cell-to-cell communication. Results: The results of two-sample MR indicated that platelets were negatively correlated with sepsis (OR = 0.976, 95% CI 0.959-0.993, P = 0.005), suggesting that platelets have a protective effect against sepsis. Additionally, reverse MR demonstrated that sepsis had no significant effect on platelets (OR = 0.909, 95% CI 0.156-5.296, P = 0.916). The mediating effect analysis revealed that monocytes and B cells were important mediators in the relationship between platelets and sepsis. Notably, the correlation between platelets and sepsis shifted from negative to positive with the involvement of monocytes and B cells. The number and strength of cell-cell interactions were decreased in sepsis. Monocytes and B cells primarily regulate platelets through the CLEC signaling pathway, contributing to the pathogenesis of sepsis. Conclusion: This study confirmed the protective role of platelets in sepsis. Monocytes and B cells mediate changes in the genetic association between platelets and sepsis. Monocytes and B cells primarily interact with platelets via the CLEC pathway, thereby modulating the genetic association between platelets and sepsis. These findings indicate that thrombocytopenia, especially when accompanied by elevated monocytes and B cells, may serve as a potential marker for sepsis.

Humans↗

PubMed2026

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma↗

PubMed2025

Learning directed acyclic graphs for ligands and receptors based on spatially resolved transcriptomic data of ovarian cancer.

To unravel the mechanism of immune activation and suppression within tumors, a critical step is to identify transcriptional signals governing cell-cell communication between tumor and immune/stromal cells in the tumor microenvironment. Central to this communication are interactions between secreted ligands and cell-surface receptors, creating a highly connected signaling network among cells. Recent advancements in in situ-omics profiling, particularly spatial transcriptomic (ST) technology, provide unique opportunities to directly characterize ligand-receptor signaling networks that power cell-cell communication. In this paper, we propose a novel statistical method, LRnetST, to characterize the ligand-receptor interaction networks between adjacent tumor and immune/stroma cells based on ST data. LRnetST utilizes a directed acyclic graph model with a novel approach to handle the zero-inflated distributions of ST data. It also leverages existing ligand-receptor regulation databases as prior information, and employs a bootstrap aggregation strategy to achieve robust network estimation. Application of LRnetST to ST data of high-grade serous ovarian tumor samples revealed both common and distinct ligand-receptor regulations across different tumors. Some of these interactions were validated through both a MERFISH dataset and a CosMx SMI dataset of independent ovarian tumor samples. These results cast light on biological processes relating to the communication between tumor and immune/stromal cells in ovarian tumors. An open-source R package of LRnetST is available on GitHub at https://github.com/jie108/LRnetST.

Humans↗

PubMed2026

Identification of Drug-resistant Cell Subpopulations in Colorectal Cancer Through Single-cell Analysis and Exploration of Potential Therapeutic Strategies.

INTRODUCTION: The therapeutic efficacy of Colorectal Cancer (CRC) is often compromised by resistance to the standard chemotherapy agent oxaliplatin. METHODS: This study obtained single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database. Differentially Expressed Genes (DEGs) between resistant and sensitive epithelial subpopulations were identified, followed by enrichment analysis. Pseudotemporal trajectory and cell-cell communication were analyzed using Monocle2 and CellChat, respectively. The candidate drug was predicted by Connectivity Map (cMAP) analysis. External validation included assessment of the EpC2 signature in an oxaliplatin-resistant cell line dataset (GSE76092), survival analysis using The Cancer Genome Atlas (TCGA) cohorts, and re-analysis of the GSE179784 dataset to assess the reproducibility of EpC2-like subpopulations and their DNA Damage Repair (DDR) scores. RESULTS: Cell subpopulations were divided into 10 clusters. Among them, epithelial cells comprised 5 subpopulations, with EPC2 identified as a potential oxaliplatin-resistant subset. DEGs were enriched in the TNF and IL-17 pathways. External validation confirmed the enrichment of EpC2 in resistant cell lines and its association with poor survival. Pseudotemporal trajectory revealed that epithelial cells underwent state transitions, forming two distinct branches. The resistant group exhibited enrichment in RNA splicing and NF-κB pathways. Cell-cell communication analysis revealed interactions involving MDK- NCL and PPIA-BSG. Dasatinib was predicted as a candidate drug. DISCUSSION: We identified an oxaliplatin-resistant subpopulation of Epithelial Cells (EpC2) in CRC, elucidated its multi-layered resistance mechanisms, and integrated multi- omics and cMAP database analyses to predict a potential intervention drug. CONCLUSION: This study provided potential therapeutic possibilities for oxaliplatin resistance, contributing to CRC treatment.

Humans↗

PubMed2026

MetaCCI: meta cell-cell interaction inference and its application to CCIs characteristics of MDS.

MOTIVATION: Cell-cell interactions (CCIs) are fundamental to multicellular organisms and play crucial roles in diverse biological processes and disease mechanisms. Understanding CCIs is vital for deciphering disease pathogenesis and developing therapeutic strategies. Although numerous computational methods have been developed to infer CCIs from complex biological data, most existing approaches rely primarily on single-gene expression levels and ligand-receptor databases, often failing to capture the nuanced network-wide changes characteristic of disease states. RESULT: We propose MetaCCI, a novel computational strategy that integrates meta-information into CCI inference by extending the traditional gene expression-based analysis to a gene regulatory network framework. MetaCCI meticulously combines established ligand-receptor pairs with quantitative insights into gene behavior within complex gene networks, enabling the precise extraction of relevant targets for CCI inference. Subsequently, CCI inference was performed using an eigen cell co-expression network, providing a more holistic view of cell-cell communication. Monte Carlo simulations demonstrated that MetaCCI consistently outperforms existing methods in CCI inference. We applied MetaCCI to characterize cell-cell communication in Myelodysplastic Syndromes (MDS). Our results identified distinct interaction patterns in MDS compared with normal cell populations, specifically highlighting the loss of CCIs between "Dendritic cells and Hematopoietic precursor cells" and between "Dendritic cells and Hematopoietic multipotent progenitor cells" as characteristic features of MDS. Furthermore, FABP5, CD63, and HMGB1 were identified as MDS-specific markers. These findings suggest that diminished CCIs involving dendritic cells, hematopoietic precursor cells, and multipotent progenitor cells are pivotal to MDS pathogenesis. AVAILABILITY AND IMPLEMENTATION: The MetaCCI software is freely available at https://github.com/HeewonGitHub/MetaCCI. An archived version of the software and example datasets used in this study is available at Zenodo: https://doi.org/10.5281/zenodo.20101527.

Myelodysplastic Syndromes↗

PubMed2026

Immune-Like Malignant Epithelial Programs Shape Tumor-Immune Interactions and Inform Prognostic Stratification in Lung Adenocarcinoma.

Lung adenocarcinoma (LUAD) is characterized by marked cellular heterogeneity, yet how malignant epithelial states contribute to immune regulation and clinical outcomes remains incompletely defined. We integrated single-cell RNA-sequencing data to map the cellular landscape of LUAD and identify malignant epithelial cells based on inferred copy-number alterations. Epithelial states were further examined through trajectory inference, transcription factor analysis, and cell-cell communication profiling. Single-cell-derived genes were subsequently integrated with TCGA and independent GEO cohorts to construct and validate a machine learning-based prognostic signature. Malignant epithelial cells displayed distinct functional programs, including an immune-like state associated with genomic instability, immune-related transcriptional activity, tumor-immune communication, and patient outcomes. The resulting immune-like malignant epithelial cell signature (IMEC-Sig) consistently stratified survival across multiple cohorts. Low IMEC-Sig scores were accompanied by greater immune infiltration, higher immune checkpoint expression, and increased immunophenoscore, whereas high scores were linked to a comparatively immunosuppressive phenotype. Pan-cancer analyses further identified KRT8 as a gene associated with unfavorable prognosis, and functional experiments showed that KRT8 silencing suppressed proliferation, migration, invasion, and colony formation in LUAD cells. Together, these findings connect malignant epithelial heterogeneity with the immune context and clinical outcomes, support IMEC-Sig as a biologically informed prognostic tool, and nominate KRT8 as a potential therapeutic target in LUAD.

Humans↗

PubMed2026

Senescent fibroblasts drive CD8+ T cell dysfunction in colorectal cancer via CD36-mediated lipid transfer and peroxidation.

BACKGROUND: Functional exhaustion of tumor-infiltrating CD8+ T cells represents a hallmark of colorectal cancer (CRC) immunosuppression, though its mechanistic drivers remain elusive. Given the established correlation between CRC progression and stromal senescence characterized by pathological lipid accumulation and impaired immunity, we investigated whether and how senescent fibroblasts actively regulate CD8+ T cell dysfunction. METHODS: Single-cell RNA sequencing (scRNA-seq) analysis was conducted to unveil the diverse fibroblast populations and the significant lipid metabolism changes between senescent fibroblasts and non-senescent fibroblasts in human CRC specimens and adjacent normal mucosa. Machine-learning identified senescent fibroblasts with a distinct gene signature. Cell-cell communication analysis was used to evaluate the interactions between senescent fibroblasts and CD8+ T cells in colorectal cancer. Co-culture experiments were conducted among senescent fibroblasts, CD8+ T cells and patient-derived organoids of CRC (CRC-PDOs), with the results evaluated with high-content imaging and propidium iodide/Hoechst 33,342 staining. Flow cytometry, ELISA and lipid pulse-chase with BODIPY FL C16 were performed to detect the alterations of CD8+ T cell cytotoxic function and metabolic status. AOM/DSS-induced CRC mouse model was used to conduct in vivo validation to evaluate whether senolytics could suppress CRC progression. Patients from the Cancer Genome Atlas colorectal cancer cohort were stratified into CD36-high and CD36-low groups by median expression, and drug sensitivity for GDSC2 compounds was predicted computationally using the oncoPredict R package. RESULTS: ScRNA-seq demonstrated the specific cell population presence and divergence of senescent fibroblasts between neoplastic and histologically normal adjacent cell clusters in CRC. Random Forest was employed for cell senescence classification. Feature importance analysis identified five genes as key contributors to the model’s decision process. Cell-cell communication analysis revealed enhanced interactions between senescent fibroblasts and CD8+ T cells in CRC. Co-culture of senescent fibroblasts significantly impaired the cytotoxic functions of CD8+ T cells on CRC-PDOs, which was reflected by the declined proportions of granzyme B (GZMB) + and interferon gamma (IFNγ) + CD8+ T cells and enhanced viability of CRC-PDOs. Mechanistically, the co-culture with senescent fibroblasts promoted the lipid shuttling into CD8+ T cells to induce lipid peroxidation and downstream impairment of cytotoxicity. Furthermore, the inhibition of CD36, the specific scavenger receptor for lipid uptake of CD8+ T cells, effectively suppressed lipid transfer and peroxidation thereby preserving the effector functions of CD8+ T cells and ultimately promoting tumor apoptosis. Complementarily, in vivo senolytic treatment significantly suppressed CRC progression in AOM-DSS CRC mouse models. Top 12 therapeutic agents were identified significantly enhanced predicted efficacy in CD36-high tumors. CONCLUSIONS: Our study identified a substantial population of senescent fibroblasts in human CRC through single cell transcriptomics, machine-learning and clinical biopsies. These senescent fibroblasts impair CD8+ T cell-mediated killing of CRC-PDOs via CD36-dependent lipid transfer, suggesting senolytic targeting of stromal cells as a promising immunotherapeutic strategy for CRC.

Colorectal Neoplasms↗

PubMed2025

A single-cell transcriptomic atlas of the pigtail macaque placenta in late gestation.

The placenta is a complex organ with multiple immune and non-immune cell types that promote fetal tolerance and facilitate the transfer of nutrients and oxygen. The nonhuman primate (NHP) is a key experimental model for studying human pregnancy complications, in part due to similarities in placental structure, which makes it essential to understand how single-cell populations compare across the human and NHP maternal-fetal interface. We constructed a single-cell RNA-Seq (scRNA-Seq) atlas of the placenta from the pigtail macaque ( Macaca nemestrina ) in the third trimester, comprising three different tissues at the maternal-fetal interface: the chorionic villi (placental disc), chorioamniotic membranes, and the maternal decidua. Each tissue was separately dissociated into single cells and processed through the 10X Genomics and Seurat pipeline, followed by aggregation, unsupervised clustering, and cluster annotation. Next, we determined the maternal-fetal origins of cell populations and analyzed single-cell RNA trajectory, Gene Ontology enrichment, and cell-cell communication. Single-cell populations in the pigtail macaque were strikingly similar in their identity and frequency to those found in the human placenta, including cells from trophoblast, stromal cell, immune, and macrophage lineages. An advantage of our approach was the deep sequencing of three tissues at the maternal-fetal interface, which yielded a rich diversity of common and rare single-cell populations. The third-trimester pigtail macaque single-cell atlas enables the identification of cellular subclusters analogous to those in humans and provides a powerful resource for understanding experimental perturbations on the NHP placenta.

Journal Article↗

PubMed2026

Cooperative care influences genome-wide levels of DNA methylation in nestling chestnut-crowned babblers.

Carers in cooperatively breeding vertebrates increase food acquisition for offspring; however, they also impact the developmental social environment. One means of linking early-life environments, such as nutrition and social structure, with later-life phenotypes is DNA methylation. Here, using whole-genome methylation sequencing, we measured how additional carers influence DNA methylation in nestlings of the cooperatively breeding chestnut-crowned babbler (Pomatostomus ruficeps). A comparison of nestlings raised by their parents (two carers) and those raised by parents plus additional helpers ('three plus' carers; mean = 4.3 ± 1.4 s.d.) revealed that additional care is associated with genome-wide differences in DNA methylation. Overall, 570 cytosine-phosphate-guanine sites from the regulatory regions of 487 genes were differentially methylated, with 85% being more methylated in nestlings reared by groups. Specifially, sites associated with genes that are integral for metabolism, growth, the regulation and promotion of sociality, the ability to cope with stressors, and cell communication were differentially methylated between the groups. Furthermore, gene ontology-term analysis revealed that differentially methylated sites were over-represented in multiple pathways, including those important for protein binding, metabolism and cell-to-cell and environment-to-cell communication. Our results suggest that the effects of being reared by groups as opposed to pairs in cooperative breeders can extend beyond those typically attributed to nutritional benefits and that these effects are molecularly mediated. This article is part of the theme issue 'Ecological epigenetics at the intersection of behaviour and life history variation in non-model animals'.

Animals↗

PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗

PubMed2026

Integrated single-cell and spatial transcriptomic analyses reveal malignant epithelial glycolytic heterogeneity and spatial niche remodeling during colorectal cancer progression.

Colorectal cancer (CRC) progression is shaped by metabolic reprogramming and complex interactions within the tumor microenvironment. However, the cellular heterogeneity, spatial organization, and clinical relevance of glycolytic activity in CRC remain incompletely understood. In this study, we integrated single-cell RNA sequencing, bulk transcriptomics, and spatial transcriptomics data to systematically characterize glycolytic heterogeneity in CRC. Glycolytic activity was quantified using five independent scoring methods, consistently showing that epithelial cells exhibited the highest glycolytic activity across the two single-cell cohorts. Stratification of CopyKAT-verified aneuploid malignant epithelial cells into high-glycolysis (HG) and low-glycolysis (LG) subgroups by glycolysis scores revealed that HG cells exhibited higher stemness scores and chromosomal copy number variations. Cell-cell communication analysis revealed that, compared with LG cells, HG cells exhibited increased interaction frequency and strength with immune and stromal populations, indicating enhanced malignant epithelial-microenvironment crosstalk. Spatial transcriptomics analyses further revealed that glycolytic activity varied across normal colorectal tissue, primary CRC, and colorectal liver metastases, accompanied by progressive remodeling of epithelial-associated spatial niches and MIF-mediated intercellular communication. Bulk transcriptomic analysis identified a glycolysis-related prognostic signature with robust predictive performance, which served as an independent prognostic factor for overall survival in CRC cohorts. Collectively, these findings indicate that glycolytic heterogeneity is a key feature of CRC malignant epithelial cells and is closely associated with tumor progression, microenvironmental remodeling, and clinical outcomes.

Humans↗

PubMed2026

Integrative multi-omics profiling of insomnia-related molecular features reveals microbiome, immune, and therapy-relevant heterogeneity in colorectal cancer.

Emerging evidence implicates insomnia as a potential risk factor in carcinogenesis, potentially involving systemic inflammation, circadian disruption, and microbiome alterations. However, the molecular associations linking insomnia-related features to colorectal cancer (CRC), particularly with respect to tumor biology, immune microenvironmental states, and therapy-relevant phenotypes, remain largely unexplored. Multi-omics integration of genomic, transcriptomic, and microbiome data from 3,026 CRC patients across seven independent cohorts, including a large, well-annotated Clinical Omics study of Colorectal Cancer in China (COCC) cohort, enabled insomnia-based molecular classification through unsupervised non-negative matrix factorization (NMF) clustering. The insomnia subtype (IS) was biologically characterized via pathway enrichment, immune deconvolution, microbial profiling, and single-cell transcriptomics. Furthermore, an insomnia score (ISscore) was developed and validated in multiple cohorts for risk stratification and assessment of treatment-response-related indicators in CRC. Unsupervised clustering revealed two distinct molecular subtypes (IS1/IS2), with IS2 demonstrating significantly poorer survival. IS2 exhibited marked activation of EMT/angiogenesis pathways versus cell cycle activation in IS1. The IS2 microenvironment showed increased immunosuppression-related infiltration and exhausted T cell signatures, together with intratumoral microbiome variation characterized by depletion of Ruminococcaceae UCG-002 and enrichment of Hungatella/Selenomonas. The ISscore system stratified survival risk and was associated with computational indicators of immunotherapy response. Single-cell analysis nominated PPIA-BSG as a potential cell-cell communication signal involving high-ISscore tumor cells, CXCL12+ endothelial cells, and CLEC9A+ dendritic cell subsets. This multi-omics characterization of insomnia-CRC interplay suggests that insomnia-related molecular features are associated with an immunologically distinct and microbiome-altered tumor ecosystem. The ISscore provides a reproducible framework for capturing insomnia-related molecular heterogeneity, supporting risk stratification and future evaluation of therapy-relevant phenotypes.IMPORTANCEChronic insomnia affects millions, but it is not typically considered a cancer risk factor. Our study, analyzing vast biological data from over 3,000 colorectal cancer patients, uncovers a potential link between a person's predisposition to insomnia and their risk of developing this disease. This suggests that the biological pathways related to sleep may play a role in cancer development. Understanding this connection opens up new avenues for identifying individuals at higher risk and developing novel prevention strategies for colorectal cancer.

colorectal cancer↗

PubMed2026

Methylation-Associated Differentiation Features Define Biological and Prognostic Heterogeneity in CMS4 Colorectal Cancer.

Consensus molecular subtype 4 (CMS4) colorectal cancer (CRC) is associated with an aggressive clinical course and poor survival, yet the biological basis of heterogeneity within this subtype remains incompletely understood. DNA methylation is an epigenetic mechanism involved in transcriptional regulation, cellular differentiation, and colorectal tumorigenesis. Here, we integrated single-cell RNA sequencing (scRNA-seq), bulk data, and promoter DNA methylation data to characterize CMS4-associated cancer cell states and methylation-related features. Using the scAB algorithm, we integrated scRNA-seq with bulk CMS4 data and identified CMS4-related cells distributed across multiple patients. Single-cell analyses of cell-cell communication and transcriptional regulation revealed a CMS4-related cancer cell population characterized by macrophage migration inhibitory factor (MIF)-centered intercellular communication, enhanced caudal type homeobox 1 (CDX1) and Kruppel-like factor 5 (KLF5) regulon activity, and gene modules enriched in differentiation-related pathways. CytoTRACE analysis further stratified CMS4 cancer cells into poorly and well-differentiated states, yielding 802 differentially expressed genes (DEGs). Linking these differentiation-associated DEGs with bulk expression and promoter methylation data identified 218 methylation-associated DEGs showing significant inverse methylation expression correlations, suggesting a link between differentiation-related heterogeneity and promoter methylation. Univariable Cox regression followed by LASSO regression further prioritized eight genes for construction of the methylation and differentiation-related prognostic model (MeDiff-PM). MeDiff-PM consistently stratified overall survival in the TCGA CMS4 cohort and two independent validation cohorts, with cutoff-independent continuous Cox analyses further supporting its prognostic association across cohorts. And MeDiff-PM remained prognostically significant after adjustment for available clinical variables. High MeDiff-PM risk scores were associated with activation of P53, WNT, and ubiquitin-mediated proteolysis pathways and with consistent predicted drug response differences for compounds across three CMS4 cohorts. While individual in silico knockout analysis suggested links between MeDiff-PM genes and metallothionein-related and immune-associated transcriptional responses. Collectively, these findings indicate that methylation-associated differentiation features represent a molecular dimension of intra-CMS4 heterogeneity and provide a biologically informed framework for prognostic stratification within CMS4 CRC.

Humans↗

PubMed2026

Whole-transcriptome-scale isoform-resolved spatial imaging of single cells in tissues.

Cell and tissue functions arise from complex interactions among numerous genes, and a systematic understanding of these functions requires isoform-resolved transcriptomic analysis of single cells with high spatial resolution. Here, we introduce an in situ RNA amplification method and its integration with multiplexed error-robust fluorescence in situ hybridization (MERFISH) to detect short RNA sequences and enable whole-transcriptome-scale, isoform-resolved spatial transcriptomics of individual cells in intact tissues. Using this approach, we imaged ∼33,000 distinct RNAs-including ∼23,000 genes and ∼10,000 isoforms-in the mouse brain. Our data enabled systematic analyses of region- and cell-type-specific gene programs and ligand-receptor-based cell-cell communications. These data further revealed rich spatial diversity and cell-type specificity in isoform usage across numerous genes, as well as brain structures particularly rich in isoform specificity. We anticipate broad application of this method for characterizing the molecular and cellular basis of tissue functions, unlocking previously inaccessible discoveries in cell and organismal biology.

Animals↗