原发部位差异对肺外低分化神经内分泌癌转移模式及生存关联的影响

《Frontiers in Endocrinology》：Primary-site variation in recorded metastatic patterns and survival associations in extrapulmonary poorly differentiated neuroendocrine carcinoma

　　

背景：肺外低分化神经内分泌癌（extrapulmonary poorly differentiated neuroendocrine carcinoma，EP-PD-NEC）具有侵袭性和异质性。本研究的主要目标是评估不同原发部位在诊断时记录的骨、脑、肝和肺转移方面的差异。次要分析检验了预后效应的修饰作用以及这些变量中的增量预测信息。 方法：研究人员分析了2010–2023年监测、流行病学和最终结果（Surveillance, Epidemiology, and End Results，SEER）数据。严格定义的EP-PD-NEC采用ICD-O-3编码8013和8041–8045。利用带有个体水平异方差一致性标准误的修正泊松回归评估诊断时记录的骨、脑、肝和肺转移。采用含原发部位与器官交互项的多变量Cox模型分析总生存期（overall survival，OS）。2010–2017年识别的关联在2018–2022年重新评估。开发了四个预设的岭惩罚Cox模型，并在具有完整四器官状态的远处分期病例中进行时间评估。主要患病率分析纳入所分析器官状态已知的患者；敏感性分析将未知状态重新编码为阴性或采用多重插补。 结果：在6044例符合条件的患者中，5820例（96.3%）拥有完整的四器官信息，2356例至少有一个记录到的器官受累。每个器官的记录转移患病率在六个原发部位组间均存在差异（所有经错误发现率校正后的P<0.001）。与泌尿生殖系原发相比，胰腺原发的肝脏调整后患病率比（prevalence ratio，PR）为3.80（95%置信区间[confidence interval，CI] 3.39–4.25），脑转移PR为2.86（95% CI 1.75–4.66）。记录的四器官计数增加与更短的OS相关。全局原发部位与器官交互作用显著（似然比χ2=43.49，自由度20，P = 0.00176），但完整的交互矩阵在时间预测中增加的预测价值极小（ΔC指数=0.00169）。在1374例后期患者中，谱系模型的C指数为0.607（95% CI 0.588–0.624），12个月曲线下面积（area under the curve，AUC）为0.672（95% CI 0.639–0.706），Brier评分为0.196（95% CI 0.185–0.207）。与四器官计数模型相比，配对bootstrap差异对于C指数为0.0085（95% CI 0.0035–0.0139），对于12个月AUC为0.0119（95% CI 0.0033–0.0231），对于Brier评分为?0.0020（95% CI ?0.0038至?0.0004）。缺失状态分析产生了相似的患病率估计值。 结论：在EP-PD-NEC中，记录到的转移模式因原发部位而异。器官身份在记录的四器官计数之外提供的预测信息极少，且原发部位与器官的交互作用带来的改善可忽略不计。总体区分能力仍然有限。

论文解读：肺外低分化神经内分泌癌的原发部位差异及其对转移模式和生存预后的影响

一、研究背景与研究目的

神经内分泌肿瘤是一类生物学行为迥异的肿瘤，其形态学、分子改变、临床特征和治疗反应差异显著。当代分类将其分为高分化神经内分泌肿瘤和低分化神经内分泌癌（poorly differentiated neuroendocrine carcinoma，NEC），后者包括小细胞和大细胞表型，通常进展迅速且预后极差。肺外低分化神经内分泌癌（extrapulmonary poorly differentiated neuroendocrine carcinoma，EP-PD-NEC）较为罕见，几乎可发生于任何器官，但由于其发病率低，前瞻性证据极为匮乏，临床管理多借鉴于小细胞肺癌的经验。既往基于人群的研究已发现不同肺外原发部位的肿瘤在分期分布和生存方面存在差异，这引发了一个重要问题：不同原发部位的转移模式及其与生存的关联究竟有何不同？此外，转移器官的身份是否提供了超越受累器官数量的额外预后信息，目前仍不明确。因此，本研究旨在系统评估EP-PD-NEC中不同原发部位在诊断时的骨、脑、肝和肺转移模式的差异，并进一步分析这些转移的生存关联是否存在原发部位特异性效应修饰，以及转移信息能否提升预后预测能力。

二、研究方法概述

研究人员开展了一项回顾性基于人群的队列研究，数据来源于美国国家癌症研究所的监测、流行病学和最终结果（Surveillance, Epidemiology, and End Results，SEER）数据库，覆盖2010年至2023年的病例。严格定义的EP-PD-NEC采用国际肿瘤疾病分类第三版（ICD-O-3）形态学编码8013和8041–8045。原发部位被预先划分为六个类别：泌尿生殖系、妇科、结直肠、胰腺、上消化道及其他肺外部位。研究首先利用修正泊松回归评估各器官转移的调整后患病率比（prevalence ratio，PR），随后采用含原发部位与器官交互项的Cox比例风险模型检验预后效应修饰。为进一步量化转移信息的预测价值，研究人员开发了四个预设的岭惩罚Cox模型：基础模型（M0）、加入原发部位的模型（M1）、加入记录的四器官计数的模型（M2）以及加入各器官单独指示变量的谱系模型（M3），并在非重叠的时间段（2010–2017年作为训练集，2018–2022年作为验证集）内进行时间外部验证。

三、研究结果

3.1 队列构成与基线特征

