KRAS与ARID1A共突变与子宫内膜癌MSI/LINE-1高甲基化表型的相关性研究

《Neoplasia》：Co-mutation of KRAS and ARID1A correlates with the MSI/LINE-1 hypermethylation profile in endometrial carcinoma

　　

子宫内膜癌（EC）是全球最常见的妇科恶性肿瘤，其中子宫内膜样亚型（EEC）占多数病例，总体预后良好。尽管如此，相当一部分EEC表现出侵袭性行为和疾病复发，反映其分子异质性。这一临床现实推动了从传统组织病理学分类向由癌症基因组图谱（TCGA）定义、并经Leon-Castillo等细化的整合分子分类的转变，该分类将EC分为四个预后亚组：POLE超突变型（POLEmut）、DNA错配修复缺陷（MMRd）或微卫星不稳定（MSI）超突变型、p53异常型（p53abn）或拷贝数高型（CN-H），以及无特定分子谱型（NSMP）或拷贝数低型（CN-L）。虽然TCGA分类提供了广泛的风险分层，但NSMP组在临床上仍具有异质性。近期证据表明，空间肿瘤-免疫表型可预测该亚型内较短的无复发生存期，凸显了识别稳健生物标志物以进一步细分MMRd/MSI和NSMP/CN-L亚组的必要性。聚焦这些亚型，KRAS突变发生于约15-30%的EEC中，界定了一个主要的分子亚群，驱动增殖信号，并与子宫内膜样组织学显著相关。近期，针对KRASG12C的共价抑制剂开发显示出选择性抑制致癌RAS信号的潜力，从而提出了靶向、精准治疗干预的可能。关键在于，EEC的发病机制还涉及显著的表观遗传改变。约20-40%的EC属于MMRd亚型，在散发病例中最常见的原因是MLH1启动子高甲基化。同时，全球DNA甲基化通常通过长散在核元件-1（LINE-1）序列甲基化评估，可作为基因组稳定性的替代生物标志物。在机制上，实验模型显示抑制TET1依赖的DNA去甲基化、驱动抑癌基因高甲基化，对KRAS介导的转化是必需的。然而，KRAS与这些关键表观遗传机制（MLH1和整体甲基化状态）之间的具体相互作用在EEC中仍不明确。越来越多关于EEC中KRAS与染色质重塑因子ARID1A共突变的证据提示RAS/MAPK信号与表观遗传失调之间存在复杂的相互作用。此外，ARID1A缺陷已被直接与MMR受损相关联，促进突变性并增强免疫检查点阻断所释放的治疗性抗肿瘤免疫。将KRAS激活与MSI表型联系起来的精确机制——不同于结直肠癌（CRC）中所见的互斥性——尚未完全阐明。在此背景下，研究人员旨在整合分子与表观遗传谱分析以刻画EEC的KRAS突变亚组。具体而言，研究人员评估了KRAS突变、MSI和LINE-1甲基化的频率与共现情况，并研究其临床病理学关联。通过结合大型公共数据集的计算分析（in silico analysis）与连续单中心EEC队列的分子谱分析，研究人员旨在识别甲基化的潜在分子驱动因素，并界定一个具有预后和治疗意义的独特KRAS相关MMRd亚组，以供进一步的生物学分类。

研究背景与问题

子宫内膜癌（EC）是全球最常见的妇科恶性肿瘤，子宫内膜样癌（EEC）占多数且总体预后良好，但相当一部分EEC表现出侵袭性行为和复发，反映其分子异质性。由癌症基因组图谱（TCGA）定义、Leon-Castillo等细化的分子分类将EC分为四个预后亚组：POLE超突变型（POLEmut）、DNA错配修复缺陷（MMRd）/微卫星不稳定（MSI）超突变型、p53异常型（p53abn）/拷贝数高型（CN-H）以及无特定分子谱型（NSMP）/拷贝数低型（CN-L）。然而NSMP组在临床上仍高度异质，亟需可靠的生物标志物进一步细分MMRd/MSI与NSMP/CN-L亚组。

