整合单细胞与批量RNA-seq分析鉴定头颈部鳞状细胞癌中与支链氨基酸（BCAA）代谢相关的预后特征

《Immunobiology》：Integrated single-cell and bulk RNA-seq analysis identifies BCAA metabolism-related prognostic signatures in head and neck squamous cell carcinoma

摘要 背景：支链氨基酸（Branched-Chain Amino Acids, BCAAs）代谢异常与癌症进展相关，但BCAA代谢在头颈部鳞状细胞癌（Head and Neck Squamous Cell Carcinoma, HNSCC）中的具体机制尚不清楚。本研究旨在鉴定HNSCC中与BCAA代谢相关的预后基因，并阐明其功能性作用。 方法：研究人员分析了HNSCC相关数据集（TCGA-HNSCC、GSE65858、GSE41613、GSE103322及GSE140042），通过差异表达分析获得的差异表达基因（DEGs）与加权基因共表达网络分析（WGCNA）中关联BCAA代谢相关基因（BCAA-MRGs）评分的关键模块基因取交集，筛选候选基因。随后采用单变量Cox回归分析、比例风险（PH）检验及最小绝对收缩和选择算子（LASSO）回归构建风险模型。研究人员进行了独立预后分析、富集分析及免疫微环境评估。此外，利用单细胞RNA测序（scRNA-seq）及拟时序（Pseudotime）分析在细胞水平评估预后基因表达。通过逆转录定量聚合酶链反应（RT-qPCR）验证HNSCC组织中基因表达水平。 结果：研究人员鉴定出SMS、PRDX6、GSTO1及ADA为风险模型的预后基因。HNSCC样本被分层为高风险组（HRG）与低风险组（LRG），LRG显示出显著更好的生存结局。基于风险评分和年龄的列线图（Nomogram）模型对HNSCC患者表现出优异的预测性能。此外，富集分析揭示戊糖磷酸途径及果糖和甘露糖代谢与HNSCC进展显著相关。HRG与LRG在20种免疫细胞类型（包括浆细胞样树突状细胞、肥大细胞及滤泡辅助T细胞）的浸润水平及10个免疫检查点（如CD276、CD27及CD40）的表达上存在差异。上皮细胞因高表达预后基因而被强调为关键细胞。值得注意的是，这些基因的表达趋势在不同细胞分化阶段存在差异。RT-qPCR验证显示SMS、GSTO1及ADA在肿瘤组中高表达，而PRDX6在肿瘤与正常组织间无显著差异。 结论：研究确定了一个四基因BCAA代谢相关特征（SMS、PRDX6、GSTO1、ADA），可将HNSCC患者划分为具有不同生存及免疫微环境的危险分组。该模型提示BCAA-免疫相互作用，特别是在上皮细胞中，为预后和治疗提供了初步见解。然而，这些结论主要基于计算预测，且RT-qPCR队列样本量小限制了实验验证强度，因此需进一步研究确认这些发现。

论文解读：整合单细胞与批量RNA-seq分析鉴定头颈部鳞状细胞癌中与支链氨基酸代谢相关的预后特征

研究背景与意义

头颈部鳞状细胞癌（Head and Neck Squamous Cell Carcinoma, HNSCC）是头颈部最常见的病理亚型，尽管治疗手段进步，其5年生存率仍仅约50%。HNSCC的异质性限制了传统分期对预后的预测精度，因此亟需新的预后标志物。异常能量代谢是肿瘤标志，支链氨基酸（Branched-Chain Amino Acids, BCAAs；包括亮氨酸、缬氨酸和异亮氨酸）代谢在肿瘤进展中起关键作用。BCAA耗竭可损害肿瘤浸润淋巴细胞功能，促进免疫逃逸。然而，BCAA代谢在HNSCC中的具体机制尚不清楚。本研究发表于《Immunobiology》，通过整合单细胞与批量RNA测序数据，旨在构建BCAA代谢相关的HNSCC预后模型，探讨其临床价值及与免疫微环境的关联。

主要关键技术方法

研究人员使用了TCGA-HNSCC（训练集，501例有完整生存信息）及GEO数据集（GSE65858、GSE41613为验证集；GSE140042、GSE103322为单细胞验证集）。首先从分子特征数据库（MSigDB）筛选27个BCAA代谢相关基因（BCAA-MRGs）。通过DESeq2进行差异表达分析（DEGs），利用加权基因共表达网络分析（WGCNA）计算BCAA-MRG评分（ssGSEA）并识别关键模块基因，取交集获得候选基因。采用单变量Cox回归、比例风险检验及LASSO回归筛选预后基因并构建风险模型（Risk Score）。通过Kaplan-Meier分析、ROC曲线、列线图（Nomogram）进行评估。利用ESTIMATE和ssGSEA分析肿瘤免疫微环境（TME），进行GSEA/GSVA富集分析。应用Seurat进行单细胞数据处理与细胞注释，Monocle2进行拟时序分析。最后收集临床样本通过RT-qPCR实验验证。