从121,264条肿瘤记录中筛选出6044例符合严格标准的EP-PD-NEC患者。其中，79.8%为小细胞形态，20.2%为大细胞形态。中位年龄67岁，男性占57.0%。原发部位分布为泌尿生殖系35.3%、妇科15.2%、结直肠15.0%、胰腺6.8%、上消化道9.6%、其他部位18.1%。52.0%的患者在诊断时为SEER远处分期。5820例（96.3%）患者拥有完整的骨、脑、肝和肺四器官转移状态信息。

3.2 原发部位相关的记录转移模式差异

在所有六个原发部位组中，每个记录器官的转移患病率均存在显著差异（经错误发现率校正后P<0.001）。最显著的对比体现在肝脏转移：胰腺原发患者肝脏转移记录率为71.1%，结直肠为52.9%，上消化道为43.9%，而泌尿生殖系仅为19.3%。经多变量调整后，与泌尿生殖系原发相比，胰腺原发的肝脏转移PR高达3.80（95%置信区间[confidence interval，CI] 3.39–4.25），结直肠为2.88（95% CI 2.56–3.23），上消化道为2.31（95% CI 2.03–2.63）。脑转移虽总体少见，但胰腺原发记录率最高（6.5%），调整后PR为2.86（95% CI 1.75–4.66）。相反，骨转移在泌尿生殖系原发中最常见（22.1%），其他部位组的调整后骨转移PR约为泌尿生殖系的0.52至0.65倍。肺部转移的差异相对较小。缺失状态敏感性分析未显著改变主要对比结果。

3.3 记录的四器官计数与转移模式

在2356例至少有一个记录器官阳性的核心转移谱系队列中，记录器官数量为1、2、3和4个的患者分别有1543、611、185和17例。随着受累器官数量增加，中位总生存期（overall survival，OS）从单个器官的6.70个月逐步缩短至2个器官的5.32个月、3个器官的4.63个月和4个器官的4.40个月。最常见的单一器官转移模式为单纯肝脏转移（1038例，中位OS 5.82个月），其次为单纯骨转移（302例，中位OS 9.99个月）。

3.4 原发部位特异性生存关联

包含所有原发部位与器官交互项的Cox模型相较于主效应模型拟合显著改善（似然比χ2=43.49，自由度20，P = 0.00176）。具体而言，肝脏转移在泌尿生殖系原发患者中的风险比（hazard ratio，HR）为1.65（95% CI 1.40–1.95），在结直肠原发中为2.35（95% CI 1.59–3.47），在胰腺原发中为3.02（95% CI 1.64–5.58），而在上消化道原发中仅为1.03。聚焦于肝脏与原发部位交互作用的简化模型显示，该交互项具有统计学显著性（似然比χ2=18.84，自由度5，P = 0.00205）。然而，当将完整的交互矩阵加入模型后，其在验证队列中的预测性能提升微乎其微，C指数仅增加0.00169，12个月受试者工作特征曲线下面积（area under the curve，AUC）甚至下降0.0131。

3.5 跨时间段关联的一致性

在2010–2017年期间识别出的七个满足错误发现率校正标准的调整后部位-器官对比，在2018–2022年期间均保持了相同的方向，且所有后期对比的95% CI均排除1。例如，胰腺与泌尿生殖系相比的肝脏转移PR在两个时期分别为3.77和3.97，表明这些转移模式差异随时间保持稳定。

3.6 简约模型的开发与时间评估

在四个预设模型中，谱系模型（M3）在验证队列中表现最佳，其C指数为0.607（95% CI 0.588–0.624），12个月AUC为0.672（95% CI 0.639–0.706），12个月Brier评分为0.196（95% CI 0.185–0.207）。与四器官计数模型（M2）相比，M3在C指数上的配对bootstrap差异为0.0085（95% CI 0.0035–0.0139），在12个月AUC上的差异为0.0119（95% CI 0.0033–0.0231），表明器官身份在计数基础上提供了微小但有统计学意义的增量预测信息。然而，整体区分能力仍然有限。

3.7 风险分层与亚组表现

应用训练集中预测的12个月死亡率三分位数切点（0.641和0.775），可将验证队列分为低、中、高风险三组，其中位OS分别为10.9、6.8和4.3个月，12个月OS分别为47.2%、28.2%和17.7%，组间差异显著。亚组分析显示，模型在上消化道原发患者中的区分能力最差（12个月AUC为0.532），而在其他部位组中最好（0.725）。

四、讨论与结论

本研究揭示了EP-PD-NEC中不同原发部位在诊断时记录转移模式的显著差异，尤其是肝脏转移在胃肠道和胰腺原发中的高发生率，这可能部分反映了门静脉引流的作用，同时也可能与不同NEC亚型的分子特征有关。然而，尽管存在这种转移异质性，原发部位与器官之间的交互作用对生存预测的贡献却非常有限，记录的四器官计数已捕获了转移变量的大部分预测增益。这一发现提示，虽然不同原发部位的转移分布存在差异，但这种差异并未转化为临床上显著的预后预测能力提升。值得注意的是，模型整体区分能力仍然偏低，特别是在上消化道原发等亚组中，这限制了其用于个体化风险估计的精确性。研究的优势在于采用了基于人群的大样本队列、严格的形态学定义以及跨时间段的验证设计，但也存在若干局限性，包括缺乏中心病理复核、无Ki-67指数或分子数据、缺少影像学强度信息、无法获取后续治疗细节等。结论部分指出：在EP-PD-NEC中，原发部位与诊断时记录的转移模式显著相关；记录的四器官计数解释了转移变量的大部分预测增益；器官身份仅提供小幅改善，而原发部位与器官的交互作用对后期预测的贡献极小；整体模型区分能力仍然有限。