《Journal of Radiation Research and Applied Sciences》：Identification of ionizing-radiation-exposure-associated hub genes and a five-gene prognostic signature for pediatric medulloblastoma: An integrated bioinformatics study
编辑推荐：
背景：髓母细胞瘤（medulloblastoma）是最常见的儿童恶性脑肿瘤，放射治疗（radiotherapy）是治疗的重要组成部分。与辐射暴露反应及预后相关的分子特征尚未完全阐明。方法：研究人员分析了来自GSE197236的配对照射与未照射人淋巴母细胞样本（
背景：髓母细胞瘤（medulloblastoma）是最常见的儿童恶性脑肿瘤，放射治疗（radiotherapy）是治疗的重要组成部分。与辐射暴露反应及预后相关的分子特征尚未完全阐明。方法：研究人员分析了来自GSE197236的配对照射与未照射人淋巴母细胞样本（6对），以鉴定与电离辐射暴露相关的差异表达基因（differentially expressed genes, DEGs），随后在来自GSE85217的103例小儿髓母细胞瘤中评估其预后关联。采用Cox回归和LASSO回归构建内部评估的预后标志物，并利用蛋白质-蛋白质相互作用（protein-protein interaction, PPI）分析进行基于网络的探索性候选基因优先级排序。结果：共鉴定出947个差异表达基因，主要与p53信号通路、细胞周期调控、DNA损伤修复及凋亡通路相关。由CCNB1、CDKN1A、RRM2、TOP2A和BIRC5组成的五基因预后模型可区分总生存期（overall survival, OS）不同的患者，在1年、3年和5年的内部估计时间依赖性AUC分别为0.834、0.812和0.786。Bootstrap 95%置信区间分别为0.764-0.901、0.730-0.879和0.681-0.862。表观C-index为0.781，经乐观校正的C-index为0.748。这些数值均来自同一开发队列，并非独立验证估计值。在探索性分析中，六个网络优先排序的枢纽基因候选与分子亚型或生存相关。结论：本研究鉴定了与电离辐射暴露相关的差异表达基因，并开发了一个五基因预后标志物，该标志物可稳健地将小儿髓母细胞瘤患者按总生存期进行风险分层，且具有较好的内部预测性能。该标志物连同六个网络优先排序的枢纽基因候选，可能为未来的风险适应性治疗策略和治疗靶点发现提供参考。
**论文解读：电离辐射暴露相关基因与小儿髓母细胞瘤预后标志物的整合生物信息学研究**
**一、研究背景与问题**
髓母细胞瘤（medulloblastoma）是儿童中枢神经系统中最常见的恶性肿瘤，约占所有儿童中枢神经系统肿瘤的20%，是儿童癌症相关死亡的重要原因。当前标准治疗包括手术切除、风险分层的颅脊髓照射（craniospinal irradiation, CSI）联合辅助化疗，5年总生存率可达60%–80%。然而，幸存者常因放疗暴露出现神经认知障碍、内分泌功能障碍、感音神经性耳聋及继发恶性肿瘤等严重远期并发症。因此，迫切需要鉴定与辐射暴露反应及预后相关的分子特征，以指导风险适应性治疗策略的优化。电离辐射的主要生物学效应是诱导DNA损伤，尤其是双链断裂，进而激活细胞周期检查点、DNA损伤修复、程序性细胞死亡和细胞衰老等一系列防御机制。p53肿瘤抑制因子在这一过程中发挥核心调控作用。既往转录组研究显示，电离辐射可在人类细胞中引起数百至数千个基因的表达改变。然而，这些辐射相关转录本在儿童髓母细胞瘤中的预后价值尚未充分阐明。现有髓母细胞瘤分子分型（WNT、SHH、Group 3和Group 4）虽已用于预后判断，但辐射暴露相关基因能否提供额外预后信息仍不清楚。此外，既往研究多采用传统生存分析，缺乏对模型内部验证和过拟合风险的控制。因此，本研究旨在通过整合公共转录组数据，系统鉴定电离辐射暴露相关差异表达基因（differentially expressed genes, DEGs），构建并内部评估多基因预后模型，同时利用蛋白质-蛋白质相互作用（protein-protein interaction, PPI）网络优先排序枢纽基因候选。
**二、研究设计与技术方法**
研究人员采用两阶段生物信息学策略。第一阶段使用GEO数据库中的GSE197236数据集（包含6对配对的人淋巴母细胞样本，接受2 Gy X射线照射与未照射对照），通过limma配对差异表达分析鉴定电离辐射暴露相关DEGs，筛选标准为|log2 fold change|>1.0且Benjamini-Hochberg校正后P<0.05。第二阶段在GSE85217数据集（包含103例具有临床预后数据的小儿髓母细胞瘤患者，其中89例具有分子亚型信息）中提取这些DEGs的表达谱，采用单因素Cox回归筛选与总生存期（overall survival, OS）相关的候选基因，再通过LASSO惩罚Cox回归构建预后风险模型。模型内部评估采用1000次嵌套bootstrap重抽样，计算时间依赖性AUC和C-index。功能富集分析采用GO和KEGG通路注释，并使用基于秩的基因集富集分析（GSEA）。PPI网络使用STRING数据库构建，MCODE模块分析和CytoHubba五种中心性算法（Degree、Betweenness、Closeness、Stress、Bottleneck）用于枢纽基因优先排序。所有计算均在R 4.5.2中完成。
**三、研究结果**
**1. 电离辐射暴露相关DEGs的鉴定**
对GSE197236中6对配对样本的分析共鉴定出947个DEGs，其中528个上调、419个下调。已知的辐射响应转录本如CDKN1A（log2FC=3.42）、FDXR（log2FC=2.98）、GADD45A（log2FC=2.87）和BAX（log2FC=2.15）均显著富集，验证了暴露相关转录程序的生物学合理性。
**2. 功能富集与通路分析**
GO注释显示DEGs显著富集于DNA损伤反应（FDR=2.34×10
-12）、细胞周期检查点调控（FDR=5.67×10
-10）、DNA修复机制（FDR=1.23×10
-9）和内在凋亡信号级联（FDR=3.45×10
-8）。KEGG通路分析发现18条显著富集通路，包括p53信号通路（FDR=1.45×10
-10）、细胞周期（FDR=3.21×10
-8）和细胞衰老（FDR=6.78×10
-7）。基于完整排序基因列表的GSEA进一步证实了p53信号（NES=2.31）、G2/M检查点（NES=2.18）和DNA修复（NES=1.95）的富集。
**3. 髓母细胞瘤中的预后景观**
在GSE85217队列中（103例，中位随访58个月，38例死亡），单因素Cox回归筛选出23个与OS显著相关的暴露相关基因（名义P<0.01）。分子亚型分布为WNT 14例、SHH 28例、Group 3 23例、Group 4 24例，另有14例亚型信息缺失。多因素分析显示，风险评分、转移状态、Group 3亚型（与WNT相比）和Group 4亚型（与WNT相比）均与OS独立相关。
**4. 预后风险模型的构建与内部评估**
LASSO回归选出CCNB1、CDKN1A、RRM2、TOP2A和BIRC5五个基因构建风险评分。以队列中位数将患者分为高风险组（n=52）和低风险组（n=51），高风险组OS显著较差（log-rank P<0.0001）。内部时间依赖性AUC在1、3、5年分别为0.834、0.812和0.786，bootstrap 95%置信区间分别为0.764-0.901、0.730-0.879和0.681-0.862。表观C-index为0.781，经乐观校正后为0.748。1、3、5年校准显示预测与观察生存率接近（0.86 vs 0.84，0.71 vs 0.73，0.61 vs 0.63）。临床/亚型模型C-index为0.721，加入风险评分后升至0.781，提示该评分具有增量价值。
**5. 枢纽基因鉴定与临床关联**
基于23个预后基因构建的PPI网络含23个节点和87条边（置信度阈值0.4）。MCODE识别出两个功能模块，中心性分析优先排序出六个候选枢纽基因：CCNB1、CDKN1A、TOP2A、BIRC5、CHEK1和MCM2。单因素Cox分析显示这些基因的每SD风险比分别为1.42、0.69、1.36、1.31、1.24和1.19。多因素分析中，前四个基因仍与OS显著相关，而CHEK1和MCM2在调整后失去显著性。STRING阈值稳健性分析在0.4和0.7置信度下保留六个候选基因，在0.9下保留五个。
**四、讨论与结论**
本研究的核心发现是鉴定出947个电离辐射暴露相关DEGs，并构建了由CCNB1、CDKN1A、RRM2、TOP2A和BIRC5组成的五基因预后模型。该模型内部评估表现良好，表明辐射暴露相关转录程序中的特定基因与儿童髓母细胞瘤患者总生存期存在关联。CCNB1、RRM2、TOP2A和BIRC5的不利预后系数可能部分反映肿瘤增殖状态，而CDKN1A的保护性关联可能反映检查点能力或分子亚型生物学。六个枢纽基因候选应视为网络优先排序结果，其中心性反映网络拓扑特征而非直接生物学因果性。BIRC5和CHEK1在其他背景下具有药理可操作性，可作为未来验证的候选靶点而非已确认的治疗靶点。
研究局限性包括：发现数据集来自淋巴母细胞而非髓母细胞组织，因此鉴定到的基因代表一般辐射暴露转录程序而非髓母细胞特异性放射敏感性特征；五基因模型在同一回顾性队列中开发和内部评估（仅38例死亡），存在过拟合风险；14例患者缺乏分子亚型信息，临床协变量不完整；标志物可能部分反映分子亚型或一般增殖生物学；缺乏直接的放疗反应、局部控制或克隆源存活终点；使用不同芯片平台和批量转录组测量可能影响可重复性。
结论：本研究提供了一个基于辐射反应来源的预后模型和一组候选枢纽基因，可能为未来的风险适应性放疗和靶向治疗设计提供参考。然而，鉴于内部评估和探索性研究性质，任何临床转化前必须进行独立的、前瞻性的验证研究。