研究结果

3.1 识别8533个DEGs和1369个关键模块基因

研究人员在TCGA-HNSCC中识别出8533个DEGs（4360上调，4173下调）。WGCNA分析确定软阈值β=4，识别出11个模块，其中MEyellow模块与BCAA-MRG评分呈最强负相关（cor=-0.7288），筛选出1369个关键模块基因（|MM|>0.2, |GS|>0.2）。

3.2 257个候选基因的多功能表征与互作

将DEGs与关键模块基因取交集获得257个候选基因。GO/KEGG富集显示这些基因显著富集于有氧呼吸、ATP合成偶联电子传递、化学致癌-活性氧、氧化磷酸化等条目。蛋白互作（PPI）网络提示UQCRFS1、UQCRC2等为枢纽节点，表明BCAA代谢与核心能量代谢重编程耦合。

3.3 确认SMS、PRDX6、GSTO1和ADA为风险模型构建的预后基因

单变量Cox回归从257个候选基因中初筛出6个生存相关基因，经PH检验后纳入LASSO回归，最终确定SMS、PRDX6、GSTO1、ADA为预后基因。风险评分公式：Risk Score = 0.2138×SMS + 0.2144×PRDX6 + 0.0512×GSTO1 + 0.1510×ADA。以中位数划分HRG/LRG，HRG生存更差（p=0.0001）。模型在训练集及GSE65858、GSE41613验证集中AUC均>0.6。批量数据显示PRDX6在肿瘤组下调，其余三者上调；RT-qPCR验证SMS、GSTO1、ADA在肿瘤组显著高表达，PRDX6无显著差异。

3.4 列线图与临床特征的预测性能探索

单变量及多变量Cox分析确定风险评分和年龄为独立预后因素。基于此构建的列线图预测1/3/5年生存，校准曲线贴合理论线，AUC分别为0.894、0.678、0.678。敏感性分析显示风险评分在不同临床亚组（如Stage IV、T2、T3/4、N2/3）中保持一致预后趋势。HRG与LRG在免疫亚型上无差异，但临床分期有差异（p=0.013）。

3.5 鉴定与HNSCC预后显著相关的信号通路

GSEA显示HRG与LRG的差异基因富集于核糖体、亨廷顿病、扩张型心肌病等通路。GSVA显示果糖和甘露糖代谢、戊糖磷酸途径、自然杀伤（NK）细胞介导的细胞毒性等通路有差异。单变量Cox证实这些通路（除JAK/STAT）多促进HNSCC进展。

3.6 风险评分与HNSCC患者免疫微环境强相关

ESTIMATE显示LRG的基质、免疫及ESTIMATE评分更高。28种免疫细胞中20种有差异，LRG普遍高浸润（除CD56^brightNK细胞）。浆细胞样树突状细胞、滤泡辅助T细胞、CD56^brightNK细胞、17型辅助T细胞、肥大细胞为前5差异细胞。HRG高表达CD276、CD40等免疫检查点，LRG高表达CD27等。TIDE分析显示LRG对功能障碍评分更敏感，HRG对排斥评分更敏感。

3.7 识别与风险评分和TMB评估相关的药物

药物敏感性分析（oncoPredict）发现多拉马匹（doramapimod）等与风险评分正相关，多西他赛（docetaxel）等负相关。体细胞突变分析显示TP53、TTN为HRG/LRG共同高频突变，HRG的肿瘤突变负荷（TMB）更高，TMB与生存相关（cor=0.14）。

3.8 预后基因随不同细胞分化阶段变化

scRNA-seq（GSE140042、GSE103322）质控聚类注释出上皮细胞、T细胞等。上皮细胞因高表达预后基因被定为关键细胞。上皮细胞再分9个亚型，拟时序分析显示ADA和PRDX6表达先降后升，GSTO1先升后降，SMS变化不大，对应尿路上皮、减数分裂前生殖细胞等分化状态。

讨论与结论总结

讨论部分指出，SMS（精胺合酶）高表达通过多胺-mTOR轴促进HNSCC；PRDX6作为抗氧化酶在其预后角色中存在一定异质性；GSTO1通过ROS-BCAA-mTOR网络发挥作用；ADA通过腺苷-ATP-BCAA代谢联系免疫抑制。富集的戊糖磷酸途径与NK细胞毒性提示代谢-免疫串扰。HRG特征为低免疫浸润、高免疫检查点，CD56^brightNK细胞减少可能与BCAA竞争有关。药物预测提供了多拉马匹、紫杉类潜在方案。scRNA-seq将关键事件定位至上皮细胞分化轨迹。

结论（翻译）

本研究基于BCAAs构建了一个相对可靠的预后模型，全面探讨了这些基因在HNSCC中的潜在作用。研究结果为BCAA-MRGs在HNSCC中的作用机制提供了新见解。然而，本研究仅建立了相关性，未证实预后基因与HNSCC的因果关系。需进一步的实验研究验证所提出的机制。