《Frontiers in Immunology》：Dissecting tumor heterogeneity in colorectal cancer: uncovering the role of BCL2L1+ cells through single-cell analysis
编辑推荐：
为阐明结直肠癌（CRC）肝转移的肿瘤异质性机制并寻找新的治疗靶点，研究人员通过单细胞RNA测序（scRNA-seq）分析了CRC原发灶及肝转移灶的细胞图谱。研究发现了一个高表达BCL2L1的C4肿瘤细胞亚型，该亚型在肝转移灶中富集，并表现出增强的增殖、代谢重编程、抗凋亡和免疫逃逸特性。研究进一步鉴定出转录因子CEBPG是关键调控因子，其沉默可抑制CRC细胞增殖与侵袭。基于此，研究者构建了BTRS预后风险评分模型，为CRC的临床风险分层和个体化治疗提供了新框架。
结直肠癌是全球范围内发病率和死亡率均位居前列的消化道恶性肿瘤。尽管诊疗技术在不断进步，但一旦发生转移，患者的预后会急剧恶化，其中肝脏是最常见的转移部位，是导致患者死亡的主要原因。当前治疗面临的巨大挑战在于肿瘤的异质性（即肿瘤内部的细胞类型和功能存在巨大差异）以及复杂的肿瘤微环境，这使得部分患者对现有疗法（包括化疗和新兴的免疫治疗）反应不佳，容易复发。为了深入“解剖”这种复杂性，揭示驱动转移和耐药的核心细胞群体，研究人员利用了一项强大的前沿技术——单细胞RNA测序，对结直肠癌的原发灶和肝转移灶进行了系统性分析，旨在发现关键的促转移细胞亚型及其调控机制，为改善患者预后寻找新的突破口。相关研究成果已发表于《Frontiers in Immunology》期刊。
为了回答上述问题，研究者主要采用了以下关键技术方法：首先，从公共数据库（GEO， 编号GSE231559）获取了结直肠癌原发灶和肝转移灶的单细胞RNA测序（scRNA-seq）数据。利用Seurat软件包进行数据质控、降维、聚类和细胞类型注释，并通过Harmony校正批次效应。通过inferCNV分析区分恶性与非恶性细胞。采用CytoTRACE和Slingshot进行细胞分化轨迹分析。利用CellChat和pySCENIC分别构建细胞间通讯网络和转录因子调控网络。此外，通过体外功能实验（包括siRNA转染、CCK-8、Transwell、伤口愈合和EdU实验）验证关键转录因子CEBPG的功能。基于TCGA数据库的转录组数据，通过Cox回归和LASSO回归构建了预后风险评分（BTRS）模型，并利用CIBERSORT等算法分析了免疫微环境特征和药物敏感性。
研究结果
scRNA-seq分析揭示CRC细胞异质性
通过单细胞转录组分析，研究者系统解析了结直肠癌组织中的细胞组成，鉴定出包括T/NK细胞、上皮细胞（EPCs）、巨噬细胞、成纤维细胞等在内的多种细胞类型。比较发现，原发灶与肝转移灶的细胞比例存在显著差异，且上皮细胞表现出较高的氧化磷酸化（OXPHOS）代谢活性，提示其与肿瘤进展密切相关。
CRC肿瘤细胞亚型的分子特征及其生物学功能富集分析
研究者进一步聚焦于上皮细胞，将其划分为五个不同的肿瘤细胞亚型：C0 ZFC3H1+TCs、C1 EEF1G+TCs、C2 PRAC1+TCs、C3 REG4+TCs和C4 BCL2L1+TCs。其中，C4 BCL2L1+TCs亚型在肝转移样本中显著富集，并且具有最高的干性评分，提示其可能具有更强的致瘤和转移潜能。
代谢谱和细胞干性揭示CRC肿瘤细胞亚型的功能异质性
代谢分析显示，C4 BCL2L1+TCs亚型特异性地高表达OXPHOS相关基因。轨迹分析（CytoTRACE和Slingshot）证实该亚型分化程度最低，处于更原始的细胞状态，这与其高干性和促转移特性相符。
CellChat分析揭示CRC微环境中复杂的细胞通讯网络和关键信号通路
细胞通讯分析表明，C4 BCL2L1+TCs亚型是肿瘤微环境中的通讯枢纽，与多种免疫细胞和基质细胞存在广泛而强烈的相互作用，特别是通过ADGRE5和MK（Midkine）信号通路。这提示该亚型可能通过调节免疫微环境来促进免疫逃逸和肿瘤进展。
系统分析CRC细胞中的转录因子调控模式
通过转录因子调控网络（pySCENIC）分析，研究者发现CEBPG是C4 BCL2L1+TCs亚型中活性最高的关键转录因子之一。其表达在该亚型中尤为突出，可能通过调控BCL2L1等下游基因来发挥作用。
敲低CEBPG抑制CRC细胞的增殖、迁移和侵袭
体外功能实验证实了CEBPG的关键作用。在SW48和SNU-81两种结直肠癌细胞系中，利用小干扰RNA（siRNA）敲低CEBPG后，细胞的增殖、迁移和侵袭能力均被显著抑制。这为CEBPG作为潜在治疗靶点提供了实验依据。
构建CRC的BTRS风险模型
基于C4 BCL2L1+TCs亚型的特征基因，研究者通过Cox回归和LASSO回归分析，构建了一个包含11个基因的预后风险评分模型（BTRS）。该模型能有效将患者分为高风险组和低风险组，高风险组患者总生存期更短，且模型在1年、3年和5年生存期预测中表现出良好的效能（AUC分别为0.72, 0.71, 0.67）。
转录组分析揭示不同BTRS亚组中预后相关分子及免疫异质性
对高低风险组的转录组比较分析发现，两组间存在显著的差异表达基因和免疫微环境特征差异。高风险组表现出免疫激活但功能可能受限的状态，例如活化的CD4记忆T细胞和M0型巨噬细胞比例升高，而自然杀伤（NK）细胞和树突状细胞比例降低。
利用风险评分全面评估药物敏感性和免疫浸润
进一步分析显示，高风险组患者对某些靶向药物（如吉非替尼、尼洛替尼等）可能更敏感。同时，高风险组中多个免疫检查点基因（如TNFSF18, HHLA2, CD44）表达上调，提示这部分患者可能从免疫检查点抑制剂治疗中获益。
研究结论与意义
本研究通过高分辨率的单细胞测序技术，成功“锁定”了一个在结直肠癌肝转移中扮演关键角色的肿瘤细胞亚型——C4 BCL2L1+TCs。该亚型具备高干性、活跃的氧化磷酸化代谢、强大的抗凋亡能力，并能作为枢纽广泛调控肿瘤微环境中的细胞通讯，从而协同促进肿瘤的进展、转移和免疫逃逸。研究进一步揭示了转录因子CEBPG是该亚型的关键上游调控因子，其表达能驱动结直肠癌细胞的恶性表型。
基于这些发现，本研究不仅从机制上深化了对结直肠癌肝转移异质性的理解，更实现了向临床应用的转化探索。一方面，CEBPG被确立为一个有潜力的治疗靶点，为开发针对该促转移亚型的新疗法提供了方向。另一方面，研究所构建的BTRS预后模型，整合了该特异亚型的分子特征与免疫微环境信息，能够超越传统的临床病理分期，更精准地预测患者预后，并初步提示了不同风险组患者对靶向治疗和免疫治疗的可能反应差异。
总之，这项工作为理解结直肠癌的进展机制提供了一套连接代谢适应、免疫逃逸与转移过程的完整框架，所鉴定的关键亚型、调控因子和预测模型，为未来实现结直肠癌的精准风险分层和个体化治疗策略制定奠定了重要的理论基础和资源储备。
