BRAF突变影响乳头状甲状腺癌瘤内异质性

摘要：瘤内异质性（intratumor heterogeneity, ITH）对患者临床结局具有重要影响，其在乳头状甲状腺癌（papillary thyroid carcinoma, PTC）中的预后价值及影响因素尚不明确，值得深入研究。研究人员首先评估了ITH对患者预后的影响，采用突变等位基因肿瘤异质性（Mutant Allele Tumor Heterogeneity, MATH）评分量化ITH水平。随后探索了与ITH相关的潜在因素，并预测可能参与ITH调控的信号通路。结果显示，在4种预后指标中，较高的ITH主要与较差的疾病无病间期（disease-free interval, DFI）相关（风险比HR = 2.64, p = 0.01），且ITH对DFI具有潜在的预测价值。进一步分析确定BRAF突变与甲状腺分化评分（thyroid differentiation score, TDS）是独立影响ITH的关键因素（均p < 0.05），其中TDS对ITH的预测能力尤为突出，并在预测中获得了良好的净获益。TDS对ITH始终表现出稳定的负向效应，而BRAF突变在单变量回归中与ITH呈正相关（β = 0.196）。通过敏感性分析、回归分析、亚组分析、交互效应分析及中介效应分析，研究人员发现TDS在BRAF突变对ITH的影响中发挥抑制作用。最后，研究人员揭示丙酸代谢通路与ITH的相关性最强。结论：ITH与PTC患者的DFI密切相关，应予以更多关注。TDS与BRAF突变是ITH的关键影响因素，但TDS对BRAF突变影响ITH的过程发挥抑制作用。

本研究针对乳头状甲状腺癌（papillary thyroid carcinoma, PTC）预后个体差异大、分子机制不清的临床难题展开。PTC占所有甲状腺癌的80%以上，尽管总体5年生存率超过95%，但部分患者仍出现局部侵犯、淋巴结转移及远处转移，导致不良临床结局。瘤内异质性（intratumor heterogeneity, ITH）指同一肿瘤内癌细胞在基因组、转录组和表观遗传水平的差异，是驱动肿瘤进展、耐药和复发的关键因素。然而，ITH在PTC中的预后作用及其调控机制尚未明确。该研究基于癌症基因组图谱（The Cancer Genome Atlas, TCGA）数据库，旨在阐明ITH对PTC患者预后的影响，识别其关键调控因素，并探索潜在分子通路，为PTC的精准治疗提供新方向。研究最终证实ITH是PTC患者疾病无病间期（disease-free interval, DFI）的独立危险因素，BRAF突变与甲状腺分化评分（thyroid differentiation score, TDS）共同调控ITH，且TDS对BRAF突变的促ITH效应具有抑制作用，丙酸代谢通路异常是ITH的核心代谢特征。该研究成果发表于《Analytical Cellular Pathology》。

研究人员采用的主要关键技术方法如下：研究数据来源于TCGA-THCA队列，纳入502例具备体细胞突变注释文件的PTC患者。采用MaFTools包计算突变等位基因肿瘤异质性（Mutant Allele Tumor Heterogeneity, MATH）评分以量化ITH。利用最大选择秩统计确定ITH最佳截断值，通过Kaplan-Meier生存分析和Cox回归模型评估ITH的预后价值，采用时间依赖受试者工作特征（time-dependent receiver operating characteristic, ROC）曲线和决策曲线分析（decision curve analysis, DCA）验证其预测效能。通过单因素和多因素线性回归筛选ITH的独立影响因素，结合敏感性分析、亚组分析、交互效应分析及中介效应分析，探究BRAF突变、TDS与ITH之间的复杂关联。采用单样本基因集富集分析（single-sample Gene Set Enrichment Analysis, ssGSEA）计算京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路活性，识别与ITH显著相关的信号通路。