KRAS突变见于约15-30%的EEC，界定了一个主要分子亚群，驱动增殖信号并与子宫内膜样组织学显著相关；针对KRAS^G12C的共价抑制剂显示出选择性抑制致癌RAS信号的潜力。同时，EEC发病涉及显著表观遗传改变：约20-40%的EC属MMRd亚型，散发病例最常见原因为MLH1启动子高甲基化；整体DNA甲基化常以长散在核元件-1（LINE-1）甲基化作为基因组稳定性的替代标志物。实验模型显示，抑制TET1依赖的DNA去甲基化、驱动抑癌基因高甲基化，对KRAS介导的转化是必需的。但KRAS与MLH1甲基化及整体甲基化状态之间的具体相互作用在EEC中尚不明确。此外，KRAS与染色质重塑因子ARID1A共突变的证据提示RAS/MAPK信号与表观遗传失调存在复杂交互；ARID1A缺陷与MMR受损相关，可促进突变性并增强免疫检查点阻断所释放的抗肿瘤免疫。KRAS激活与MSI表型的联系机制（不同于结直肠癌中的互斥性）尚未阐明。因此研究人员整合分子与表观遗传谱分析，刻画EEC的KRAS突变亚组，评估KRAS突变、MSI与LINE-1甲基化的频率及共现，并探讨其临床病理学关联，以期识别甲基化的潜在分子驱动因素，界定具有预后和治疗意义的KRAS相关MMRd亚组。

研究方法

研究纳入意大利Varese的ASST-Sette Laghi/Insubria大学解剖病理科2018-2019年连续收集的83例EEC队列（属一项更大回顾性研究，伦理批件号51759），并对cBioPortal公共数据集中筛选的1300例EEC原发肿瘤进行计算机模拟分析（in silico analysis）。关键技术方法包括：免疫组化（IHC）检测四种MMR蛋白（hMLH1、hMSH2、hMSH6、hPMS2）表达；基于Archer? VariantPlex? Solid Tumour panel（覆盖70个基因，含POLE、POLD1及ARID1A）的新一代测序（NGS）；OncoMate? MSI Dx分析系统检测五个单核苷酸标志物的MSI状态；甲基化特异性多重连接依赖探针扩增（MS-MLPA）检测MLH1启动子甲基化；焦磷酸测序定量LINE-1四个CpG位点甲基化（以模型聚类确定59%为高/低甲基化界值）；统计学采用Fisher精确检验、Benjamini-Hochberg校正、Student t检验、Kaplan-Meier/log-rank及Cox回归，并在RStudio中进行多因素logistic回归。

研究结果

队列概况：cBioPortal队列含320例（24.6%）KRAS突变与980例（75.4%）KRAS野生型（WT）肿瘤；KRAS突变组低级别（G1/2）病例显著占优（p=0.029）；MSI表型在KRAS突变组占优（53.8%），MSS在KRAS-WT组占优（65.3%，p=0.003）；MLH1甲基化在KRAS突变组为56.2%，显著高于KRAS-WT组的30.8%（p=0.0001）。机构队列83例中17例（20.5%）KRAS突变，多数为低级别（84.3%）和早期（97.6%）；因样本量有限多数关联未达统计学显著，但呈现KRAS突变见于低级别、早期及高甲基化肿瘤的趋势。

KRAS突变与MLH1甲基化的关联：在cBioPortal队列中，MLH1甲基化MSI肿瘤中37/97（38.1%）为KRAS突变，而MLH1未甲基化肿瘤中突变频率为32/180（17.8%），OR为2.9（p=0.0004）。机构队列中相应为6/19（31.6%）对11/64（17.2%），OR为2.2但未达显著（p=0.2），因MLH1甲基化检测例数有限需谨慎解读。

KRAS突变与LINE-1甲基化的关联：LINE-1甲基化范围为48-70%（均值61.9±5.1%）；KRAS突变肿瘤（n=17）LINE-1甲基化水平（64.4±3.6%）高于WT肿瘤（n=66，61.3±5.3%，p=0.0079），提示KRAS激活与EEC整体高甲基化谱相关；按59%界值分组后KRAS突变组高甲基化比例更高（88.2%对72.7%），但未达显著（p=0.2）。

KRAS突变谱：cBioPortal队列328个KRAS变异以G12D（27.4%）和G12V（25%）为主，其次为G13D/C/V（18.3%）、G12A（11.3%）、G12C（6.4%）等；机构队列17个变异中G12D、G12V、G12C各4例（23.5%），G13D/C 3例（17.6%）。

