《Frontiers in Genetics》：Longitudinal single-cell analysis identifies a chemotherapy-suppressed proliferative program associated with adverse prognosis in neuroblastoma
编辑推荐：
背景： 高危神经母细胞瘤在强化治疗后常复发，但治疗相关的恶性神经母细胞程序对基线肿瘤风险的提示作用仍未完全明确。研究人员利用纵向单细胞/单核RNA测序（single-cell/single-nucleus RNA-seq）数据，识别治疗后恶性神经母细胞的转录变
背景： 高危神经母细胞瘤在强化治疗后常复发，但治疗相关的恶性神经母细胞程序对基线肿瘤风险的提示作用仍未完全明确。研究人员利用纵向单细胞/单核RNA测序（single-cell/single-nucleus RNA-seq）数据，识别治疗后恶性神经母细胞的转录变化，并在诊断性批量RNA测序（bulk RNA-seq）队列中评估其预后相关性。
方法： 研究人员分析了诊断（DX）和治疗后（PTX）神经母细胞瘤样本的纵向单细胞/单核RNA-seq图谱。将神经母细胞按供体和治疗状态水平聚合为伪批量（pseudo-bulk）表达谱，并进行配对差异表达分析，以识别PTX调控的恶性神经母细胞程序。通过Hallmark过表达分析对PTX下调（PTX-down）和PTX上调（PTX-up）基因集进行功能表征，并在批量队列中通过基因集变异分析（GSVA）进一步估计Hallmark通路活性。随后，研究人员通过计算PTX-down基因的队列内平均z-score（PTX-down增殖评分），将该程序投射到诊断性批量RNA-seq队列中。在一个主要诊断性批量RNA-seq验证队列和三个独立外部诊断性批量RNA-seq队列中评估其预后价值，并进行了临床模型增量、跨队列荟萃分析、临床特征关联及Hallmark相关性分析。
结果： 配对伪批量分析显示，PTX恶性神经母细胞谱中富集E2F靶点、G2/M检查点、MYC靶点、DNA修复和有丝分裂纺锤体通路的增殖-复制程序活性降低，同时相对富集与缺氧、炎症信号、凋亡和上皮-间质转化相关的适应状态特征。在主要验证队列中，高PTX-down增殖评分与较差的总生存期（OS）和无事件生存期（EFS）相关，且在调整年龄、COG定义的高危分类、MYCN状态和INSS分期后仍具有预后价值。将PTX-down评分加入模型后，C-index及1年、3年和5年时间依赖性AUC均优于COG模型和临床模型。在四个批量RNA-seq队列中，较高的PTX-down评分一致地与不良OS相关（随机效应合并HR = 1.83，95% CI 1.20–2.80；P = 0.005）。该评分在MYCN扩增、4期疾病和年龄≥18个月的肿瘤中也较高。核心基因和Hallmark相关性分析将PTX-down评分定义为E2F/G2M/MYC驱动的增殖轴。
结论： 纵向单细胞分析识别出恶性神经母细胞中化疗抑制的增殖程序。将该程序投射到诊断性批量肿瘤中可捕获高危生物学特征，改善预后判别能力，并预测神经母细胞瘤的不良结局。
神经母细胞瘤化疗抑制性增殖程序：从纵向单细胞分析到诊断性预后标志物的转化研究解读
一、研究背景与科学问题
神经母细胞瘤（neuroblastoma）是儿童最常见的颅外实体肿瘤，也是儿科肿瘤学中癌症相关死亡的主要原因之一。尽管多模式治疗改善了部分患者的预后，高危神经母细胞瘤仍以复发和治疗耐药为特征。肿瘤在强化治疗后的复发提示，恶性神经母细胞经历了深刻的治疗相关重塑，而治疗塑造的转录程序可能包含关于肿瘤侵袭性的临床相关信息。神经母细胞瘤的生物学异质性已被广泛联系到神经嵴发育、发育细胞状态、MYCN激活、增殖信号、染色质调控、免疫微环境特征和谱系可塑性等方面。基因组和转录组研究已确立，侵袭性神经母细胞瘤由复发性高危改变和转录状态多样性塑造，包括MYCN扩增、结构性染色体改变、端粒维持机制和发育可塑性。近年来的单细胞、单核、多组学和空间转录组研究进一步表明，神经母细胞瘤包含多样的恶性和非恶性区室，恶性神经母细胞占据与肾上腺能（ADRN）身份、间充质（MES）或应激相关程序、增殖活性和发育可塑性相关的状态。然而，大多数预后转录组分析仍基于诊断性批量肿瘤谱。这些批量分析可以识别与结局相关的程序，但往往缺乏关于治疗后肿瘤细胞状态如何变化的直接信息。因此，化疗相关的恶性神经母细胞重塑与基线预后生物学之间的关系仍未完全阐明。
纵向单细胞或单核RNA测序数据集提供了直接比较治疗前后恶性神经母细胞的独特机会。近期对高危神经母细胞瘤诱导化疗前后的纵向分析已证明，肿瘤和微环境区室均发生显著的治疗相关重塑，不良结局与增殖性和代谢活跃的肿瘤状态相关，而更分化的神经元样状态与较有利的生物学相关。通过聚焦配对的诊断和治疗后肿瘤细胞谱，可以识别在诱导治疗后受抑制或保留的转录程序。此类程序可能反映治疗敏感的增殖状态、残留适应状态或应激相关表型。重要的是，如果一个在治疗后受抑制的程序在诊断时高度活跃，它可能标记具有高基线增殖能力和侵袭性临床行为的肿瘤。因此，治疗相关的单细胞程序不仅可作为治疗后重塑的描述符，还可作为诊断性肿瘤风险的生物学可解释标志物。
基于上述背景，本研究分析了一个纵向单细胞/单核RNA-seq神经母细胞瘤图谱，以定义恶性神经母细胞中化疗相关的转录程序。研究人员将神经母细胞按供体和治疗状态聚合成配对伪批量谱，实现供体感知的诊断（DX）与治疗后（PTX）样本比较。该分析识别出一个PTX下调程序，富集E2F靶点、G2/M检查点、MYC靶点、DNA修复和有丝分裂纺锤体通路，共同代表化疗抑制的增殖-复制轴。同时，PTX上调基因富集缺氧、炎症信号、凋亡和上皮-间质转化相关程序，与治疗相关的适应状态一致。随后，研究人员将该PTX下调程序作为PTX-down增殖评分投射到诊断性批量RNA-seq队列中，并在多个神经母细胞瘤数据集中评估其临床相关性。
二、主要关键技术方法
本研究采用的核心技术方法包括：基于纵向单细胞/单核RNA测序图谱的恶性神经母细胞伪批量生成，利用来源图谱的肿瘤细胞注释，按供体和治疗状态聚合为44个伪批量谱（来自22个配对供体）；使用edgeR进行配对伪批量差异表达分析，识别出1,478个PTX-down和1,601个PTX-up基因（FDR < 0.05）；通过MSigDB Hallmark基因集进行过度代表分析表征基因集功能；在批量RNA-seq队列中，采用队列内基因-wise z-score标准化策略计算PTX-down增殖评分，即PTX-down基因的平均z-score；在主要验证队列GSE62564（n = 498）中开展Kaplan-Meier生存分析、Cox比例风险回归、C-index比较和时间依赖性AUC分析（基于IPCW的timeROC框架）；跨GSE181559（n = 97）、TARGET-NBL（n = 161）和fwr144（n = 144）三个外部队列进行交叉验证和随机效应荟萃分析；以及通过GSVA估计Hallmark通路活性并进行Spearman相关性分析。此外还进行了随机基因集零分布检验以评估评分的特异性。
三、研究结果
纵向单细胞分析定义恶性神经母细胞区室用于配对伪批量分析： 通过UMAP可视化展示了完整图谱中的多种细胞群体，包括神经母细胞、免疫细胞、基质细胞、内皮细胞、施万细胞、肾上腺皮质细胞、肝细胞和肾脏相关细胞。MDS分析显示DX和PTX谱之间存在治疗相关的分离，支持使用供体感知的配对伪批量差异表达框架。
PTX恶性神经母细胞谱显示增殖-复制程序降低和适应状态特征相对富集： 配对差异表达分析识别出1,478个PTX-down和1,601个PTX-up基因。PTX-down基因包括MKI67、TOP2A、PCNA、MCM2、MCM4、MCM6、MCM10、CDC45、E2F1、UHRF1、BIRC5和AURKB等典型增殖、细胞周期和DNA复制基因；PTX-up基因包括CALCA、ADCYAP1、VIP、DDIT3、HBEGF、NR4A1和NR4A3等神经内分泌、应激反应和适应状态相关基因。Hallmark富集分析显示，PTX-down基因富集于E2F靶点、G2/M检查点、MYC靶点、DNA修复、有丝分裂纺锤体、mTORC1信号和糖酵解通路；PTX-up基因富集于TNFA信号通过NFKB、缺氧、炎症反应、IL6/JAK/STAT3信号、凋亡、血红素代谢、TGF-beta信号、上皮-间质转化、KRAS信号和补体通路。
PTX-down增殖评分预测GSE62564中的不良OS和EFS： 在GSE62564中，高PTX-down评分组OS的HR为6.84（95% CI 4.06–11.51，log-rank P < 0.001），EFS的HR为3.62（95% CI 2.61–5.02，log-rank P < 0.001）。多变量Cox回归在调整年龄、COG高危分类、MYCN状态和INSS分期后，PTX-down评分仍与OS和EFS风险增加相关。敏感性分析显示，使用top 200、top 300或top 400基因集以及不同排序策略时，评分与不良OS和EFS的关联保持稳健。随机基因集零分布检验的实证单侧P值对OS和EFS均<0.001，支持PTX-down程序的特异性。
PTX-down增殖评分改善超越临床风险变量的预后判别能力： 在GSE62564中，将PTX-down评分加入COG模型和常规临床模型（年龄、MYCN状态、INSS分期）后，OS和EFS的C-index均增加。时间依赖性AUC分析显示，在1年、3年和5年预测区间，加入PTX-down评分一致地改善OS和EFS的AUC。COG模型中，OS的AUC改善范围为+0.035至+0.058，EFS为+0.059至+0.070；临床模型中，OS的AUC改善范围为+0.020至+0.025，EFS为+0.026至+0.029。
跨队列验证将PTX-down增殖评分与不良OS和高危临床生物学联系起来： 在四个批量RNA-seq队列中，较高的PTX-down评分对OS均显示一致的不良方向。随机效应荟萃分析确认了显著合并关联，合并HR为1.83（95% CI 1.20–2.80，P = 0.005），但队列间异质性较大（I2 = 87.7%）。在所有四个队列中，该评分在MYCN扩增、4期疾病和年龄≥18个月的肿瘤中均高于相应参考组。
PTX-down评分反映诊断性肿瘤中的E2F/G2M/MYC增殖轴： 核心基因热图显示，高PTX-down评分肿瘤中E2F1、E2F2、MYBL2、UHRF1、PCNA、MCM家族、CDC45、EXO1、MKI67、TOP2A、AURKB、BIRC5、CDK1、CCNB1、CCNB2、CENPF、KIF11和KIF18B等增殖相关基因协调上调。通路相关性分析显示，PTX-down评分与E2F靶点、G2/M检查点、MYC靶点、DNA修复和有丝分裂纺锤体的Hallmark活性强正相关，而PTX-up评分与这些增殖通路呈负相关，与凋亡、缺氧、TNFA信号和上皮-间质转化通路正相关。
四、讨论总结与结论翻译
讨论部分指出，本研究将治疗相关的单细胞肿瘤状态与诊断性风险分层联系起来，识别出一个以E2F靶点、G2/M检查点、MYC靶点、DNA修复和有丝分裂纺锤体等增殖相关通路为特征的PTX-down恶性神经母细胞程序。该评分与既往基于批量转录组的预后标志物（如铁死亡相关、免疫或微环境相关、癌症相关成纤维细胞相关模型）在概念上不同，其核心创新在于首先从配对纵向恶性神经母细胞伪批量分析中定义为诱导治疗后降低的程序，随后投射到诊断性批量转录组，提供治疗信息化的肿瘤细胞固有增殖-复制状态读数。PTX-down与PTX-up程序的相互对立关系与药物耐受持久细胞概念一致，即残留恶性细胞抑制MYC/MYCN连接的增殖程序并激活炎症或应激反应通路。研究还提出，PTX-down高肿瘤可能值得未来功能评估沿增殖-复制和复制应激轴的脆弱性。局限性包括：单细胞发现数据集受样本可用性、治疗时机和配对标本结构限制；批量验证队列在表达单位、临床变量、终点定义、治疗背景和随访时长方面存在差异；GSE62564同时作为主要验证队列和荟萃分析队列之一可能向上偏倚合并估计；多变量Cox模型调整的变量存在部分共线性；评分基于回顾性公共数据集评估，需前瞻性验证后方可临床转化。
结论部分翻译：纵向单细胞分析识别出恶性神经母细胞中化疗抑制的增殖-复制程序。当投射到诊断性批量RNA-seq队列时，该PTX-down程序的基线活性预测了不良生存，改善了超越临床风险变量的预后判别能力，并追踪了已确立的高危生物学特征。这些发现将PTX-down增殖评分定义为生物学可解释的E2F/G2M/MYC驱动轴，将化疗抑制的肿瘤细胞重塑与高危神经母细胞瘤生物学和不良预后联系起来。
