《Frontiers in Genetics》：Risk factors and time-to-event prediction of brain metastasis in patients with lung adenocarcinoma: a retrospective cohort study
编辑推荐：
背景：脑转移（brain metastasis, BM）是肺腺癌（lung adenocarcinoma, LUAD）的一种临床重要并发症，然而其预测因随访差异、删失、治疗效果及竞争性死亡而变得复杂。本研究评估了基线无BM的LUAD患者的时间依赖性BM风险，并
背景：脑转移（brain metastasis, BM）是肺腺癌（lung adenocarcinoma, LUAD）的一种临床重要并发症，然而其预测因随访差异、删失、治疗效果及竞争性死亡而变得复杂。本研究评估了基线无BM的LUAD患者的时间依赖性BM风险，并开发了一个探索性预测模型。方法：研究人员对2018年至2024年间接受治疗的连续经病理确诊的LUAD患者进行了回顾性筛选。患者于基线阴性增强脑磁共振成像（contrast-enhanced brain MRI）之日进入队列。使用Cox比例风险回归分析新发BM。候选预测因子被预先指定，评估了共线性，并使用LASSO惩罚Cox回归进行预测因子缩减。内部验证使用了1000次Bootstrap重抽样。采用Fine-Gray分析评估竞争性死亡，并将表皮生长因子受体-酪氨酸激酶抑制剂（EGFR-TKI）暴露作为时间依赖性协变量进行检验。结果：在94名筛选患者中，66名被纳入；在中位随访27.8个月期间，26名患者发生了新发BM。LASSO保留了T分期、N分期、EGFR突变状态和血清癌胚抗原（CEA）。在简约Cox模型中，T3-4期疾病、N1-3期受累和EGFR突变与更高的BM风险相关，而CEA的估计值精确度较低。表观C指数为0.781，乐观校正后C指数为0.742。时间依赖性受试者工作特征曲线下面积（AUC）在12个月时为0.773，在24个月时为0.758。竞争风险和经治疗调整的分析显示方向一致的结果。结论：一个基于Cox的模型整合了T分期、N分期、EGFR状态和CEA，在估计12个月和24个月的无BM概率方面显示出中等程度的内部验证性能。在临床应用前需进行外部验证。
研究背景与目的
肺癌是全球癌症相关死亡的主要原因之一，其中非小细胞肺癌（non-small cell lung cancer, NSCLC）约占80%-85%，而肺腺癌（lung adenocarcinoma, LUAD）是其最常见的组织学亚型。尽管分子靶向治疗和免疫治疗显著改变了LUAD的治疗格局，但远处转移进展仍是发病率和死亡率的主要决定因素。脑转移（brain metastasis, BM）是NSCLC中尤为重要的转移部位，其发生与神经功能损伤、生活质量下降、治疗复杂性增加以及不良生存结局相关。因此，在基线脑部评估阴性后，识别后续BM高风险患者对于风险分层和纵向监测具有临床意义。
目前已有一些针对NSCLC和LUAD的BM预测模型，但存在重要的方法论局限性。部分研究将BM视为二元结局而未考虑异质性随访时长，另一些则聚焦于高度选择的人群。此外，在小样本回顾性队列中，若使用广泛的单变量筛选和逐步变量选择而未进行充分的内部验证，预测性能易受过度拟合影响。由于BM本质上是时间-事件结局，且BM前的死亡是一个潜在的竞争事件，因此需要更严谨的方法。同时，不同人群的分子流行病学存在差异，例如亚洲人群中EGFR突变频率远高于欧洲人群，这提示在不同临床或分子背景下开发的模型可能无法直接适用。
因此，本研究旨在评估一个中国LUAD队列中基线无BM证据的患者发生新发BM的时间依赖性风险，并开发一个探索性的基于Cox的比例风险预测模型，用于估计12个月和24个月的无脑转移生存概率。
关键技术方法
这是一项回顾性队列研究，数据来源于聊城市第二人民医院胸肿瘤科。研究人员从2018年1月至2024年12月期间连续筛查经病理确诊的LUAD患者，最终纳入66名符合条件的患者。主要统计方法包括：使用Cox比例风险回归分析新发BM；采用LASSO惩罚Cox回归进行预测因子缩减以避免过度拟合；使用1,000次Bootstrap重抽样进行内部验证以校正乐观偏差；通过Harrell's C指数和时间依赖性受试者工作特征曲线下面积（AUC）评估区分度，通过校准图和Brier评分评估校准度；使用Fine-Gray竞争风险模型处理BM前的死亡事件；并将EGFR-TKI暴露作为时间依赖性协变量进行敏感性分析。
研究结果
3.1 队列构建、随访与脑转移事件
在94名筛查患者中，28名因各种原因被排除，最终66名患者构成分析队列。中位随访27.8个月期间，26名患者发生新发BM，中位发生时间为15.7个月。9名患者在BM前死亡，31名患者存活至末次随访且无BM记录。
3.2 基线临床与治疗特征
发生BM的患者相较于未发生者，更常表现为T3~4期（61.5% vs. 30.0%）、N1~3期淋巴结受累（80.8% vs. 52.5%）、EGFR突变阳性（73.1% vs. 42.5%）以及血清癌胚抗原（CEA）>5 ng/mL（65.4% vs. 37.5%）。
3.3 脑转移的时间-事件分析
单变量Cox回归显示，T3~4期、N1~3期、EGFR突变和CEA>5 ng/mL均与更高的BM风险相关。LASSO惩罚Cox回归保留了T分期、N分期、EGFR突变状态和血清CEA这四个变量。在相应的简约Cox模型中，T3~4期（HR 2.31）、N1~3期（HR 2.42）和EGFR突变（HR 2.49）与更高的BM风险显著相关，而CEA>5 ng/mL的关联则不够精确（HR 2.04）。Kaplan-Meier分析显示，EGFR突变型肿瘤患者的无脑转移生存期（BMFS）显著短于野生型患者（中位22.4个月 vs. 未达到，log-rank P=0.009）。
3.4 探索性时间-事件预测模型与内部验证
基于上述四个变量构建的Cox列线图，用于估计个体化的12个月和24个月无BM概率。模型的表观C指数为0.781，经Bootstrap校正后的乐观校正C指数为0.742。时间依赖性ROC分析显示，12个月和24个月的AUC分别为0.773和0.758。校正后的校准图显示预测值与观测值之间具有合理的一致性，校正后的校准斜率为0.87，表明存在适度的残余过度拟合。
3.5 竞争性死亡与敏感性分析
Fine-Gray竞争风险分析的结果与主Cox模型方向一致：T3~4期（sHR 2.18）、N1~3期（sHR 2.27）和EGFR突变（sHR 2.36）仍与更高的BM风险相关。在将EGFR-TKI暴露作为时间依赖性协变量的额外Cox敏感性分析中，各预测因子的关联方向保持一致，但效应估计值有所减弱。
3.6 按EGFR突变状态的探索性分析
在EGFR突变型患者中，BM发生率更高。携带外显子19缺失和L858R突变的患者中BM发生数量相近，但由于样本量过小，未对不同突变亚型的转移风险进行正式推断。
讨论与结论
本研究发现，晚期T分期、区域淋巴结受累、EGFR突变和血清CEA升高与LUAD患者后续发生BM的风险增高相关。这些关联与既往证据一致，具有生物学合理性。研究的一个重要方法是采用时间-事件分析框架，而非简单的二分类终点，从而有效处理了异质性随访时间和右删失数据。同时，通过竞争风险分析和时间依赖性治疗暴露分析，增强了结果的稳健性。
然而，该预测模型应被视为探索性工具，而非经过充分验证的临床决策手段。其主要价值在于提出假设和进行初步风险分层，尚需在更大的独立队列中进行外部验证。研究的局限性包括单中心小样本设计、非统一的影像监测方案、异质性的治疗方案以及缺乏外部验证。
研究结论：在该回顾性LUAD队列中，晚期T分期、淋巴结受累、EGFR突变和血清CEA升高与后续BM的时间依赖性风险增加相关。一个整合了这些常规可用因素的探索性Cox模型显示出中等水平的内部验证区分度和校准度。这些发现支持进一步研究个体化BM风险分层，但不能作为直接改变治疗方案的依据。需要进行更大规模的前瞻性多中心研究，并进行外部验证、标准化影像监测、详细的治疗暴露评估和竞争风险分析，以确定该方法是否能对临床管理产生有意义的影响。
