Integrating Bulk and Single-Cell RNA Sequencing to Identify Prognostic Genes Associated with Adaptive Immunity and the Tumor Microenvironment in Cholangiocarcinoma for Risk Model Construction
Wang Jie, Zhou Jixiang, Bao Yang, Wang Yanting, Liu Ling
{"title":"Integrating Bulk and Single-Cell RNA Sequencing to Identify Prognostic Genes Associated with Adaptive Immunity and the Tumor Microenvironment in Cholangiocarcinoma for Risk Model Construction","authors":"Wang Jie, Zhou Jixiang, Bao Yang, Wang Yanting, Liu Ling","doi":"10.1016/j.imbio.2026.153213","DOIUrl":null,"url":null,"abstract":"<div><h3>Background</h3><div>Cholangiocarcinoma (CCA) has a complex tumor microenvironment (TME) and poor immunotherapy response. Understanding how the adaptive immune response (AIR) affects the TME and prognosis is critical.</div></div><div><h3>Methods</h3><div>All data were obtained from publicly available resources. Differentially expressed genes (DEGs) were identified by comparing tumor vs. normal and stromal/immune score-based samples, and were then intersected with an AIR gene set to obtain candidate genes. Prognostic genes identified by univariate Cox regression were used to build a risk model. Analyses included functional enrichment, immune landscape, drug sensitivity, and single-cell RNA sequencing (scRNA-seq) to identify key cell types, cell–cell communication, and pseudotime trajectories.</div></div><div><h3>Results</h3><div>Intersection analysis yielded 68 candidate genes, from which eight prognostic genes (C1QA, CD55, ZAP70, CLEC4M, FCAMR, LILRA1, CD80, C17orf99) were identified. The risk model stratified patients into high-risk (HRG) and low-risk (LRG) groups with significantly different overall survival. LRG was enriched for immune activation pathways, whereas HRG exhibited an immune-excluded phenotype (lower CD8+ T cell infiltration, higher exclusion score). LRG showed higher sensitivity to gemcitabine and ABT-737. Macrophages were identified as key cells, with enhanced malignant cell–macrophage communication via the MIF–(CD74 + CXCR4) axis in CCA. Macrophages comprised seven heterogeneous subtypes and five distinct differentiation states. Firmicutes abundance was lower in the HRG. A nomogram integrating risk score and Firmicutes demonstrated favorable predictive performance.</div></div><div><h3>Conclusion</h3><div>An AIR-related gene signature stratifies CCA prognosis by reflecting TME composition rather than tumor-intrinsic drivers, highlighting immune cell–centric mechanisms and intratumoral microbiota. The nomogram may aid prognostic stratification in research settings.</div></div>","PeriodicalId":13270,"journal":{"name":"Immunobiology","volume":"231 4","pages":"Article 153213"},"PeriodicalIF":3.1000,"publicationDate":"2026-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Immunobiology","FirstCategoryId":"3","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0171298526000598","RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2026/7/6 0:00:00","PubModel":"Epub","JCR":"Q3","JCRName":"IMMUNOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
Background
Cholangiocarcinoma (CCA) has a complex tumor microenvironment (TME) and poor immunotherapy response. Understanding how the adaptive immune response (AIR) affects the TME and prognosis is critical.
Methods
All data were obtained from publicly available resources. Differentially expressed genes (DEGs) were identified by comparing tumor vs. normal and stromal/immune score-based samples, and were then intersected with an AIR gene set to obtain candidate genes. Prognostic genes identified by univariate Cox regression were used to build a risk model. Analyses included functional enrichment, immune landscape, drug sensitivity, and single-cell RNA sequencing (scRNA-seq) to identify key cell types, cell–cell communication, and pseudotime trajectories.
Results
Intersection analysis yielded 68 candidate genes, from which eight prognostic genes (C1QA, CD55, ZAP70, CLEC4M, FCAMR, LILRA1, CD80, C17orf99) were identified. The risk model stratified patients into high-risk (HRG) and low-risk (LRG) groups with significantly different overall survival. LRG was enriched for immune activation pathways, whereas HRG exhibited an immune-excluded phenotype (lower CD8+ T cell infiltration, higher exclusion score). LRG showed higher sensitivity to gemcitabine and ABT-737. Macrophages were identified as key cells, with enhanced malignant cell–macrophage communication via the MIF–(CD74 + CXCR4) axis in CCA. Macrophages comprised seven heterogeneous subtypes and five distinct differentiation states. Firmicutes abundance was lower in the HRG. A nomogram integrating risk score and Firmicutes demonstrated favorable predictive performance.
Conclusion
An AIR-related gene signature stratifies CCA prognosis by reflecting TME composition rather than tumor-intrinsic drivers, highlighting immune cell–centric mechanisms and intratumoral microbiota. The nomogram may aid prognostic stratification in research settings.
期刊介绍:
Immunobiology is a peer-reviewed journal that publishes highly innovative research approaches for a wide range of immunological subjects, including
• Innate Immunity,
• Adaptive Immunity,
• Complement Biology,
• Macrophage and Dendritic Cell Biology,
• Parasite Immunology,
• Tumour Immunology,
• Clinical Immunology,
• Immunogenetics,
• Immunotherapy and
• Immunopathology of infectious, allergic and autoimmune disease.