研究结果部分，首先探讨ITH与PTC患者预后的关系。研究人员通过生存分析发现，高ITH组患者的DFI和无进展生存期（progression-free survival, PFS）显著低于低ITH组（p < 0.05），而ITH对总生存期（overall survival, OS）和疾病特异性生存期（disease-specific survival, DSS）无显著影响。时间依赖ROC分析显示，ITH预测1年、3年和5年DFI的曲线下面积（area under the curve, AUC）分别为0.72、0.61和0.59，表明其对短期DFI具有中等预测效能。DCA分析进一步证实，当阈值概率在0.10至0.16范围内时，ITH对全期DFI具有临床净获益。

其次，筛选与ITH相关的临床特征。通过比较高低ITH组的临床病理特征，研究人员发现肿瘤大小、T分期、腺外侵犯、淋巴结转移情况、TDS、BRAF突变状态、组织学类型、调节性T细胞（regulatory T cells, Tregs）浸润水平、术后放疗情况及个人病史等存在差异（p < 0.05）。排除共线性变量后，多因素线性回归分析确定TDS和BRAF突变为ITH的独立影响因素。ROC和DCA分析显示，TDS预测ITH的AUC达0.86，优于BRAF突变的0.62，且在更广的阈值范围内具有临床净获益。

第三，解析关键临床特征对ITH的作用模式。研究人员发现，TDS在单因素和多因素分析中均与ITH呈稳定负相关，而BRAF突变在单因素分析中呈正相关，在多因素分析中转为负相关。敏感性分析显示，仅在调整TDS后BRAF突变与ITH的负相关关系出现，提示TDS是导致这种关联转变的核心因素。亚组分析将患者按BRAF突变状态和TDS水平分为四组，结果显示无论BRAF突变状态如何，低TDS均显著增加ITH；在低TDS状态下，无BRAF突变组的ITH反而高于BRAF突变组。逻辑回归和线性回归进一步验证了上述发现。交互效应分析未发现TDS与BRAF突变对ITH存在交互作用（交互项p > 0.05）。中介效应分析表明，BRAF突变显著降低TDS（效应值-1.425, p < 0.001），而TDS降低显著升高ITH（效应值-0.188, p < 0.001），TDS在BRAF突变影响ITH的过程中发挥显著的抑制性中介效应。

第四，探索ITH预后影响的潜在机制。研究人员通过ssGSEA分析比较高低ITH组的186条KEGG通路活性，发现126条通路存在差异。相关性分析和多因素线性回归最终确定丙酸代谢通路是唯一独立于TDS和BRAF突变、与ITH显著相关的通路（β = -2.123, p < 0.001）。该通路在高ITH组中活性下调，其预测ITH的AUC为0.77，且在阈值概率0.02至0.61范围内具有临床净获益。

讨论部分指出，BRAF V600E突变通过持续激活MAPK通路促进细胞增殖、诱导DNA损伤和基因组不稳定、重塑肿瘤微环境，从而增加ITH。TDS反映肿瘤细胞保留正常甲状腺滤泡细胞特征的程度，低分化状态伴随基因组不稳定性增加、放射性碘抵抗及代谢重编程，导致ITH升高。TDS通过维持基因组稳定性、限制代谢可塑性及调控BRAF下游效应，掩盖了BRAF突变对ITH的促进作用。丙酸代谢通路异常通过组蛋白丙酰化修饰、线粒体功能障碍及活性氧（reactive oxygen species, ROS）积累，诱导复制应激和DNA损伤，驱动基因突变和染色体不稳定，从而促进ITH形成。

研究结论部分总结：ITH可作为PTC患者的重要预后标志物。BRAF突变和TDS与ITH显著相关，但TDS对BRAF突变影响ITH的过程具有抑制作用。此外，丙酸代谢通路功能异常与ITH相关，提示该通路为肿瘤异质性提供能量支持和代谢适应性。这些发现为PTC的个性化治疗和精准医疗提供了新的研究方向。