利用单细胞测序和转录组测序建立结直肠癌免疫细胞的9基因预后模型

唐楠, 黄慧霞, 刘晓健

  1. 1.复旦大学附属肿瘤医院肿瘤内科,复旦大学上海医学院肿瘤学系,上海 200032
    2.复旦大学附属肿瘤医院病理科,复旦大学上海医学院肿瘤学系,上海 200032
    3.复旦大学附属肿瘤医院肿瘤研究所,上海 200032
    4.复旦大学病理学研究所,上海 200032
  • 收稿日期:2024-02-06 修回日期:2024-06-05 出版日期:2024-06-30 发布日期:2024-07-16
  • 通信作者: 刘晓健 E-mail:liuxiaojian3000@hotmail.com
  • 作者简介:唐 楠(ORCID:0000-0002-4771-5324),硕士研究生;
    黄慧霞(ORCID:0000-0002-6967-6872),硕士研究生。

摘要/Abstract

摘要:

背景与目的:结直肠癌(colorectal carcinoma,CRC)是一种常见的恶性肿瘤,其发病率和死亡率仅次于胃癌和食管癌,在消化系统恶性肿瘤中居第二位。越来越多的证据表明,免疫细胞在CRC的发生、发展中发挥着重要作用。本研究旨在构建一个与免疫细胞相关的基因的预后模型,以预测CRC患者的预后并进行精准管理。方法:从基因表达综合(Gene Expression Omnibus,GEO)数据库和癌症基因图谱(The Cancer Genome Atlas,TCGA)数据库下载结直肠癌的单细胞RNA测序(single-cell RNA sequencing,scRNA-seq)和普通转录组测序(RNA sequencing,RNA-seq)数据以及临床信息,提取免疫细胞亚型的差异基因,通过Cox回归和LASSO回归分析在TCGA数据中筛选出与预后相关的基因,同时使用GSE39582和GSE41258进行外部验证。基于预后模型进行化疗药物敏感性分析、免疫治疗效果分析、风险评分相关的通路分析以及临床相关性分析。采用实时荧光定量聚合酶链反应(real-time fluorescence quantitative polymerase chain reaction,RTFQ-PCR)和免疫组织化学检测验证模型基因在10例经术后病理学检查的新鲜冷冻结直肠癌样本和细胞系中的表达水平,样本采集均已获得复旦大学附属肿瘤医院伦理委员会的批准(伦理编号:050432-4-2108*)。结果:本研究根据scRNA-seq数据集(GSE161277)定义了16个细胞群,并通过R包celldex将这些群标记为不同的细胞类型。然后对免疫细胞亚型进行差异分析,共获得了374个差异表达基因。通过单因素Cox和LASSO回归分析,本研究在CRC中构建了一个9基因风险预后模型。该风险模型表现出对预后的可靠预测效果,对预测抗肿瘤药物敏感性、免疫治疗效果、潜在分子机制及临床特征等方面起重要作用,高风险分数的患者从免疫疗法中受益概率较低。结论:我们基于CRC肿瘤微环境中免疫细胞的异质性构建了一个包含9个基因的风险预后模型,并据此预测了CRC患者的生存和治疗效果。

关键词: 结直肠癌, 预后模型, 肿瘤免疫与免疫治疗, 药物筛选, scRNA-seq

Abstract:

Background and purpose: Colorectal carcinoma (CRC) is a common malignant tumor with incidence and mortality rates second only to gastric cancer and esophageal cancer among digestive system malignancies. Increasing evidence suggests that immune cells play a significant role in the occurrence and development of CRC. The aim of this study was to construct a prognostic model related to immune cell-associated genes to predict the prognosis of CRC patients and enable precise management. Methods: Single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (RNA-seq) data, along with clinical information for colorectal cancer, were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Differential genes of immune cell subtypes were extracted, and genes related to prognosis were screened in TCGA data using Cox regression and LASSO regression, with external validation using GSE39582 and GSE41258. The prognostic model was used to analyze chemotherapy drug sensitivity, assess immunotherapy efficacy, analyze pathways related to risk scores, and perform clinical correlation analysis. Finally, the expression levels of model genes were validated in 10 fresh frozen CRC tissues with post-operative pathological examination and cell lines using real-time fluorescence quantitative polymerase chain reaction (RTFQ-PCR) and immunohistochemistry. All samples were approved by the Fudan University Shanghai Cancer Center Ethics Committee (No. 050432-4-2108*). Results: We defined 16 cell clusters based on the scRNA-seq dataset (GSE161277) and labeled these clusters as different cell types using the R package celldex. Differential analysis of immune cell subtypes yielded 374 differentially expressed genes. Using univariate Cox and LASSO analyses, we constructed a 9-gene risk prognostic model in CRC. This risk model exhibited reliable prognostic prediction performance and played an important role in predicting anti-tumor drug sensitivity, immunotherapy efficacy, potential molecular mechanisms, and clinical characteristics. Patients with high-risk scores had a lower probability of benefiting from immunotherapy. Conclusion: We constructed a 9-gene risk prognostic model based on the heterogeneity of immune cells in the CRC tumor microenvironment, which predicted the survival and treatment outcomes of CRC patients.