KRAS与其他基因突变的共现：cBioPortal队列常见突变为PTEN（80.5%）、ARID1A（59.3%）、PIK3CA（51.2%）等；KRAS突变病例中ARID1A突变强烈共现（p=0.001），而CTNNB1、MLH1、MSH2、MSH6、POLE突变与KRAS显著互斥；经Benjamini-Hochberg校正后ARID1A、CTNNB1、MLH1、MSH2和POLE的关联仍显著。机构队列中ARID1A在76%的KRAS突变肿瘤中共现（p=0.168），CTNNB1显示显著互斥（p=0.0336），校正后均不显著。

KRAS/ARID1A突变作为表观遗传调节因子：在cBioPortal队列四个互斥组中，双突变组27/52（52%）为MSI，而双WT组90/122（74%）为MSS（p=0.0006）；双突变组27/48（56%）为MLH1甲基化，双WT组86/116（74%）为MLH1未甲基化（p=0.0005）。机构队列呈现相似趋势（双突变5/13为MSI；4例受检突变肿瘤全部MLH1甲基化）。LINE-1甲基化在双突变组（n=13，中位数64.8）显著高于双WT组（n=29，中位数61.8，p=0.0474）；仅KRAS突变组（n=4）中位数最高为65.5，与仅ARID1A突变组（n=37，中位数63.1）差异显著（p=0.0235）。

多因素分析：在283例数据完整的cBioPortal病例中，KRAS与ARID1A突变关联在未校正模型中显著（OR=2.32，p=0.0029），校正MSI状态（OR=2.15，p=0.0079）、MLH1甲基化（OR=2.13，p=0.009）、POLE状态及年龄、分级、分期后仍显著；在262例有分子分型数据者中校正分子亚型后仍显著（OR=2.06，p=0.0201）。

生存分析：cBioPortal队列中KRAS突变与总生存（OS）无关，但显示无病生存（DFS）延长的趋势（log-rank p=0.0594；Cox HR=0.51，p=0.0461）；ARID1A突变者OS显著改善（log-rank p=0.0188，HR=0.5）且DFS呈边缘延长趋势（p=0.0651）。机构队列事件数少，KRAS突变组DFS呈较好趋势（p=0.1771），未行Cox回归。

讨论与结论

研究人员认为，KRAS突变主要见于低级别、早期EEC，与ARID1A突变、MLH1启动子甲基化共关联，并以LINE-1高甲基化为特征，可能预后较好；多因素分析表明KRAS-ARID1A共现独立于既有分子亚型、MSI状态、MLH1甲基化及临床病理变量，提示其具有超出标准分类框架的独特生物学作用。KRAS与TP53或CTNNB1的互斥性符合这些驱动事件激活替代分子路径的观点；ARID1A与KRAS共突变与SWI/SNF复合物破坏和RAS/MAPK通路激活在子宫内膜肿瘤发生中的协同一致。研究人员提出假说：KRAS激活可能促成EEC的高甲基化表型，类似CRC中BRAF突变与CpG岛甲基化表型（CIMP）的关联；但研究为观察性，应解释为潜在生物学关联而非确证因果关系，需功能验证。KRAS^G12D与KRAS^G12V占优（通常导致强烈持久的MAPK通路激活）为持续致癌信号驱动慢性表观遗传改变提供合理性，也使KRAS成为前瞻性治疗靶点。KRAS相关高甲基化表型不仅限于局部基因沉默，还扩展至整体DNA甲基化，这与癌症进展中KRAS信号导致整体低甲基化的经典模型相反。临床上KRAS突变EEC呈较好预后趋势，与CRC、胰腺癌和肺癌中的侵袭性作用形成对比，反映组织限制性趋向性和谱系特异性依赖。研究优势在于整合设计与独立机构队列的全面谱分析；局限包括机构队列规模有限、MLH1甲基化检测例数少、缺乏功能验证、单区块测序无法评估瘤内异质性、早期病例为主限制生存事件数。从转化角度看，ARID1A突变与KRAS频繁共现提示潜在治疗脆弱性，ARID1A缺陷可诱导甲羟戊酸通路激活，形成可用他汀类药物靶向的代谢依赖。研究结论：KRAS突变EEC可能代表一个分子和表观遗传上独特的亚群，以ARID1A共突变、MSI/MLH1甲基化及LINE-1高甲基化为特征；这些发现对细分异质性NSMP组的潜在价值值得在更大独立队列中进一步评估，以确定其临床和治疗适用性。