《Frontiers in Immunology》：Molecular subtypes associated with lymph node metastasis in colorectal cancer reveal the tumor-promoting role of GDI1 and the preclinical activity of deslanoside
编辑推荐：
背景： 淋巴结转移（lymph node metastasis, LNM）是结直肠癌分期和治疗决策的核心环节，标志着肿瘤从局部生长向区域播散的关键转变。与淋巴结转移相关的分子异质性及其相应的干预靶点需要系统性研究。
方法： 本研究基于淋巴结转移相关表达特征建立
背景： 淋巴结转移（lymph node metastasis, LNM）是结直肠癌分期和治疗决策的核心环节，标志着肿瘤从局部生长向区域播散的关键转变。与淋巴结转移相关的分子异质性及其相应的干预靶点需要系统性研究。
方法： 本研究基于淋巴结转移相关表达特征建立了C1至C4分子亚型。整合多个转录组队列、临床构成、通路特征以及TIDE、IPS和ESTIMATE评分来表征亚型差异。通过多个队列的分析筛选与进展相关的候选基因。开展细胞实验和皮下异种移植实验评估GDI1功能。利用GDI1结构进行虚拟筛选和分子动力学分析，随后在细胞和异种移植模型中对地高辛苷元（Deslanoside）进行表型评估。
结果： 四种分子亚型展现出不同的表达谱、临床构成和肿瘤微环境特征。C4亚型协同富集了侵袭、间充质、缺氧和基质相关程序，并伴随与免疫效应活性受限相关的算法特征。多队列分析鉴定GDI1为与不良预后和进展相关的候选基因。两条独立的shRNA在两种结直肠癌细胞系中均降低了GDI1表达、细胞生长和迁移相关表型。使用shGDI1 #1敲低GDI1也降低了皮下异种移植瘤负荷。地高辛苷元在计算筛选中位列领先候选化合物。地高辛苷元处理在细胞和异种移植模型中产生的抗肿瘤表型与GDI1耗竭产生的表型方向一致。
结论： 本研究 delineates 了结直肠癌淋巴结转移相关的分子异质性，支持GDI1的促肿瘤作用，并为地高辛苷元作为候选化合物的进一步研究提供了基础。
论文解读文章
一、研究背景与科学问题
结直肠癌（colorectal cancer, CRC）在全球癌症发病率和死亡率中占据重要份额。淋巴结转移（lymph node metastasis, LNM）是结直肠癌分期和治疗决策的核心依据，标志着肿瘤从局部生长向区域播散的关键转变。这一过程受到肿瘤细胞状态、基质重塑、淋巴管通路和免疫生态的共同塑造。单一标志物通常只能捕捉这一过程的一个维度，难以充分代表肿瘤间在转移能力和微环境依赖性方面的差异。因此，从淋巴结转移相关表达特征出发定义分子亚型，有助于在统一框架内识别与进展相关的疾病状态，并为后续功能筛选提供更清晰的生物学背景。
既往结直肠癌共识分子亚型（consensus molecular subtypes, CMS）研究揭示了上皮程序、代谢状态、免疫反应和基质浸润的系统性差异，并定义了以转化生长因子-β（transforming growth factor-β, TGFβ）信号和基质重塑为特征的特殊类型。多区域转录组学进一步表明，瘤内转录异质性和微环境组成共同影响CMS分型，其中基质成分对CMS4分类具有重要贡献。整合单细胞和 bulk 转录组研究进一步区分了肿瘤上皮细胞内在状态与纤维化微环境的贡献，并显示结直肠癌细胞状态和微环境组成沿连续谱变化。近期在晚期结直肠癌肝转移模型和患者样本中的研究还表明，TGFβ可同时限制T细胞募集并诱导免疫抑制性SPP1+巨噬细胞状态，从而将纤维化基质与免疫逃逸联系起来。这些发现提示，聚焦淋巴结转移的分子分型应同时解析肿瘤内在程序、基质状态和免疫环境。
从疾病状态转向干预的分子候选靶点，需要将跨队列关联与功能扰动相结合。Rab依赖性囊泡运输调控受体循环、膜蛋白定位和分泌，其失调可促进肿瘤细胞迁移和侵袭。GDP解离抑制因子1（GDP dissociation inhibitor 1, GDI1）是Rab循环的重要调控因子。既往结直肠癌队列研究将其高表达与侵袭性临床特征和不良预后相关联。近期在HCT116和HT29细胞中的siRNA研究也发现，GDI1沉默后增殖和迁移减少、凋亡增加，为促肿瘤作用提供了独立的体外证据。在此基础上，本研究建立了基于淋巴结转移相关表达特征的结直肠癌分子亚型，跨队列筛选候选基因，使用两条独立短发夹RNA（short hairpin RNA, shRNA）在HCT116和SW480细胞中评估GDI1功能，并将验证扩展到皮下异种移植模型。随后利用GDI1结构进行虚拟筛选和分子动力学分析，在细胞和异种移植模型中对地高辛苷元（Deslanoside）进行表型评估。这一设计将分子亚型鉴定、功能候选验证和药理学假说的形成连接起来。
二、主要技术方法概述
本研究整合了多层次的生物信息学分析与实验验证手段。数据来源包括TCGA-CRC队列（584例样本）及GSE12945、GSE13294、GSE14333、GSE17536、GSE26682、GSE39582、GSE41258、GSE87211、GSE92921、GSE143985、GSE161158等多个GEO转录组队列。研究采用limma差异表达分析、单变量Cox回归、共识聚类（consensus clustering）建立分子亚型，使用最近模板预测（nearest template prediction）进行跨队列亚型投射。通路分析采用单样本基因集富集分析（ssGSEA）、PROGENy和CancerSEA功能状态评分。免疫微环境分析整合了ESTIMATE、TIDE和IPS评分。功能实验采用两条独立shRNA在HCT116和SW480细胞中敲低GDI1，通过CCK-8、Transwell迁移和划痕愈合实验评估细胞表型，并建立皮下异种移植模型进行体内验证。药物发现方面，利用AlphaFold预测的GDI1三维结构（UniProt P31150）进行AutoDock Vina 1.2.5虚拟筛选（2,990个化合物），并对地高辛苷元-GDI1复合物进行100 ns分子动力学模拟。
三、研究结果
淋巴结转移相关表达特征定义结直肠癌分子亚型。 研究人员首先评估淋巴结转移相关基因的生存关联。森林图将GABRD、MPP2、PIANP、GDI1和ATP6V1B1与较高风险相关联，而OXSM、NUDT6、LARS2、PSMA5和GSR显示相反方向的关联，其中GDI1位列领先的风险相关基因。基于该特征集的共识聚类在四簇解下显示出高内部一致性，将样本分为C1、C2、C3和C4四种分子亚型。主成分分析和UMAP证实四种亚型在表达空间中占据相对独立的区域，且四种亚型在发现队列中呈现不同的总生存曲线。
分子亚型展现不同的表达谱和临床构成。 每种亚型相对于其余样本呈现独特的差异表达谱。C1以SPOCK1、TNS1、FLT1、HIP1和KCNE4上调为特征，涉及细胞外基质、细胞黏附和血管过程；C2以INAFM1、ZNF316、MOSPD3和YPEL3为代表；C3以PPA2、COX5A、NAT1、UQCRFS1和GSR上调为主，指向线粒体代谢和氧化还原过程；C4以ZFHX4、FAM43B、NRK、PNMA8B和NALCN高表达为特征。将亚型标志基因模板投射至TCGA-CRC和多个GEO队列，均识别出对应的表达模式。临床构成分析显示，与无淋巴结转移（nLNM）样本相比，LNM样本中C2和C4比例更高、C3比例更低；晚期临床分期、M1状态和更高T分期样本中同样呈现C2和C4增加、C3减少的趋势，其中C4的变化最为集中。
C4富集侵袭、基质和缺氧相关分子程序。 Hallmark基因集分析显示C4富集上皮间质转化（EMT）、缺氧、Notch信号、Wnt/β-catenin信号、KRAS信号、肌生成、补体和顶侧连接程序，而C3以氧化磷酸化、细胞周期、MYC靶点、DNA修复、脂肪酸代谢和p53通路为主。KEGG通路分析同样突出C4中的基质相关程序，包括黏着斑、PI3K-Akt信号、Rap1信号、细胞外基质受体相互作用和细胞黏附分子通路。PROGENy分析显示TGFβ和缺氧信号在C4中突出。CancerSEA功能状态评分中，C4在EMT、侵袭、转移、血管生成、炎症和缺氧方面得分较高。体细胞突变谱显示APC、TP53、KRAS、PIK3CA、FAT4和ZFHX4等高频突变基因在不同亚型中分布不同，TP53突变在C2中更为集中。亚型相关基因形成包含FLT4、NGF、COMP、THBS3、THBS4、PRELP、PRG4和BGN等节点的蛋白互作模块，涉及细胞外基质组织和组织重塑。转录因子活性分析显示PRDM14、CTCF和SOX2在C4中活性较高。
C4呈现免疫抑制特征和低预测免疫治疗应答。 TIDE分析显示预测应答者在C3中比例较高，而预测非应答者在C4中占主导。C4的TIDE、Dysfunction、Exclusion和癌症相关成纤维细胞（CAF）评分普遍较高，C3较低。IPS分析显示C4的综合IPS较低，C3较高。ESTIMATE分析显示C4具有较高的基质、免疫和综合评分及较低的肿瘤纯度。免疫调节分子表达谱显示，HLA-DPA1、HLA-DPB1、HLA-DRA等抗原呈递分子和TIGIT、CTLA4、PDCD1、LAG3、HAVCR2等共抑制检查点在C4中高表达，共刺激受体和细胞因子也呈现协调变化。这些特征将C4界定为基质和免疫成分富集、免疫调节信号增强且预测免疫治疗应答低的微环境状态。
多队列筛选鉴定GDI1为与不良预后相关的候选基因。 跨队列风险图谱显示GDI1在GSE17536、GSE39582、GSE41258、GSE71187和TCGA-CRC中风险方向一致，按平均风险排名第一。GDI1在全部五个队列中均显示风险比大于1且具有统计学显著性。Kaplan-Meier分析显示GDI1高表达组在TCGA-CRC和多个GEO队列中总生存较差，且在TCGA-CRC中疾病特异性生存、无病生存和无进展生存均呈不利方向。临床表型层面，GDI1表达随淋巴结分期进展而升高，其表达可区分LNM与nLNM状态，且晚期临床分期样本中GDI1表达更高。基因集富集分析将GDI1高表达与EMT、肌生成、顶侧连接、TNFα/NF-κB信号、凝血和血管生成相关联。
GDI1敲低抑制结直肠癌细胞生长和迁移并降低皮下异种移植瘤生长。 两条独立shRNA在HCT116和SW480细胞中均有效降低GDI1转录水平。CCK-8实验显示GDI1敲低后两种细胞的生长曲线均降低，且差异随时间扩大。Transwell实验显示任一shRNA处理后穿过膜的细胞均显著减少。划痕愈合实验显示GDI1敲低组伤口闭合率低于NC组。在HCT116皮下异种移植模型中，shGDI1 #1组形成更小的终点肿瘤，肿瘤体积曲线更低，终点肿瘤重量相应减少。SW480异种移植模型呈现相同方向的发现。两种细胞系的体外一致发现及其相应的体内致瘤性发现表明GDI1促进结直肠癌细胞生长、迁移和异种移植瘤生长。
虚拟筛选和分子动力学模拟鉴定地高辛苷元为候选GDI1配体。 基于GDI1三维结构的虚拟筛选在2,990个成功对接的化合物中，地高辛苷元和两性霉素B共享最优预测评分?10.40 kcal/mol，苏拉明（?10.37）、佐他莫司（?10.26）、芦丁（?10.25）和来迪派韦（?10.23）紧随其后。预测的地高辛苷元-GDI1复合物将地高辛苷元置于GDI1表面的预测结合口袋内，与GLU233、GLN236、SER45、TYR38、ARG35和ARG55等残基形成预测接触。100 ns分子动力学模拟显示蛋白骨架均方根偏差在初始快速调整后进入相对稳定的波动范围。
地高辛苷元和GDI1敲低产生一致的抗肿瘤表型。 研究人员在HCT116和SW480细胞中建立NC、NC加地高辛苷元、shGDI1 #1加DMSO、shGDI1 #1加地高辛苷元四组。CCK-8实验显示地高辛苷元暴露和GDI1敲低均降低两种细胞的生长曲线，联合组同样维持在低水平。Transwell实验显示地高辛苷元暴露或GDI1敲低组的迁移细胞均少于NC组，GDI1敲低后加用地高辛苷元的迁移细胞数与单独GDI1敲低组相似。划痕愈合实验获得一致结果。在HCT116皮下异种移植瘤中，地高辛苷元给药和GDI1敲低均降低肿瘤生长曲线和终点肿瘤重量，联合组的终点肿瘤重量与GDI1敲低组相似。SW480异种移植瘤重现相同方向的改变。
四、讨论与结论总结
讨论部分指出，C4亚型最突出的特征是侵袭、间充质、缺氧和基质相关程序的协同激活，并在淋巴结转移和晚期临床特征样本中富集。C4与经典CMS4的相似性指共享的基质和间充质程序，而非分类一致性或等价性。C4的另一个defining特征是基质富集与免疫效应活性受限并存，提示基质富集可能与抗肿瘤免疫活性受限共存，为在该亚型中研究TGFβ和CAF相关机制提供了依据。
GDI1将分子亚型所代表的进展相关状态与可重复的功能表型联系起来。本研究在多个队列中发现一致的关联，在HCT116和SW480细胞中用两条独立shRNA观察到一致的生长和迁移相关表型，并在两种皮下异种移植模型中观察到较低的肿瘤负荷。这些发现通过独立扰动复制、额外结直肠癌细胞系和体内肿瘤生长证据，扩展了既往队列关联和体外功能证据。
药理学研究引入地高辛苷元作为候选干预。既往研究在前列腺癌细胞和异种移植模型中观察到地高辛苷元处理后生长和迁移相关表型的变化。相关强心苷Oleandrin也在结直肠癌细胞和患者来源结肠癌类器官中引发抗肿瘤表型，涉及Na+/K+-ATP酶和GRP78应激调控。强心苷的经典药理学涉及Na+/K+-ATP酶抑制和下游信号调控。本研究中，对接和分子动力学提供了预测相互作用，地高辛苷元处理和GDI1耗竭产生方向一致的表型，但直接结合和GDI1依赖性药物作用在缺乏靶点参与和拯救实验的情况下仍未证实。强心苷药理学也引发全身抗癌治疗的安全性担忧，瘤内异种移植实验不能定义肿瘤与心脏暴露、临床可达到的全身浓度或心脏安全性。
研究结论：本研究 delineates 了结直肠癌淋巴结转移相关的分子异质性，支持GDI1的促肿瘤作用，并为地高辛苷元作为候选化合物的进一步研究提供了基础。该研究将表达关联推进至功能验证和药理学假说，为结直肠癌淋巴结转移相关亚型的干预策略提供了新的研究方向。
