阐明AML核糖体生物发生重塑预后评估与免疫微环境：整合单细胞和批量RNA分析

《Frontiers in Immunology》：Elucidating AML ribosome biogenesis reshaping prognostic assessment and immune microenvironment integrated single-cell and bulk RNA analysis

　　

背景：急性髓系白血病（AML）是一种侵袭性血液恶性肿瘤，因其异质性和高复发率面临重大临床挑战。失调的核糖体生物发生（ribosome biogenesis）是公认的肿瘤发生和治疗耐药驱动因素，但其对AML预后和免疫微环境的全面影响尚待充分阐明。方法：研究人员构建了一个AML单细胞图谱并量化了核糖体生物发生活性。通过整合机器学习开发了一种预后核糖体生物发生标志（RBS），并使用SHapley Additive exPlanations（SHAP）框架进行解释。研究人员评估了RBS与预后、肿瘤免疫微环境、免疫治疗应答和药物敏感性的关联。关键发现通过qPCR、Western blot和分子对接验证。结果：单细胞分析显示，在AML干细胞和祖细胞中，核糖体生物发生活性评分（RAS，一种量化细胞核糖体合成能力的复合指标）显著升高。使用随机生存森林算法建立的RBS模型显示出强大的预后能力。SHAP解释将EXOSC2确定为最重要的预测因子。低RBS患者表现出免疫激活的微环境和增强的免疫治疗应答。分子对接表明哇巴因和地高辛与EXOSC2具有高亲和力结合。结论：本研究描绘了核糖体生物发生在AML进展中的关键作用，并建立了一个可解释的RBS预后标志。该工具有效评估免疫治疗应答性并揭示新的靶向机会，为改善AML预后的临床决策提供了有价值的见解。

**论文解读：急性髓系白血病核糖体生物发生重塑预后评估与免疫微环境——整合单细胞与批量RNA分析**

**研究背景与问题**
急性髓系白血病（AML）是一种源自造血干/祖细胞的侵袭性血液恶性肿瘤，以分化阻滞、异常增殖和凋亡失调为特征，全球年新发病例超过20万，5年总生存率低于30%。现有治疗（如靶向疗法、免疫治疗）仍受限于肿瘤异质性、治疗耐药和免疫逃逸。核糖体生物发生（RiboSis）失调是肿瘤发生和耐药的公认驱动因素，但在AML中，其对预后和免疫微环境的综合影响尚不明确。此前研究多聚焦单基因标志物，缺乏基于单细胞分辨率的多基因整合分析。因此，本研究旨在通过单细胞转录组学解析RiboSis在AML中的动态调控，并构建可解释的预后标志（RBS），以指导风险分层和免疫治疗决策。

**为何开展研究**
鉴于AML的异质性和高复发率，亟需识别精准预后生物标志物和治疗靶点。RiboSis失调在肿瘤进化中起关键作用，但现有证据多来自批量分析，无法解析细胞亚群特异性。单细胞RNA测序（scRNA-seq）技术可精细描绘细胞异质性，为揭示RiboSis在AML干细胞和免疫微环境中的功能提供新视角。本研究旨在填补这一空白，开发临床可用的预后工具，并探索潜在治疗靶点。

**研究内容与结论**
研究人员通过整合Zenodo数据库（5例AML骨髓穿刺液和4例健康对照）及GSE116256队列（16例初诊AML）的scRNA-seq数据，构建了AML单细胞图谱。利用331个RiboSis相关基因（RRGs）计算核糖体生物发生活性评分（RAS），发现AML干细胞（HSCs）、粒细胞-巨噬细胞祖细胞（GMPs）和巨核-红系祖细胞（MEPs）中RAS显著升高，并伴随基因组不稳定性增强和代谢重编程（如糖酵解、氧化磷酸化、嘌呤代谢激活）。CellChat分析显示，高RAS亚群通过MIF-(CD74+CD44)轴驱动免疫抑制微环境。基于40个预后相关RRGs的共识聚类将AML分为两个亚型，其中Cluster A预后更差，富集代谢-免疫共失调通路。随后，通过整合10种机器学习算法，以随机生存森林（RSF）算法构建了包含9个核心基因（EXOSC2、LYAR、NSUN5、RPUSD2、NSUN3、RIOK3、ANG、RRP7A、RPUSD1）的RBS模型，并在GSE37642、GSE12417、BeatAML、TCGA-LAML、OHSU等独立队列中验证其稳健预后能力（C-index最高达0.75）。SHAP解释性分析将EXOSC2、LYAR和NSUN5识别为关键风险驱动因子，而RPUSD2和NSUN3为保护因子。高RBS患者呈现免疫抑制微环境（M2巨噬细胞增多、CD8⁺ T细胞耗竭、免疫排斥评分升高），低RBS患者则具有免疫激活状态（效应T细胞扩增、趋化因子CCL5/CXCL13上调）。分子对接显示哇巴因和地高辛与EXOSC2高亲和力结合（结合能< -8 kcal/mol），提示其潜在治疗价值。qRT-PCR和Western blot验证了EXOSC2在AML中显著上调，NSUN5下调。

