《Signal Transduction and Targeted Therapy》：Mechanisms and therapeutic targets in distant metastasis of lung adenocarcinoma
编辑推荐：
远处转移以器官趋向性为特征，是肺腺癌（LUAD）死亡的主要原因。在本研究中，研究人员整合了来自52例LUAD患者的数字空间分析（DSP）、多重免疫荧光（mIF）和临床数据，以开发器官特异性转移风险模型，并研究了转移性器官趋向性的分子机制。基于原发肿瘤空间转录组
远处转移以器官趋向性为特征，是肺腺癌（LUAD）死亡的主要原因。在本研究中，研究人员整合了来自52例LUAD患者的数字空间分析（DSP）、多重免疫荧光（mIF）和临床数据，以开发器官特异性转移风险模型，并研究了转移性器官趋向性的分子机制。基于原发肿瘤空间转录组学的随机森林模型准确预测了向脑（AUC=0.974）、肝（AUC=0.975）、肾上腺（AUC=0.929）和骨（AUC=0.907）的转移。研究人员鉴定了与器官趋向性转移相关的关键区室特异性基因表达特征，包括在肿瘤中表达的（例如，脑转移的FKBP1A和肝转移的MOCOS）、免疫细胞中表达的（例如，肝转移的ADAMTSL2）和基质细胞中表达的（例如，脑转移的CKAP2）。通路分析揭示了与器官趋向性相关的不同生物学过程，例如脑转移中细胞死亡通路的富集和肝转移中细胞外基质（ECM）重塑。转移后生存模型强调基质基因表达（例如，总生存期（OS）的PKM和无进展生存期（PFS）的VCAM1）和免疫抑制微环境（例如，M2巨噬细胞浸润）作为关键预后因素。本研究中鉴定的高精度预测模型和关键分子特征增强了对“种子-土壤”相互作用动态的理解，并为未来的临床应用提供了有前景的生物标志物和治疗靶点。
**研究背景与问题**
肺腺癌（LUAD）远处转移是导致患者死亡的主要原因，其过程呈现明显的器官趋向性（organotropism），即肿瘤细胞倾向于向特定器官（如脑、肝、肾上腺、骨）转移。尽管“种子-土壤”假说已提出多年，但现有研究仍无法准确基于原发肿瘤特征预测器官特异性转移风险，且对转移后生存的预测模型十分匮乏。此外，不同转移灶对治疗反应的异质性及其分子机制尚未阐明。为此，本研究旨在通过整合空间多组学数据，系统解析LUAD器官趋向性转移的分子机制，并构建高精度预测模型，为临床提供新的生物标志物和治疗靶点。该论文发表在《Signal Transduction and Targeted Therapy》。
**主要关键技术方法**
研究人员从52例晚期LUAD患者中收集了匹配的原发肿瘤与转移灶样本（包括23对脑转移、26对肝转移、3对肾上腺转移及相应原发肿瘤，同时获取邻近正常组织作为对照）。主要采用以下技术：（1）数字空间分析（DSP）技术，对肿瘤、免疫、基质三个区室进行高分辨率空间转录组学分析；（2）多重免疫荧光（mIF）染色，评估不同免疫细胞亚型（如M1/M2巨噬细胞、树突状细胞、T细胞等）的浸润比例；（3）基于H&E染色计算HEindex（每个肿瘤细胞周围非肿瘤细胞的平均比例），量化肿瘤微环境的空间结构。基于这些数据，研究人员构建了随机森林模型用于器官特异性转移风险预测和转移后生存预测。
**研究结果**
**1. 脑转移风险模型及关键特征**
通过随机森林模型，研究人员鉴定出与LUAD脑转移风险最相关的基因，包括肿瘤细胞表达的FKBP1A和MRPS33，免疫细胞表达的AFAP1L1、ELP6等，以及基质细胞表达的CKAP2。该模型在测试队列中AUC达0.974。基因集富集分析（GSEA）显示，脑转移高危原发肿瘤的肿瘤、免疫和基质区室中细胞死亡通路显著富集。多重免疫荧光显示，肿瘤相关成纤维细胞I型（CAFI）和正常成纤维细胞（NF）比例与模型风险评分正相关。通过多层次筛选，研究人员鉴定出78个核心差异表达基因（DEGs），并构建了基因调控网络，发现脑转移后肿瘤细胞中HNRNPA2B1、H3C13、CYCS上调，BST2、COL3A1下调；免疫细胞中HSPA8、HSP90AA1等上调；基质细胞中RPS2、COL3A1等下调，提示蛋白折叠、应激反应及细胞外基质（ECM）相互作用的重编程。
**2. 肝转移风险模型及关键特征**
肝转移风险模型的关键基因包括肿瘤细胞表达的MOCOS、PDGFRB等，以及免疫细胞表达的ADAMTSL2、SPCS1等，模型AUC达0.975。GSEA显示，肝转移高危原发肿瘤的肿瘤区室中细胞死亡、代谢与能量通路下调，而免疫区室中细胞死亡和免疫通路下调。mIF显示，M1巨噬细胞和髓系树突状细胞（mDCs）比例在肝转移组显著升高，且与风险评分正相关。通过DEGs交集，鉴定出72个核心基因，富集于DNA复制依赖性染色质组装和核小体通路。基因调控网络显示，肝转移后肿瘤细胞中组蛋白基因（H3C15、H2AC19等）和核糖体蛋白基因（RPL36A、RPL38等）上调，而LCN12、SERPINA5下调；免疫和基质区室中也出现类似上调模式，提示转录抑制和癌细胞休眠是肝转移的关键适应机制。
**3. 肾上腺转移风险模型**
肾上腺转移风险模型的关键基因包括肿瘤细胞表达的PTGR2、IGKC，免疫细胞表达的ADIG、ASB3等，以及基质细胞表达的IRF3、MSLN等，模型AUC为0.929。GSEA显示，高危原发肿瘤的基质区室中细胞死亡通路富集。mIF发现，mDCs比例在肾上腺转移组显著升高，而浆细胞样树突状细胞（pDCs）比例与风险评分负相关。HEindex与风险评分正相关。通过DEGs交集（296个基因），富集到凋亡和氧化还原酶通路。由于样本量小，未进行加权调控网络分析。
**4. 骨转移风险模型**
骨转移风险模型的关键基因包括肿瘤细胞表达的ARMCX6、PLPP1，免疫细胞表达的DDIAS、POLE等，以及基质细胞表达的B3GNT5，模型AUC为0.907。GSEA未发现显著富集的通路，提示骨转移倾向可能源于更广泛的生物学过程。mIF显示M1巨噬细胞比例在骨转移组升高，但与风险评分无显著相关性。由于缺乏骨转移灶组织，无法进行配对分析。
**5. 转移后生存预测模型**
通过配对比较原发与转移灶的DSP数据，研究人员鉴定出转移过程中显著改变的DEGs，并基于其中FDR<0.01的基因构建随机森林生存模型。总生存期（OS）相关关键基因包括ATP8A1_Immune、EPC1_Stroma、PKM_Stroma等；无进展生存期（PFS）相关基因包括MDH1_Tumor、VCAM1_Stroma等。GSEA显示，OS不良与肿瘤和基质区室中细胞周期、遗传及表观遗传信息通路上调相关；PFS不良则与细胞死亡、代谢等通路相关。mIF分析显示，转移灶中M2巨噬细胞和效应T细胞浸润与OS风险评分负相关，M2巨噬细胞、效应T细胞和调节性T细胞与PFS风险评分负相关，提示免疫微环境对生存的复杂影响。HEindex与OS和PFS风险评分均负相关。
**讨论与结论**
讨论部分指出，本研究的主要发现是原发肿瘤的空间转录组特征可高精度预测器官特异性转移风险（AUC>0.9），支持了“种子-土壤”假说中转移程序已编码于原发肿瘤分子蓝图的观点。不同器官转移的分子机制存在异质性：脑转移涉及细胞死亡通路和蛋白折叠重编程，肝转移涉及组蛋白和核糖体蛋白上调及转录抑制，肾上腺转移涉及凋亡和氧化还原通路，而骨转移未发现特定通路。转移后生存模型显示基质基因（如PKM、VCAM1）和免疫抑制微环境（M2巨噬细胞浸润与较好生存相关，但需谨慎解释）是重要预后因素。研究还提出了一个两步预测框架：初始诊断时基于原发肿瘤进行器官特异性风险分层，指导监测和预防措施；转移后基于转移灶构建生存预测模型，指导治疗决策。HEindex作为新参数，可量化肿瘤微环境中的非肿瘤细胞浸润，并与肾上腺转移风险正相关，与转移后生存负相关。研究局限性包括样本量较小（尤其肾上腺和骨转移）、缺乏多组学数据（如基因组学、表观基因组学），以及关键基因的功能验证不足。未来需要更大规模队列和体外/体内模型验证。
**研究结论**（翻译原文结论部分）
总之，从独特的空间转录组学视角，本研究揭示了原发LUAD病灶的分子特征与器官特异性转移潜能相关，而转移灶的分子和免疫特征影响转移后生存。本工作中鉴定的高精度预测模型和关键分子特征不仅增强了对“种子-土壤”相互作用时空动态的理解，也为未来的临床应用提供了有前景的生物标志物和治疗靶点。
