《HUMAN MUTATION》：Identification of Immuno-Diagnostic and Therapeutic Markers in Tumors and Inflammatory Diseases Based on Single-Cell and Machine Learning Methods
编辑推荐：
背景：化疗是结直肠癌（CRC）全身治疗的基石，但耐药性和肿瘤异质性显著限制了患者生存期。线粒体代谢日益被认为是肿瘤进展、免疫格局和治疗反应的调节因子。然而，线粒体代谢相关转录模式在CRC中的预后意义及其与化疗反应特征和免疫特征的关系仍不完全清楚。
方法：分析了
背景：化疗是结直肠癌（CRC）全身治疗的基石，但耐药性和肿瘤异质性显著限制了患者生存期。线粒体代谢日益被认为是肿瘤进展、免疫格局和治疗反应的调节因子。然而，线粒体代谢相关转录模式在CRC中的预后意义及其与化疗反应特征和免疫特征的关系仍不完全清楚。
方法：分析了癌症基因组图谱结直肠癌队列（TCGA-CRC）及基因表达综合数据库队列GSE17537和GSE104645的临床和转录组数据。利用一组精选的线粒体代谢相关基因（MMRGs）计算MMRG评分，同时采用单样本基因集富集分析（ssGSEA），将GSE104645中化疗敏感病例上调的基因用于推导药物敏感性基因集评分。通过加权基因共表达网络分析（WGCNA）、差异表达分析、单变量Cox回归、最小绝对收缩和选择算子回归以及多变量Cox回归，构建了一个五基因风险评分（RiskScore）。采用Kaplan-Meier曲线、时间依赖性受试者工作特征曲线、外部验证和亚组分析评估预后关联性。评估了免疫细胞富集和pRRophetic预测的半最大抑制浓度值，并通过体外表型实验检测了CD36。
结果：肿瘤组织的MMRG评分低于正常组织，较高的MMRG评分与更好的总生存期（OS）相关。药物敏感性基因集评分在肿瘤中更高，但与OS无显著关联。由五个基因（CD36、CPA3、FCRL2、POU2F2和S100B）组成的特征生成的RiskScore在TCGA训练集、测试集、整体队列以及GSE17537队列中均能区分OS结局，但其时间依赖性辨别能力在不同队列和随访期间存在差异。在按分期、年龄和性别划分的亚组中也观察到类似的生存关联；但未评估其独立或增量预后价值。高风险肿瘤表现出改变的免疫细胞富集，经Benjamini-Hochberg校正后，RiskScore与pRRophetic预测的药物IC50值相关。CD36在SW1116细胞中的表达高于NCM460细胞，敲低CD36可抑制增殖、迁移和侵袭。在这些实验中未直接测量线粒体功能。
结论：五基因RiskScore与CRC的OS、免疫细胞富集及计算机模拟药物敏感性估计相关。它可能作为预后分层的一个候选工具，但其增量临床价值和治疗预测效用需前瞻性验证。CD36敲低实验支持其在SW1116细胞中促进肿瘤的表型，但并未确立其对线粒体代谢重编程的依赖性。
论文解读文章
一、研究背景与科学问题
结直肠癌（CRC）是全球癌症相关死亡的主要原因之一。尽管筛查技术、内镜手段和围手术期管理取得了显著进步，但许多患者在确诊时已处于局部晚期或转移阶段，且相当一部分早期患者最终会出现复发和远处转移。氟尿嘧啶类为基础的化疗（常联合奥沙利铂或伊立替康，如FOLFOX和FOLFIRI方案）是CRC系统性治疗的基石，并逐渐辅以抗血管内皮生长因子（VEGF）或抗表皮生长因子受体（EGFR）单克隆抗体等靶向治疗，以及在特定亚群中使用免疫检查点抑制剂。然而，原发性耐药和获得性耐药的快速出现是临床实践中面临的严峻挑战。一旦发生化疗耐药，治疗选择变得极为有限，患者常出现快速疾病进展和不良长期生存。此外，传统的临床病理参数，包括肿瘤-淋巴结-转移（TNM）分期和基本实验室指标，不足以准确捕捉CRC的生物学异质性，也无法可靠预测哪些患者将从化疗中获益以及哪些会早期复发。因此，识别稳健的、具有生物学基础的分子标志物，能够同时根据预后和化疗敏感性对CRC患者进行分层，对于实现更精细的风险分类和个体化治疗策略至关重要。
代谢重编程是多种癌症特征的基础，包括对应激的适应、转移和治疗抵抗。尽管历史上大量研究聚焦于“瓦伯格效应”（Warburg effect）和有氧糖酵解，但新兴证据强调线粒体代谢在许多癌症类型（包括CRC）中仍然高度相关。线粒体控制三羧酸（TCA）循环、脂肪酸氧化、氧化磷酸化（OXPHOS）、代谢物生成和氧化还原稳态，并密切参与程序性细胞死亡和细胞应激反应的调控。近年来，多项研究表明特定的线粒体代谢状态（如升高的OXPHOS能力或增强的脂肪酸氧化）可通过在应激下维持三磷酸腺苷（ATP）产生、缓冲活性氧（ROS）、调节线粒体外膜通透性和重塑凋亡阈值来赋予对细胞毒性药物或靶向治疗的抵抗。在CRC中，线粒体生物合成、线粒体DNA含量、OXPHOS相关基因表达或线粒体融合/分裂动力学的改变与更具侵袭性的表型、更高的转移潜能和更差的预后相关。然而，现有的多数CRC预后或预测模型集中于免疫或基质基因特征、糖酵解/脂质代谢基因集或上皮-间充质转化（EMT）标志物，而系统的、以线粒体为中心的分析——尤其是将线粒体代谢与化疗耐药和肿瘤免疫联系起来的研究——仍然稀缺。
近期CRC研究已开发出线粒体相关基因特征用于评估预后、表征免疫微环境和计算预测治疗反应。这些模型大多主要源自肿瘤转录组，并针对生存分类或免疫治疗相关推断进行了优化。相对较少的研究将化疗反应注释数据集与线粒体基因集评分、共表达网络分析、外部生存验证和免疫谱分析相结合。因此，本研究被设计为一项化疗反应信息指导下的预后探索性研究，而非直接的化疗获益临床预测工具。
这些见解共同表明，线粒体代谢可能作为连接化疗反应、恶性进展和免疫微环境调节的核心整合器。线粒体活性影响内在肿瘤细胞行为（如增殖、凋亡耐受、运动和EMT），也塑造肿瘤微环境（TME）的特征，如营养可用性、耗氧量、氧化还原平衡以及影响免疫细胞功能的代谢物可用性。例如，线粒体衍生的代谢物、ROS和线粒体动力学可以调节抗原呈递细胞活性、T细胞活化/耗竭、髓系细胞极化和免疫抑制性龛形成。同时，化疗对CRC细胞的暴露施加了代谢和氧化应激，可能选择性偏好具有更灵活或高OXPHOS线粒体表型的细胞亚群，从而促进耐药克隆的出现。因此，反映线粒体代谢的转录程序——特别是与化疗耐药和免疫格局相关的程序——有可能被用来推导一个与化疗反应相关转录特征、免疫特征和长期预后相关的基因特征。尽管已有数个线粒体特征被报道，但对这些特征的协调评估仍然有限。
本研究旨在开发和验证一个基于与CRC总生存期（OS）相关的线粒体代谢相关基因（MMRGs）的预后特征，并表征其与化疗反应衍生转录特征、免疫细胞富集和计算机模拟药物敏感性的关系。来自癌症基因组图谱（TCGA）的临床和RNA测序数据与包含一线化疗反应注释的GSE104645以及独立生存队列GSE17537进行了整合。在内部、外部和临床亚组分析中评估了预后性能。由于未进行与传统临床病理模型的直接比较或前瞻性治疗反应验证，该分析旨在为未来验证产生假设，而非建立临床独立性或指导治疗选择。
二、关键技术方法概述
研究人员从TCGA数据库获取结肠腺癌（COAD）和直肠腺癌（READ）的RNA-seq表达谱及临床数据，并从基因表达综合数据库（GEO）下载GSE17537（用于外部OS验证）和GSE104645（用于鉴定化疗反应相关基因）微阵列数据集。采用单样本基因集富集分析（ssGSEA）计算每个样本的MMRG评分和药物敏感性基因集评分。运用加权基因共表达网络分析（WGCNA）识别与这两个评分相关的基因模块。通过差异表达分析、单变量Cox回归、最小绝对收缩和选择算子（LASSO）回归以及多变量Cox回归，构建了一个五基因风险评分（RiskScore）。利用Kaplan-Meier（KM）生存分析和时间依赖性受试者工作特征（ROC）曲线评估预后表现。通过ssGSEA评估28种免疫细胞亚群的富集分数，并使用pRRophetic包预测药物的半最大抑制浓度（IC50）值。最后，在SW1116人结肠腺癌细胞系中通过小干扰RNA（siRNA）敲低CD36，进行细胞活力（CCK-8）、伤口愈合和Transwell侵袭实验。
三、主要研究结果
3.1 GSE104645中化疗反应相关基因的鉴定
通过对GSE104645数据集中化疗敏感组与耐药组进行差异表达分析，鉴定出118个基因（p < 0.05），其中51个在敏感组中上调，67个在耐药组中上调。这51个敏感组上调基因被用作特征集，通过ssGSEA计算药物敏感性基因集评分。
3.2 MMRG评分和药物敏感性评分的分布及预后价值分析
在TCGA-CRC队列中，肿瘤组织的MMRG评分显著低于邻近正常组织，而药物敏感性评分在肿瘤组中显著更高。KM生存分析显示，高MMRG评分患者的OS显著优于低评分患者，而高、低药物敏感性评分组之间的OS无显著差异（log-rank p = 0.71）。
3.3 通过WGCNA鉴定与MMRG和药物敏感性评分相关的关键基因模块
利用WGCNA构建共表达网络，鉴定出八个不同的共表达模块。其中，绿色模块与MMRG评分（r = -0.39, p = 1.64 × 10-22）和药物敏感性评分（r = -0.30, p = 1.98 × 10-13）均呈显著负相关，因此被选为与两个性状联合相关的模块。
3.4 差异表达基因与绿色模块的交集及功能富集分析
将TCGA-CRC肿瘤与正常组织的差异表达基因（DEGs）与绿色模块基因取交集，得到287个重叠基因。基因本体论（GO）生物学过程（BP）富集分析显示这些基因与免疫过程（如T细胞活化和白细胞增殖）密切相关，京都基因与基因组百科全书（KEGG）通路富集分析则涉及细胞粘附分子和趋化因子信号通路等。
3.5 风险模型构建与验证
在TCGA训练集中，通过顺序Cox和LASSO回归推导出一个五基因（CD36, CPA3, FCRL2, POU2F2, S100B）预后特征，计算公式为：RiskScore = (0.682 × CD36) + (-0.460 × CPA3) + (0.632 × FCRL2) + (0.542 × POU2F2) + (-0.448 × S100B)。根据最佳截断值将患者分为高低风险组。KM分析显示，在训练集、测试集、整个TCGA-CRC队列和GSE17537队列中，高风险组的OS均较短。时间依赖性AUC值在不同队列和随访时间点间有所变化（训练集0.68-0.77，测试集0.45-0.67，整个TCGA-CRC队列0.55-0.70，GSE17537队列0.61-0.67），表明预后辨别能力中等且因队列而异。
3.6 RiskScore的亚组生存关联
分层生存分析显示，RiskScore在不同临床亚组（早期/晚期分期、年龄≤67岁/>67岁、女性/男性）中均能将高风险组与较差的OS关联起来。但这些分析并未证明RiskScore相对于传统临床病理变量的独立或增量预后价值。
3.7 风险评分与免疫浸润特征及药物敏感性的关系分析
高低风险组在多种免疫细胞亚群的ssGSEA富集分数上存在差异，包括活化及记忆CD8+ T细胞、γδ T细胞、中性粒细胞、肥大细胞、巨噬细胞和髓源性抑制细胞（MDSCs）。经Benjamini-Hochberg校正后，RiskScore与乙胺嘧啶、MS.275、AZ628和硼替佐米的pRRophetic预测IC50值呈负相关，而与顺铂的预测IC50值呈正相关。
3.8 CD36在SW1116细胞中的表达及功能作用的实验验证
RT-qPCR结果显示，CD36、POU2F2和S100B在SW1116细胞中的mRNA表达显著高于正常结肠上皮细胞NCM460，而CPA3和FCRL2无显著差异。通过siRNA敲低CD36后，CCK-8实验表明SW1116细胞增殖能力受损，伤口愈合实验显示细胞迁移受到抑制，Transwell侵袭实验发现侵袭细胞数量减少。这些结果表明CD36在SW1116细胞中发挥促癌作用，但未直接测量线粒体功能，因此不能确定观察到的表型依赖于线粒体代谢重编程。
四、讨论与结论总结
讨论部分总结：
该研究通过整合MMRGs、化疗反应衍生转录评分及多个CRC转录组队列，构建了一个五基因RiskScore。该评分与OS、免疫细胞富集改变及计算机模拟药物敏感性估计相关联。与其他近期描述的线粒体特征相比，本模型的特点在于构建前将来自GSE104645的药物反应衍生ssGSEA性状纳入了WGCNA，并同时评估了免疫细胞富集和预测IC50相关性以及外部生存验证。然而，分析未直接比较RiskScore与临床病理参考模型或测试一致性、校准度或决策曲线获益。亚组分析显示了分期、年龄和性别分层内的生存关联，但未确立独立或增量预后价值。pRRophetic分析反映的是计算机模拟药物敏感性估计，而非观察到的化疗获益。
在生物学层面，富集和免疫细胞分析暗示了线粒体代谢-免疫微环境联系。287个交集基因富集于T细胞活化、白细胞增殖等免疫相关通路。高风险肿瘤中多种免疫细胞（包括活化/记忆CD8+ T细胞、γδ T细胞等）富集分数降低，但由于MDSCs等具有上下文依赖功能的免疫群体也降低，因此不支持将高风险组简单地归类为均匀免疫抑制或“冷”肿瘤。
在五个模型基因中，CD36作为连接线粒体代谢、脂质处理和恶性行为的代表候选基因被重点讨论。CD36是一种脂肪酸转位酶，介导长链脂肪酸摄取，支持能量生产和膜生物合成。CD36在SW1116细胞中高表达，敲低后抑制增殖、迁移和侵袭，这与促癌作用一致，但未直接测试其机制是否依赖线粒体。其他四个基因（CPA3、FCRL2、POU2F2、S100B）可能更多地将RiskScore与免疫细胞组成和肿瘤细胞状态联系起来，而非直接代表线粒体调节因子。
研究局限性包括：回顾性分析、基于批量转录组数据、缺乏多变量临床病理比较、外部验证队列样本量小（GSE17537仅55例）、免疫和药物敏感性为计算推断、功能验证仅限于单个CRC细胞系中的CD36且未进行直接线粒体功能测量。
研究结论翻译：
由CD36、CPA3、FCRL2、POU2F2和S100B组成的五基因线粒体代谢相关RiskScore与内部和外部CRC队列的OS相关，并在多个临床亚组中显示出相似的生存关联。该RiskScore还与免疫细胞富集改变和探索性计算机模拟药物敏感性估计相关。CD36敲低降低了SW1116细胞的增殖、迁移和侵袭，但未检验其线粒体依赖性机制。因此，该特征应被视为一个候选预后框架，需要前瞻性临床验证和进一步的机制研究，而非一个已确立的治疗选择工具。