**主要技术方法**
研究人员利用Zenodo和GSE116256队列的scRNA-seq数据（共5例AML、4例健康对照及16例AML）进行单细胞分析，通过Seurat处理、Harmony校正批次效应、InferCNV/CopyKAT分析拷贝数变异。RAS计算采用irGSEA包中的三种算法（singscore、ssgsea、viper）并整合为加权主成分分析（WPCA）得分。RBS开发基于GSE37642（训练集）和GSE12417（内部验证），采用10种机器学习算法组合，最终选择RSF模型。SHAP框架用于模型解释。免疫微环境通过IOBR包整合8种反卷积算法评估。药物敏感性使用oncoPredict基于GDSC数据库预测，分子对接通过AutoDock Vina进行。关键基因表达由qRT-PCR和Western blot验证（样本来自温州医科大学第三附属医院，8例AML和8例健康对照）。

**研究结果**
- **3.1 scRNA-seq分析AML**：通过UMAP降维和聚类，鉴定出15个细胞簇，注释为10种细胞类型，包括HSCs、GMPs、恶性细胞等。AML与健康样本细胞比例差异显著，且不同亚群富集不同生物学功能（如GO和KEGG通路）。
- **3.2 scRNA-seq中RRGs特征**：RAS在HSCs、GMPs、MEPs中显著升高，高RAS细胞显示基因组不稳定性（如del(5q)、-7/del(7q)）、HALLMARK通路激活（如G2M检查点、DNA修复、MYC靶点）和代谢重编程（糖酵解、氧化磷酸化）。
- **3.3 高/低RAS组间细胞通讯差异**：CellChat分析表明，高RAS组通过MIF-(CD74+CD44)轴增强免疫抑制信号，HSCs和GMPs成为主要通讯枢纽，而低RAS组以单核/巨噬细胞和T/NK细胞为信号节点。
- **3.4 基于RRGs的AML亚型及功能/免疫微环境特征**：共识聚类将AML分为Cluster A（预后差）和Cluster B（预后好）。Cluster A富集代谢和线粒体应激通路，免疫抑制细胞（MDSC、Tregs）浸润高；Cluster B适应性免疫增强，肿瘤纯度更高。
- **3.5 整合机器学习开发和验证RBS**：RSF模型在6个独立队列中表现稳健，优于8个已发表标志。多因素分析确认风险评分（HR=1.072）、年龄和RUNX1突变为独立预后因素。列线图整合风险评分和年龄，校准曲线和决策曲线分析显示良好临床净获益。
- **3.6 RSF模型特征重要性评估**：SHAP分析识别EXOSC2、LYAR、NSUN5为主要风险基因，RPUSD2和NSUN3为保护基因。个体样本分析显示，高风险基因表达高时风险评分升高。
- **3.7 RBS相关生物学机制**：高RBS样本富集能量代谢、核苷酸合成和DNA修复通路，低RBS样本富集自噬、MHC I类抗原呈递和表观遗传沉默通路。
- **3.8 免疫微环境、细胞因子和免疫治疗应答分析**：高RBS组显示免疫抑制（M2巨噬细胞、Tregs、IDO1升高，CD8⁺ T/NK细胞降低），免疫排斥评分和TAM M2评分高，CD8评分和Merck18评分低，提示免疫治疗应答差。
- **3.9 药物敏感性和分子对接**：高RBS患者对CDK4/6抑制剂、BCL-2抑制剂、HDAC抑制剂等敏感。EXOSC2与哇巴因（结合能-8.3 kcal/mol）和地高辛（-8.8 kcal/mol）高亲和力结合。
- **3.10 qRT-PCR和Western blot验证**：EXOSC2在AML中mRNA和蛋白水平上调，NSUN5下调，LYAR和RPUSD2的mRNA变化趋势与SHAP预测一致。

**讨论总结**
本研究通过单细胞转录组学揭示了RiboSis在AML中的动态调控，证明其通过促进代谢重编程和免疫逃逸驱动疾病进展。RBS模型作为多基因预后标志，在独立队列中稳健地将患者分为免疫激活和免疫抑制亚型，为免疫治疗分层提供依据。EXOSC2被鉴定为关键可靶向蛋白，心脏糖苷类化合物（哇巴因/地高辛）可能通过抑制EXOSC2发挥抗白血病作用。研究局限性包括回顾性设计、单细胞数据可能遗漏稀有相互作用、机制推断主要基于关联和计算，需未来功能验证。

**研究结论翻译**
本研究通过单细胞转录组学探讨了核糖体生物发生（RiboSis）在AML进展中的关键作用。研究人员发现，恶性祖细胞（HSCs/GMPs/MEPs）表现出特异性升高的RiboSis活性，并与基因组不稳定性、代谢重编程以及MIF-CD74诱导的免疫抑制微环境密切相关。研究人员建立了一个多算法整合的RBS预后模型，并应用SHAP框架阐明其特征基因的预后贡献，将预测性能与生物学机制联系起来。RBS将患者分为免疫激活（低RBS）和免疫抑制（高RBS）亚型，揭示了肿瘤微环境、免疫治疗应答和临床结局的显著差异。实验验证证实了核心基因的病理失调。关键的是，分子对接预测了哇巴因/地高辛与EXOSC2的高亲和力结合（结合能< -8 kcal/mol），证明了其可靶向性，并提出了新的临床转化途径。总之，本研究将核糖体生物发生定位为连接AML恶性转化与免疫逃逸的核心调控因子，为精准预后分层、免疫治疗获益人群筛选和个性化靶向策略开发提供了理论基础和实践指导。