Key words: Colorectal carcinoma, Prognostic model, Cancer immunity and immunotherapy, Drug screening, scRNA-seq

中图分类号: 

相关文章

[1] 翁俊勇, 叶紫岚, 张若昕, 刘琪, 李心翔. 探究不良病理学特征数量对Ⅰ~Ⅲ期结直肠癌复发风险分层的指导作用:对9 875例病例的回顾性队例研究[J]. 中国癌症杂志, 2024, 34(6): 527-536.
[2] 陈虹, 曹治云. 人源胰腺癌类器官模型的构建及应用新进展[J]. 中国癌症杂志, 2024, 34(6): 590-597.
[3] 张若昕, 叶紫岚, 翁俊勇, 李心翔. 高龄与Ⅱ期结直肠癌患者预后不良的相关性研究[J]. 中国癌症杂志, 2024, 34(5): 485-492.
[4] 陆悦, 卢仁泉, 张杰, 郑慧. 联合凝血功能指标在结直肠癌患者化疗后高凝状态监测中的应用价值[J]. 中国癌症杂志, 2024, 34(3): 278-285.
[5] 上海市抗癌协会大肠癌专业委员会. 结直肠癌早筛、早诊、早治上海方案(2023年版)[J]. 中国癌症杂志, 2024, 34(1): 13-66.
[6] 周聪, 何丽娜, 成小姣, 黄廷磊, 涂水平. RSPO3抑制体内结直肠癌移植瘤生长和增加NK细胞浸润的作用研究[J]. 中国癌症杂志, 2023, 33(7): 664-672.
[7] 冶俊玲, 郑小影, 郭新建, 陈瑞慧, 杨柳, 苟笑丹, 蒋汉梅. LncRNA SNHG5/miR-26a-5p/MTDH信号轴促进结直肠癌转移的机制研究[J]. 中国癌症杂志, 2023, 33(7): 673-685.
[8] 蔡家洛, 朱锐秋, 李森, 曹亦军, 黄坊. 炎症癌相关成纤维细胞介导结直肠癌细胞耐药的机制研究[J]. 中国癌症杂志, 2023, 33(12): 1065-1072.
[9] 陈玺, 曾晓颖, 陈佳艳, 刘菲, 唐曦. 错配修复蛋白联合血清肿瘤标志物与Ki-67增殖指数对结直肠癌预后的临床价值[J]. 中国癌症杂志, 2022, 32(3): 243-250.
[10] 王乾, 张静, 郁成礼, 包龙龙, 蔡旭, 蒋文华, 黄丹, 盛伟琪, 朱晓丽, 周晓燕, 柏乾明. 全自动IdyllaTM系统在中国结直肠癌微卫星不稳定性检测中的可行性研究[J]. 中国癌症杂志, 2022, 32(11): 1084-1090.
[11] 王若彤, 王欣, 沈波. 类器官在肿瘤转化医学中的应用和进展[J]. 中国癌症杂志, 2022, 32(11): 1105-1114.
[12] 贾真真, 何双, 李扬扬, 温菲菲, 许晓阳, 郭宁杰, 吴淑华. 结直肠癌中DPD与LC3、P62表达的相关性及其临床意义[J]. 中国癌症杂志, 2022, 32(1): 24-33.
[13] 马祎菲, 梁新军, 魏少忠 . 炎症与免疫指标在可切除性结直肠癌中的预后价值[J]. 中国癌症杂志, 2021, 31(9): 845-851.
[14] 张迎强, 林岩松 . 90 Y微球选择性内放射治疗在结直肠癌肝转移中的应用及研究进展[J]. 中国癌症杂志, 2021, 31(8): 762-768.
[15] 花 晴, 申雪芳, 许平波. KLK8通过调节EGF对结直肠癌细胞凋亡的影响及其机制[J]. 中国癌症杂志, 2021, 31(5): 383-389.