《Oxygen》：Identification and Characterization of Oxidative Stress-Related Prognostic Biomarkers in Acute Myeloid Leukemia Using Integrated Bioinformatics Analysis
编辑推荐：
氧化应激（Oxidative Stress）促进急性髓系白血病（Acute Myeloid Leukemia, AML）的生物学进程和治疗耐药性，但其预后相关性尚未完全阐明。研究人员分析了来自GSE37642数据集的553例患者，将其中417例患者的配对GPL
氧化应激（Oxidative Stress）促进急性髓系白血病（Acute Myeloid Leukemia, AML）的生物学进程和治疗耐药性，但其预后相关性尚未完全阐明。研究人员分析了来自GSE37642数据集的553例患者，将其中417例患者的配对GPL96/GPL97测量值作为发现队列，而GPL570（n=136）则作为平台定义的内部-外部验证队列。从Enrichr基因本体论生物学过程2023库中检索了注释为细胞对氧化应激反应（GO:0034599）的人类基因；在两个队列定义中共有110个基因可被检测。在验证前，发现队列独有的预处理步骤被固定。17个名义单变量Cox p<0.05的基因进入交叉验证的LASSO-Cox建模，最终保留了包含13个基因的特征集（ATP13A2, AKT1, DAPK1, HSPA1A, LRRK2, CPEB2, CAMKK2, CCS, PRKAA1, RBM11, TP53, GPX1, and TPM1）。使用固定的发现队列衍生系数和中位截断值（0.00807），该评分在发现队列中区分了高风险和低风险患者（对数秩p=3.51×10-12；连续评分HR=3.24，95%CI 2.46-4.29；C指数=0.627，95%CI 0.591-0.661），并在平台定义的验证队列中仍与生存期相关（对数秩p=0.019；HR=1.90，95%CI 1.20-3.00；C指数=0.580，95%CI 0.519-0.641）。验证队列的曲线下面积（AUC）在1年、3年和5年分别为0.590、0.621和0.655。完全嵌套的五折诊断分析得出平均留出C指数为0.571。在验证队列中，将风险评分加入年龄改善了嵌套模型拟合度（似然比p=0.028），但判别能力的改善很小（ΔC指数=0.019），其配对袋外自助法95%置信区间包含零（-0.070至0.056）；相应的ΔAUC区间也包含零。Schoenfeld残差诊断未提示比例风险假设违反。这些发现支持了一个生物学上可解释的、与氧化应激相关的预后信号，具有适度的跨平台判别能力，同时强调了在临床应用前需要进行独立的当代验证。
论文解读文章
一、研究背景与目的
急性髓系白血病是一种侵袭性极强的克隆性髓系前体细胞疾病，尽管欧洲白血病网（European LeukemiaNet, ELN）分类利用细胞遗传学和分子特征改进了风险分层，并且近年来涌现了大量新型靶向和表观遗传药物，但仍有大量患者，尤其是不良遗传风险类别者，预后不佳。现有遗传风险分类未能完全捕捉影响治疗应答的全部因素。氧化还原失调作为一种跨越遗传亚型的生物学过程，日益受到关注。恶性细胞常利用过量的活性氧（Reactive Oxygen Species, ROS）支持增殖、逃避凋亡和抵抗治疗。在AML中，母细胞处于持续升高的氧化应激状态，这反过来又加剧基因组不稳定性并维持白血病克隆的存活信号。然而，系统地评估氧化应激相关基因在AML中的预后价值的研究尚属空白。因此，本研究旨在利用一项可重复的、基于基因本体论（Gene Ontology, GO）的策略，鉴定并表征与AML总体生存期相关的氧化应激相关转录组特征，并量化其在不同芯片平台间的可移植性和增量预后信息。
二、主要研究方法
研究人员采用回顾性生物信息学研究框架。数据来源于NCBI GEO数据库的GSE37642数据集。研究队列分为两部分：发现队列由417例患者组成，其基因表达谱来自配对的GPL96和GPL97平台；验证队列由136例患者组成，其数据来自GPL570平台，被称为平台定义的内部-外部验证队列。氧化应激相关基因通过程序化检索Enrichr GO Biological Process 2023库中注释为“细胞对氧化应激反应”（GO:0034599）的基因来定义，最终纳入110个可检测基因。在发现队列中进行单变量Cox比例风险回归筛选，将满足名义p<0.05标准的17个基因引入LASSO-Cox回归模型进行特征选择和模型开发。模型开发后，所有发现队列衍生的预处理参数（包括中位数填补值、均值和标准差）、模型系数以及风险评分截断值（中位数0.008067）被“冻结”，并原封不动地应用于验证队列，无任何重新校准。研究还进行了多项稳健性分析，包括完全嵌套的五折交叉验证、200次自助法重采样稳定性评估、年龄模型比较、Schoenfeld比例风险假设检验、替代惩罚模型敏感性分析以及平台偏移评估。
三、主要研究结果
1. 患者特征
最终生存分析数据集包含553例AML患者：417例发现队列患者（308例死亡，占73.9%）和136例GPL570验证队列患者（98例死亡，占72.1%）。患者身份审计确认了GPL96和GPL97的完全配对，且发现队列与GPL570验证队列之间无患者重叠。发现队列中位年龄为57岁，验证队列为59.5岁；中位总生存期分别为316天和440.5天。
2. 单变量Cox回归揭示与总生存期相关的氧化应激相关基因
在发现队列中对110个可测的GO:0034599基因进行单变量Cox回归分析。结果显示，17个基因满足预设的名义p<0.05筛选标准，并被带入后续的LASSO-Cox回归模型。
3. 使用LASSO-Cox回归开发氧化应激相关预后特征
通过十折交叉验证选择的最优正则化参数alpha为0.01669，最终保留了13个非零系数的基因：ATP13A2, AKT1, DAPK1, HSPA1A, LRRK2, CPEB2, CAMKK2, CCS, PRKAA1, RBM11, TP53, GPX1和TPM1。其中，ATP13A2的正系数最大（0.203），AKT1的负系数绝对值最大（-0.184）。基于此，计算了个体化的加权风险评分。
4. 氧化应激相关预后特征的验证
使用固定的发现队列衍生系数、预处理参数和截断值，在发现队列中将209例患者划分为高风险，208例为低风险；在验证队列中，66例为高风险，70例为低风险。Kaplan-Meier分析显示，两组在发现队列（对数秩p=3.51×10-12）和验证队列（对数秩p=0.019）中的生存差异均有统计学意义。连续风险评分与死亡率显著相关，在发现队列中风险比（Hazard Ratio, HR）为3.24，在验证队列中为1.90。Harrell's C指数在发现队列为0.627，在验证队列为0.580，表明跨平台判别能力适中。
5. 时间依赖性预测性能
时间依赖性受试者工作特征曲线下面积（Area Under the Curve, AUC）分析显示，在发现队列中，1年、3年和5年的AUC分别为0.687、0.750和0.726。在验证队列中，相应的AUC分别为0.590、0.621和0.655，进一步证实了跨平台转移后适度的判别能力。
6. 总体预后表现
固定模型的总体表现总结为：发现队列C指数0.627，验证队列0.580；发现队列连续评分HR为3.24，验证队列为1.90。这些结果支持一个可复现的生存关联，但跨平台转移后的个体水平判别能力有限。
7. 增量和多变量预后评估
通过比较仅年龄、仅风险和年龄加风险模型来评估增量预后价值。在发现队列中，年龄加风险模型相比仅年龄模型在嵌套模型似然比检验中显著更优（p=8.57×10-11），但配对袋外自助法95%置信区间包含零。在验证队列中，似然比检验同样显示年龄加风险模型拟合更佳（p=0.028），但ΔC指数和ΔAUC的95%置信区间均包含零。在包含年龄、RUNX1-RUNX1T1融合状态、RUNX1突变状态和风险评分的完整病例多变量Cox模型中，风险评分在发现队列中仍与总生存期独立相关，但在验证队列中方向一致但未达到传统显著性水平（p=0.065）。Schoenfeld残差诊断未发现比例风险假设的违反证据。
8. 稳健性和敏感性分析
完全嵌套的五折管道诊断的平均留出C指数为0.571。在200次自助法重采样中，13个最终基因的选择频率从0.675到1.000不等，其中12个基因在至少75%的重采样中被选中。弹性网络（Elastic Net）产生了与LASSO相似的C指数，而岭回归样（ridge-like）惩罚模型表现较差。平台偏移诊断显示，最终基因的平均z分数偏移较小，其中TPM1的偏移最大。这些分析支持主要信号的稳定性，同时强化了验证队列中适度的判别能力。
四、讨论与结论总结
讨论部分总结： 该13基因氧化应激特征能够在发现队列中按生存期区分AML患者，并在固定转移到GPL570平台后仍保持统计学上可检测的关联。研究认为，该模型应被解释为一个具有适度判别能力的生物学知情预后信号，而非临床即用型分类器。该风险评分最可能捕获的是白血病克隆协调应激抵抗、损伤修复、代谢适应和调节性细胞死亡的能力，而非直接的ROS读数。研究深入讨论了各个基因（如ATP13A2, DAPK1, HSPA1A, AKT1等）在AML背景下可能的生物学角色及其与系数方向的潜在关系，强调了转录组水平的关联不能直接等同于功能或蛋白水平的因果关系。研究指出，临床上相关的氧化应激并非单一抗氧化通路的作用，而是一个涵盖代谢、细胞器质量控制、RNA调控、蛋白质稳态、凋亡和微环境适应的协调细胞状态。
研究结论： 该13基因组合标记了一种与AML生存期相关的氧化应激相关转录表型。这种关联在发现队列中最强，并在固定转移到GPL570后仍可检测到，但跨平台判别能力适中，且超越年龄的增量判别能力在配对自助法区间内并不稳健。因此，该评分不能作为独立的临床测试，也不应被解释为经过独立的外部验证。其当前的价值在于产生假说：它确定了一个可复现的氧化还原适应轴，需要在具有完整ELN变量的独立当代AML队列和正交实验验证中进行检验。
