《Journal of Radiation Research and Applied Sciences》：A radiotherapy-related five-gene signature predicts prognosis and regulates radiosensitivity in neuroblastoma
编辑推荐：
本研究旨在通过整合转录组分析、机器学习、免疫分析、调控网络推断和体外验证，确定神经母细胞瘤中放疗相关的预后生物标志物和放射增敏候选靶点。从照射相关数据集中选择放疗相关基因，并在两个神经母细胞瘤患者队列中构建和验证预后模型。进行了功能和免疫关联的计算分析，并评估
本研究旨在通过整合转录组分析、机器学习、免疫分析、调控网络推断和体外验证，确定神经母细胞瘤中放疗相关的预后生物标志物和放射增敏候选靶点。从照射相关数据集中选择放疗相关基因，并在两个神经母细胞瘤患者队列中构建和验证预后模型。进行了功能和免疫关联的计算分析，并评估了候选基因沉默对神经母细胞瘤细胞系照射反应的影响。研究人员确定了一种放疗相关的五基因特征（CDCA8、CDCA2、FAM83D、ARHGAP11A和GAS2L3），可用于将神经母细胞瘤患者分为不同的生存风险组，并在独立队列中进行了外部验证。在GSE49710（HR=3.80, P<0.001）和GSE85047（HR=4.33, P<0.001）中，高风险组患者的总生存期均显著较差。功能分析将五基因特征与细胞周期调控和DNA损伤修复联系起来，而探索性免疫分析和预测的lncRNA–miRNA–mRNA网络分析提示了与免疫背景和上游非编码RNA调控的潜在关联。分级照射实验确定SK-N-BE(2)和SK-N-AS在五个神经母细胞瘤细胞系中相对放射抵抗。在这两种细胞系中，siRNA介导的五个基因沉默增强了照射诱导的生长抑制并降低了克隆存活。进一步的γ-H2AX分析显示，在6 Gy照射后24小时，CDCA8或CDCA2的沉默与两种细胞系中残留的γ-H2AX灶数量增加相关。本研究提出了一种与放疗相关的神经母细胞瘤五基因预后特征。这五个基因与预后和增殖相关通路相关。它们的沉默增强了照射诱导的生长抑制，并降低了体外克隆存活，但该特征可能部分归因于一般的增殖性肿瘤生物学，而不一定是放疗特异性机制。
神经母细胞瘤(neuroblastoma, NBL)是儿童最常见的颅外实体肿瘤，起源于神经嵴细胞衍生的交感神经系统前体细胞，临床行为高度异质，高危患者预后差。放疗是局部控制的重要手段，但肿瘤对放射的反应存在显著异质性，这与DNA修复能力、细胞周期检查点、死亡信号通路及微环境等因素有关。目前多数NBL预后模型基于肿瘤初始生物学特征，未充分整合放疗相关转录变化和功能反应信息。因此，建立放疗相关预后特征并识别放射增敏靶点具有重要临床意义。
研究人员利用基因表达综合数据库(GEO)的转录组数据，结合机器学习、免疫分析和体外实验，系统开发了放疗相关五基因预后特征。研究首先从GSE197124（含24例放疗组和24例非放疗组）中筛选出32个放疗相关差异表达基因；在GSE49710训练集（483例）中通过Cox回归和117种算法组合构建模型，并在GSE85047验证集（269例）中外部验证。最终确定由CDCA8、CDCA2、FAM83D、ARHGAP11A和GAS2L3组成的五基因特征。该特征能有效区分高低风险患者，高危组总生存显著较差。体外实验表明，沉默这些基因可增强放射诱导的生长抑制和克隆死亡。研究发表在《Journal of Radiation Research and Applied Sciences》上，为NBL放疗反应预测和放射增敏提供了候选靶点。
主要技术方法如下：研究人员利用GSE197124、GSE49710和GSE85047三个数据集；采用limma筛选差异表达基因，通过GO/KEGG和STRING/PPI分析功能；Cox回归筛选预后基因，并应用117种机器学习算法组合（包括LASSO、RSF、SuperPC等）构建模型；使用ssGSEA和ESTIMATE评估免疫微环境；利用miRDB和StarBase预测ceRNA网络；体外使用SK-N-BE(2)、SK-N-AS、IMR-32、SK-N-MC、SH-SY5Y及hiPSC-NCCs细胞，通过CCK-8、RT-qPCR、siRNA转染、克隆形成和γ-H2AX免疫荧光实验验证。
结果如下：
3.1 放疗相关基因的鉴定和功能表征：在GSE197124中，与未照射组相比，照射组有32个基因差异表达，其中11个上调、21个下调。GO和KEGG分析显示，这些基因富集于神经元凋亡、神经营养因子-TRK信号、微管相关结构和凋亡代谢通路。STRING构建的PPI网络包含21个节点和78条边，提示基因间存在潜在相互作用。
3.2 基于Cox回归和机器学习的预后特征构建：在GSE49710训练集中，经比例风险假说检验和单变量Cox回归，6个基因与总生存显著相关。通过比较117种算法组合的C-index，最终选择SuperPC模型，保留CDCA8、CDCA2、FAM83D、ARHGAP11A和GAS2L3五个基因。这些基因的HR均大于1，提示高表达与死亡风险增加有关。
3.3 预后特征的外部验证和效能评估：按各队列中位风险评分分组，训练集和验证集中高危组患者总生存均显著较差，HR分别为3.80和4.33。时间依赖ROC分析显示，该模型在训练集1、3、5年AUC分别为0.781、0.684、0.653，验证集为0.704、0.761、0.782，表明具有跨队列的稳定预测能力。
3.4 核心风险基因的预后意义和功能注释：五个核心基因各自的高表达均与较差总生存相关。GSEA分析显示，这些基因富集于DNA复制、细胞周期、错配修复、剪接体、核苷酸切除修复、细胞粘附分子以及免疫相关通路，如移植物排斥反应、移植物抗宿主病等。
3.5 免疫微环境表征和调控网络分析：ssGSEA显示28种免疫细胞中有27种在高低风险组间存在显著差异。ESTIMATE分析显示低风险组具有更高的基质评分、免疫评分和ESTIMATE评分。通过miRDB和StarBase预测到90个miRNA和53个lncRNA，构建了lncRNA-miRNA-mRNA调控网络。
3.6 相对放射抵抗细胞模型选择和核心基因表达验证：经0-8 Gy梯度照射，SK-N-BE(2)和SK-N-AS细胞活力相对较高，被选为相对放射抵抗模型。RT-qPCR显示五个基因在NBL细胞系中的表达高于hiPSC-NCCs，且在所选两种细胞中表达最高。
3.7 沉默核心基因增强放射诱导的生长抑制：siRNA敲低任一基因后，6 Gy照射下的细胞活力较单独照射或单独敲低显著下降，联合处理组活力最低。
3.8 沉默核心基因增强放射诱导的克隆存活抑制：克隆形成实验显示，敲低任一基因与照射联合均显著降低克隆存活率。γ-H2AX分析表明，CDCA8或CDCA2沉默增加照射后24小时残留灶数量，提示DNA双链断裂修复可能受损。
讨论部分：研究人员指出，五基因特征与细胞周期和增殖通路密切相关，可能部分反映一般不良预后生物学而非放疗特异性机制。该特征与免疫细胞浸润和ESTIMATE评分的关联提示肿瘤内在基因程序与微环境存在相互作用。然而，研究存在若干局限性：数据来自回顾性队列，缺乏放疗剂量和分割信息；未充分校正年龄、分期、MYCN状态等临床因素；风险阈值基于队列中位数，缺乏固定截断值；ceRNA网络和免疫分析基于计算预测，尚需实验验证。
结论：本研究开发并外部验证了放疗相关五基因预后特征，包括CDCA8、CDCA2、FAM83D、ARHGAP11A和GAS2L3。该特征可对患者生存进行分层，并与细胞周期、DNA复制修复、免疫微环境及预测的ceRNA特征相关。沉默这些基因增强了放射诱导的生长抑制并降低了克隆存活，其中CDCA8或CDCA2沉默还增加了γ-H2AX残留灶，提示放射增敏效应可能与持续DNA损伤有关。在确立其独立预后价值或临床转化应用之前，仍需进一步的临床和机制验证。
